From 897fd186990518f7c77aa2e73423f4a3a109a3d1 Mon Sep 17 00:00:00 2001
From: Martin Schupfner <k204212@levante2.atos.local>
Date: Thu, 1 Dec 2022 19:19:07 +0100
Subject: [PATCH] Updated mpiesm cmor scripts, cmor tables and CV

- CV:
  + added interactive transient deglaciation experiments
  + added AWI-ESM source_id
- cmor tables: added Icen* and Iyr* tables
- cmor mpiesm:
  + cdocmorinfo: moved calendar to submodel specific config files
  + aggr/runpp: moved jsbach input files definition to runpp script
  + runpp: ncatted - manually changing frequency to 'cen' since CMOR does not support centennial frequency
---
 cmor/mpiesm/scripts/cdocmorinfo_echam6_CR     |   2 +
 cmor/mpiesm/scripts/cdocmorinfo_jsbach_CR     |   2 +
 cmor/mpiesm/scripts/cdocmorinfo_mpiom_CR      |   2 +
 cmor/mpiesm/scripts/cdocmorinfo_mpism_CR      |   3 +-
 .../conf/PalMod2_requested_vars_AllExp.conf   |  87 +----------
 .../palmod2_aggregation_MPI-ESM1-2_jsbach.h   |   9 --
 cmor/mpiesm/scripts/pismtest.runpp            |   5 +-
 cmor/mpiesm/scripts/pismtest_cdocmorinfo      |   1 -
 .../ssp119-future-methane_r1i1p1f1-CR.runpp   |  15 +-
 ...119-future-methane_r1i1p1f1-CR_cdocmorinfo |   2 +-
 ...rescribed-glac1d-methane_r1i1p1f1-CR.runpp |  15 +-
 ...bed-glac1d-methane_r1i1p1f1-CR_cdocmorinfo |   2 +-
 ...iation-prescribed-glac1d_r1i1p1f1-CR.runpp |  15 +-
 ...-prescribed-glac1d_r1i1p1f1-CR_cdocmorinfo |   2 +-
 ...iation-prescribed-glac1d_r1i1p2f2-CR.runpp |  15 +-
 ...-prescribed-glac1d_r1i1p2f2-CR_cdocmorinfo |   2 +-
 ...iation-prescribed-glac1d_r1i1p3f2-CR.runpp |  15 +-
 ...-prescribed-glac1d_r1i1p3f2-CR_cdocmorinfo |   2 +-
 ...ciation-prescribed-ice6g_r1i1p1f1-CR.runpp |  15 +-
 ...n-prescribed-ice6g_r1i1p1f1-CR_cdocmorinfo |   2 +-
 ...ciation-prescribed-ice6g_r1i1p2f2-CR.runpp |  15 +-
 ...n-prescribed-ice6g_r1i1p2f2-CR_cdocmorinfo |   2 +-
 ...ciation-prescribed-ice6g_r1i1p3f2-CR.runpp |  15 +-
 ...n-prescribed-ice6g_r1i1p3f2-CR_cdocmorinfo |   2 +-
 cmor_tables/PalMod2_CV.json                   |  76 +++++++++
 cmor_tables/PalMod2_IcenAnt.json              | 144 ++++++++++++++++++
 cmor_tables/PalMod2_IcenGre.json              | 144 ++++++++++++++++++
 cmor_tables/PalMod2_IyrAnt.json               | 144 ++++++++++++++++++
 cmor_tables/PalMod2_IyrGre.json               | 144 ++++++++++++++++++
 29 files changed, 791 insertions(+), 108 deletions(-)
 create mode 100644 cmor_tables/PalMod2_IcenAnt.json
 create mode 100644 cmor_tables/PalMod2_IcenGre.json
 create mode 100644 cmor_tables/PalMod2_IyrAnt.json
 create mode 100644 cmor_tables/PalMod2_IyrGre.json

diff --git a/cmor/mpiesm/scripts/cdocmorinfo_echam6_CR b/cmor/mpiesm/scripts/cdocmorinfo_echam6_CR
index 5f8fb26..9721adb 100644
--- a/cmor/mpiesm/scripts/cdocmorinfo_echam6_CR
+++ b/cmor/mpiesm/scripts/cdocmorinfo_echam6_CR
@@ -1,3 +1,5 @@
 #NOMINAL_RESOLUTION read from external file
 #For Atm/Srf (ECHAM6/JSBACH)
 NOMINAL_RESOLUTION="500 km"
+# Calendar
+CALENDAR=proleptic_gregorian
diff --git a/cmor/mpiesm/scripts/cdocmorinfo_jsbach_CR b/cmor/mpiesm/scripts/cdocmorinfo_jsbach_CR
index 10c944a..6b464bc 100644
--- a/cmor/mpiesm/scripts/cdocmorinfo_jsbach_CR
+++ b/cmor/mpiesm/scripts/cdocmorinfo_jsbach_CR
@@ -5,3 +5,5 @@ NOMINAL_RESOLUTION="500 km"
 grid_info="/work/bm0021/PalMod2/cmor/mpiesm/scripts/gridinfo_JSBACH_L5.nc"
 switch_z=y
 switch_xy=n
+# Calendar
+CALENDAR=proleptic_gregorian
diff --git a/cmor/mpiesm/scripts/cdocmorinfo_mpiom_CR b/cmor/mpiesm/scripts/cdocmorinfo_mpiom_CR
index e99ee39..75885ad 100644
--- a/cmor/mpiesm/scripts/cdocmorinfo_mpiom_CR
+++ b/cmor/mpiesm/scripts/cdocmorinfo_mpiom_CR
@@ -5,3 +5,5 @@ NOMINAL_RESOLUTION="500 km"
 grid_info="/work/bm0021/PalMod2/cmor/mpiesm/scripts/gridinfo_MPIOM_L40.nc"
 switch_z=y
 switch_xy=n
+# Calendar
+CALENDAR=proleptic_gregorian
diff --git a/cmor/mpiesm/scripts/cdocmorinfo_mpism_CR b/cmor/mpiesm/scripts/cdocmorinfo_mpism_CR
index 2fe30f8..6f8cbb8 100644
--- a/cmor/mpiesm/scripts/cdocmorinfo_mpism_CR
+++ b/cmor/mpiesm/scripts/cdocmorinfo_mpism_CR
@@ -1,8 +1,9 @@
 #NOMINAL_RESOLUTION read from external file
 #For ice (mPISM)
-CALENDAR="365day"
 NOMINAL_RESOLUTION="25 km"
 #Level/Depth Bounds
 grid_info="/work/bm0021/PalMod2/cmor/mpiesm/scripts/gridinfo_mPISM_L121.nc"
 switch_z=y
 switch_xy=n
+# Calendar
+CALENDAR=proleptic_gregorian
diff --git a/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExp.conf b/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExp.conf
index f8f6e3a..5d8c839 100755
--- a/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExp.conf
+++ b/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExp.conf
@@ -33,97 +33,24 @@ TOTAL    : TOTAL  = False # only produce variables explicitly set to True
 ################
 # Aggregation
 ################
