Updates and fixes for MPI-M transient-deglaciation-interactive scripts.

5795b4e8 · Martin Schupfner · 51431a94 · 5795b4e8 · 5795b4e8 · 5795b4e8
Commit 5795b4e8 authored 3 months ago by Martin Schupfner
--- a/cmor/mpiesm/scripts/ncatted_CS.sh
+++ b/cmor/mpiesm/scripts/ncatted_CS.sh
@@ -6,14 +6,14 @@ module load nco

 # Attributes
 TITLE="MPI-ESM1-2-1 output prepared for PalMod2"
-REFERENCES="The paper ist currently in open discussion under the following citation:\nMikolajewicz, U., Kapsch, M.-L., Schannwell, C., Six, K. D., Ziemen, F. A., Bagge, M., Baudouin, J.-P., Erokhina, O., Gayler, V., Klemann, V., Meccia, V. L., Mouchet, A., and Riddick, T.: Deglaciation and abrupt events in a coupled comprehensive atmosphere–ocean–ice sheet–solid earth model, Clim. Past Discuss. [preprint], https://doi.org/10.5194/cp-2024-55, in review, 2024"
+REFERENCES="The paper is currently in open discussion under the following citation:\nMikolajewicz, U., Kapsch, M.-L., Schannwell, C., Six, K. D., Ziemen, F. A., Bagge, M., Baudouin, J.-P., Erokhina, O., Gayler, V., Klemann, V., Meccia, V. L., Mouchet, A., and Riddick, T.: Deglaciation and abrupt events in a coupled comprehensive atmosphere–ocean–ice sheet–solid earth model, Clim. Past Discuss. [preprint], https://doi.org/10.5194/cp-2024-55, in review, 2024"
 #NOMINAL_RESOLUTION="500 km" # echam6, jsbach, mpiom, hamocc
 #NOMINAL_RESOLUTION="km"    # vilma
 #NOMINAL_RESOLUTION="25 km" # mPISM
-SOURCE="MPI-ESM1.2.1-CR (2022): \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="MPI-ESM1.2.1-CR (2022): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none, prescribed\nland: JSBACH3.20, River Transport Model\nlandIce: none / mPISM 0.7 (10 km x 10 km (NH), 15 km x 15 km (SH), 121 levels)\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none, prescribed\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none / VILMA-1D"
 CONVENTIONS="CF-1.7 PalMod-2.0"
-CREATION_DATE="2024-09-20T18:00:00Z"
-HISTORY="2024-09-20T18:00:00Z ; CMOR rewrote data to be consistent with PalMod2, CF-1.7 PalMod-2.0 and CF standards."
+CREATION_DATE="2025-01-05T18:00:00Z"
+HISTORY="2025-01-05T18:00:00Z ; CMOR rewrote data to be consistent with PalMod2, CF-1.7 PalMod-2.0 and CF standards."
 PARENTEXPID="no parent"
 declare -A TABLE_INFO
 # 00.00.02
@@ -25,27 +25,27 @@ declare -A TABLE_INFO
 #TABLE_INFO[fx]="Creation Date:(14 February 2022) MD5:26bac4b397da7babac4e36561826e6b0"
 #TABLE_INFO[Ofx]="Creation Date:(14 February 2022) MD5:26bac4b397da7babac4e36561826e6b0"
 # 00.00.03