-echammonplev19  : interpolation = True
-echammonafter   : afterburner   = True
 jsbachdecfx     : input         = True
 echamdecfx      : input         = True
-Omon            : areacello     = True
 ################
 # Ocean
 ################
-#basin    : Ofx    = True
-#sltbasin : EmonZ  = True
-areacello: Ofx    = True
-#volcello : Ofx    = True
-#volcello : Omon   = True
-sos      : Omon   = True
-hfds     : Omon   = True
-tauvo    : Omon   = True
-tos      : Omon   = True
-zos      : Omon   = True
-siconc   : SImon  = True
-simass   : SImon  = True
-sisnmass : SImon  = True
-#sisnconc : SImon  = True
-tauuo    : Omon   = True
-sisnthick: SImon  = True
-sithick  : SImon  = True
-siu      : SImon  = True
-siv      : SImon  = True
-hfx      : Omon   = True
-hfy      : Omon   = True
-sftof    : Ofx    = True
 sftof    : Odec   = True
 ################
 # Atmosphere
 ################
-#hurs     : Amon   = True
-#huss     : Amon   = True
-cl       : Amon   = True
-clt      : Amon   = True
-hfls     : Amon   = True
-hfss     : Amon   = True
-hur      : Amon   = True
-hus      : Amon   = True
-prc      : Amon   = True
-prsn     : Amon   = True
+orog     : dec    = True
 ps       : Amon   = True
-psl      : Amon   = True
-rlus     : Amon   = True
-rlut     : Amon   = True
-rsdt     : Amon   = True
-rsus     : Amon   = True
-rsut     : Amon   = True
-rlds     : Amon   = True
-rsds     : Amon   = True
-sfcWind  : Amon   = True
-ta       : Amon   = True
-tasmax   : Amon   = True
-tasmin   : Amon   = True
-uas      : Amon   = True
-vas      : Amon   = True
-wap      : Amon   = True
-zg       : Amon   = True
-tdps     : Emon   = True
-snm      : LImon  = True
-snd      : LImon  = True
-mrro     : Lmon   = True
-sitemptop : SImon = True
-evspsbl  : Amon   = True
-snw      : LImon  = True
-mrros    : Lmon   = True
-ua       : AERmon = True
-va       : AERmon = True
 #################
 # Land
 #################
-snc           : LImon  = True
-sbl           : LImon  = True
-tran          : Lmon   = True
-mrsofc        : fx     = True
-mrsofc        : dec    = True
-rootd         : fx     = True
-rootd         : dec    = True
-sftgif        : fx     = True
-sftgif        : dec    = True
-sftlf         : fx     = True
-sftlf         : dec    = True
-lai           : Lmon   = True
-tsl           : Lmon   = True
-vegFrac       : Emon   = True
-grassFrac     : Lmon   = True
-landCoverFrac : Lmon   = True
-shrubFrac     : Lmon   = True
-treeFrac      : Lmon   = True
-mrso          : Lmon   = True
+sftlf    : dec    = True
+mrsofc   : dec    = True
+rootd    : dec    = True
+sftgif   : dec    = True
+landCoverFrac : Lmon = True
+vegFrac       : Emon = True
 #<---- Specify your settings for Experiment AllExp here
diff --git a/cmor/mpiesm/scripts/palmod2_aggregation_MPI-ESM1-2_jsbach.h b/cmor/mpiesm/scripts/palmod2_aggregation_MPI-ESM1-2_jsbach.h
index a3cbaf3..c5acc53 100644
--- a/cmor/mpiesm/scripts/palmod2_aggregation_MPI-ESM1-2_jsbach.h
+++ b/cmor/mpiesm/scripts/palmod2_aggregation_MPI-ESM1-2_jsbach.h
@@ -4,15 +4,6 @@
 suf=grb
 
 { (if_requested $member $srfmod jsbachdecfx input $chunk '???1123124' && {
-  # Read input file
-  if [[ "${iniyeararr[$y]}" == "$period" ]]
-  then
-    echo "Reading from parent ${RAW_EXP_ID_PARENT[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
-    ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_ID_PARENT[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
-  else
-    echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
-    ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
-  fi
 
   # Create decadal "fx"-variables
 
diff --git a/cmor/mpiesm/scripts/pismtest.runpp b/cmor/mpiesm/scripts/pismtest.runpp
index 1be7bbb..3955a57 100755
--- a/cmor/mpiesm/scripts/pismtest.runpp
+++ b/cmor/mpiesm/scripts/pismtest.runpp
@@ -603,7 +603,10 @@ echo "-> Looping over ${#flist[@]} netCDF files!"
 for f in ${flist[@]}
 do
     echo "... editing $f"
-    ncatted -O -h -a further_info_url,global,d,, $f && echo "... successful!" || echo "... failed (ERROR)"
+    ncatted -O -h -a further_info_url,global,d,, $f && echo "... deleted 'further_info_url'." || echo "... failed (ERROR)"
+    [[ "$(basename $f)" == *"cen"* ]] && {
+        ncatted -O -h -a frequency,global,m,c,"cen" $f && echo "... changed frequency to 'cen'." || echo "... failed (ERROR)"
+    }
 done
 
 print "$(date +%Y-%m-%dT%H:%M:%S): NCO cleanup ended for $dr"
diff --git a/cmor/mpiesm/scripts/pismtest_cdocmorinfo b/cmor/mpiesm/scripts/pismtest_cdocmorinfo
index 5849d13..89b47fe 100644
--- a/cmor/mpiesm/scripts/pismtest_cdocmorinfo
+++ b/cmor/mpiesm/scripts/pismtest_cdocmorinfo
@@ -48,7 +48,6 @@ REFERENCES="Kapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022)
 SOURCE="MPI-ESM1.2-CR (2017): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none\nland: JSBACH3.20\nlandIce: none/prescribed\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none/prescribed"
 SOURCE_TYPE="AOGCM"
 #CALENDAR=proleptic_gregorian
-CALENDAR="365_day"
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
 #MAPPING_TABLE=MPIESM1_CMIP6.txt
 #GRID_FILE_DIR=__GRID_FILE_DIR
diff --git a/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR.runpp b/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR.runpp
index a52b245..82519da 100755
--- a/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR.runpp
@@ -332,6 +332,16 @@ function run_agg
         err="${errtemp_agg}/err${period}"
         rm -f ${err}*
 
+        # Read input file
+        if [[ "${iniyeararr[$y]}" == "$period" ]]
+        then
+            echo "Reading from parent ${RAW_EXP_ID_PARENT[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/bm1030/from_Mistral/bm1030/m300020/mpiesm-palmod_tk_r10186_megan_next/experiments/${RAW_EXP_ID_PARENT[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        else
+            echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/bm1030/from_Mistral/bm1030/m300020/mpiesm-palmod_tk_r10186_megan_next/experiments/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        fi
+
         # Perform aggregation for each submodel
         for submodel in $atmmod $ocemod $srfmod; do # $bgcmod; do
 
@@ -528,7 +538,10 @@ echo "-> Looping over ${#flist[@]} netCDF files!"
 for f in ${flist[@]}
 do
     echo "... editing $f"
-    ncatted -O -h -a further_info_url,global,d,, $f && echo "... successful!" || echo "... failed (ERROR)"
+    ncatted -O -h -a further_info_url,global,d,, $f && echo "... deleted 'further_info_url'." || echo "... failed (ERROR)"
+    [[ "$(basename $f)" == *"cen"* ]] && {
+        ncatted -O -h -a frequency,global,m,c,"cen" $f && echo "... changed frequency to 'cen'." || echo "... failed (ERROR)"
+    }
 done
 
 print "$(date +%Y-%m-%dT%H:%M:%S): NCO cleanup ended for $dr"
diff --git a/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR_cdocmorinfo b/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR_cdocmorinfo
index ef41c34..813ef38 100644
--- a/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR_cdocmorinfo
@@ -58,7 +58,7 @@ SOURCE_ID=MPI-ESM1-2-CR
 REFERENCES="Kapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022). Ocean response in transient simulations of the last deglaciation dominated by underlying ice-sheet reconstruction and method of meltwater distribution. Geophysical Research Letters, 49, e2021GL096767. https://doi.org/10.1029/2021GL096767\nMPI-ESM: Mauritsen, T. et al. (2019), Developments in the MPI‐M Earth System Model version 1.2 (MPI‐ESM1.2) and Its Response to Increasing CO2, J. Adv. Model. Earth Syst.,11, 998-1038, doi:10.1029/2018MS001400,\nMueller, W.A. et al. (2018): A high‐resolution version of the Max Planck Institute Earth System Model MPI‐ESM1.2‐HR. J. Adv. Model. EarthSyst.,10,1383–1413, doi:10.1029/2017MS001217"
 SOURCE="MPI-ESM1.2-CR (2017): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none\nland: JSBACH3.20\nlandIce: none/prescribed\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none/prescribed"
 SOURCE_TYPE="AOGCM CHEM"
-CALENDAR=proleptic_gregorian
+#CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
 #MAPPING_TABLE=MPIESM1_CMIP6.txt
 #GRID_FILE_DIR=__GRID_FILE_DIR
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR.runpp
index 13fa383..e98209b 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR.runpp
@@ -346,6 +346,16 @@ function run_agg
         err="${errtemp_agg}/err${period}"
         rm -f ${err}*
 
+        # Read input file
+        if [[ "${iniyeararr[$y]}" == "$period" ]]
+        then
+            echo "Reading from parent ${RAW_EXP_ID_PARENT[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/bm1030/from_Mistral/bm1030/m300020/mpiesm-palmod_tk_r10081/experiments/${RAW_EXP_ID_PARENT[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        else
+            echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/bm1030/from_Mistral/bm1030/m300020/mpiesm-palmod_tk_r10081/experiments/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        fi
+
         # Perform aggregation for each submodel
         for submodel in $atmmod $ocemod $srfmod; do # $bgcmod; do
 
@@ -542,7 +552,10 @@ echo "-> Looping over ${#flist[@]} netCDF files!"
 for f in ${flist[@]}
 do
     echo "... editing $f"
-    ncatted -O -h -a further_info_url,global,d,, $f && echo "... successful!" || echo "... failed (ERROR)"
+    ncatted -O -h -a further_info_url,global,d,, $f && echo "... deleted 'further_info_url'." || echo "... failed (ERROR)"
+    [[ "$(basename $f)" == *"cen"* ]] && {
+        ncatted -O -h -a frequency,global,m,c,"cen" $f && echo "... changed frequency to 'cen'." || echo "... failed (ERROR)"
+    }
 done
 
 print "$(date +%Y-%m-%dT%H:%M:%S): NCO cleanup ended for $dr"
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR_cdocmorinfo b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR_cdocmorinfo
index 73e52ea..c416f66 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR_cdocmorinfo
@@ -52,7 +52,7 @@ SOURCE_ID=MPI-ESM1-2-CR
 REFERENCES="Kleinen, T., Mikolajewicz, U. and Brovkin, V. (2020). Terrestrial methane emissions from the Last Glacial Maximum to the preindustrial period. Climate of the Past, 16(2):575–595. doi: 10.5194/cp-16-575-2020.\nKleinen, T., Gromov, S., Steil, B. and Brovkin, V. (2021). Atmospheric methane underestimated in future climate projections. Environmental Research Letters, 16(9):094006. doi: 10.1088/1748-9326/ac1814.\nKapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022). Ocean response in transient simulations of the last deglaciation dominated by underlying ice-sheet reconstruction and method of meltwater distribution. Geophysical Research Letters, 49, e2021GL096767. https://doi.org/10.1029/2021GL096767\nMPI-ESM: Mauritsen, T. et al. (2019), Developments in the MPI‐M Earth System Model version 1.2 (MPI‐ESM1.2) and Its Response to Increasing CO2, J. Adv. Model. Earth Syst.,11, 998-1038, doi:10.1029/2018MS001400,\nMueller, W.A. et al. (2018): A high‐resolution version of the Max Planck Institute Earth System Model MPI‐ESM1.2‐HR. J. Adv. Model. EarthSyst.,10,1383–1413, doi:10.1029/2017MS001217"
 SOURCE="MPI-ESM1.2-CR (2017): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none\nland: JSBACH3.20\nlandIce: none/prescribed\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none/prescribed" # todo: replace
 SOURCE_TYPE="AOGCM CHEM"
-CALENDAR=proleptic_gregorian
+#CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
 #MAPPING_TABLE=MPIESM1_CMIP6.txt
 #GRID_FILE_DIR=__GRID_FILE_DIR
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp
index ffb9bd9..60866f5 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp
@@ -345,6 +345,16 @@ function run_agg
         err="${errtemp_agg}/err${period}"
         rm -f ${err}*
 
+        # Read input file
+        if [[ "${iniyeararr[$y]}" == "$period" ]]
+        then
+            echo "Reading from parent ${RAW_EXP_ID_PARENT[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_ID_PARENT[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        else
+            echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        fi
+
         # Perform aggregation for each submodel
         for submodel in $atmmod $ocemod $srfmod; do # $bgcmod; do
 
@@ -541,7 +551,10 @@ echo "-> Looping over ${#flist[@]} netCDF files!"
 for f in ${flist[@]}
 do
     echo "... editing $f"
-    ncatted -O -h -a further_info_url,global,d,, $f && echo "... successful!" || echo "... failed (ERROR)"
+    ncatted -O -h -a further_info_url,global,d,, $f && echo "... deleted 'further_info_url'." || echo "... failed (ERROR)"
+    [[ "$(basename $f)" == *"cen"* ]] && {
+        ncatted -O -h -a frequency,global,m,c,"cen" $f && echo "... changed frequency to 'cen'." || echo "... failed (ERROR)"
+    }
 done
 
 print "$(date +%Y-%m-%dT%H:%M:%S): NCO cleanup ended for $dr"
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR_cdocmorinfo b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR_cdocmorinfo
index eadd02b..a93a9c3 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR_cdocmorinfo
@@ -47,7 +47,7 @@ SOURCE_ID=MPI-ESM1-2-CR
 REFERENCES="Kapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022). Ocean response in transient simulations of the last deglaciation dominated by underlying ice-sheet reconstruction and method of meltwater distribution. Geophysical Research Letters, 49, e2021GL096767. https://doi.org/10.1029/2021GL096767\nMPI-ESM: Mauritsen, T. et al. (2019), Developments in the MPI‐M Earth System Model version 1.2 (MPI‐ESM1.2) and Its Response to Increasing CO2, J. Adv. Model. Earth Syst.,11, 998-1038, doi:10.1029/2018MS001400,\nMueller, W.A. et al. (2018): A high‐resolution version of the Max Planck Institute Earth System Model MPI‐ESM1.2‐HR. J. Adv. Model. EarthSyst.,10,1383–1413, doi:10.1029/2017MS001217"
 SOURCE="MPI-ESM1.2-CR (2017): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none\nland: JSBACH3.20\nlandIce: none/prescribed\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none/prescribed"
 SOURCE_TYPE="AOGCM"
-CALENDAR=proleptic_gregorian
+#CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
 #MAPPING_TABLE=MPIESM1_CMIP6.txt
 #GRID_FILE_DIR=__GRID_FILE_DIR
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp
index a7b2f57..3bbbcf9 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp
@@ -347,6 +347,16 @@ function run_agg
         err="${errtemp_agg}/err${period}"
         rm -f ${err}*
 
+        # Read input file
+        if [[ "${iniyeararr[$y]}" == "$period" ]]
+        then
+            echo "Reading from parent ${RAW_EXP_ID_PARENT[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_ID_PARENT[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        else
+            echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        fi
+
         # Perform aggregation for each submodel
         for submodel in $atmmod $ocemod $srfmod; do # $bgcmod; do
 
@@ -543,7 +553,10 @@ echo "-> Looping over ${#flist[@]} netCDF files!"
 for f in ${flist[@]}
 do
     echo "... editing $f"
-    ncatted -O -h -a further_info_url,global,d,, $f && echo "... successful!" || echo "... failed (ERROR)"
+    ncatted -O -h -a further_info_url,global,d,, $f && echo "... deleted 'further_info_url'." || echo "... failed (ERROR)"
+    [[ "$(basename $f)" == *"cen"* ]] && {
+        ncatted -O -h -a frequency,global,m,c,"cen" $f && echo "... changed frequency to 'cen'." || echo "... failed (ERROR)"
+    }
 done
 
 print "$(date +%Y-%m-%dT%H:%M:%S): NCO cleanup ended for $dr"
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR_cdocmorinfo b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR_cdocmorinfo
index 56e3791..943b864 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR_cdocmorinfo
@@ -47,7 +47,7 @@ SOURCE_ID=MPI-ESM1-2-CR
 REFERENCES="Kapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022). Ocean response in transient simulations of the last deglaciation dominated by underlying ice-sheet reconstruction and method of meltwater distribution. Geophysical Research Letters, 49, e2021GL096767. https://doi.org/10.1029/2021GL096767\nMPI-ESM: Mauritsen, T. et al. (2019), Developments in the MPI‐M Earth System Model version 1.2 (MPI‐ESM1.2) and Its Response to Increasing CO2, J. Adv. Model. Earth Syst.,11, 998-1038, doi:10.1029/2018MS001400,\nMueller, W.A. et al. (2018): A high‐resolution version of the Max Planck Institute Earth System Model MPI‐ESM1.2‐HR. J. Adv. Model. EarthSyst.,10,1383–1413, doi:10.1029/2017MS001217"
 SOURCE="MPI-ESM1.2-CR (2017): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none\nland: JSBACH3.20\nlandIce: none/prescribed\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none/prescribed"
 SOURCE_TYPE="AOGCM"
-CALENDAR=proleptic_gregorian
+#CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
 #MAPPING_TABLE=MPIESM1_CMIP6.txt
 #GRID_FILE_DIR=__GRID_FILE_DIR
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp
index ab08948..9629828 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp
@@ -347,6 +347,16 @@ function run_agg
         err="${errtemp_agg}/err${period}"
         rm -f ${err}*
 
+        # Read input file
+        if [[ "${iniyeararr[$y]}" == "$period" ]]
+        then
+            echo "Reading from parent ${RAW_EXP_ID_PARENT[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_ID_PARENT[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        else
+            echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        fi
+
         # Perform aggregation for each submodel
         for submodel in $atmmod $ocemod $srfmod; do # $bgcmod; do
 
@@ -543,7 +553,10 @@ echo "-> Looping over ${#flist[@]} netCDF files!"
 for f in ${flist[@]}
 do
     echo "... editing $f"
-    ncatted -O -h -a further_info_url,global,d,, $f && echo "... successful!" || echo "... failed (ERROR)"
+    ncatted -O -h -a further_info_url,global,d,, $f && echo "... deleted 'further_info_url'." || echo "... failed (ERROR)"
+    [[ "$(basename $f)" == *"cen"* ]] && {
+        ncatted -O -h -a frequency,global,m,c,"cen" $f && echo "... changed frequency to 'cen'." || echo "... failed (ERROR)"
+    }
 done
 
 print "$(date +%Y-%m-%dT%H:%M:%S): NCO cleanup ended for $dr"
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR_cdocmorinfo b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR_cdocmorinfo
index 400aa3e..3219c9a 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR_cdocmorinfo
@@ -47,7 +47,7 @@ SOURCE_ID=MPI-ESM1-2-CR
 REFERENCES="Kapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022). Ocean response in transient simulations of the last deglaciation dominated by underlying ice-sheet reconstruction and method of meltwater distribution. Geophysical Research Letters, 49, e2021GL096767. https://doi.org/10.1029/2021GL096767\nMPI-ESM: Mauritsen, T. et al. (2019), Developments in the MPI‐M Earth System Model version 1.2 (MPI‐ESM1.2) and Its Response to Increasing CO2, J. Adv. Model. Earth Syst.,11, 998-1038, doi:10.1029/2018MS001400,\nMueller, W.A. et al. (2018): A high‐resolution version of the Max Planck Institute Earth System Model MPI‐ESM1.2‐HR. J. Adv. Model. EarthSyst.,10,1383–1413, doi:10.1029/2017MS001217"
 SOURCE="MPI-ESM1.2-CR (2017): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none\nland: JSBACH3.20\nlandIce: none/prescribed\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none/prescribed"
 SOURCE_TYPE="AOGCM"
-CALENDAR=proleptic_gregorian
+#CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
 #MAPPING_TABLE=MPIESM1_CMIP6.txt
 #GRID_FILE_DIR=__GRID_FILE_DIR
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp
index e937563..463971e 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp
@@ -347,6 +347,16 @@ function run_agg
         err="${errtemp_agg}/err${period}"
         rm -f ${err}*
 
+        # Read input file
+        if [[ "${iniyeararr[$y]}" == "$period" ]]
+        then
+            echo "Reading from parent ${RAW_EXP_ID_PARENT[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_ID_PARENT[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        else
+            echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        fi
+
         # Perform aggregation for each submodel
         for submodel in $atmmod $ocemod $srfmod; do # $bgcmod; do
 
@@ -543,7 +553,10 @@ echo "-> Looping over ${#flist[@]} netCDF files!"
 for f in ${flist[@]}
 do
     echo "... editing $f"
-    ncatted -O -h -a further_info_url,global,d,, $f && echo "... successful!" || echo "... failed (ERROR)"
+    ncatted -O -h -a further_info_url,global,d,, $f && echo "... deleted 'further_info_url'." || echo "... failed (ERROR)"
+    [[ "$(basename $f)" == *"cen"* ]] && {
+        ncatted -O -h -a frequency,global,m,c,"cen" $f && echo "... changed frequency to 'cen'." || echo "... failed (ERROR)"
+    }
 done
 
 print "$(date +%Y-%m-%dT%H:%M:%S): NCO cleanup ended for $dr"
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR_cdocmorinfo b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR_cdocmorinfo
index 819123a..6b716ca 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR_cdocmorinfo
@@ -47,7 +47,7 @@ SOURCE_ID=MPI-ESM1-2-CR
 REFERENCES="Kapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022). Ocean response in transient simulations of the last deglaciation dominated by underlying ice-sheet reconstruction and method of meltwater distribution. Geophysical Research Letters, 49, e2021GL096767. https://doi.org/10.1029/2021GL096767\nMPI-ESM: Mauritsen, T. et al. (2019), Developments in the MPI‐M Earth System Model version 1.2 (MPI‐ESM1.2) and Its Response to Increasing CO2, J. Adv. Model. Earth Syst.,11, 998-1038, doi:10.1029/2018MS001400,\nMueller, W.A. et al. (2018): A high‐resolution version of the Max Planck Institute Earth System Model MPI‐ESM1.2‐HR. J. Adv. Model. EarthSyst.,10,1383–1413, doi:10.1029/2017MS001217"
 SOURCE="MPI-ESM1.2-CR (2017): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none\nland: JSBACH3.20\nlandIce: none/prescribed\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none/prescribed"
 SOURCE_TYPE="AOGCM"
-CALENDAR=proleptic_gregorian
+#CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
 #MAPPING_TABLE=MPIESM1_CMIP6.txt
 #GRID_FILE_DIR=__GRID_FILE_DIR
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp
index a6024d6..45c83a3 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp
@@ -347,6 +347,16 @@ function run_agg
         err="${errtemp_agg}/err${period}"
         rm -f ${err}*
 
+        # Read input file
+        if [[ "${iniyeararr[$y]}" == "$period" ]]
+        then
+            echo "Reading from parent ${RAW_EXP_ID_PARENT[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_ID_PARENT[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        else
+            echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        fi
+
         # Perform aggregation for each submodel
         for submodel in $atmmod $ocemod $srfmod; do # $bgcmod; do
 
@@ -543,7 +553,10 @@ echo "-> Looping over ${#flist[@]} netCDF files!"
 for f in ${flist[@]}
 do
     echo "... editing $f"
-    ncatted -O -h -a further_info_url,global,d,, $f && echo "... successful!" || echo "... failed (ERROR)"
+    ncatted -O -h -a further_info_url,global,d,, $f && echo "... deleted 'further_info_url'." || echo "... failed (ERROR)"
+    [[ "$(basename $f)" == *"cen"* ]] && {
+        ncatted -O -h -a frequency,global,m,c,"cen" $f && echo "... changed frequency to 'cen'." || echo "... failed (ERROR)"
+    }
 done
 
 print "$(date +%Y-%m-%dT%H:%M:%S): NCO cleanup ended for $dr"
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR_cdocmorinfo b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR_cdocmorinfo
index 68d1516..f2980af 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR_cdocmorinfo
@@ -47,7 +47,7 @@ SOURCE_ID=MPI-ESM1-2-CR
 REFERENCES="Kapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022). Ocean response in transient simulations of the last deglaciation dominated by underlying ice-sheet reconstruction and method of meltwater distribution. Geophysical Research Letters, 49, e2021GL096767. https://doi.org/10.1029/2021GL096767\nMPI-ESM: Mauritsen, T. et al. (2019), Developments in the MPI‐M Earth System Model version 1.2 (MPI‐ESM1.2) and Its Response to Increasing CO2, J. Adv. Model. Earth Syst.,11, 998-1038, doi:10.1029/2018MS001400,\nMueller, W.A. et al. (2018): A high‐resolution version of the Max Planck Institute Earth System Model MPI‐ESM1.2‐HR. J. Adv. Model. EarthSyst.,10,1383–1413, doi:10.1029/2017MS001217"
 SOURCE="MPI-ESM1.2-CR (2017): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none\nland: JSBACH3.20\nlandIce: none/prescribed\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none/prescribed"
 SOURCE_TYPE="AOGCM"
-CALENDAR=proleptic_gregorian
+#CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
 #MAPPING_TABLE=MPIESM1_CMIP6.txt
 #GRID_FILE_DIR=__GRID_FILE_DIR
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp
index b982c08..ae31651 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp
@@ -347,6 +347,16 @@ function run_agg
         err="${errtemp_agg}/err${period}"
         rm -f ${err}*
 
+        # Read input file
+        if [[ "${iniyeararr[$y]}" == "$period" ]]
+        then
+            echo "Reading from parent ${RAW_EXP_ID_PARENT[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_ID_PARENT[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        else
+            echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc"
+            ifile=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
+        fi
+
         # Perform aggregation for each submodel
         for submodel in $atmmod $ocemod $srfmod; do # $bgcmod; do
 
@@ -543,7 +553,10 @@ echo "-> Looping over ${#flist[@]} netCDF files!"
 for f in ${flist[@]}
 do
     echo "... editing $f"
-    ncatted -O -h -a further_info_url,global,d,, $f && echo "... successful!" || echo "... failed (ERROR)"
+    ncatted -O -h -a further_info_url,global,d,, $f && echo "... deleted 'further_info_url'." || echo "... failed (ERROR)"
+    [[ "$(basename $f)" == *"cen"* ]] && {
+        ncatted -O -h -a frequency,global,m,c,"cen" $f && echo "... changed frequency to 'cen'." || echo "... failed (ERROR)"
+    }
 done
 
 print "$(date +%Y-%m-%dT%H:%M:%S): NCO cleanup ended for $dr"
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR_cdocmorinfo b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR_cdocmorinfo
index 4fb529b..05cdf3e 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR_cdocmorinfo
@@ -47,7 +47,7 @@ SOURCE_ID=MPI-ESM1-2-CR
 REFERENCES="Kapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022). Ocean response in transient simulations of the last deglaciation dominated by underlying ice-sheet reconstruction and method of meltwater distribution. Geophysical Research Letters, 49, e2021GL096767. https://doi.org/10.1029/2021GL096767\nMPI-ESM: Mauritsen, T. et al. (2019), Developments in the MPI‐M Earth System Model version 1.2 (MPI‐ESM1.2) and Its Response to Increasing CO2, J. Adv. Model. Earth Syst.,11, 998-1038, doi:10.1029/2018MS001400,\nMueller, W.A. et al. (2018): A high‐resolution version of the Max Planck Institute Earth System Model MPI‐ESM1.2‐HR. J. Adv. Model. EarthSyst.,10,1383–1413, doi:10.1029/2017MS001217"
 SOURCE="MPI-ESM1.2-CR (2017): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none\nland: JSBACH3.20\nlandIce: none/prescribed\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none/prescribed"
 SOURCE_TYPE="AOGCM"
-CALENDAR=proleptic_gregorian
+#CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
 #MAPPING_TABLE=MPIESM1_CMIP6.txt
 #GRID_FILE_DIR=__GRID_FILE_DIR
diff --git a/cmor_tables/PalMod2_CV.json b/cmor_tables/PalMod2_CV.json
index 808e683..ea44d5c 100644
--- a/cmor_tables/PalMod2_CV.json
+++ b/cmor_tables/PalMod2_CV.json
@@ -147,6 +147,19 @@
                 "source_id":"AWI-ESM-2-1-LR",
                 "source":"AWI-ESM 2.1 LR (2019): \naerosol: none\natmos: ECHAM6.3.04p1 (T63L47 native atmosphere T63 gaussian grid; 192 x 96 longitude/latitude; 47 levels; top level 80 km)\natmosChem: none\nland: JSBACH 3.20 with dynamic vegetation\nlandIce: none\nocean: FESOM 2 (unstructured grid in the horizontal with 126858 wet nodes; 48 levels; top grid cell 0-5 m)\nocnBgchem: none\nseaIce: FESOM 2\nsolidLand: none/prescribed"
             },
+            "AWI-ESM2-1-LR-PISM":{
+               "activity_participation":[
+                    "PalMod2"
+                ],
+                "cohort":[
+                    "Registered"
+                ],
+                "institution_id":[
+                    "AWI"
+                ],
+                "source_id":"AWI-ESM2-1-LR-PISM",
+                "source":"AWI-ESM2-1-LR_PISM: \n aerosol: none, prescribed MAC-v1\n atmos: ECHAM6.3.05 (spectral T63; 192 x 96 longitude/latitude; 47 levels; top level 0.01 hPa)\n atmosChem: none\n land: JSBACH3.20p1 (dynamic vegetation)\n landIce: PISM 1.2 (10 km x 10 km, 101 levels)\n ocean: FESOM2 (CORE2 mesh, 126,858 variable resolution; 47 levels; top grid cell 0-5 m)\n ocnBgchem: none\n seaIce: FESIM (elastic-viscous-plastic sea ice solver)\nsolidLand: none"
+            },
             "EC-Earth3":{
                 "activity_participation":[
                     "CMIP",
@@ -834,6 +847,69 @@
                     "none"
                 ]
             },
+            "transient-deglaciation-interactive":{
+                "activity_id":[
+                    "PalMod2"
+                ],
+                "additional_allowed_model_components":[
+                ],
+                "experiment":"synchronously coupled transient deglaciation with interactive ice sheets and solid earth",
+                "experiment_id":"transient_deglaciation_interactive",
+                "parent_activity_id":[
+                    "no parent"
+                ],
+                "parent_experiment_id":[
+                    "no parent"
+                ],
+                "required_model_components":[
+                    "AOGCM ISM SLM"
+                ],
+                "sub_experiment_id":[
+                    "none"
+                ]
+            },
+            "transient-deglaciation-interactive-mis3":{
+                "activity_id":[
+                    "PalMod2"
+                ],
+                "additional_allowed_model_components":[
+                ],
+                "experiment":"synchronously coupled transient simulation of marine isotope stage 3 (MIS3) with interactive ice sheets and solid earth",
+                "experiment_id":"transient-deglaciation-interactive-mis3",
+                "parent_activity_id":[
+                    "no parent"
+                ],
+                "parent_experiment_id":[
+                    "no parent"
+                ],
+                "required_model_components":[
+                    "AOGCM ISM SLM"
+                ],
+                "sub_experiment_id":[
+                    "none"
+                ]
+            },
+            "transient-deglaciation-interactive-mis3-asynchronous":{
+                "activity_id":[
+                    "PalMod2"
+                ],
+                "additional_allowed_model_components":[
+                ],
+                "experiment":"asynchronously coupled transient simulation of marine isotope stage 3 (MIS3) with interactive ice sheets and solid earth",
+                "experiment_id":"transient-deglaciation-interactive-mis3-asynchronous",
+                "parent_activity_id":[
+                    "no parent"
+                ],
+                "parent_experiment_id":[
+                    "no parent"
+                ],
+                "required_model_components":[
+                    "AOGCM ISM SLM"
+                ],
+                "sub_experiment_id":[
+                    "none"
+                ]
+            },
             "esm-past1000":{
                 "activity_id":[
                     "PMIP"
diff --git a/cmor_tables/PalMod2_IcenAnt.json b/cmor_tables/PalMod2_IcenAnt.json
new file mode 100644
index 0000000..d220ce1
--- /dev/null
+++ b/cmor_tables/PalMod2_IcenAnt.json
@@ -0,0 +1,144 @@
+{
+    "Header": {
+        "data_specs_version": "00.00.03",
+        "cmor_version": "3.6.1",
+        "table_id": "Table IcenAnt",
+        "realm": "landIce land",
+        "table_date": "14 February 2022",
+        "missing_value": "1e20",
+        "int_missing_value": "-999",
+        "product": "model-output",
+        "approx_interval": "36500.00000",
+        "generic_levels": "alevel",
+        "Conventions": "CF-1.7 PalMod-2.0"
+    },
+    "variable_entry": {
+        "orogIs": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "surface_altitude",
+            "units": "m",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Altitude",
+            "comment": "The surface called 'surface' means the lower boundary of the atmosphere. Altitude is the (geometric) height above the geoid, which is the reference geopotential surface. The geoid is similar to mean sea level.",
+            "dimensions": "longitude latitude time",
+            "out_name": "orogIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "tsIs": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "surface_temperature",
+            "units": "K",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Temperature",
+            "comment": "Temperature of the lower boundary of the atmosphere",
+            "dimensions": "longitude latitude time",
+            "out_name": "tsIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "acabfIs": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_surface_specific_mass_balance_flux",
+            "units": "kg m-2 s-1",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Mass Balance Flux",
+            "comment": "Specific mass balance means the net rate at which ice is added per unit area at the land ice surface. Computed as the total surface mass balance on the land ice portion of the grid cell divided by land ice area in the grid cell. A negative value means loss of ice",
+            "dimensions": "longitude latitude time",
+            "out_name": "acabfIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "libmassbffl": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_basal_specific_mass_balance_flux",
+            "units": "kg m-2 s-1",
+            "cell_methods": "area: time: mean where floating_ice_shelf (comment: mask=sftflf)",
+            "cell_measures": "area: areacellg",
+            "long_name": "Basal Specific Mass Balance Flux of Floating Ice Shelf",
+            "comment": "Specific mass balance means the net rate at which ice is added per unit area at the land ice base.  A negative value means loss of ice. Computed as the total basal mass balance on the floating land ice (floating ice shelf) portion of the grid cell divided by floating land ice (floating ice shelf) area in the grid cell. Cell_methods: area: mean where floating_ice_shelf",
+            "dimensions": "longitude latitude time",
+            "out_name": "libmassbffl",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "tlIs": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_temperature",
+            "units": "K",
+            "cell_methods": "area: time: mean where ice_sheet (Weighted Time Mean on Ice Sheet)",
+            "cell_measures": "area: areacellg",
+            "long_name": "Temperature of Land Ice",
+            "comment": "Land ice means glaciers, ice-caps and ice-sheets resting on bedrock and also includes ice-shelves.",
+            "dimensions": "longitude latitude alevel time",
+            "out_name": "tlIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "sftgifIt": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_area_fraction",
+            "units": "%",
+            "cell_methods": "area: mean where ice_sheet over all_area_types time: mean",
+            "cell_measures": "area: areacellg",
+            "long_name": "Type of Land Ice",
+            "comment": "Mask of land ice type (ice sheet, ice shelf, ocean, ice-free land).",
+            "dimensions": "longitude latitude typelice time",
+            "out_name": "sftgifIt",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "dtb": {
+            "frequency": "dec",
+            "modeling_realm": "land",
+            "standard_name": "bedrock_altitude",
+            "units": "m",
+            "cell_methods": "area: time: mean",
+            "cell_measures": "area: areacella",
+            "long_name": "Bedrock Altitude",
+            "comment": "Altitude is the (geometric) height above the geoid, which is the reference geopotential surface. The geoid is similar to mean sea level. Bedrock is the solid Earth surface beneath land ice, ocean water or soil.",
+            "dimensions": "longitude latitude time",
+            "out_name": "dtb",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        }
+    }
+}
+
diff --git a/cmor_tables/PalMod2_IcenGre.json b/cmor_tables/PalMod2_IcenGre.json
new file mode 100644
index 0000000..bc0755f
--- /dev/null
+++ b/cmor_tables/PalMod2_IcenGre.json
@@ -0,0 +1,144 @@
+{
+    "Header": {
+        "data_specs_version": "00.00.03",
+        "cmor_version": "3.6.1",
+        "table_id": "Table IcenGre",
+        "realm": "landIce land",
+        "table_date": "14 February 2022",
+        "missing_value": "1e20",
+        "int_missing_value": "-999",
+        "product": "model-output",
+        "approx_interval": "36500.00000",
+        "generic_levels": "alevel",
+        "Conventions": "CF-1.7 PalMod-2.0"
+    },
+    "variable_entry": {
+        "orogIs": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "surface_altitude",
+            "units": "m",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Altitude",
+            "comment": "The surface called 'surface' means the lower boundary of the atmosphere. Altitude is the (geometric) height above the geoid, which is the reference geopotential surface. The geoid is similar to mean sea level.",
+            "dimensions": "longitude latitude time",
+            "out_name": "orogIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "tsIs": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "surface_temperature",
+            "units": "K",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Temperature",
+            "comment": "Temperature of the lower boundary of the atmosphere",
+            "dimensions": "longitude latitude time",
+            "out_name": "tsIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "acabfIs": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_surface_specific_mass_balance_flux",
+            "units": "kg m-2 s-1",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Mass Balance Flux",
+            "comment": "Specific mass balance means the net rate at which ice is added per unit area at the land ice surface. Computed as the total surface mass balance on the land ice portion of the grid cell divided by land ice area in the grid cell. A negative value means loss of ice",
+            "dimensions": "longitude latitude time",
+            "out_name": "acabfIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "libmassbffl": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_basal_specific_mass_balance_flux",
+            "units": "kg m-2 s-1",
+            "cell_methods": "area: time: mean where floating_ice_shelf (comment: mask=sftflf)",
+            "cell_measures": "area: areacellg",
+            "long_name": "Basal Specific Mass Balance Flux of Floating Ice Shelf",
+            "comment": "Specific mass balance means the net rate at which ice is added per unit area at the land ice base.  A negative value means loss of ice. Computed as the total basal mass balance on the floating land ice (floating ice shelf) portion of the grid cell divided by floating land ice (floating ice shelf) area in the grid cell. Cell_methods: area: mean where floating_ice_shelf",
+            "dimensions": "longitude latitude time",
+            "out_name": "libmassbffl",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "tlIs": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_temperature",
+            "units": "K",
+            "cell_methods": "area: time: mean where ice_sheet (Weighted Time Mean on Ice Sheet)",
+            "cell_measures": "area: areacellg",
+            "long_name": "Temperature of Land Ice",
+            "comment": "Land ice means glaciers, ice-caps and ice-sheets resting on bedrock and also includes ice-shelves.",
+            "dimensions": "longitude latitude alevel time",
+            "out_name": "tlIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "sftgifIt": {
+            "frequency": "dec",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_area_fraction",
+            "units": "%",
+            "cell_methods": "area: mean where ice_sheet over all_area_types time: mean",
+            "cell_measures": "area: areacellg",
+            "long_name": "Type of Land Ice",
+            "comment": "Mask of land ice type (ice sheet, ice shelf, ocean, ice-free land).",
+            "dimensions": "longitude latitude typelice time",
+            "out_name": "sftgifIt",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "dtb": {
+            "frequency": "dec",
+            "modeling_realm": "land",
+            "standard_name": "bedrock_altitude",
+            "units": "m",
+            "cell_methods": "area: time: mean",
+            "cell_measures": "area: areacella",
+            "long_name": "Bedrock Altitude",
+            "comment": "Altitude is the (geometric) height above the geoid, which is the reference geopotential surface. The geoid is similar to mean sea level. Bedrock is the solid Earth surface beneath land ice, ocean water or soil.",
+            "dimensions": "longitude latitude time",
+            "out_name": "dtb",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        }
+    }
+}
+
diff --git a/cmor_tables/PalMod2_IyrAnt.json b/cmor_tables/PalMod2_IyrAnt.json
new file mode 100644
index 0000000..06b3226
--- /dev/null
+++ b/cmor_tables/PalMod2_IyrAnt.json
@@ -0,0 +1,144 @@
+{
+    "Header": {
+        "data_specs_version": "00.00.03",
+        "cmor_version": "3.6.1",
+        "table_id": "Table IyrAnt",
+        "realm": "landIce land",
+        "table_date": "14 February 2022",
+        "missing_value": "1e20",
+        "int_missing_value": "-999",
+        "product": "model-output",
+        "approx_interval": "365.00000",
+        "generic_levels": "alevel",
+        "Conventions": "CF-1.7 PalMod-2.0"
+    },
+    "variable_entry": {
+        "orogIs": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "surface_altitude",
+            "units": "m",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Altitude",
+            "comment": "The surface called 'surface' means the lower boundary of the atmosphere. Altitude is the (geometric) height above the geoid, which is the reference geopotential surface. The geoid is similar to mean sea level.",
+            "dimensions": "longitude latitude time",
+            "out_name": "orogIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "tsIs": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "surface_temperature",
+            "units": "K",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Temperature",
+            "comment": "Temperature of the lower boundary of the atmosphere",
+            "dimensions": "longitude latitude time",
+            "out_name": "tsIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "acabfIs": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_surface_specific_mass_balance_flux",
+            "units": "kg m-2 s-1",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Mass Balance Flux",
+            "comment": "Specific mass balance means the net rate at which ice is added per unit area at the land ice surface. Computed as the total surface mass balance on the land ice portion of the grid cell divided by land ice area in the grid cell. A negative value means loss of ice",
+            "dimensions": "longitude latitude time",
+            "out_name": "acabfIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "libmassbffl": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_basal_specific_mass_balance_flux",
+            "units": "kg m-2 s-1",
+            "cell_methods": "area: time: mean where floating_ice_shelf (comment: mask=sftflf)",
+            "cell_measures": "area: areacellg",
+            "long_name": "Basal Specific Mass Balance Flux of Floating Ice Shelf",
+            "comment": "Specific mass balance means the net rate at which ice is added per unit area at the land ice base.  A negative value means loss of ice. Computed as the total basal mass balance on the floating land ice (floating ice shelf) portion of the grid cell divided by floating land ice (floating ice shelf) area in the grid cell. Cell_methods: area: mean where floating_ice_shelf",
+            "dimensions": "longitude latitude time",
+            "out_name": "libmassbffl",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "tlIs": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_temperature",
+            "units": "K",
+            "cell_methods": "area: time: mean where ice_sheet (Weighted Time Mean on Ice Sheet)",
+            "cell_measures": "area: areacellg",
+            "long_name": "Temperature of Land Ice",
+            "comment": "Land ice means glaciers, ice-caps and ice-sheets resting on bedrock and also includes ice-shelves.",
+            "dimensions": "longitude latitude alevel time",
+            "out_name": "tlIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "sftgifIt": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_area_fraction",
+            "units": "%",
+            "cell_methods": "area: mean where ice_sheet over all_area_types time: mean",
+            "cell_measures": "area: areacellg",
+            "long_name": "Type of Land Ice",
+            "comment": "Mask of land ice type (ice sheet, ice shelf, ocean, ice-free land).",
+            "dimensions": "longitude latitude typelice time",
+            "out_name": "sftgifIt",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "dtb": {
+            "frequency": "yr",
+            "modeling_realm": "land",
+            "standard_name": "bedrock_altitude",
+            "units": "m",
+            "cell_methods": "area: time: mean",
+            "cell_measures": "area: areacella",
+            "long_name": "Bedrock Altitude",
+            "comment": "Altitude is the (geometric) height above the geoid, which is the reference geopotential surface. The geoid is similar to mean sea level. Bedrock is the solid Earth surface beneath land ice, ocean water or soil.",
+            "dimensions": "longitude latitude time",
+            "out_name": "dtb",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        }
+    }
+}
+
diff --git a/cmor_tables/PalMod2_IyrGre.json b/cmor_tables/PalMod2_IyrGre.json
new file mode 100644
index 0000000..f8c67ab
--- /dev/null
+++ b/cmor_tables/PalMod2_IyrGre.json
@@ -0,0 +1,144 @@
+{
+    "Header": {
+        "data_specs_version": "00.00.03",
+        "cmor_version": "3.6.1",
+        "table_id": "Table IyrGre",
+        "realm": "landIce land",
+        "table_date": "14 February 2022",
+        "missing_value": "1e20",
+        "int_missing_value": "-999",
+        "product": "model-output",
+        "approx_interval": "365.00000",
+        "generic_levels": "alevel",
+        "Conventions": "CF-1.7 PalMod-2.0"
+    },
+    "variable_entry": {
+        "orogIs": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "surface_altitude",
+            "units": "m",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Altitude",
+            "comment": "The surface called 'surface' means the lower boundary of the atmosphere. Altitude is the (geometric) height above the geoid, which is the reference geopotential surface. The geoid is similar to mean sea level.",
+            "dimensions": "longitude latitude time",
+            "out_name": "orogIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "tsIs": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "surface_temperature",
+            "units": "K",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Temperature",
+            "comment": "Temperature of the lower boundary of the atmosphere",
+            "dimensions": "longitude latitude time",
+            "out_name": "tsIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "acabfIs": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_surface_specific_mass_balance_flux",
+            "units": "kg m-2 s-1",
+            "cell_methods": "area: time: mean where ice_sheet",
+            "cell_measures": "area: areacellg",
+            "long_name": "Ice Sheet Surface Mass Balance Flux",
+            "comment": "Specific mass balance means the net rate at which ice is added per unit area at the land ice surface. Computed as the total surface mass balance on the land ice portion of the grid cell divided by land ice area in the grid cell. A negative value means loss of ice",
+            "dimensions": "longitude latitude time",
+            "out_name": "acabfIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "libmassbffl": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_basal_specific_mass_balance_flux",
+            "units": "kg m-2 s-1",
+            "cell_methods": "area: time: mean where floating_ice_shelf (comment: mask=sftflf)",
+            "cell_measures": "area: areacellg",
+            "long_name": "Basal Specific Mass Balance Flux of Floating Ice Shelf",
+            "comment": "Specific mass balance means the net rate at which ice is added per unit area at the land ice base.  A negative value means loss of ice. Computed as the total basal mass balance on the floating land ice (floating ice shelf) portion of the grid cell divided by floating land ice (floating ice shelf) area in the grid cell. Cell_methods: area: mean where floating_ice_shelf",
+            "dimensions": "longitude latitude time",
+            "out_name": "libmassbffl",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "tlIs": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_temperature",
+            "units": "K",
+            "cell_methods": "area: time: mean where ice_sheet (Weighted Time Mean on Ice Sheet)",
+            "cell_measures": "area: areacellg",
+            "long_name": "Temperature of Land Ice",
+            "comment": "Land ice means glaciers, ice-caps and ice-sheets resting on bedrock and also includes ice-shelves.",
+            "dimensions": "longitude latitude alevel time",
+            "out_name": "tlIs",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "sftgifIt": {
+            "frequency": "yr",
+            "modeling_realm": "landIce",
+            "standard_name": "land_ice_area_fraction",
+            "units": "%",
+            "cell_methods": "area: mean where ice_sheet over all_area_types time: mean",
+            "cell_measures": "area: areacellg",
+            "long_name": "Type of Land Ice",
+            "comment": "Mask of land ice type (ice sheet, ice shelf, ocean, ice-free land).",
+            "dimensions": "longitude latitude typelice time",
+            "out_name": "sftgifIt",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        },
+        "dtb": {
+            "frequency": "yr",
+            "modeling_realm": "land",
+            "standard_name": "bedrock_altitude",
+            "units": "m",
+            "cell_methods": "area: time: mean",
+            "cell_measures": "area: areacella",
+            "long_name": "Bedrock Altitude",
+            "comment": "Altitude is the (geometric) height above the geoid, which is the reference geopotential surface. The geoid is similar to mean sea level. Bedrock is the solid Earth surface beneath land ice, ocean water or soil.",
+            "dimensions": "longitude latitude time",
+            "out_name": "dtb",
+            "type": "real",
+            "positive": "",
+            "valid_min": "",
+            "valid_max": "",
+            "ok_min_mean_abs": "",
+            "ok_max_mean_abs": ""
+        }
+    }
+}
+
-- 
GitLab