-dreqversion="00.00.03"
-TABLE_INFO[AERmon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[Amon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[centennial]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[dec]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[Emon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[EmonZ]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[fx]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[IcenAnt]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[IcenGre]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[IdecAnt]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[IdecGre]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[IyrAnt]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[IyrGre]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[LImon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[Lmon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[Odec]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[Ofx]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[Omon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[Oyr]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
-TABLE_INFO[SImon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#dreqversion="00.00.03"
+#TABLE_INFO[AERmon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[Amon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[centennial]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[dec]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[Emon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[EmonZ]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[fx]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[IcenAnt]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[IcenGre]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[IdecAnt]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[IdecGre]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[IyrAnt]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[IyrGre]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[LImon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[Lmon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[Odec]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[Ofx]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[Omon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[Oyr]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
+#TABLE_INFO[SImon]="Creation Date:(1 February 2023) MD5:437d9b7f4524e7483e04d57b09f4aaaa"
 # 00.00.04
 dreqversion="00.00.04"
 TABLE_INFO[AERmon]="Creation Date:(18 September 2024) MD5:028d52f4a33ef4da346678b341512dda"
@@ -71,8 +71,25 @@ TABLE_INFO[SImon]="Creation Date:(18 September 2024) MD5:028d52f4a33ef4da346678b
 TABLE_INFO[SLdec]="Creation Date:(18 September 2024) MD5:028d52f4a33ef4da346678b341512dda"
 #TABLE_INFO[]=
 declare -A VARIANT_INFO
-VARIANT_INFO["4"]="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. The simulation was forced with time-varying volcanic forcing (Schindlbeck-Belo et al., 2024). GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978)."
+# !!!!!!!
+# Beware: the variant_info does also vary with the forcing parameter and not only with the physics parameter!
+# !!!!!!!
+#p1f1 / pmo0007
 VARIANT_INFO["1"]="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. The simulation used climatological volcanic forcing representing PI. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978)."
+#p1f2 / pmo0012
+#VARIANT_INFO["1"]="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978). Equivalent to p1 but forced with time-varying volcanic forcing (Schindlbeck-Belo et al., 2024)."
+#p2f1 / pmo0009
+VARIANT_INFO["2"]="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. The simulation used climatological volcanic forcing representing PI. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978). Equivalent to p1 but with increased vertical mixing in the ocean. For details see Mikoljawicz et al. (2024)."
+#p3f1 / pmo0014
+VARIANT_INFO["3"]="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. The simulation used climatological volcanic forcing representing PI. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978). Equivalent to p1 but with decreased vertical mixing in the ocean. For details see Mikoljawicz et al. (2024)."
+#p4f2 / pmo0016
+VARIANT_INFO["4"]="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. Equivalent to p1 but with retuned atmosphere as well as modified sea-ice rigidity, ocean mixing, and basal sliding parameters. In addition, the simulation was forced with time-varying volcanic forcing (Schindlbeck-Belo et al., 2024). GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978)."
+#p5f2 / pmo0018
+VARIANT_INFO["5"]="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978). Equivalent to p4 but with the atmosphere tuning of p1. The simulation was forced with time-varying volcanic forcing (Schindlbeck-Belo et al., 2024)."
+#p6f2 / pmo0019
+VARIANT_INFO["6"]="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978). Equivalent to p4 but with brighter albedo parameters. The simulation was forced with time-varying volcanic forcing (Schindlbeck-Belo et al., 2024)."
+#p7f2 / pmo0020
+VARIANT_INFO["7"]="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978). Equivalent to p4 but with decreased basal sliding near the grounding line. The simulation was forced with time-varying volcanic forcing (Schindlbeck-Belo et al., 2024)."
 echo "------------------------"
 echo $TITLE
 echo $REFERENCES
@@ -84,7 +101,7 @@ echo $HISTORY
 echo ${TABLE_INFO[@]}
 echo "-------------------------"

-exp_id=transient-deglaciation-interactive_r1i1p4f2-CR
+exp_id=transient-deglaciation-interactive_r1i1p2f1-CR

 ifolder=/work/kd1292/ESGF_Buff/k204212/palmod_CS/${exp_id}/PalMod2/

@@ -108,6 +125,7 @@ echo $i $ifile
    #ncatted -O -h -a variant_info,global,m,c,"${VARIANT_INFO[$physics_index]}" $ifile || echo "ERROR $ifile"
    #ncatted -O -h -a tracking_id,global,m,c,"$uuid" $ifile || echo "ERROR $ifile"
    #ncatted -O -h -a parent_experiment_id,global,m,c,"${PARENTEXPID}" -a source,global,m,c,"$SOURCE" -a references,global,m,c,"$REFERENCES" -a variant_info,global,m,c,"${VARIANT_INFO[$physics_index]}" -a table_info,global,m,c,"${TABLE_INFO[$table]}" -a data_specs_version,global,m,c,${dreqversion} -a history,global,m,c,"$HISTORY" -a creation_date,global,m,c,"$CREATION_DATE" -a Conventions,global,m,c,"$CONVENTIONS" -a title,global,m,c,"$TITLE" -a tracking_id,global,m,c,"$uuid" $ifile || echo "ERROR $ifile"
+    #ncatted -O -h -a history,global,m,c,"$HISTORY" -a creation_date,global,m,c,"$CREATION_DATE" $ifile || echo "ERROR $ifile"

    # Adjust filename timestamp
    if [[ "$table" != *"fx"* ]]; then

--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f1-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f1-CR.cmor_ctl
@@ -12,7 +12,7 @@ n_subsim=4
 spinup=1000
 iniyeararr=( $((4000+spinup)) 2000 2000 2000 )
 finyeararr=( 8999 8999 8999 8999 )
-iniyearoffset=1000
+iniyearoffset=0
 iniyear=$(( 1 + $iniyearoffset + $spinup ))
 finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
@@ -27,7 +27,7 @@ inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10
 # Chunksize
 #  Note that there are prerequisites concerning the chunksize:
 #  - the defined decadal timeslizes may not overlap between two chunks!
-#  - decadal bounds have to be defined accordingly - "???01231" in this case 
+#  - decadal bounds have to be defined accordingly - "???01231" in this case
 #  - usually this is the case by default due to decadal model output!
 chunksize=100


--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f1-CR.runpp
@@ -43,7 +43,7 @@ alias ncatted=/sw/spack-levante/nco-5.0.6-3xkdth/bin/ncatted
 #########################################

 #Root directory of all scripts, outdata, logs etc
-EXP_ID=transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR
+EXP_ID=transient-deglaciation-interactive_r1i1p1f1-CR
 #EXP_DIR_TRUNK=/work/bk1192/k204212/palmod_CS/transient-deglaciation-interactive_r1i1p1f1-CR
 EXP_DIR_TRUNK=/work/kd1292/ESGF_Buff/k204212/palmod_CS/transient-deglaciation-interactive_r1i1p1f1-CR

@@ -135,7 +135,7 @@ ca+=([vilma]="${SCRIPT_ROOT}/cdocmorinfo_vilma_CR")
 experiment=transient-deglaciation-interactive
 member=r1i1p1f1
 mip=PalMod2
-vd=v20240920
+vd=v20250105
 #vd=v$(date '+%Y%m%d') #CMOR version directory name, should be a date in format YYYYMMDD

 #Where to store errors
@@ -388,7 +388,7 @@ function run_agg
        fi

        # Perform aggregation for each submodel
-        for submodel in $slmod; do #$icemod $slmod $atmmod $ocemod $srfmod; do
+        for submodel in $icemod $slmod $atmmod $ocemod $srfmod; do

            if [[ "$submodel" == "$icemod" ]]; then
                period=$(( inv_chunk_year - 10 ))
@@ -416,7 +416,7 @@ function run_agg
            sdir=${SDIR[$y]}/outdata/$submodel
            rawsdir=${RAWSDIR[$y]}/outdata/$submodel
            rawsdir_pism_ant=$(echo $rawsdir | rev | cut -d "/" -f 3- | rev)/PISM_ANT
-            rawsdir_pism_nh=$(echo $rawsdir | rev | cut -d "/" -f 3- | rev)/PISM_NH            
+            rawsdir_pism_nh=$(echo $rawsdir | rev | cut -d "/" -f 3- | rev)/PISM_NH
            echo "... writing output to $sdir/out_aggr"
            echo "... writing log to $errtemp_agg"

@@ -478,7 +478,7 @@ function run_diag
        rm -f ${err}*

        # Perform diagnostic for each submodel
-        for submodel in $slmod; do #$icemod $slmod $atmmod $ocemod $srfmod; do
+        for submodel in $icemod $slmod $atmmod $ocemod $srfmod; do

            if [[ "$submodel" == "$icemod" ]]; then
                period=$(( inv_chunk_year - 10 ))
@@ -585,7 +585,7 @@ function run_cmor
        rm -f ${err}*

        # Perform CMOR rewrite for each submodel
-        for submodel in $slmod; do #$icemod $slmod $atmmod $ocemod $srfmod; do
+        for submodel in $icemod $slmod $atmmod $ocemod $srfmod; do

            if [[ "$submodel" == "$icemod" ]]; then
                period=$(( inv_chunk_year - 10 ))

--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f1-CR_cdocmorinfo
@@ -19,7 +19,7 @@ FORCING_INDEX=1
 VARIANT_INFO="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. The simulation used climatological volcanic forcing representing PI. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978)."
 COMMENT="Earth viscosity configuration for the solid earth model VILMA. The vertical viscosity structure is piecewise constant for the following specified layers of the earth: (1) 0 to 3.480000e+06 m from earth core: 1.000000e+19 Pa s. (2) 3.480000e+06 to 5.701000e+06 m from earth core: 1.000000e+22 Pa s. (3) 5.701000e+06 to 6.291000e+06 m from earth core: 4.000000e+20 Pa s. (4) 6.291000e+06 to 6.371000e+06 m from earth core: 1.000000e+30 Pa s."
 ACTIVITY_ID=PalMod2
-#MIP_ERA=PalMod2                   #ONLY CMIP6  
+#MIP_ERA=PalMod2                   #ONLY CMIP6
 PROJECT_ID=PalMod2
 REQUIRED_TIME_UNITS="days since 1-1-1 00:00:00"
 #FURTHER_INFO_URL="https://furtherinfo.es-doc.org/<mip_era>.<institution_id>.<source_id>.<experiment_id>.<sub_experiment_id>.<variant_label>
@@ -28,7 +28,7 @@ REQUIRED_TIME_UNITS="days since 1-1-1 00:00:00"
 #  else use BRANCH_DATES to give the dates in YYYYMMDD,YYYYMMDD format
 #BRANCH_TIME_IN_PARENT=
 #BRANCH_TIME_IN_CHILD=
-#BRANCH_DATES=19000101,18500101 
+#BRANCH_DATES=19000101,18500101
 #BRANCH_METHOD="standard"
 #PARENT_MIP_ERA=none              #ONLY CMIP6
 #PARENT_ACTIVITY_ID=none
@@ -45,7 +45,7 @@ SUB_EXPERIMENT=none
 #
 SOURCE_ID=MPI-ESM1-2-1-CR
 #MODEL_ID=MPI-ESM1-2-1-CR
-REFERENCES="The paper ist currently in open discussion under the following citation:\nMikolajewicz, U., Kapsch, M.-L., Schannwell, C., Six, K. D., Ziemen, F. A., Bagge, M., Baudouin, J.-P., Erokhina, O., Gayler, V., Klemann, V., Meccia, V. L., Mouchet, A., and Riddick, T.: Deglaciation and abrupt events in a coupled comprehensive atmosphere–ocean–ice sheet–solid earth model, Clim. Past Discuss. [preprint], https://doi.org/10.5194/cp-2024-55, in review, 2024"
+REFERENCES="The paper is currently in open discussion under the following citation:\nMikolajewicz, U., Kapsch, M.-L., Schannwell, C., Six, K. D., Ziemen, F. A., Bagge, M., Baudouin, J.-P., Erokhina, O., Gayler, V., Klemann, V., Meccia, V. L., Mouchet, A., and Riddick, T.: Deglaciation and abrupt events in a coupled comprehensive atmosphere–ocean–ice sheet–solid earth model, Clim. Past Discuss. [preprint], https://doi.org/10.5194/cp-2024-55, in review, 2024"
 SOURCE="MPI-ESM1.2.1-CR (2022): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none, prescribed\nland: JSBACH3.20, River Transport Model\nlandIce: none / mPISM 0.7 (10 km x 10 km (NH), 15 km x 15 km (SH), 121 levels)\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none, prescribed\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none / VILMA-1D"
 SOURCE_TYPE="AOGCM ISM SLM"
 #CALENDAR=proleptic_gregorian

--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f2-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f2-CR.cmor_ctl
+#! /bin/sh -e
+#SBATCH --account=bm0021
+#SBATCH --qos=esgf
+#SBATCH --time=48:00:00
+#SBATCH --partition=shared
+#SBATCH --mail-type=FAIL
+#SBATCH --output=%x_%j.log
+
+# Number of sub-simulations
+n_subsim=4
+#initial and final year of the experiment
+spinup=1000
+iniyeararr=( $((4000+spinup)) 2000 2000 2000 )
+finyeararr=( 8999 8999 8999 8949 )
+iniyearoffset=0
+iniyear=$(( 1 + $iniyearoffset + $spinup ))
+finyear=$(( 25950 + $iniyearoffset - $spinup ))
+#offset compared to startyear
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) )
+# Time axis (26000-51 BP)
+##     exp time       total yrs   new exp time   offset
+#a - 26000-21001 BP - 5000 yrs -     1- 5000 -   - 3999
+#b - 21000-14001 BP - 7000 yrs -  5001-12000 -     3001
+#c - 14000- 7001 BP - 7000 yrs - 12001-19000 -    10001
+#d -  7000-   51 BP - 7000 yrs - 19001-25950 -    17001
+
+# Chunksize
+#  Note that there are prerequisites concerning the chunksize:
+#  - the defined decadal timeslizes may not overlap between two chunks!
+#  - decadal bounds have to be defined accordingly - "???01231" in this case
+#  - usually this is the case by default due to decadal model output!
+chunksize=100
+
+# Create list
+for i in $(seq 0 $(( n_subsim - 1 )) )
+do
+  for yyyy in $(seq ${iniyeararr[$i]} ${finyeararr[$i]} )
+  do
+    paramarr+=( "$i $yyyy" )
+  done
+done
+
+rm -f cmor_ctl_list_interactive_synch_r1i1p1f2
+in=0
+for i in $(seq $((iniyear-iniyear)) $chunksize $((finyear-iniyear)))
+do
+  [[ $i -eq $((iniyear-iniyear)) ]] && continue
+  li=$(echo ${paramarr[$in]} | cut -d " " -f1)
+  ly=$(echo ${paramarr[$in]} | cut -d " " -f2)
+  ui=$(echo ${paramarr[$((i-1))]} | cut -d " " -f1)
+  uy=$(echo ${paramarr[$((i-1))]} | cut -d " " -f2)
+  in=$i
+  echo "$li $ui $ly $uy" >> cmor_ctl_list_interactive_synch_r1i1p1f2
+done
+[[ $i -ne $((finyear-iniyear)) ]] && {
+  li=$(echo ${paramarr[$in]} | cut -d " " -f1)
+  ly=$(echo ${paramarr[$in]} | cut -d " " -f2)
+  ui=$(echo ${paramarr[$((finyear-iniyear))]} | cut -d " " -f1)
+  uy=$(echo ${paramarr[$((finyear-iniyear))]} | cut -d " " -f2)
+  echo "$li $ui $ly $uy" >> cmor_ctl_list_interactive_synch_r1i1p1f2
+}
+
+make -k -j 19 -f  transient-deglaciation-interactive_r1i1p1f2-CR.cmor_ctl.mk "$@"
+
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f2-CR.cmor_ctl.mk
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f2-CR.cmor_ctl.mk
+EXP_ID = transient-deglaciation-interactive_r1i1p1f2-CR
+
+ACCOUNT = bm0021
+PARTITION = compute
+
+SUFFIX =
+ARCHIVE_SUFFIX =
+
+TARGET_STEP=cmor
+NUMCHUNKS=$(shell cat cmor_ctl_list_interactive_synch_r1i1p1f2 | wc -l )
+
+AGGR_TIME = 00:30:00
+DIAGS_TIME = 01:00:00
+CMOR_TIME = 04:00:00
+
+CMOR_FLAGS =
+
+SBATCHFLAGS=-W
+
+MESSAGE = $$(date +%Y-%m-%dT%H:%M:%S): CMIP6 diagnostics and CMOR rewriting ($(SUFFIX:_%=%))
+
+.PHONY: all
+.PRECIOUS: targets/$(EXP_ID).aggr$(SUFFIX).% targets/$(EXP_ID).diags$(SUFFIX).% targets/$(EXP_ID).cmor$(SUFFIX).%
+
+
+
+all:
+	@echo "$(MESSAGE) started"
+	mkdir -vp targets
+	$(MAKE) $(MFLAGS) -f $(MAKEFILE_LIST) $(addprefix targets/$(EXP_ID).$(TARGET_STEP)$(SUFFIX).,$(shell seq 1 $(NUMCHUNKS)))
+	@echo "$(MESSAGE) finished"
+
+targets/$(EXP_ID).aggr$(SUFFIX).%:
+	sbatch $(SBATCHFLAGS) --time=$(AGGR_TIME) --job-name=$(EXP_ID)_runpp_aggr$(SUFFIX) --output=%x_%j.log --comment=$* $(EXP_ID).runpp -s '$(SUFFIX)' -S '$(ARCHIVE_SUFFIX)' -A $(shell sed -n '$(*)'p < cmor_ctl_list_interactive_synch_r1i1p1f2)
+	@touch $@
+
+targets/$(EXP_ID).diags$(SUFFIX).%: targets/$(EXP_ID).aggr$(SUFFIX).%
+	sbatch $(SBATCHFLAGS) --time=$(DIAGS_TIME) --job-name=$(EXP_ID)_runpp_diag$(SUFFIX) --output=%x_%j.log --comment=$* $(EXP_ID).runpp -s '$(SUFFIX)' -S '$(ARCHIVE_SUFFIX)' -d $(shell sed -n '$(*)'p < cmor_ctl_list_interactive_synch_r1i1p1f2)
+	@touch $@
+
+targets/$(EXP_ID).cmor$(SUFFIX).%: targets/$(EXP_ID).diags$(SUFFIX).%
+	sbatch $(SBATCHFLAGS) --time=$(CMOR_TIME) --job-name=$(EXP_ID)_runpp_cmor$(SUFFIX) --output=%x_%j.log --comment=$* $(EXP_ID).runpp -s '$(SUFFIX)' -S '$(ARCHIVE_SUFFIX)' -c $(CMOR_FLAGS) $(shell sed -n '$(*)'p < cmor_ctl_list_interactive_synch_r1i1p1f2)
+	@touch $@
+
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f2-CR.runpp
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f2-CR_cdocmorinfo
+#
+# experiment info:
+#
+_CONTROLLED_VOCABULARY_FILE=PalMod2_CV.json
+_FORMULA_VAR_FILE=PalMod2_formula_terms.json
+_AXIS_ENTRY_FILE=PalMod2_coordinate.json
+CONVENTIONS="CF-1.7 PalMod-2.0"
+TITLE="MPI-ESM1-2-1 output prepared for PalMod2"
+_history_template="%s ; CMOR rewrote data to be consistent with <activity_id>, <Conventions> and CF standards."
+
+EXPERIMENT_ID=transient-deglaciation-interactive
+EXPERIMENT="synchronously coupled transient deglaciation with interactive ice sheets and solid earth"
+#Specify VARIANT_LABEL explicitly or realization_index, physics_index, initialization_index and forcing_index:
+VARIANT_LABEL=r1i1p1f2
+REALIZATION_INDEX=1
+INITIALIZATION_INDEX=1
+PHYSICS_INDEX=1
+FORCING_INDEX=2
+VARIANT_INFO="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978). Equivalent to p1 but forced with time-varying volcanic forcing (Schindlbeck-Belo et al., 2024)."
+COMMENT="Earth viscosity configuration for the solid earth model VILMA. The vertical viscosity structure is piecewise constant for the following specified layers of the earth: (1) 0 to 3.480000e+06 m from earth core: 1.000000e+19 Pa s. (2) 3.480000e+06 to 5.701000e+06 m from earth core: 1.000000e+22 Pa s. (3) 5.701000e+06 to 6.291000e+06 m from earth core: 4.000000e+20 Pa s. (4) 6.291000e+06 to 6.371000e+06 m from earth core: 1.000000e+30 Pa s."
+ACTIVITY_ID=PalMod2
+#MIP_ERA=PalMod2                   #ONLY CMIP6
+PROJECT_ID=PalMod2
+REQUIRED_TIME_UNITS="days since 1-1-1 00:00:00"
+#FURTHER_INFO_URL="https://furtherinfo.es-doc.org/<mip_era>.<institution_id>.<source_id>.<experiment_id>.<sub_experiment_id>.<variant_label>
+#Branch Info: If needed use the following 2 attributes to define the branch time
+#  in Parent_Time_Units and Required_Time_Units (=child)
+#  else use BRANCH_DATES to give the dates in YYYYMMDD,YYYYMMDD format
+#BRANCH_TIME_IN_PARENT=
+#BRANCH_TIME_IN_CHILD=
+#BRANCH_DATES=19000101,18500101
+#BRANCH_METHOD="standard"
+#PARENT_MIP_ERA=none              #ONLY CMIP6
+#PARENT_ACTIVITY_ID=none
+#PARENT_EXPERIMENT=none
+#cmip6-spinup-HR
+#PARENT_EXPERIMENT_ID=none
+#PARENT_SOURCE_ID=none
+#PARENT_TIME_UNITS="days since 1850-1-1 00:00:00"
+#PARENT_VARIANT_LABEL=none
+SUB_EXPERIMENT_ID=none
+SUB_EXPERIMENT=none
+#
+#model info:
+#
+SOURCE_ID=MPI-ESM1-2-1-CR
+#MODEL_ID=MPI-ESM1-2-1-CR
+REFERENCES="The paper is currently in open discussion under the following citation:\nMikolajewicz, U., Kapsch, M.-L., Schannwell, C., Six, K. D., Ziemen, F. A., Bagge, M., Baudouin, J.-P., Erokhina, O., Gayler, V., Klemann, V., Meccia, V. L., Mouchet, A., and Riddick, T.: Deglaciation and abrupt events in a coupled comprehensive atmosphere–ocean–ice sheet–solid earth model, Clim. Past Discuss. [preprint], https://doi.org/10.5194/cp-2024-55, in review, 2024"
+SOURCE="MPI-ESM1.2.1-CR (2022): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none, prescribed\nland: JSBACH3.20, River Transport Model\nlandIce: none / mPISM 0.7 (10 km x 10 km (NH), 15 km x 15 km (SH), 121 levels)\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none, prescribed\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none / VILMA-1D"
+SOURCE_TYPE="AOGCM ISM SLM"
+#CALENDAR=proleptic_gregorian
+#MAPPING_TABLE_DIR=cmip6_mapping_tables
+#MAPPING_TABLE=MPIESM1_CMIP6.txt
+#GRID_FILE_DIR=__GRID_FILE_DIR
+#GRID_FILE=__GRID_FILE
+GRID=gn
+GRID_LABEL=gn
+#NOMINAL_RESOLUTION read from external file
+#For MPIOM/HAMOCC CR / GR30
+#NOMINAL_RESOLUTION="500 km"
+#For ECHAM6/JSBACH CR / T31
+#NOMINAL_RESOLUTION="500 km"
+CHAR_AXIS_TYPELICE="ice_free_land","ice_sheet","floating_ice","ice_free_sea"
+CHAR_AXIS_VEGTYPE=bare_land,glacier,tropical_evergreen_trees,tropical_deciduous_trees,extra-tropical_evergreen_trees,extra-tropical_deciduous_trees,raingreen_shrubs,deciduous_shrubs,C3_grass,C4_grass,C3_pasture,C4_pasture,C3_crops,C4_crops
+CHAR_AXIS_SOILPOOLS=acid-soluble_below_ground_leaf_litter,acid-soluble_below_ground_woody_litter,water-soluble_below_ground_leaf_litter,water-soluble_below_ground_woody_litter,ethanol-soluble_below_ground_leaf_litter,ethanol-soluble_below_ground_woody_litter,non-soluble_below_ground_leaf_litter,non-soluble_below_ground_woody_litter,humus_from_leaf_litter,humus_from_woody_litter
+CHAR_AXIS_LANDUSE=primary_and_secondary_land,crops,pastures,urban
+CHAR_AXIS_BASIN=global_ocean,atlantic_arctic_ocean,indian_pacific_ocean
+CHAR_AXIS_OLINE=barents_opening,bering_strait,canadian_archipelago,denmark_strait,drake_passage,english_channel,pacific_equatorial_undercurrent,faroe_scotland_channel,florida_bahamas_strait,fram_strait,iceland_faroe_channel,indonesian_throughflow,mozambique_channel,taiwan_luzon_straits,windward_passage
+CHAR_AXIS_SILINE=barents_opening,bering_strait,canadian_archipelago,fram_strait
+#spectband for the 3 aeropt variables
+CHAR_AXIS_SPECTBAND_aeroptbnd=1818182
+CHAR_AXIS_SPECTBAND_aeroptbnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_aeroptbnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_aerssabnd=1818182
+CHAR_AXIS_SPECTBAND_aerssabnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_aerssabnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_aerasymbnd=1818182
+CHAR_AXIS_SPECTBAND_aerasymbnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_aerasymbnd_UNITS="m-1"
+# 2d-spectband for the albedos
+CHAR_AXIS_SPECTBAND_albdirbnd=3225689,766690
+CHAR_AXIS_SPECTBAND_albdirbnd_BOUNDS=5000000,1451379,1451379,82001
+CHAR_AXIS_SPECTBAND_albdirbnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_albdiffbnd=3225689,766690
+CHAR_AXIS_SPECTBAND_albdiffbnd_BOUNDS=5000000,1451379,1451379,82001
+CHAR_AXIS_SPECTBAND_albdiffbnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_solbnd=1818182
+CHAR_AXIS_SPECTBAND_solbnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_solbnd_UNITS="m-1"
+#T_AXIS="cmip" causes the time axis to be rewritten according to the requested time stamps
+T_AXIS="cmip"
+#OUTPUT_MODE="r"
+OUTPUT_MODE="a"
+MAX_SIZE=0
+SAVE_CHUNK="n"
+#For files without time axis, if one is requested:
+#firsttimeval=1 (1850-01-02 00:00:00 wenn days since 1850-01-01 00:00:00)
+FIRSTTIMEVAL=1
+#
+# institution and contact info:
+#
+# none of these can be set in the command line
+#
+DEFLATE_LEVEL=1
+CONTACT=palmod2-mpi-esm@dkrz.de
+#
+# MPI-M
+#
+INSTITUTION_ID=MPI-M
+INSTITUTION="Max Planck Institute for Meteorology, Hamburg 20146, Germany"
+LICENSE="PalMod2 model data produced by MPI-M is licensed under a Creative Commons Attribution ShareAlike 4.0 International License (https://creativecommons.org/licenses). The data producers and data providers make no warranty, either express or implied, including, but not limited to, warranties of merchantability and fitness for a particular purpose. All liabilities arising from the supply of the information (including any liability arising in negligence) are excluded to the fullest extent permitted by law."
+#MIP_TABLE_DIR="/pool/data/CMIP6/cmip6-cmor-tables/Tables/"
+MIP_TABLE_DIR="/work/bm0021/PalMod2/cmor_tables/"
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p2f1-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p2f1-CR.cmor_ctl
+#! /bin/sh -e
+#SBATCH --account=bm0021
+#SBATCH --qos=esgf
+#SBATCH --time=48:00:00
+#SBATCH --partition=shared
+#SBATCH --mail-type=FAIL
+#SBATCH --output=%x_%j.log
+
+# Number of sub-simulations
+n_subsim=4
+#initial and final year of the experiment
+spinup=1000
+iniyeararr=( $((4000+spinup)) 2000 2000 2000 )
+finyeararr=( 8999 8999 8999 8999 )
+iniyearoffset=0
+iniyear=$(( 1 + $iniyearoffset + $spinup ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
+#offset compared to startyear
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) )
+# Time axis (26000-1 BP)
+##     exp time       total yrs   new exp time   offset
+#a - 26000-21001 BP - 5000 yrs -     1- 5000 -   - 3999
+#b - 21000-14001 BP - 7000 yrs -  5001-12000 -     3001
+#c - 14000- 7001 BP - 7000 yrs - 12001-19000 -    10001
+#d -  7000-    1 BP - 7000 yrs - 19001-26000 -    17001
+
+# Chunksize
+#  Note that there are prerequisites concerning the chunksize:
+#  - the defined decadal timeslizes may not overlap between two chunks!
+#  - decadal bounds have to be defined accordingly - "???01231" in this case
+#  - usually this is the case by default due to decadal model output!
+chunksize=100
+
+# Create list
+for i in $(seq 0 $(( n_subsim - 1 )) )
+do
+  for yyyy in $(seq ${iniyeararr[$i]} ${finyeararr[$i]} )
+  do
+    paramarr+=( "$i $yyyy" )
+  done
+done
+
+rm -f cmor_ctl_list_interactive_synch_r1i1p2f1
+in=0
+for i in $(seq $((iniyear-iniyear)) $chunksize $((finyear-iniyear)))
+do
+  [[ $i -eq $((iniyear-iniyear)) ]] && continue
+  li=$(echo ${paramarr[$in]} | cut -d " " -f1)
+  ly=$(echo ${paramarr[$in]} | cut -d " " -f2)
+  ui=$(echo ${paramarr[$((i-1))]} | cut -d " " -f1)
+  uy=$(echo ${paramarr[$((i-1))]} | cut -d " " -f2)
+  in=$i
+  echo "$li $ui $ly $uy" >> cmor_ctl_list_interactive_synch_r1i1p2f1
+done
+[[ $i -ne $((finyear-iniyear)) ]] && {
+  li=$(echo ${paramarr[$in]} | cut -d " " -f1)
+  ly=$(echo ${paramarr[$in]} | cut -d " " -f2)
+  ui=$(echo ${paramarr[$((finyear-iniyear))]} | cut -d " " -f1)
+  uy=$(echo ${paramarr[$((finyear-iniyear))]} | cut -d " " -f2)
+  echo "$li $ui $ly $uy" >> cmor_ctl_list_interactive_synch_r1i1p2f1
+}
+
+make -k -j 19 -f  transient-deglaciation-interactive_r1i1p2f1-CR.cmor_ctl.mk "$@"
+
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p2f1-CR.cmor_ctl.mk
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p2f1-CR.cmor_ctl.mk
+EXP_ID = transient-deglaciation-interactive_r1i1p2f1-CR
+
+ACCOUNT = bm0021
+PARTITION = compute
+
+SUFFIX =
+ARCHIVE_SUFFIX =
+
+TARGET_STEP=cmor
+NUMCHUNKS=$(shell cat cmor_ctl_list_interactive_synch_r1i1p2f1 | wc -l )
+
+AGGR_TIME = 00:30:00
+DIAGS_TIME = 01:00:00
+CMOR_TIME = 04:00:00
+
+CMOR_FLAGS =
+
+SBATCHFLAGS=-W
+
+MESSAGE = $$(date +%Y-%m-%dT%H:%M:%S): CMIP6 diagnostics and CMOR rewriting ($(SUFFIX:_%=%))
+
+.PHONY: all
+.PRECIOUS: targets/$(EXP_ID).aggr$(SUFFIX).% targets/$(EXP_ID).diags$(SUFFIX).% targets/$(EXP_ID).cmor$(SUFFIX).%
+
+
+
+all:
+	@echo "$(MESSAGE) started"
+	mkdir -vp targets
+	$(MAKE) $(MFLAGS) -f $(MAKEFILE_LIST) $(addprefix targets/$(EXP_ID).$(TARGET_STEP)$(SUFFIX).,$(shell seq 1 $(NUMCHUNKS)))
+	@echo "$(MESSAGE) finished"
+
+targets/$(EXP_ID).aggr$(SUFFIX).%:
+	sbatch $(SBATCHFLAGS) --time=$(AGGR_TIME) --job-name=$(EXP_ID)_runpp_aggr$(SUFFIX) --output=%x_%j.log --comment=$* $(EXP_ID).runpp -s '$(SUFFIX)' -S '$(ARCHIVE_SUFFIX)' -A $(shell sed -n '$(*)'p < cmor_ctl_list_interactive_synch_r1i1p2f1)
+	@touch $@
+
+targets/$(EXP_ID).diags$(SUFFIX).%: targets/$(EXP_ID).aggr$(SUFFIX).%
+	sbatch $(SBATCHFLAGS) --time=$(DIAGS_TIME) --job-name=$(EXP_ID)_runpp_diag$(SUFFIX) --output=%x_%j.log --comment=$* $(EXP_ID).runpp -s '$(SUFFIX)' -S '$(ARCHIVE_SUFFIX)' -d $(shell sed -n '$(*)'p < cmor_ctl_list_interactive_synch_r1i1p2f1)
+	@touch $@
+
+targets/$(EXP_ID).cmor$(SUFFIX).%: targets/$(EXP_ID).diags$(SUFFIX).%
+	sbatch $(SBATCHFLAGS) --time=$(CMOR_TIME) --job-name=$(EXP_ID)_runpp_cmor$(SUFFIX) --output=%x_%j.log --comment=$* $(EXP_ID).runpp -s '$(SUFFIX)' -S '$(ARCHIVE_SUFFIX)' -c $(CMOR_FLAGS) $(shell sed -n '$(*)'p < cmor_ctl_list_interactive_synch_r1i1p2f1)
+	@touch $@
+
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p2f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p2f1-CR.runpp
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p2f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p2f1-CR_cdocmorinfo
+#
+# experiment info:
+#
+_CONTROLLED_VOCABULARY_FILE=PalMod2_CV.json
+_FORMULA_VAR_FILE=PalMod2_formula_terms.json
+_AXIS_ENTRY_FILE=PalMod2_coordinate.json
+CONVENTIONS="CF-1.7 PalMod-2.0"
+TITLE="MPI-ESM1-2-1 output prepared for PalMod2"
+_history_template="%s ; CMOR rewrote data to be consistent with <activity_id>, <Conventions> and CF standards."
+
+EXPERIMENT_ID=transient-deglaciation-interactive
+EXPERIMENT="synchronously coupled transient deglaciation with interactive ice sheets and solid earth"
+#Specify VARIANT_LABEL explicitly or realization_index, physics_index, initialization_index and forcing_index:
+VARIANT_LABEL=r1i1p2f1
+REALIZATION_INDEX=1
+INITIALIZATION_INDEX=1
+PHYSICS_INDEX=2
+FORCING_INDEX=1
+VARIANT_INFO="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. The simulation used climatological volcanic forcing representing PI. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978). Equivalent to p1 but with increased vertical mixing in the ocean. For details see Mikoljawicz et al. (2024)." 
+COMMENT="Earth viscosity configuration for the solid earth model VILMA. The vertical viscosity structure is piecewise constant for the following specified layers of the earth: (1) 0 to 3.480000e+06 m from earth core: 1.000000e+19 Pa s. (2) 3.480000e+06 to 5.701000e+06 m from earth core: 1.000000e+22 Pa s. (3) 5.701000e+06 to 6.291000e+06 m from earth core: 4.000000e+20 Pa s. (4) 6.291000e+06 to 6.371000e+06 m from earth core: 1.000000e+30 Pa s."
+ACTIVITY_ID=PalMod2
+#MIP_ERA=PalMod2                   #ONLY CMIP6
+PROJECT_ID=PalMod2
+REQUIRED_TIME_UNITS="days since 1-1-1 00:00:00"
+#FURTHER_INFO_URL="https://furtherinfo.es-doc.org/<mip_era>.<institution_id>.<source_id>.<experiment_id>.<sub_experiment_id>.<variant_label>
+#Branch Info: If needed use the following 2 attributes to define the branch time
+#  in Parent_Time_Units and Required_Time_Units (=child)
+#  else use BRANCH_DATES to give the dates in YYYYMMDD,YYYYMMDD format
+#BRANCH_TIME_IN_PARENT=
+#BRANCH_TIME_IN_CHILD=
+#BRANCH_DATES=19000101,18500101
+#BRANCH_METHOD="standard"
+#PARENT_MIP_ERA=none              #ONLY CMIP6
+#PARENT_ACTIVITY_ID=none
+#PARENT_EXPERIMENT=none
+#cmip6-spinup-HR
+#PARENT_EXPERIMENT_ID=none
+#PARENT_SOURCE_ID=none
+#PARENT_TIME_UNITS="days since 1850-1-1 00:00:00"
+#PARENT_VARIANT_LABEL=none
+SUB_EXPERIMENT_ID=none
+SUB_EXPERIMENT=none
+#
+#model info:
+#
+SOURCE_ID=MPI-ESM1-2-1-CR
+#MODEL_ID=MPI-ESM1-2-1-CR
+REFERENCES="The paper is currently in open discussion under the following citation:\nMikolajewicz, U., Kapsch, M.-L., Schannwell, C., Six, K. D., Ziemen, F. A., Bagge, M., Baudouin, J.-P., Erokhina, O., Gayler, V., Klemann, V., Meccia, V. L., Mouchet, A., and Riddick, T.: Deglaciation and abrupt events in a coupled comprehensive atmosphere–ocean–ice sheet–solid earth model, Clim. Past Discuss. [preprint], https://doi.org/10.5194/cp-2024-55, in review, 2024"
+SOURCE="MPI-ESM1.2.1-CR (2022): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none, prescribed\nland: JSBACH3.20, River Transport Model\nlandIce: none / mPISM 0.7 (10 km x 10 km (NH), 15 km x 15 km (SH), 121 levels)\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none, prescribed\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none / VILMA-1D"
+SOURCE_TYPE="AOGCM ISM SLM"
+#CALENDAR=proleptic_gregorian
+#MAPPING_TABLE_DIR=cmip6_mapping_tables
+#MAPPING_TABLE=MPIESM1_CMIP6.txt
+#GRID_FILE_DIR=__GRID_FILE_DIR
+#GRID_FILE=__GRID_FILE
+GRID=gn
+GRID_LABEL=gn
+#NOMINAL_RESOLUTION read from external file
+#For MPIOM/HAMOCC CR / GR30
+#NOMINAL_RESOLUTION="500 km"
+#For ECHAM6/JSBACH CR / T31
+#NOMINAL_RESOLUTION="500 km"
+CHAR_AXIS_TYPELICE="ice_free_land","ice_sheet","floating_ice","ice_free_sea"
+CHAR_AXIS_VEGTYPE=bare_land,glacier,tropical_evergreen_trees,tropical_deciduous_trees,extra-tropical_evergreen_trees,extra-tropical_deciduous_trees,raingreen_shrubs,deciduous_shrubs,C3_grass,C4_grass,C3_pasture,C4_pasture,C3_crops,C4_crops
+CHAR_AXIS_SOILPOOLS=acid-soluble_below_ground_leaf_litter,acid-soluble_below_ground_woody_litter,water-soluble_below_ground_leaf_litter,water-soluble_below_ground_woody_litter,ethanol-soluble_below_ground_leaf_litter,ethanol-soluble_below_ground_woody_litter,non-soluble_below_ground_leaf_litter,non-soluble_below_ground_woody_litter,humus_from_leaf_litter,humus_from_woody_litter
+CHAR_AXIS_LANDUSE=primary_and_secondary_land,crops,pastures,urban
+CHAR_AXIS_BASIN=global_ocean,atlantic_arctic_ocean,indian_pacific_ocean
+CHAR_AXIS_OLINE=barents_opening,bering_strait,canadian_archipelago,denmark_strait,drake_passage,english_channel,pacific_equatorial_undercurrent,faroe_scotland_channel,florida_bahamas_strait,fram_strait,iceland_faroe_channel,indonesian_throughflow,mozambique_channel,taiwan_luzon_straits,windward_passage
+CHAR_AXIS_SILINE=barents_opening,bering_strait,canadian_archipelago,fram_strait
+#spectband for the 3 aeropt variables
+CHAR_AXIS_SPECTBAND_aeroptbnd=1818182
+CHAR_AXIS_SPECTBAND_aeroptbnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_aeroptbnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_aerssabnd=1818182
+CHAR_AXIS_SPECTBAND_aerssabnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_aerssabnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_aerasymbnd=1818182
+CHAR_AXIS_SPECTBAND_aerasymbnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_aerasymbnd_UNITS="m-1"
+# 2d-spectband for the albedos
+CHAR_AXIS_SPECTBAND_albdirbnd=3225689,766690
+CHAR_AXIS_SPECTBAND_albdirbnd_BOUNDS=5000000,1451379,1451379,82001
+CHAR_AXIS_SPECTBAND_albdirbnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_albdiffbnd=3225689,766690
+CHAR_AXIS_SPECTBAND_albdiffbnd_BOUNDS=5000000,1451379,1451379,82001
+CHAR_AXIS_SPECTBAND_albdiffbnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_solbnd=1818182
+CHAR_AXIS_SPECTBAND_solbnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_solbnd_UNITS="m-1"
+#T_AXIS="cmip" causes the time axis to be rewritten according to the requested time stamps
+T_AXIS="cmip"
+#OUTPUT_MODE="r"
+OUTPUT_MODE="a"
+MAX_SIZE=0
+SAVE_CHUNK="n"
+#For files without time axis, if one is requested:
+#firsttimeval=1 (1850-01-02 00:00:00 wenn days since 1850-01-01 00:00:00)
+FIRSTTIMEVAL=1
+#
+# institution and contact info:
+#
+# none of these can be set in the command line
+#
+DEFLATE_LEVEL=1
+CONTACT=palmod2-mpi-esm@dkrz.de
+#
+# MPI-M
+#
+INSTITUTION_ID=MPI-M
+INSTITUTION="Max Planck Institute for Meteorology, Hamburg 20146, Germany"
+LICENSE="PalMod2 model data produced by MPI-M is licensed under a Creative Commons Attribution ShareAlike 4.0 International License (https://creativecommons.org/licenses). The data producers and data providers make no warranty, either express or implied, including, but not limited to, warranties of merchantability and fitness for a particular purpose. All liabilities arising from the supply of the information (including any liability arising in negligence) are excluded to the fullest extent permitted by law."
+#MIP_TABLE_DIR="/pool/data/CMIP6/cmip6-cmor-tables/Tables/"
+MIP_TABLE_DIR="/work/bm0021/PalMod2/cmor_tables/"
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p3f1-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p3f1-CR.cmor_ctl
+#! /bin/sh -e
+#SBATCH --account=bm0021
+#SBATCH --qos=esgf
+#SBATCH --time=48:00:00
+#SBATCH --partition=shared
+#SBATCH --mail-type=FAIL
+#SBATCH --output=%x_%j.log
+
+# Number of sub-simulations
+n_subsim=4
+#initial and final year of the experiment
+spinup=1000
+iniyeararr=( $((4000+spinup)) 2000 2000 2000 )
+finyeararr=( 8999 8999 8999 8999 )
+iniyearoffset=0
+iniyear=$(( 1 + $iniyearoffset + $spinup ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
+#offset compared to startyear
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) )
+# Time axis (26000-1 BP)
+##     exp time       total yrs   new exp time   offset
+#a - 26000-21001 BP - 5000 yrs -     1- 5000 -   - 3999
+#b - 21000-14001 BP - 7000 yrs -  5001-12000 -     3001
+#c - 14000- 7001 BP - 7000 yrs - 12001-19000 -    10001
+#d -  7000-    1 BP - 7000 yrs - 19001-26000 -    17001
+
+# Chunksize
+#  Note that there are prerequisites concerning the chunksize:
+#  - the defined decadal timeslizes may not overlap between two chunks!
+#  - decadal bounds have to be defined accordingly - "???01231" in this case
+#  - usually this is the case by default due to decadal model output!
+chunksize=100
+
+# Create list
+for i in $(seq 0 $(( n_subsim - 1 )) )
+do
+  for yyyy in $(seq ${iniyeararr[$i]} ${finyeararr[$i]} )
+  do
+    paramarr+=( "$i $yyyy" )
+  done
+done
+
+rm -f cmor_ctl_list_interactive_synch_r1i1p3f1
+in=0
+for i in $(seq $((iniyear-iniyear)) $chunksize $((finyear-iniyear)))
+do
+  [[ $i -eq $((iniyear-iniyear)) ]] && continue
+  li=$(echo ${paramarr[$in]} | cut -d " " -f1)
+  ly=$(echo ${paramarr[$in]} | cut -d " " -f2)
+  ui=$(echo ${paramarr[$((i-1))]} | cut -d " " -f1)
+  uy=$(echo ${paramarr[$((i-1))]} | cut -d " " -f2)
+  in=$i
+  echo "$li $ui $ly $uy" >> cmor_ctl_list_interactive_synch_r1i1p3f1
+done
+[[ $i -ne $((finyear-iniyear)) ]] && {
+  li=$(echo ${paramarr[$in]} | cut -d " " -f1)
+  ly=$(echo ${paramarr[$in]} | cut -d " " -f2)
+  ui=$(echo ${paramarr[$((finyear-iniyear))]} | cut -d " " -f1)
+  uy=$(echo ${paramarr[$((finyear-iniyear))]} | cut -d " " -f2)
+  echo "$li $ui $ly $uy" >> cmor_ctl_list_interactive_synch_r1i1p3f1
+}
+
+make -k -j 19 -f  transient-deglaciation-interactive_r1i1p3f1-CR.cmor_ctl.mk "$@"
+
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p3f1-CR.cmor_ctl.mk
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p3f1-CR.cmor_ctl.mk
+EXP_ID = transient-deglaciation-interactive_r1i1p3f1-CR
+
+ACCOUNT = bm0021
+PARTITION = compute
+
+SUFFIX =
+ARCHIVE_SUFFIX =
+
+TARGET_STEP=cmor
+NUMCHUNKS=$(shell cat cmor_ctl_list_interactive_synch_r1i1p3f1 | wc -l )
+
+AGGR_TIME = 00:30:00
+DIAGS_TIME = 01:00:00
+CMOR_TIME = 04:00:00
+
+CMOR_FLAGS =
+
+SBATCHFLAGS=-W
+
+MESSAGE = $$(date +%Y-%m-%dT%H:%M:%S): CMIP6 diagnostics and CMOR rewriting ($(SUFFIX:_%=%))
+
+.PHONY: all
+.PRECIOUS: targets/$(EXP_ID).aggr$(SUFFIX).% targets/$(EXP_ID).diags$(SUFFIX).% targets/$(EXP_ID).cmor$(SUFFIX).%
+
+
+
+all:
+	@echo "$(MESSAGE) started"
+	mkdir -vp targets
+	$(MAKE) $(MFLAGS) -f $(MAKEFILE_LIST) $(addprefix targets/$(EXP_ID).$(TARGET_STEP)$(SUFFIX).,$(shell seq 1 $(NUMCHUNKS)))
+	@echo "$(MESSAGE) finished"
+
+targets/$(EXP_ID).aggr$(SUFFIX).%:
+	sbatch $(SBATCHFLAGS) --time=$(AGGR_TIME) --job-name=$(EXP_ID)_runpp_aggr$(SUFFIX) --output=%x_%j.log --comment=$* $(EXP_ID).runpp -s '$(SUFFIX)' -S '$(ARCHIVE_SUFFIX)' -A $(shell sed -n '$(*)'p < cmor_ctl_list_interactive_synch_r1i1p3f1)
+	@touch $@
+
+targets/$(EXP_ID).diags$(SUFFIX).%: targets/$(EXP_ID).aggr$(SUFFIX).%
+	sbatch $(SBATCHFLAGS) --time=$(DIAGS_TIME) --job-name=$(EXP_ID)_runpp_diag$(SUFFIX) --output=%x_%j.log --comment=$* $(EXP_ID).runpp -s '$(SUFFIX)' -S '$(ARCHIVE_SUFFIX)' -d $(shell sed -n '$(*)'p < cmor_ctl_list_interactive_synch_r1i1p3f1)
+	@touch $@
+
+targets/$(EXP_ID).cmor$(SUFFIX).%: targets/$(EXP_ID).diags$(SUFFIX).%
+	sbatch $(SBATCHFLAGS) --time=$(CMOR_TIME) --job-name=$(EXP_ID)_runpp_cmor$(SUFFIX) --output=%x_%j.log --comment=$* $(EXP_ID).runpp -s '$(SUFFIX)' -S '$(ARCHIVE_SUFFIX)' -c $(CMOR_FLAGS) $(shell sed -n '$(*)'p < cmor_ctl_list_interactive_synch_r1i1p3f1)
+	@touch $@
+
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p3f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p3f1-CR.runpp
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p3f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p3f1-CR_cdocmorinfo
+#
+# experiment info:
+#
+_CONTROLLED_VOCABULARY_FILE=PalMod2_CV.json
+_FORMULA_VAR_FILE=PalMod2_formula_terms.json
+_AXIS_ENTRY_FILE=PalMod2_coordinate.json
+CONVENTIONS="CF-1.7 PalMod-2.0"
+TITLE="MPI-ESM1-2-1 output prepared for PalMod2"
+_history_template="%s ; CMOR rewrote data to be consistent with <activity_id>, <Conventions> and CF standards."
+
+EXPERIMENT_ID=transient-deglaciation-interactive
+EXPERIMENT="synchronously coupled transient deglaciation with interactive ice sheets and solid earth"
+#Specify VARIANT_LABEL explicitly or realization_index, physics_index, initialization_index and forcing_index:
+VARIANT_LABEL=r1i1p3f1
+REALIZATION_INDEX=1
+INITIALIZATION_INDEX=1
+PHYSICS_INDEX=3
+FORCING_INDEX=1
+VARIANT_INFO="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. The simulation used climatological volcanic forcing representing PI. GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978). Equivalent to p1 but with decreased vertical mixing in the ocean. For details see Mikoljawicz et al. (2024)."
+COMMENT="Earth viscosity configuration for the solid earth model VILMA. The vertical viscosity structure is piecewise constant for the following specified layers of the earth: (1) 0 to 3.480000e+06 m from earth core: 1.000000e+19 Pa s. (2) 3.480000e+06 to 5.701000e+06 m from earth core: 1.000000e+22 Pa s. (3) 5.701000e+06 to 6.291000e+06 m from earth core: 4.000000e+20 Pa s. (4) 6.291000e+06 to 6.371000e+06 m from earth core: 1.000000e+30 Pa s."
+ACTIVITY_ID=PalMod2
+#MIP_ERA=PalMod2                   #ONLY CMIP6
+PROJECT_ID=PalMod2
+REQUIRED_TIME_UNITS="days since 1-1-1 00:00:00"
+#FURTHER_INFO_URL="https://furtherinfo.es-doc.org/<mip_era>.<institution_id>.<source_id>.<experiment_id>.<sub_experiment_id>.<variant_label>
+#Branch Info: If needed use the following 2 attributes to define the branch time
+#  in Parent_Time_Units and Required_Time_Units (=child)
+#  else use BRANCH_DATES to give the dates in YYYYMMDD,YYYYMMDD format
+#BRANCH_TIME_IN_PARENT=
+#BRANCH_TIME_IN_CHILD=
+#BRANCH_DATES=19000101,18500101
+#BRANCH_METHOD="standard"
+#PARENT_MIP_ERA=none              #ONLY CMIP6
+#PARENT_ACTIVITY_ID=none
+#PARENT_EXPERIMENT=none
+#cmip6-spinup-HR
+#PARENT_EXPERIMENT_ID=none
+#PARENT_SOURCE_ID=none
+#PARENT_TIME_UNITS="days since 1850-1-1 00:00:00"
+#PARENT_VARIANT_LABEL=none
+SUB_EXPERIMENT_ID=none
+SUB_EXPERIMENT=none
+#
+#model info:
+#
+SOURCE_ID=MPI-ESM1-2-1-CR
+#MODEL_ID=MPI-ESM1-2-1-CR
+REFERENCES="The paper is currently in open discussion under the following citation:\nMikolajewicz, U., Kapsch, M.-L., Schannwell, C., Six, K. D., Ziemen, F. A., Bagge, M., Baudouin, J.-P., Erokhina, O., Gayler, V., Klemann, V., Meccia, V. L., Mouchet, A., and Riddick, T.: Deglaciation and abrupt events in a coupled comprehensive atmosphere–ocean–ice sheet–solid earth model, Clim. Past Discuss. [preprint], https://doi.org/10.5194/cp-2024-55, in review, 2024"
+SOURCE="MPI-ESM1.2.1-CR (2022): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none, prescribed\nland: JSBACH3.20, River Transport Model\nlandIce: none / mPISM 0.7 (10 km x 10 km (NH), 15 km x 15 km (SH), 121 levels)\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none, prescribed\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none / VILMA-1D"
+SOURCE_TYPE="AOGCM ISM SLM"
+#CALENDAR=proleptic_gregorian
+#MAPPING_TABLE_DIR=cmip6_mapping_tables
+#MAPPING_TABLE=MPIESM1_CMIP6.txt
+#GRID_FILE_DIR=__GRID_FILE_DIR
+#GRID_FILE=__GRID_FILE
+GRID=gn
+GRID_LABEL=gn
+#NOMINAL_RESOLUTION read from external file
+#For MPIOM/HAMOCC CR / GR30
+#NOMINAL_RESOLUTION="500 km"
+#For ECHAM6/JSBACH CR / T31
+#NOMINAL_RESOLUTION="500 km"
+CHAR_AXIS_TYPELICE="ice_free_land","ice_sheet","floating_ice","ice_free_sea"
+CHAR_AXIS_VEGTYPE=bare_land,glacier,tropical_evergreen_trees,tropical_deciduous_trees,extra-tropical_evergreen_trees,extra-tropical_deciduous_trees,raingreen_shrubs,deciduous_shrubs,C3_grass,C4_grass,C3_pasture,C4_pasture,C3_crops,C4_crops
+CHAR_AXIS_SOILPOOLS=acid-soluble_below_ground_leaf_litter,acid-soluble_below_ground_woody_litter,water-soluble_below_ground_leaf_litter,water-soluble_below_ground_woody_litter,ethanol-soluble_below_ground_leaf_litter,ethanol-soluble_below_ground_woody_litter,non-soluble_below_ground_leaf_litter,non-soluble_below_ground_woody_litter,humus_from_leaf_litter,humus_from_woody_litter
+CHAR_AXIS_LANDUSE=primary_and_secondary_land,crops,pastures,urban
+CHAR_AXIS_BASIN=global_ocean,atlantic_arctic_ocean,indian_pacific_ocean
+CHAR_AXIS_OLINE=barents_opening,bering_strait,canadian_archipelago,denmark_strait,drake_passage,english_channel,pacific_equatorial_undercurrent,faroe_scotland_channel,florida_bahamas_strait,fram_strait,iceland_faroe_channel,indonesian_throughflow,mozambique_channel,taiwan_luzon_straits,windward_passage
+CHAR_AXIS_SILINE=barents_opening,bering_strait,canadian_archipelago,fram_strait
+#spectband for the 3 aeropt variables
+CHAR_AXIS_SPECTBAND_aeroptbnd=1818182
+CHAR_AXIS_SPECTBAND_aeroptbnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_aeroptbnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_aerssabnd=1818182
+CHAR_AXIS_SPECTBAND_aerssabnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_aerssabnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_aerasymbnd=1818182
+CHAR_AXIS_SPECTBAND_aerasymbnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_aerasymbnd_UNITS="m-1"
+# 2d-spectband for the albedos
+CHAR_AXIS_SPECTBAND_albdirbnd=3225689,766690
+CHAR_AXIS_SPECTBAND_albdirbnd_BOUNDS=5000000,1451379,1451379,82001
+CHAR_AXIS_SPECTBAND_albdirbnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_albdiffbnd=3225689,766690
+CHAR_AXIS_SPECTBAND_albdiffbnd_BOUNDS=5000000,1451379,1451379,82001
+CHAR_AXIS_SPECTBAND_albdiffbnd_UNITS="m-1"
+CHAR_AXIS_SPECTBAND_solbnd=1818182
+CHAR_AXIS_SPECTBAND_solbnd_BOUNDS=2262443,1600000
+CHAR_AXIS_SPECTBAND_solbnd_UNITS="m-1"
+#T_AXIS="cmip" causes the time axis to be rewritten according to the requested time stamps
+T_AXIS="cmip"
+#OUTPUT_MODE="r"
+OUTPUT_MODE="a"
+MAX_SIZE=0
+SAVE_CHUNK="n"
+#For files without time axis, if one is requested:
+#firsttimeval=1 (1850-01-02 00:00:00 wenn days since 1850-01-01 00:00:00)
+FIRSTTIMEVAL=1
+#
+# institution and contact info:
+#
+# none of these can be set in the command line
+#
+DEFLATE_LEVEL=1
+CONTACT=palmod2-mpi-esm@dkrz.de
+#
+# MPI-M
+#
+INSTITUTION_ID=MPI-M
+INSTITUTION="Max Planck Institute for Meteorology, Hamburg 20146, Germany"
+LICENSE="PalMod2 model data produced by MPI-M is licensed under a Creative Commons Attribution ShareAlike 4.0 International License (https://creativecommons.org/licenses). The data producers and data providers make no warranty, either express or implied, including, but not limited to, warranties of merchantability and fitness for a particular purpose. All liabilities arising from the supply of the information (including any liability arising in negligence) are excluded to the fullest extent permitted by law."
+#MIP_TABLE_DIR="/pool/data/CMIP6/cmip6-cmor-tables/Tables/"
+MIP_TABLE_DIR="/work/bm0021/PalMod2/cmor_tables/"
--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p4f2-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p4f2-CR.cmor_ctl
@@ -12,7 +12,7 @@ n_subsim=4
 spinup=1000
 iniyeararr=( $((4000+spinup)) 2000 2000 2000 )
 finyeararr=( 8999 8999 8999 8999 )
-iniyearoffset=1000
+iniyearoffset=0
 iniyear=$(( 1 + $iniyearoffset + $spinup ))
 finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
@@ -27,7 +27,7 @@ inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10
 # Chunksize
 #  Note that there are prerequisites concerning the chunksize:
 #  - the defined decadal timeslizes may not overlap between two chunks!
-#  - decadal bounds have to be defined accordingly - "???01231" in this case 
+#  - decadal bounds have to be defined accordingly - "???01231" in this case
 #  - usually this is the case by default due to decadal model output!
 chunksize=100


--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p4f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p4f2-CR.runpp
@@ -43,7 +43,7 @@ alias ncatted=/sw/spack-levante/nco-5.0.6-3xkdth/bin/ncatted
 #########################################

 #Root directory of all scripts, outdata, logs etc
-EXP_ID=transient-deglaciation-prescribed-glac1d_r1i1p4f2-CR
+EXP_ID=transient-deglaciation-interactive_r1i1p4f2-CR
 #EXP_DIR_TRUNK=/work/bk1192/k204212/palmod_CS/transient-deglaciation-interactive_r1i1p4f2-CR
 EXP_DIR_TRUNK=/work/kd1292/ESGF_Buff/k204212/palmod_CS/transient-deglaciation-interactive_r1i1p4f2-CR

@@ -135,7 +135,7 @@ ca+=([vilma]="${SCRIPT_ROOT}/cdocmorinfo_vilma_CR")
 experiment=transient-deglaciation-interactive
 member=r1i1p4f2
 mip=PalMod2
-vd=v20240920
+vd=v20250105
 #vd=v$(date '+%Y%m%d') #CMOR version directory name, should be a date in format YYYYMMDD

 #Where to store errors
@@ -388,7 +388,7 @@ function run_agg
        fi

        # Perform aggregation for each submodel
-        for submodel in $slmod; do #$icemod $slmod $atmmod $ocemod $srfmod; do
+        for submodel in $icemod $slmod $atmmod $ocemod $srfmod; do

            if [[ "$submodel" == "$icemod" ]]; then
                period=$(( inv_chunk_year - 10 ))
@@ -416,7 +416,7 @@ function run_agg
            sdir=${SDIR[$y]}/outdata/$submodel
            rawsdir=${RAWSDIR[$y]}/outdata/$submodel
            rawsdir_pism_ant=$(echo $rawsdir | rev | cut -d "/" -f 3- | rev)/PISM_ANT
-            rawsdir_pism_nh=$(echo $rawsdir | rev | cut -d "/" -f 3- | rev)/PISM_NH            
+            rawsdir_pism_nh=$(echo $rawsdir | rev | cut -d "/" -f 3- | rev)/PISM_NH
            echo "... writing output to $sdir/out_aggr"
            echo "... writing log to $errtemp_agg"

@@ -478,7 +478,7 @@ function run_diag
        rm -f ${err}*

        # Perform diagnostic for each submodel
-        for submodel in $slmod; do #$icemod $slmod $atmmod $ocemod $srfmod; do
+        for submodel in $icemod $slmod $atmmod $ocemod $srfmod; do

            if [[ "$submodel" == "$icemod" ]]; then
                period=$(( inv_chunk_year - 10 ))
@@ -585,7 +585,7 @@ function run_cmor
        rm -f ${err}*

        # Perform CMOR rewrite for each submodel
-        for submodel in $slmod; do #$icemod $slmod $atmmod $ocemod $srfmod; do
+        for submodel in $icemod $slmod $atmmod $ocemod $srfmod; do

            if [[ "$submodel" == "$icemod" ]]; then
                period=$(( inv_chunk_year - 10 ))

--- a/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p4f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p4f2-CR_cdocmorinfo
@@ -16,10 +16,10 @@ REALIZATION_INDEX=1
 INITIALIZATION_INDEX=1
 PHYSICS_INDEX=4
 FORCING_INDEX=2
-VARIANT_INFO="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. The simulation was forced with time-varying volcanic forcing (Schindlbeck-Belo et al., 2024). GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978)."
+VARIANT_INFO="The data are from a synchronously coupled model simulation of the last deglaciation with the MPI-ESM1.2.1-CR-mPISM-VILMA setup, hence, including interactive ice sheets and solid earth components. Equivalent to p1 but with retuned atmosphere as well as modified sea-ice rigidity, ocean mixing, and basal sliding parameters. In addition, the simulation was forced with time-varying volcanic forcing (Schindlbeck-Belo et al., 2024). GHG forcing is after Koehler et al. (2017), Orbit after Berger (1978)."
 COMMENT="Earth viscosity configuration for the solid earth model VILMA. The vertical viscosity structure is piecewise constant for the following specified layers of the earth: (1) 0 to 3.480000e+06 m from earth core: 1.000000e+19 Pa s. (2) 3.480000e+06 to 5.701000e+06 m from earth core: 1.000000e+22 Pa s. (3) 5.701000e+06 to 6.291000e+06 m from earth core: 4.000000e+20 Pa s. (4) 6.291000e+06 to 6.371000e+06 m from earth core: 1.000000e+30 Pa s."
 ACTIVITY_ID=PalMod2
-#MIP_ERA=PalMod2                   #ONLY CMIP6  
+#MIP_ERA=PalMod2                   #ONLY CMIP6
 PROJECT_ID=PalMod2
 REQUIRED_TIME_UNITS="days since 1-1-1 00:00:00"
 #FURTHER_INFO_URL="https://furtherinfo.es-doc.org/<mip_era>.<institution_id>.<source_id>.<experiment_id>.<sub_experiment_id>.<variant_label>
@@ -28,7 +28,7 @@ REQUIRED_TIME_UNITS="days since 1-1-1 00:00:00"
 #  else use BRANCH_DATES to give the dates in YYYYMMDD,YYYYMMDD format
 #BRANCH_TIME_IN_PARENT=
 #BRANCH_TIME_IN_CHILD=
-#BRANCH_DATES=19000101,18500101 
+#BRANCH_DATES=19000101,18500101
 #BRANCH_METHOD="standard"
 #PARENT_MIP_ERA=none              #ONLY CMIP6
 #PARENT_ACTIVITY_ID=none
@@ -45,7 +45,7 @@ SUB_EXPERIMENT=none
 #
 SOURCE_ID=MPI-ESM1-2-1-CR
 #MODEL_ID=MPI-ESM1-2-1-CR
-REFERENCES="The paper ist currently in open discussion under the following citation:\nMikolajewicz, U., Kapsch, M.-L., Schannwell, C., Six, K. D., Ziemen, F. A., Bagge, M., Baudouin, J.-P., Erokhina, O., Gayler, V., Klemann, V., Meccia, V. L., Mouchet, A., and Riddick, T.: Deglaciation and abrupt events in a coupled comprehensive atmosphere–ocean–ice sheet–solid earth model, Clim. Past Discuss. [preprint], https://doi.org/10.5194/cp-2024-55, in review, 2024"
+REFERENCES="The paper is currently in open discussion under the following citation:\nMikolajewicz, U., Kapsch, M.-L., Schannwell, C., Six, K. D., Ziemen, F. A., Bagge, M., Baudouin, J.-P., Erokhina, O., Gayler, V., Klemann, V., Meccia, V. L., Mouchet, A., and Riddick, T.: Deglaciation and abrupt events in a coupled comprehensive atmosphere–ocean–ice sheet–solid earth model, Clim. Past Discuss. [preprint], https://doi.org/10.5194/cp-2024-55, in review, 2024"
 SOURCE="MPI-ESM1.2.1-CR (2022): \naerosol: none, prescribed Kinne (2010)\natmos: ECHAM6.3 (spectral T31; 96 x 48 longitude/latitude; 31 levels; top level 10 hPa)\natmosChem: none, prescribed\nland: JSBACH3.20, River Transport Model\nlandIce: none / mPISM 0.7 (10 km x 10 km (NH), 15 km x 15 km (SH), 121 levels)\nocean: MPIOM1.63 (bipolar GR3.0, approximately 300km; 122 x 101 longitude/latitude; 40 levels; top grid cell 0-15 m)\nocnBgchem: none, prescribed\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none / VILMA-1D"
 SOURCE_TYPE="AOGCM ISM SLM"
 #CALENDAR=proleptic_gregorian