diff --git a/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExpTK.conf b/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExpTK.conf
index 41f1887423485a699c25ee8c844b08713cad1268..aacd79de895f08d6659df3b21d6bb4c595c10bc5 100755
--- a/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExpTK.conf
+++ b/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExpTK.conf
@@ -29,7 +29,7 @@ fx            : fx            = slice: first_step
 
 USERSETTINGS
 # ---> Specify your settings for Experiment AllExp here
-TOTAL    : TOTAL  = False # only produce variables explicitly set to True
+#TOTAL    : TOTAL  = False # only produce variables explicitly set to True
 ################
 # Aggregation
 ################
@@ -37,30 +37,20 @@ echammonplev19  : interpolation = True
 echammonafter   : afterburner   = True
 jsbachdecfx     : input         = True
 echamdecfx      : input         = True
+mpiom           : gridfile      = False # basin file does not exist for GR30
 Omon            : areacello     = True
 Ofx             : areacello     = True
+#TOTAL:TOTAL=False
+#Ofx:Ofx=True
+#fx:fx=True
 ################
 # Ocean
 ################
-hfx      : Omon   = True
-msftmz   : Omon   = True
 ################
 # Atmosphere
 ################
-tas      : Amon   = True
-psl      : Amon   = True
-lossch4  : Emon   = True
-ch4      : Emon   = True
-ch4brdn  : Emon   = True
-emilnox  : AERmon = True
-fSoilCH4 : Emon   = True
-emich4   : Emon   = True
 #################
 # Land
 #################
-wetlandCH4 : Emon = True
-wetlandFrac: Emon = True
-termiteCH4 : Emon = True
-fVegFireCH4: Emon = True
-herbivoreCH4:Emon = True
+#herbivoreCH4:Emon = False # for SSP* experiments
 #<---- Specify your settings for Experiment AllExp here
diff --git a/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExpTKCH4.conf b/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExpTKCH4.conf
new file mode 100755
index 0000000000000000000000000000000000000000..6ce3a1112f4744e85c85c45328513cda9f1a7772
--- /dev/null
+++ b/cmor/mpiesm/scripts/conf/PalMod2_requested_vars_AllExpTKCH4.conf
@@ -0,0 +1,67 @@
+# ##################################################
+# ##################################################
+# CMIP6 PostProcessing Workflow Configuration File
+# ##################################################
+# This configuration file specifies which variables will be
+#  processed for which time intervals.
+# It contains + DREQSETTINGS (computed from the CMIP6-DataRequest) and
+#             + USERSETTINGS (specified by the user)
+# for each experiment.
+#
+# There are 3 types of settings: Standard, TimeSlice, GridAndTimeSlice.
+#
+# -Use '#' for comments (also inline).
+# ##################################################
+
+
+# Data Request Version: 00.00.02
+# Supported MIPs (0):
+# Supported Experiments (0):
+
+
+##################################################
+EXP=AllExp
+##################################################
+
+DREQSETTINGS
+Ofx           : Ofx           = slice: first_step
+fx            : fx            = slice: first_step
+
+USERSETTINGS
+# ---> Specify your settings for Experiment AllExp here
+TOTAL    : TOTAL  = False # only produce variables explicitly set to True
+################
+# Aggregation
+################
+echammonplev19  : interpolation = True
+echammonafter   : afterburner   = True
+jsbachdecfx     : input         = True
+echamdecfx      : input         = True
+mpiom           : gridfile      = False # basin file does not exist for GR30
+Omon            : areacello     = True
+Ofx             : areacello     = True
+################
+# Ocean
+################
+#hfx      : Omon   = True
+#msftmz   : Omon   = True
+################
+# Atmosphere
+################
+#tas      : Amon   = True
+#psl      : Amon   = True
+lossch4  : Emon   = True
+ch4      : Emon   = True
+ch4brdn  : Emon   = True
+emilnox  : AERmon = True
+fSoilCH4 : Emon   = True
+emich4   : Emon   = True
+#################
+# Land
+#################
+wetlandCH4 : Emon = True
+wetlandFrac: Emon = True
+termiteCH4 : Emon = True
+fVegFireCH4: Emon = True
+herbivoreCH4:Emon = True # Not available for SSP* simulations
+#<---- 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 c5acc532e5176bc0e936bad4e5e43a7f71718d57..f2040997bf32e6692428448a091779afc31071e0 100644
--- a/cmor/mpiesm/scripts/palmod2_aggregation_MPI-ESM1-2_jsbach.h
+++ b/cmor/mpiesm/scripts/palmod2_aggregation_MPI-ESM1-2_jsbach.h
@@ -3,9 +3,10 @@
 ################################
 suf=grb
 
-{ (if_requested $member $srfmod jsbachdecfx input $chunk '???1123124' && {
+{ (if_requested $member $srfmod jsbachdecfx input $chunk "$srfmod_input_dec_chunk" && {
 
   # Create decadal "fx"-variables
+  ifile=$input_file_jsbach
 
   #-- slm
   $cdo selname,slm $ifile $sdir/tmp_aggr/tmp_slm_$period || echo ERROR
diff --git a/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_echam6_auto.h b/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_echam6_auto.h
index 5a65c0d5b89935aa138114c8e1268febc4e07abf..d8d9881d5b21813e76056130a88fe7b0ff997c54 100644
--- a/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_echam6_auto.h
+++ b/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_echam6_auto.h
@@ -65,7 +65,7 @@ for var in $cn; do
 done
 
 #-- CMOR-rewrite for echam6 (ESM: MPI-ESM1-2) Emon
-cn='ch4 lossch4'
+cn='ch4'
 for var in $cn; do
   { (if_requested $member $atmmod Emon $var $chunk || continue
   mkdir -p $dr/$submodel/Emon_${var}
@@ -76,7 +76,7 @@ for var in $cn; do
 done
 
 #-- CMOR-rewrite for echam6 (ESM: MPI-ESM1-2) Emon
-cn='ch4brdn emich4 fSoilCH4'
+cn='ch4brdn emich4 fSoilCH4 lossch4'
 find_file "$sdir" "*_echam6_ch4_mm_${period}*" ifile >> $err.find_file.Emon 2>&1
 for var in $cn; do
   { (if_requested $member $atmmod Emon $var $chunk || continue
@@ -166,7 +166,7 @@ done
 #-- CMOR-rewrite for echam6 (ESM: MPI-ESM1-2) dec
 cn='orog'
 for var in $cn; do
-  { (if_requested $member $atmmod dec $var $chunk '???0123124' || continue
+  { (if_requested $member $atmmod dec $var $chunk '${atmmod_dec_chunk}' || continue
   mkdir -p $dr/$submodel/dec_${var}
   ifile=${sdir}/out_diag/dec_${var}_$period.nc
   echo $cdo cmor,dec,i=$it,mt=$mt,dr=$dr/$submodel/dec_${var},vd=$vd,cn=$var,ta=cmip,di=$di ${cdochain-} $ifile
diff --git a/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_jsbach_auto.h b/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_jsbach_auto.h
index 55651fc066a037e07eca9c64bf621ef918a05b7f..0e6410b6ea889152c1f8b625b30e263792edafc8 100644
--- a/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_jsbach_auto.h
+++ b/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_jsbach_auto.h
@@ -1,5 +1,5 @@
 #-- CMOR-rewrite for jsbach (ESM: MPI-ESM1-2) Emon
-cn='fVegFireCH4 termiteCH4 vegFrac wetlandCH4'
+cn='fVegFireCH4 herbivoreCH4 termiteCH4 vegFrac wetlandCH4'
 for var in $cn; do
   { (if_requested $member $srfmod Emon $var $chunk || continue
   mkdir -p $dr/$submodel/Emon_${var}
@@ -91,7 +91,7 @@ done
 #-- CMOR-rewrite for jsbach (ESM: MPI-ESM1-2) dec
 cn='mrsofc rootd sftgif sftlf'
 for var in $cn; do
-  { (if_requested $member $srfmod dec $var $chunk '???0123124' || continue
+  { (if_requested $member $srfmod dec $var $chunk '${srfmod_dec_chunk}' || continue
   mkdir -p $dr/$submodel/dec_${var}
   ifile=${sdir}/out_diag/dec_${var}_$period.nc
   echo $cdo cmor,dec,i=$it,mt=$mt,dr=$dr/$submodel/dec_${var},vd=$vd,cn=$var,ta=cmip,di=$di ${cdochain-} $ifile
diff --git a/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_mpiom_auto.h b/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_mpiom_auto.h
index 29c36a9fbd95bfda60f03ffb9ab616d6df04dade..ce06818f37b26cffeec44a2252046c2d3e38d1ce 100644
--- a/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_mpiom_auto.h
+++ b/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_mpiom_auto.h
@@ -12,7 +12,7 @@ done
 #-- CMOR-rewrite for mpiom (ESM: MPI-ESM1-2) Odec
 cn='sftof'
 for var in $cn; do
-  { (if_requested $member $ocemod Odec $var $chunk '???1123124' || continue
+  { (if_requested $member $ocemod Odec $var $chunk '${ocemod_dec_chunk}' || continue
   mkdir -p $dr/$submodel/Odec_${var}
   ifile=${sdir}/out_diag/Odec_${var}_$period.nc
   echo $cdo cmor,Odec,i=$it,mt=$mt,dr=$dr/$submodel/Odec_${var},vd=$vd,cn=$var,ta=cmip,di=$dioce ${cdochain-} $ifile
diff --git a/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_mpism_auto.h b/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_mpism_auto.h
index 2c63cb541cd4825a306c7267e1a290dc63490b51..8222d200c900fa142293a7d11d3368a0186cb23f 100644
--- a/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_mpism_auto.h
+++ b/cmor/mpiesm/scripts/palmod2_cmor-rewrite_MPI-ESM1-2_mpism_auto.h
@@ -1,7 +1,7 @@
 #-- CMOR-rewrite for mpism (ESM: MPI-ESM1-2) IdecAnt
 cn='sftgifIt'
 for var in $cn; do
-  { (if_requested $member $icemod IdecAnt $var $chunk '???1123124' || continue
+  { (if_requested $member $icemod IdecAnt $var $chunk '${icemod_dec_chunk}' || continue
   mkdir -p $dr/$submodel/IdecAnt_${var}
   ifile=${sdir}/out_diag/IdecAnt_${var}_$period.nc
   echo $cdo cmor,IdecAnt,i=$it,mt=$mt,dr=$dr/$submodel/IdecAnt_${var},vd=$vd,cn=$var,ta=cmip,di=$di $ifile
@@ -12,7 +12,7 @@ done
 #-- CMOR-rewrite for mpism (ESM: MPI-ESM1-2) IdecAnt
 cn='libmassbffl'
 for var in $cn; do
-  { (if_requested $member $icemod IdecAnt $var $chunk '???1123124' || continue
+  { (if_requested $member $icemod IdecAnt $var $chunk '${icemod_dec_chunk}' || continue
   mkdir -p $dr/$submodel/IdecAnt_${var}
   ifile=${sdir}/out_diag/IdecAnt_${var}_$period.nc
   echo $cdo cmor,IdecAnt,i=$it,mt=$mt,dr=$dr/$submodel/IdecAnt_${var},vd=$vd,cn=$var,ta=cmip,di=$di ${cdochain-} $ifile
@@ -23,7 +23,7 @@ done
 #-- CMOR-rewrite for mpism (ESM: MPI-ESM1-2) IdecGre
 cn='sftgifIt'
 for var in $cn; do
-  { (if_requested $member $icemod IdecGre $var $chunk '???1123124' || continue
+  { (if_requested $member $icemod IdecGre $var $chunk '${icemod_dec_chunk}' || continue
   mkdir -p $dr/$submodel/IdecGre_${var}
   ifile=${sdir}/out_diag/IdecGre_${var}_$period.nc
   echo $cdo cmor,IdecGre,i=$it,mt=$mt,dr=$dr/$submodel/IdecGre_${var},vd=$vd,cn=$var,ta=cmip,di=$di $ifile
@@ -34,7 +34,7 @@ done
 #-- CMOR-rewrite for mpism (ESM: MPI-ESM1-2) IdecGre
 cn='libmassbffl'
 for var in $cn; do
-  { (if_requested $member $icemod IdecGre $var $chunk '???1123124' || continue
+  { (if_requested $member $icemod IdecGre $var $chunk '${icemod_dec_chunk}' || continue
   mkdir -p $dr/$submodel/IdecGre_${var}
   ifile=${sdir}/out_diag/IdecGre_${var}_$period.nc
   echo $cdo cmor,IdecGre,i=$it,mt=$mt,dr=$dr/$submodel/IdecGre_${var},vd=$vd,cn=$var,ta=cmip,di=$di ${cdochain-} $ifile
@@ -45,7 +45,7 @@ done
 #-- CMOR-rewrite for mpism (ESM: MPI-ESM1-2) IdecAnt
 cn='acabfIs dtb orogIs tlIs tsIs'
 for var in $cn; do
-  { (if_requested $member $icemod IdecAnt $var $chunk '???1123124' || continue
+  { (if_requested $member $icemod IdecAnt $var $chunk '${icemod_dec_chunk}' || continue
   find_file "$sdir" "pism_ant_${period}*.nc" ifile >> $err.find_file.IdecAnt 2>&1
   mkdir -p $dr/$submodel/IdecAnt_${var}
   echo $cdo cmor,IdecAnt,i=$it,mt=$mt,dr=$dr/$submodel/IdecAnt_${var},vd=$vd,cn=$var,ta=cmip,di=$di ${cdochain-} $ifile
@@ -56,7 +56,7 @@ done
 #-- CMOR-rewrite for mpism (ESM: MPI-ESM1-2) IdecGre
 cn='acabfIs dtb orogIs tlIs tsIs'
 for var in $cn; do
-  { (if_requested $member $icemod IdecGre $var $chunk '???1123124' || continue
+  { (if_requested $member $icemod IdecGre $var $chunk '${icemod_dec_chunk}' || continue
   find_file "$sdir" "pism_nh_${period}*.nc" ifile >> $err.find_file.IdecGre 2>&1
   mkdir -p $dr/$submodel/IdecGre_${var}
   echo $cdo cmor,IdecGre,i=$it,mt=$mt,dr=$dr/$submodel/IdecGre_${var},vd=$vd,cn=$var,ta=cmip,di=$di ${cdochain-} $ifile
diff --git a/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_echam6_auto.h b/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_echam6_auto.h
index d7a3591cb82cd3ed848565b67dc4ade0b74dd9fd..5fb66aa4e82a8885fa86d38764d4d569a6d888f0 100644
--- a/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_echam6_auto.h
+++ b/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_echam6_auto.h
@@ -50,15 +50,6 @@
     ${sdir}/out_diag/Emon_ch4_$period.nc || echo ERROR
 }; )&; }>$err.ch4.Emon 2>&1
 
-#-- Diagnostic for echam6 (ESM: MPI-ESM1-2) variable lossch4 / table Emon
-# Editor Note: ch4_dec_mass
-{ (if_requested $member $atmmod Emon lossch4 $chunk && {
-  find_file -e            "$sdir" "*_echam6_ch4_mm_${period}*" ifile
-  $cdo -f nc -O \
-    expr,'lossch4=(var8)/(var45);' \
-    $ifile ${sdir}/out_diag/Emon_lossch4_$period.nc || echo ERROR
-}; )&; }>$err.lossch4.Emon 2>&1
-
 #-- Diagnostic for echam6 (ESM: MPI-ESM1-2) variable snw / table LImon
 { (if_requested $member $atmmod LImon snw $chunk && {
   find_file -e            "$sdir" "*_echam_${period}*" ifile
@@ -83,7 +74,7 @@
     yearmean \
     -expr,'orog=var129/9.81;' \
     $ifile ${sdir}/tmp_diag/dec_orog_$period.nc || echo ERROR
-  if_requested $member $atmmod dec orog $chunk '???0123124' && {
+  if_requested $member $atmmod dec orog $chunk '${atmmod_dec_chunk}' && {
     eval $cdo -O ensmean \
       ${sdir}/tmp_diag/dec_orog_{$((period-9))..${period}}.nc ${sdir}/out_diag/dec_orog_${period}.nc && \
       eval rm ${sdir}/tmp_diag/dec_orog_{$((period-9))..${period}}.nc || echo ERROR
diff --git a/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_jsbach_auto.h b/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_jsbach_auto.h
index b8e2e19a7899a57d67df6925eeb8da17df0b2440..6d96b511c758b10815b5f2c5b73e854cdcc6d11a 100644
--- a/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_jsbach_auto.h
+++ b/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_jsbach_auto.h
@@ -7,6 +7,15 @@
     $ifile ${sdir}/out_diag/Emon_fVegFireCH4_$period.nc || echo ERROR
 }; )&; }>$err.fVegFireCH4.Emon 2>&1
 
+#-- Diagnostic for jsbach (ESM: MPI-ESM1-2) variable herbivoreCH4 / table Emon
+# Editor Note: *_jsbach_methane_mm_DATE*|box_CH4_out_herbivore_acc
+{ (if_requested $member $srfmod Emon herbivoreCH4 $chunk && {
+  find_file -e            "$sdir" "*_jsbach_methane_mm_${period}*" ifile
+  $cdo -f nc -O \
+    expr,'herbivoreCH4=-1e-3*16.04*vertsum(var172);' \
+    $ifile ${sdir}/out_diag/Emon_herbivoreCH4_$period.nc || echo ERROR
+}; )&; }>$err.herbivoreCH4.Emon 2>&1
+
 #-- Diagnostic for jsbach (ESM: MPI-ESM1-2) variable termiteCH4 / table Emon
 # Editor Note: box_CH4_out_termite_acc
 { (if_requested $member $srfmod Emon termiteCH4 $chunk && {
@@ -197,7 +206,7 @@
     yearmean \
     -selname,mrsofc $ifile \
     ${sdir}/tmp_diag/dec_mrsofc_$period.nc || echo ERROR
-  if_requested $member $srfmod dec mrsofc $chunk '???0123124' && {
+  if_requested $member $srfmod dec mrsofc $chunk '${srfmod_dec_chunk}' && {
     eval $cdo -O ensmean \
       ${sdir}/tmp_diag/dec_mrsofc_{$((period-9))..${period}}.nc ${sdir}/out_diag/dec_mrsofc_${period}.nc && \
       eval rm ${sdir}/tmp_diag/dec_mrsofc_{$((period-9))..${period}}.nc || echo ERROR
@@ -210,7 +219,7 @@
     yearmean \
     -selname,root_depth $ifile \
     ${sdir}/tmp_diag/dec_rootd_$period.nc || echo ERROR
-  if_requested $member $srfmod dec rootd $chunk '???0123124' && {
+  if_requested $member $srfmod dec rootd $chunk '${srfmod_dec_chunk}' && {
     eval $cdo -O ensmean \
       ${sdir}/tmp_diag/dec_rootd_{$((period-9))..${period}}.nc ${sdir}/out_diag/dec_rootd_${period}.nc && \
       eval rm ${sdir}/tmp_diag/dec_rootd_{$((period-9))..${period}}.nc || echo ERROR
@@ -223,7 +232,7 @@
     yearmean \
     -selname,glac $ifile \
     ${sdir}/tmp_diag/dec_sftgif_$period.nc || echo ERROR
-  if_requested $member $srfmod dec sftgif $chunk '???0123124' && {
+  if_requested $member $srfmod dec sftgif $chunk '${srfmod_dec_chunk}' && {
     eval $cdo -O ensmean \
       ${sdir}/tmp_diag/dec_sftgif_{$((period-9))..${period}}.nc ${sdir}/out_diag/dec_sftgif_${period}.nc && \
       eval rm ${sdir}/tmp_diag/dec_sftgif_{$((period-9))..${period}}.nc || echo ERROR
@@ -236,7 +245,7 @@
     yearmean \
     -selname,slm $ifile \
     ${sdir}/tmp_diag/dec_sftlf_$period.nc || echo ERROR
-  if_requested $member $srfmod dec sftlf $chunk '???0123124' && {
+  if_requested $member $srfmod dec sftlf $chunk '${srfmod_dec_chunk}' && {
     eval $cdo -O ensmean \
       ${sdir}/tmp_diag/dec_sftlf_{$((period-9))..${period}}.nc ${sdir}/out_diag/dec_sftlf_${period}.nc && \
       eval rm ${sdir}/tmp_diag/dec_sftlf_{$((period-9))..${period}}.nc || echo ERROR
diff --git a/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_mpiom_auto.h b/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_mpiom_auto.h
index 87b5d34a80e7fc458c2a8353526456fd07951eda..666b41f8ac11e043030e9a9b9b04c9a8671a1f10 100644
--- a/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_mpiom_auto.h
+++ b/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_mpiom_auto.h
@@ -1,6 +1,6 @@
 #-- Diagnostic for mpiom (ESM: MPI-ESM1-2) variable sftof / table Odec
 # Editor Note: Calculated by div sao sao or div sos sos
-{ (if_requested $member $ocemod Odec sftof $chunk '???1123124' && {
+{ (if_requested $member $ocemod Odec sftof $chunk '${ocemod_dec_chunk}' && {
   find_file -e            "$sdir" "weto_${period}.nc" ifile
   $cdo -f nc -O \
     timmean \
diff --git a/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_mpism_auto.h b/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_mpism_auto.h
index 8d2517c379ed9a25cd30f142b450e003baf57946..ba8065a27c530bc4eaf32cfa75c35de51f8ad2c1 100644
--- a/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_mpism_auto.h
+++ b/cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_mpism_auto.h
@@ -1,6 +1,6 @@
 #-- Diagnostic for mpism (ESM: MPI-ESM1-2) variable libmassbffl / table IdecAnt
 # Editor Note: This variable needs to be converted to the correct units. So it should bmelt*ice_density (917)*secondsperyear. Requires aggregation to account for leap years. Added recipe just for testing.
-{ (if_requested $member $icemod IdecAnt libmassbffl $chunk '???1123124' && {
+{ (if_requested $member $icemod IdecAnt libmassbffl $chunk '${icemod_dec_chunk}' && {
   find_file -e            "$sdir" "pism_ant_${period}*.nc" ifile
   echo "$cdo -f nc -O \
     -expr,'libmassbffl=bmelt*917*86400*365;' \
@@ -12,7 +12,7 @@
 
 #-- Diagnostic for mpism (ESM: MPI-ESM1-2) variable libmassbffl / table IdecGre
 # Editor Note: This variable needs to be converted to the correct units. So it should bmelt*ice_density (917)*secondsperyear. Requires aggregation to account for leap years. Added recipe just for testing.
-{ (if_requested $member $icemod IdecGre libmassbffl $chunk '???1123124' && {
+{ (if_requested $member $icemod IdecGre libmassbffl $chunk '${icemod_dec_chunk}' && {
   find_file -e            "$sdir" "pism_nh_${period}*.nc" ifile
   echo "$cdo -f nc -O \
     -expr,'libmassbffl=bmelt*917*86400*365;' \
@@ -23,7 +23,7 @@
 }; )&; }>$err.libmassbffl.IdecGre 2>&1
 
 #-- Diagnostic for mpism (ESM: MPI-ESM1-2) variable sftgifIt / table IdecAnt
-{ (if_requested $member $icemod IdecAnt sftgifIt $chunk '???1123124' && {
+{ (if_requested $member $icemod IdecAnt sftgifIt $chunk '${icemod_dec_chunk}' && {
   find_file -e            "$sdir" "pism_ant_${period}*.nc" ifile
   echo "$cdo -f nc -O \
     ${cdochain-} -expr,'ifl=mask==0;gi=mask==2;fi=mask==3;o=mask==4;' \
@@ -34,7 +34,7 @@
 }; )&; }>$err.sftgifIt.IdecAnt 2>&1
 
 #-- Diagnostic for mpism (ESM: MPI-ESM1-2) variable sftgifIt / table IdecGre
-{ (if_requested $member $icemod IdecGre sftgifIt $chunk '???1123124' && {
+{ (if_requested $member $icemod IdecGre sftgifIt $chunk '${icemod_dec_chunk}' && {
   find_file -e            "$sdir" "pism_nh_${period}*.nc" ifile
   echo "$cdo -f nc -O \
     ${cdochain-} -expr,'ifl=mask==0;gi=mask==2;fi=mask==3;o=mask==4;' \
diff --git a/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR.runpp b/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR.runpp
index 82519da9bd51ae4be59051560c570228e0c5a722..5e117a37ed742aee1354a5b0cb6fccb2294f26df 100755
--- a/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR.runpp
@@ -45,11 +45,13 @@ alias ncatted=/sw/spack-levante/nco-5.0.6-3xkdth/bin/ncatted
 #Root directory of all scripts, outdata, logs etc
 EXP_ID=ssp119-future-methane_r1i1p1f1-CR
 EXP_DIR_TRUNK=/work/kd1292/k204212/palmod_TK/ssp119-future-methane_r1i1p1f1-CR
+EXP_DIR_TRUNK=/work/bk1192/k204212/palmod_TK/test_ssp119-future-methane_r1i1p1f1-CR
 
-#TODO todo initial files for aggregation
+#initial files for aggregation
+INPUT_DIR_TRUNK=/work/bm1030/from_Mistral/bm1030/m300020/mpiesm-palmod_tk_r10186_megan_next/experiments
 
 # RAW Data - path and sub-simulations
-RAW_EXP_DIR_TRUNK=/work/bk1192/WG2/WP2.3/future_scenarios_2021/SSP1-1.9
+RAW_EXP_DIR_TRUNK=/work/bk1192/from_Mistral/bk1192/WG2/WP2.3/future_scenarios_2021/SSP1-1.9
 # As usually the Experiment is composed of several sub-simulations, define them as array
 RAW_EXP_IDS=( pmt0550_20RN_SSP119_r2 )
 # Parent experiment that the first RAW_EXP_ID branched from
@@ -57,11 +59,13 @@ RAW_EXP_IDS=( pmt0550_20RN_SSP119_r2 )
 # !! spinup - the aggregation will fail if the spinup period
 # !!          is as long or longer than
 # !!          the length of the first sub simulation!
-spinup=300
+spinup=0
 RAW_EXP_ID_PARENT=( pmt0531_d )
-if [[ $spinup -ne 0 ]]; then
-    RAW_EXP_ID_PARENT[0]=${RAW_EXP_IDS[0]}
-fi
+# actually, there was a 300y spinup after the branch from pmt0531_d@5600
+#  but for simplicity this is ignored here
+#if [[ $spinup -ne 0 ]]; then
+RAW_EXP_ID_PARENT[0]=${RAW_EXP_IDS[0]}
+#fi
 
 echo ${RAW_EXP_IDS[@]}
 echo ${RAW_EXP_ID_PARENT[@]}
@@ -102,6 +106,11 @@ chunks+=([echam6]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([jsbach]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([mpiom]=" 0 ")
 chunks+=([jsbachdec]=" 0 ")
+atmmod_dec_chunk="???9123124"
+srfmod_dec_chunk="???9123124"
+srfmod_input_dec_chunk="???0123124"
+ocemod_dec_chunk="???0123124"
+icemod_dec_chunk="???0123124"
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -116,7 +125,7 @@ ca+=([jsbach]="${SCRIPT_ROOT}/cdocmorinfo_jsbach_CR")
 experiment=ssp119-future-methane
 member=r1i1p1f1
 mip=PalMod2
-vd=v20220222
+vd=v20230105
 #vd=v$(date '+%Y%m%d') #CMOR version directory name, should be a date in format YYYYMMDD
 
 #Where to store errors
@@ -126,12 +135,12 @@ errdir_diag=${errdir}_diag
 errdir_agg=${errdir}_agg
 
 #initial and final year of the experiment
-iniyeararr=($((5600+spinup))) # 1850-3049
+iniyeararr=($((5600+spinup))) # 1850-3049, after the 300y spinup
 finyeararr=(7099)
 iniyear=1850
 finyear=3049
 #offset compared to startyear
-inioffsetarr=(-4050) # 1850-3049 
+inioffsetarr=(-4050) # 1850-3049
 #  -  5900-7099 - 1200 yrs - 1850-3049 - -4050
 
 # Command line options
@@ -195,7 +204,7 @@ if [[ ${#iniyeararr[@]} != ${#inioffsetarr[@]} ]] || [[ ${#iniyeararr[@]} != ${#
 else
   for y in $(seq 0 $(( n_arr - 1 )) ); do
     SDIR+=( ${EXP_DIR_TRUNK}/${RAW_EXP_IDS[$y]} )
-    RAWSDIR+=( ${RAW_EXP_DIR_TRUNK}/${RAW_EXP_IDS[$y]} )
+    RAWSDIR+=( ${RAW_EXP_DIR_TRUNK} )
     WORK_DIR+=( ${WORK_DIR_TRUNK}/${RAW_EXP_IDS[$y]} )
   done
   [[ -z "$2" ]] && [[ "$PROCESS_CLIM" == "false" ]] && die 'invalid number of parameters; need index of start and end simulations, start year (of first index) and end year (of second index) as YYYY'
@@ -332,15 +341,9 @@ 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
+        # Read input file - for SSP* runs the jsbach input is time invariant
+        echo "Reading from ${RAW_EXP_IDS[$y]} jsbach_T31GR30_11tiles_5layers_natural-veg_100k.nc"
+        input_file_jsbach=$INPUT_DIR_TRUNK/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_100k.nc
 
         # Perform aggregation for each submodel
         for submodel in $atmmod $ocemod $srfmod; do # $bgcmod; do
@@ -354,7 +357,7 @@ function run_agg
 
             # Location of input/output
             sdir=${SDIR[$y]}/outdata/$submodel
-            rawsdir=${RAWSDIR[$y]}/outdata/$submodel
+            rawsdir=${RAWSDIR[$y]}/$submodel
             echo "... writing output to $sdir/out_aggr"
             echo "... writing log to $errtemp_agg"
 
@@ -363,6 +366,7 @@ function run_agg
             mkdir -p ${sdir}/out_aggr
 
             # Create softlinks to raw model output
+            echo "... creating links from '$rawsdir' for $period, $submodel"
             ln -sf $rawsdir/*_${period}* $sdir/ 2>/dev/null 1>&2 || { echo "ERROR creating links for $period ${RAW_EXP_IDS[$y]} $submodel" && exit 1 ; }
 
             # Load and run the aggr ScriptFragment
@@ -405,7 +409,7 @@ function run_diag
 
             #Location of input/output
             sdir=${SDIR[$y]}/outdata/$submodel
-            rawsdir=${RAWSDIR[$y]}/outdata/$submodel
+            rawsdir=${RAWSDIR[$y]}/$submodel
             echo "... writing output to $sdir/out_diag"
             echo "... writing log to $errtemp_diag"
 
@@ -491,7 +495,7 @@ function run_cmor
 
             #Location of input/output
             sdir=${SDIR[$y]}/outdata/$submodel
-            rawsdir=${RAWSDIR[$y]}/outdata/$submodel
+            rawsdir=${RAWSDIR[$y]}/$submodel
             echo "... writing output to $dr/$submodel"
             echo "... writing log to $errtemp_cmor"
 
diff --git a/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR_cdocmorinfo b/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR_cdocmorinfo
index 813ef387f3a6adaccef37e95d6e95467b2e59af2..82d3f5d7fc1c1752c7aaac20ad9c7b124de6bd5c 100644
--- a/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/ssp119-future-methane_r1i1p1f1-CR_cdocmorinfo
@@ -16,12 +16,9 @@ REALIZATION_INDEX=1
 INITIALIZATION_INDEX=1
 PHYSICS_INDEX=1
 FORCING_INDEX=1
-
-# To be set:
-#VARIANT_INFO=""
-
+VARIANT_INFO="PalMod MPI-ESM1-2-CR setup with full methane cycle. Model physics identical to MPI-ESM1-2-CR simulation transient-deglaciation-prescribed_r1i1p1f1."
 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>
@@ -33,13 +30,13 @@ REQUIRED_TIME_UNITS="days since 1-1-1 00:00:00"
 #
 # todo: set BRANCH_DATES
 #
-#BRANCH_DATES=19000101,18500101 
-#BRANCH_METHOD="standard"
+BRANCH_DATES=239010101,18500101
+BRANCH_METHOD="standard"
 #PARENT_MIP_ERA=none              #ONLY CMIP6
 PARENT_ACTIVITY_ID=PalMod2
 PARENT_EXPERIMENT="transient deglaciation with GLAC-1D ice sheets, methane preliminary"
 PARENT_EXPERIMENT_ID=transient-deglaciation-prescribed-glac1d-methane
-PARENT_SOURCE_ID=MPI-ESM1-2-CR
+PARENT_SOURCE_ID=MPI-ESM1-2-CR-CH4
 PARENT_TIME_UNITS="days since 1-1-1 00:00:00"
 #
 # todo: set PARENT_VARIANT_LABEL
@@ -50,13 +47,11 @@ SUB_EXPERIMENT=none
 #
 #model info:
 #
-SOURCE_ID=MPI-ESM1-2-CR
-#MODEL_ID=MPI-ESM1-2-CR
-#
-# todo: update REFERENCES
-#
-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_ID=MPI-ESM1-2-CR-CH4
+#MODEL_ID=MPI-ESM1-2-CR-CH4
+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: unnamed (accelerated methane chemistry)\nland: JSBACH3.20, River Transport Model\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, prescribed\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none, prescribed"
 SOURCE_TYPE="AOGCM CHEM"
 #CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
@@ -97,6 +92,8 @@ 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"
+# PalMod
+CHAR_AXIS_TYPELICE="ice sheet","ice shelf","ocean","ice-free land"
 #T_AXIS="cmip" causes the time axis to be rewritten according to the requested time stamps
 T_AXIS="cmip"
 #OUTPUT_MODE="r"
diff --git a/cmor/mpiesm/scripts/tables/MPI-ESM1-2_echam6_PalMod2_mapping.txt b/cmor/mpiesm/scripts/tables/MPI-ESM1-2_echam6_PalMod2_mapping.txt
index e7a1795c97dfeb70adcf4b8bc160674fb3bb037b..2c486c56e793ded9f16a5299f9fd96a39226f042 100644
--- a/cmor/mpiesm/scripts/tables/MPI-ESM1-2_echam6_PalMod2_mapping.txt
+++ b/cmor/mpiesm/scripts/tables/MPI-ESM1-2_echam6_PalMod2_mapping.txt
@@ -1,7 +1,7 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Automatically created Variable-Mapping-Table for ECHAM6
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-! Creation Date: 2022-09-08 08:11:47.861865 UTC
+! Creation Date: 2023-01-12 12:50:26.815773 UTC
 ! Data Request Version: 00.00.03
 ! $URL$
 ! $Rev$  $Date$
@@ -38,7 +38,6 @@
 ! OtherAdditional keywords are allowed; they will be ignored by the cmor operator.
 !
 !----------------------------------------------------------------------------------------
-&parameter cmor_name="emilnox"                      code="3"   project_mip_table="AERmon"       units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="ua"                           code="131" project_mip_table="AERmon"       units="m s-1"         cell_methods="m" z_axis="alevel"      /
 &parameter cmor_name="va"                           code="132" project_mip_table="AERmon"       units="m s-1"         cell_methods="m" z_axis="alevel"      /
 &parameter cmor_name="cl"                           code="223" project_mip_table="Amon"         units="1"             cell_methods="m" z_axis="alevel"      /
@@ -73,6 +72,7 @@
 &parameter cmor_name="ch4brdn"                      code="8"   project_mip_table="Emon"         units="kg m-2"        cell_methods="m"  /
 &parameter cmor_name="emich4"                       code="11"  project_mip_table="Emon"         units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="fSoilCH4"                     code="22"  project_mip_table="Emon"         units="kg m-2 s-1"    cell_methods="m"  /
+&parameter cmor_name="lossch4"                      code="45"  project_mip_table="Emon"         units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="tdps"                         code="168" project_mip_table="Emon"         units="K"             cell_methods="m"  /
 &parameter cmor_name="snd"                          code="141" project_mip_table="LImon"        units="m"             cell_methods="m"  /
 &parameter cmor_name="snm"                          code="218" project_mip_table="LImon"        units="kg m-2 s-1"    cell_methods="m"  /
@@ -81,12 +81,12 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Diagnosed Variables listed in the following: !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+&parameter cmor_name="emilnox"                                 project_mip_table="AERmon"       units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="evspsbl"                                 project_mip_table="Amon"         units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="pr"                                      project_mip_table="Amon"         units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="rlds"                                    project_mip_table="Amon"         units="W m-2"         cell_methods="m" positive="d"  /
 &parameter cmor_name="rsds"                                    project_mip_table="Amon"         units="W m-2"         cell_methods="m" positive="d"  /
 &parameter cmor_name="ch4"                                     project_mip_table="Emon"         units="kg kg-1"       cell_methods="m" z_axis="alevel"      /
-&parameter cmor_name="lossch4"                                 project_mip_table="Emon"         units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="snw"                                     project_mip_table="LImon"        units="kg m-2"        cell_methods="m"  /
 &parameter cmor_name="mrros"                                   project_mip_table="Lmon"         units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="orog"                                    project_mip_table="dec"          units="m"             cell_methods="m"  /
diff --git a/cmor/mpiesm/scripts/tables/MPI-ESM1-2_jsbach_PalMod2_mapping.txt b/cmor/mpiesm/scripts/tables/MPI-ESM1-2_jsbach_PalMod2_mapping.txt
index 1f45f498c2c534d0fe9613a6ff7fe493ee47d785..4b74c264b5c658e6ecefb5dc19eff1f86f951867 100644
--- a/cmor/mpiesm/scripts/tables/MPI-ESM1-2_jsbach_PalMod2_mapping.txt
+++ b/cmor/mpiesm/scripts/tables/MPI-ESM1-2_jsbach_PalMod2_mapping.txt
@@ -1,7 +1,7 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Automatically created Variable-Mapping-Table for JSBACH
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-! Creation Date: 2022-09-07 18:03:29.765270 UTC
+! Creation Date: 2023-01-11 15:00:46.272866 UTC
 ! Data Request Version: 00.00.03
 ! $URL$
 ! $Rev$  $Date$
@@ -69,6 +69,7 @@
 ! Diagnosed Variables listed in the following: !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 &parameter cmor_name="fVegFireCH4"                             project_mip_table="Emon"         units="kg m-2 s-1"    cell_methods="m"  /
+&parameter cmor_name="herbivoreCH4"                            project_mip_table="Emon"         units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="termiteCH4"                              project_mip_table="Emon"         units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="wetlandCH4"                              project_mip_table="Emon"         units="kg m-2 s-1"    cell_methods="m"  /
 &parameter cmor_name="sbl"                                     project_mip_table="LImon"        units="kg m-2 s-1"    cell_methods="m"  /
diff --git a/cmor/mpiesm/scripts/test_methane b/cmor/mpiesm/scripts/test_methane
index 980cb8ff1005033be181ea3ad8552a7efccfabb5..56338563691e04feb9f07495ad2ef19307388f72 100755
--- a/cmor/mpiesm/scripts/test_methane
+++ b/cmor/mpiesm/scripts/test_methane
@@ -10,19 +10,19 @@
 n_subsim=4
 #initial and final year of the experiment
 spinup=1000
-iniyeararr=( $(( 2100 + spinup )) 2000 2000 2000 )
+iniyeararr=( $(( 3000 + spinup )) 2000 2000 2000 )
 finyeararr=( 8999 8999 8999 5999 )
 iniyear=$(( 1 + spinup ))
-finyear=24900
+finyear=24000
 #offset compared to startyear
-inioffsetarr=(-2099 4901 11901 18901)
+inioffsetarr=(-2999 4001 11001 18001)
 
-# Time axis (24900-1 BP)
+# Time axis (24000-1 BP)
 ##     exp time       total yrs   new exp time   offset
-#a - 24900-18001 BP - 6900 yrs -     1- 6900 -   - 2099
-#b - 18000-11001 BP - 7000 yrs -  6901-13900 -     4901 
-#c - 11000- 4001 BP - 7000 yrs - 13901-20900 -    11901
-#d -  4000-    1 BP - 4000 yrs - 20901-24900 -    18901
+#a - 24000-18001 BP - 6000 yrs -     1- 6000 -   - 2999
+#b - 18000-11001 BP - 7000 yrs -  6001-13000 -     4001
+#c - 11000- 4001 BP - 7000 yrs - 13001-20000 -    11001
+#d -  4000-    1 BP - 4000 yrs - 20001-24000 -    18001
 
 # Chunksize
 #  Note that there are prerequisites concerning the chunksize:
@@ -36,7 +36,7 @@ chunksize=100
 for i in 0
 do
   #for yyyy in $(seq ${iniyeararr[$i]} ${finyeararr[$i]} )
-  for yyyy in $(seq 6000 6999)
+  for yyyy in $(seq 4000 6999)
   do
     paramarr+=( "$i $yyyy" )
   done
diff --git a/cmor/mpiesm/scripts/test_ssp119 b/cmor/mpiesm/scripts/test_ssp119
index f1019e08d110ee2efad48ee778c7e3c3f5d0a643..873db8cfdc0de38cd335592ea810f6b6e92abd78 100755
--- a/cmor/mpiesm/scripts/test_ssp119
+++ b/cmor/mpiesm/scripts/test_ssp119
@@ -30,7 +30,7 @@ chunksize=100
 for i in 0
 do
   #for yyyy in $(seq ${iniyeararr[$i]} ${finyeararr[$i]} )
-  for yyyy in $(seq 6000 6999)
+  for yyyy in $(seq 5900 6999)
   do
     paramarr+=( "$i $yyyy" )
   done
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR.cmor_ctl b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR.cmor_ctl
index ce5dbcdbe8e6dc4262ffba0869d02835fdefaeeb..10e36464cb7eb3f1574295937347e54fd0ae2c15 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR.cmor_ctl
@@ -10,19 +10,19 @@
 n_subsim=4
 #initial and final year of the experiment
 spinup=1000
-iniyeararr=( $(( 2100 + spinup )) 2000 2000 2000 )
+iniyeararr=( $(( 3000 + spinup )) 2000 2000 2000 )
 finyeararr=( 8999 8999 8999 5999 )
 iniyear=$(( 1 + spinup ))
-finyear=24900
+finyear=24000
 #offset compared to startyear
-inioffsetarr=(-2099 4901 11901 18901)
+inioffsetarr=(-2999 4001 11001 18001)
 
-# Time axis (24900-1 BP)
+# Time axis (24000-1 BP)
 ##     exp time       total yrs   new exp time   offset
-#a - 24900-18001 BP - 6900 yrs -     1- 6900 -   - 2099
-#b - 18000-11001 BP - 7000 yrs -  6901-13900 -     4901 
-#c - 11000- 4001 BP - 7000 yrs - 13901-20900 -    11901
-#d -  4000-    1 BP - 4000 yrs - 20901-24900 -    18901
+#a - 24000-18001 BP - 6000 yrs -     1- 6000 -   - 2999
+#b - 18000-11001 BP - 7000 yrs -  6001-13000 -     4001
+#c - 11000- 4001 BP - 7000 yrs - 13001-20000 -    11001
+#d -  4000-    1 BP - 4000 yrs - 20001-24000 -    18001
 
 # Chunksize
 #  Note that there are prerequisites concerning the chunksize:
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 e98209bb1bed1e2f8f0ae9f30af174f099db5028..9c212c5f4d1c408ce8cdc33ff8b559fdf25dd1eb 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
@@ -44,23 +44,25 @@ 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-methane_r1i1p1f1-CR
-EXP_DIR_TRUNK=/work/kd1292/k204212/palmod_TK/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR
+#EXP_DIR_TRUNK=/work/kd1292/k204212/palmod_TK/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR
+EXP_DIR_TRUNK=/work/bk1192/k204212/palmod_TK/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR
 
-#TODO todo initial files for aggregation
+#initial files for aggregation
+INPUT_DIR_TRUNK=/work/bk1192/WG2/WP2.3/transient_deglaciation_2022/
 
 # RAW Data - path and sub-simulations
-RAW_EXP_DIR_TRUNK=/work/bk1192/WG2/WP2.3/transient_deglaciation_2020/
+RAW_EXP_DIR_TRUNK=/work/bk1192/WG2/WP2.3/transient_deglaciation_2022/
 # As usually the Experiment is composed of several sub-simulations, define them as array
-RAW_EXP_IDS=( pmt0531_a  pmt0531_b  pmt0531_c  pmt0531_d )
+RAW_EXP_IDS=( pmt0585_a  pmt0585_b  pmt0585_c  pmt0585_d )
 # Parent experiment that the first RAW_EXP_ID branched from
 #  (used for input data of the parent that is required for aggregation)
 # !! spinup - the aggregation will fail if the spinup period
 # !!          is as long or longer than
 # !!          the length of the first sub simulation!
 spinup=1000
-RAW_EXP_ID_PARENT=( None pmt0531_a  pmt0531_b  pmt0531_c ) # cmorizing 24900-1 BP as 1-24900 AD
+RAW_EXP_ID_PARENT=( None pmt0585_a  pmt0585_b  pmt0585_c )
 if [[ $spinup -ne 0 ]]; then
-    RAW_EXP_ID_PARENT[0]=${RAW_EXP_IDS[0]} # cmorizing (24900-spinup)-1 BP as 1-(24900-spinup) AD
+    RAW_EXP_ID_PARENT[0]=${RAW_EXP_IDS[0]}
 fi
 
 echo ${RAW_EXP_IDS[@]}
@@ -102,12 +104,17 @@ chunks+=([echam6]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([jsbach]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([mpiom]=" 0 ")
 chunks+=([jsbachdec]=" 0 ")
+atmmod_dec_chunk="???0123124"
+srfmod_dec_chunk="???0123124"
+srfmod_input_dec_chunk="???1123124"
+ocemod_dec_chunk="???1123124"
+icemod_dec_chunk="???1123124"
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
 # Multiple files possible with "," as delimiter, eg:
 #  it=expinfo.txt,userinfo.txt,modelinfo.txt
-it_temp="${SCRIPT_ROOT}/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR_cdocmorinfo"
+it_temp="${SCRIPT_ROOT}/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR_cdocmorinfo"
 ca+=([mpiom]="${SCRIPT_ROOT}/cdocmorinfo_mpiom_CR")
 ca+=([echam6]="${SCRIPT_ROOT}/cdocmorinfo_echam6_CR")
 ca+=([jsbach]="${SCRIPT_ROOT}/cdocmorinfo_jsbach_CR")
@@ -116,7 +123,7 @@ ca+=([jsbach]="${SCRIPT_ROOT}/cdocmorinfo_jsbach_CR")
 experiment=transient-deglaciation-prescribed-glac1d-methane
 member=r1i1p1f1
 mip=PalMod2
-vd=v20220222
+vd=v20230105
 #vd=v$(date '+%Y%m%d') #CMOR version directory name, should be a date in format YYYYMMDD
 
 #Where to store errors
@@ -126,26 +133,26 @@ errdir_diag=${errdir}_diag
 errdir_agg=${errdir}_agg
 
 #initial and final year of the experiment
-iniyeararr=($((3000+spinup)) 2000 2000 2000) # 26000-1 BP
-finyeararr=(7999 8999 8999 8999)
+iniyeararr=($((3000+spinup)) 2000 2000 2000) # 24000-1 BP
+finyeararr=(8999 8999 8999 5999)
 iniyear=1
-finyear=$(( 26000 - $spinup ))
+finyear=$(( 24000 - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-2999 3001 10001 17001) # 26000-1 BP
-# Time axis (24900-1 BP)
+inioffsetarr=(-2999 4001 11001 18001) # 24000-1 BP
+# Time axis (24000-1 BP)
 ##     exp time       total yrs   new exp time   offset
-#a - 24900-18001 BP - 6900 yrs -     1- 6900 -   - 2099
-#b - 18000-11001 BP - 7000 yrs -  6901-13900 -     4901 
-#c - 11000- 4001 BP - 7000 yrs - 13901-20900 -    11901
-#d -  4000-    1 BP - 4000 yrs - 20901-24900 -    18901
+#a - 24000-18001 BP - 6000 yrs -     1- 6000 -   - 2999
+#b - 18000-11001 BP - 7000 yrs -  6001-13000 -     4001
+#c - 11000- 4001 BP - 7000 yrs - 13001-20000 -    11001
+#d -  4000-    1 BP - 4000 yrs - 20001-24000 -    18001
 
-# /work/bk1192/WG2/WP2.3/transient_deglaciation_2020 contains a transient deglaciation experiment, 
+# /work/bk1192/WG2/WP2.3/transient_deglaciation_2020 contains a transient deglaciation experiment,
 # from 24900 yrs BP to 0 BP (present = 1950CE). For technical reasons it is split into 4 parts:
 #
-# 1) pmt0531_a from 24900 yrs BP (model year 2100) to 18001 yrs BP (model year 8999)
-# 2) pmt0531_b from 18000 yrs BP (model year 2000) to 11001 yrs BP (model year 8999)
-# 3) pmt0531_c from 11000 yrs BP (model year 2000) to 4001 yrs BP  (model year 8999)
-# 4) pmt0531_d from 4000 yrs BP (model year 2000) to 1 yr BP (model year 5999)
+# pmt0585_a (24000-18001 BP = 3000-8999);
+# pmt0585_b (18000-11001 BP = 2000-8999);
+# pmt0585_c (11000-4001 BP = 2000-8999);
+# pmt0585_d (4000-1 BP = 2000-5999)
 #
 
 # Command line options
@@ -350,10 +357,10 @@ function run_agg
         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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${RAW_EXP_ID_PARENT[$y]}/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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${RAW_EXP_IDS[$y]}/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
         fi
 
         # Perform aggregation for each submodel
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
old mode 100755
new mode 100644
index c416f66c9d0ab0ed572dab6e4348408c9720ff76..6c78e003f7c138ab86c8ac75408c485ddee8e4c0
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d-methane_r1i1p1f1-CR_cdocmorinfo
@@ -16,12 +16,9 @@ REALIZATION_INDEX=1
 INITIALIZATION_INDEX=1
 PHYSICS_INDEX=1
 FORCING_INDEX=1
-#
-# todo: set VARIANT_INFO
-#
-VARIANT_INFO="tbd" # we can set it in relation to Marie's definition of model version
+VARIANT_INFO="PalMod MPI-ESM1-2-CR setup with full methane cycle. Model physics identical to MPI-ESM1-2-CR simulation transient-deglaciation-prescribed_r1i1p1f1."
 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>
@@ -30,7 +27,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
@@ -44,13 +41,11 @@ SUB_EXPERIMENT=none
 #
 #model info:
 #
-SOURCE_ID=MPI-ESM1-2-CR
-#MODEL_ID=MPI-ESM1-2-CR
-#
-# todo: Update References
-#
-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_ID=MPI-ESM1-2-CR-CH4
+#MODEL_ID=MPI-ESM1-2-CR-CH4
+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: unnamed (accelerated methane chemistry)\nland: JSBACH3.20, River Transport Model\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, prescribed\nseaIce: unnamed (thermodynamic (Semtner zero-layer) dynamic (Hibler 79) sea ice model)\nsolidLand: none, prescribed"
 SOURCE_TYPE="AOGCM CHEM"
 #CALENDAR=proleptic_gregorian
 #MAPPING_TABLE_DIR=cmip6_mapping_tables
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 60866f56646f9d4b93a4b21bc448f0e7bb940515..9d78a7001428cd1fffe00ca86004d6d77b054a9f 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp
@@ -10,8 +10,8 @@
 #SBATCH --mem=100G
 #SBATCH --time=01:30:00
 #SBATCH --mail-type=none
-#SBATCH --account=bm0021
-#SBATCH --qos=esgf
+#SBATCH --account=bk1192
+##SBATCH --qos=esgf
 #SBATCH --exclusive
 #SBATCH --output=transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp_%j.log
 
@@ -46,6 +46,9 @@ alias ncatted=/sw/spack-levante/nco-5.0.6-3xkdth/bin/ncatted
 EXP_ID=transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR
 EXP_DIR_TRUNK=/work/kd1292/k204212/palmod_errataBatch/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR
 
+#initial files for aggregation
+INPUT_DIR_TRUNK=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments
+
 # RAW Data - path and sub-simulations
 RAW_EXP_DIR_TRUNK=/work/bk1192/from_Mistral/bk1192/WG1/WP1.1/MPIM_transient_deglaciation_prescribed/
 # As usually the Experiment is composed of several sub-simulations, define them as array
@@ -100,6 +103,11 @@ chunks+=([echam6]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([jsbach]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([mpiom]=" 0 ")
 chunks+=([jsbachdec]=" 0 ")
+atmmod_dec_chunk="???0123124"
+srfmod_dec_chunk="???0123124"
+srfmod_input_dec_chunk="???1123124"
+ocemod_dec_chunk="???1123124"
+icemod_dec_chunk="???1123124"
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -349,10 +357,10 @@ function run_agg
         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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
         fi
 
         # Perform aggregation for each submodel
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 a93a9c36b3e6af283f2a1d3875ad3b10d9a66ed0..a67d648c65e77e6ccc3a4e2e95295c56a8afdf6e 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR_cdocmorinfo
@@ -18,7 +18,7 @@ PHYSICS_INDEX=1
 FORCING_INDEX=1
 VARIANT_INFO="Predecessor model version to p2. p2 contains a range of bug fixes and an updated model tuning: it comprises the introduction of a brine plume parameterization in MPI-OM and a modification of sea-ice lead closure parameters, improvement of the automated topography generation scripts (e.g. baroclinic exchange in shallow fjords), and a replacement of the direct use of topographic roughness in ECHAM6 by a parameterized treatment like the one in the standard ECHAM6 model. Tuning parameters are taken from Mauritsen and Roeckner (2020) and comprise the surface value of the critical humidity profile for cloud formation (crs), the threshold between cloud water and ice (csecfrl) and the fall speed of cloud ice (cvtfall). Additionally, a mean volcanic aerosol forcing was taken into account. For details see Supporting Information of Kapsch et al., 2022."
 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>
@@ -27,7 +27,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
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 3bbbcf9383fde82a7a8f0f52d7cc6ef574a133f4..2ed21f33c7dcef000ffc123bf2f82001699aebf5 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp
@@ -10,8 +10,8 @@
 #SBATCH --mem=100G
 #SBATCH --time=01:30:00
 #SBATCH --mail-type=none
-#SBATCH --account=bm0021
-#SBATCH --qos=esgf
+#SBATCH --account=bk1192
+##SBATCH --qos=esgf
 #SBATCH --exclusive
 #SBATCH --output=transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp_%j.log
 
@@ -46,6 +46,9 @@ alias ncatted=/sw/spack-levante/nco-5.0.6-3xkdth/bin/ncatted
 EXP_ID=transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR
 EXP_DIR_TRUNK=/work/kd1292/k204212/palmod_errataBatch/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR
 
+#initial files for aggregation
+INPUT_DIR_TRUNK=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments
+
 # RAW Data - path and sub-simulations
 RAW_EXP_DIR_TRUNK=/work/bk1192/from_Mistral/bk1192/WG1/WP1.1/MPIM_transient_deglaciation_prescribed/
 # As usually the Experiment is composed of several sub-simulations, define them as array
@@ -78,9 +81,7 @@ SCRIPT_DIR=/work/bm0021/PalMod2/cmor/mpiesm/scripts
 ############################
 
 #Where to find cdo (incl. CMOR operator), eg.:
-#cdo="/home/k/k204210/cdo-git/bin/cdo -v"  # PalMod cdo cmor version
-#cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-05-11_cmor3.6.0_gcc/bin/cdo -v" # hopefully the correct version
-cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v"
+cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v" # latest version
 
 # Base directory for DataRequest related headers/scripts
 SCRIPT_ROOT=${SCRIPT_DIR}
@@ -102,6 +103,11 @@ chunks+=([echam6]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([jsbach]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([mpiom]=" 0 ")
 chunks+=([jsbachdec]=" 0 ")
+atmmod_dec_chunk="???0123124"
+srfmod_dec_chunk="???0123124"
+srfmod_input_dec_chunk="???1123124"
+ocemod_dec_chunk="???1123124"
+icemod_dec_chunk="???1123124"
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -351,10 +357,10 @@ function run_agg
         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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
         fi
 
         # Perform aggregation for each submodel
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 943b864a4d8943517786f6992eb7dae1824f56d3..12ae0fc03e3700373e3935edf675729f45de0435 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR_cdocmorinfo
@@ -18,7 +18,7 @@ PHYSICS_INDEX=2
 FORCING_INDEX=2
 VARIANT_INFO="For details see Supporting Information of Kapsch et al., 2022."
 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>
@@ -27,7 +27,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
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 9629828f6ffbeb04db8dd2c50c7e6e2d8bae7e90..6db525f31d3d89c4edf09016bd2bcd88ee90b83f 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp
@@ -10,8 +10,8 @@
 #SBATCH --mem=100G
 #SBATCH --time=01:30:00
 #SBATCH --mail-type=none
-#SBATCH --account=bm0021
-#SBATCH --qos=esgf
+#SBATCH --account=bk1192
+##SBATCH --qos=esgf
 #SBATCH --exclusive
 #SBATCH --output=transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp_%j.log
 
@@ -46,6 +46,9 @@ alias ncatted=/sw/spack-levante/nco-5.0.6-3xkdth/bin/ncatted
 EXP_ID=transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR
 EXP_DIR_TRUNK=/work/kd1292/k204212/palmod_errataBatch/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR
 
+#initial files for aggregation
+INPUT_DIR_TRUNK=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments
+
 # RAW Data - path and sub-simulations
 RAW_EXP_DIR_TRUNK=/work/bk1192/from_Mistral/bk1192/WG1/WP1.1/MPIM_transient_deglaciation_prescribed/
 # As usually the Experiment is composed of several sub-simulations, define them as array
@@ -78,9 +81,7 @@ SCRIPT_DIR=/work/bm0021/PalMod2/cmor/mpiesm/scripts
 ############################
 
 #Where to find cdo (incl. CMOR operator), eg.:
-#cdo="/home/k/k204210/cdo-git/bin/cdo -v"  # PalMod cdo cmor version
-#cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-05-11_cmor3.6.0_gcc/bin/cdo -v" # hopefully the correct version
-cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v"
+cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v" # latest version
 
 # Base directory for DataRequest related headers/scripts
 SCRIPT_ROOT=${SCRIPT_DIR}
@@ -102,6 +103,11 @@ chunks+=([echam6]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([jsbach]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([mpiom]=" 0 ")
 chunks+=([jsbachdec]=" 0 ")
+atmmod_dec_chunk="???0123124"
+srfmod_dec_chunk="???0123124"
+srfmod_input_dec_chunk="???1123124"
+ocemod_dec_chunk="???1123124"
+icemod_dec_chunk="???1123124"
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -351,10 +357,10 @@ function run_agg
         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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
         fi
 
         # Perform aggregation for each submodel
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 3219c9a5532bb5f485642d2498757d497fd60694..3f7922fde061c3ed6cd09a5ceeaf01f957d912ee 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR_cdocmorinfo
@@ -18,7 +18,7 @@ PHYSICS_INDEX=3
 FORCING_INDEX=2
 VARIANT_INFO="Equivalent to p2 but with a different model tuning. For details see Supporting Information of Kapsch et al., 2022."
 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>
@@ -27,7 +27,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
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 463971e628bd48d1eef779f6c86a93d8b9c2372f..b32ded7fc37306f57814d3421673a7710e728588 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp
@@ -10,8 +10,8 @@
 #SBATCH --mem=100G
 #SBATCH --time=01:30:00
 #SBATCH --mail-type=none
-#SBATCH --account=bm0021
-#SBATCH --qos=esgf
+#SBATCH --account=bk1192
+##SBATCH --qos=esgf
 #SBATCH --exclusive
 #SBATCH --output=transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp_%j.log
 
@@ -46,6 +46,9 @@ alias ncatted=/sw/spack-levante/nco-5.0.6-3xkdth/bin/ncatted
 EXP_ID=transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR
 EXP_DIR_TRUNK=/work/kd1292/k204212/palmod_errataBatch/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR
 
+#initial files for aggregation
+INPUT_DIR_TRUNK=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments
+
 # RAW Data - path and sub-simulations
 RAW_EXP_DIR_TRUNK=/work/bk1192/from_Mistral/bk1192/WG1/WP1.1/MPIM_transient_deglaciation_prescribed/
 # As usually the Experiment is composed of several sub-simulations, define them as array
@@ -78,9 +81,7 @@ SCRIPT_DIR=/work/bm0021/PalMod2/cmor/mpiesm/scripts
 ############################
 
 #Where to find cdo (incl. CMOR operator), eg.:
-#cdo="/home/k/k204210/cdo-git/bin/cdo -v"  # PalMod cdo cmor version
-#cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-05-11_cmor3.6.0_gcc/bin/cdo -v" # hopefully the correct version
-cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v"
+cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v" # latest version
 
 # Base directory for DataRequest related headers/scripts
 SCRIPT_ROOT=${SCRIPT_DIR}
@@ -102,6 +103,11 @@ chunks+=([echam6]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([jsbach]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([mpiom]=" 0 ")
 chunks+=([jsbachdec]=" 0 ")
+atmmod_dec_chunk="???0123124"
+srfmod_dec_chunk="???0123124"
+srfmod_input_dec_chunk="???1123124"
+ocemod_dec_chunk="???1123124"
+icemod_dec_chunk="???1123124"
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -351,10 +357,10 @@ function run_agg
         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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
         fi
 
         # Perform aggregation for each submodel
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 6b716cac422ddc569bcade844bc3e4bd98ed319f..00443059f50337567ed188ec0cfc7f358c25e09d 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR_cdocmorinfo
@@ -18,7 +18,7 @@ PHYSICS_INDEX=1
 FORCING_INDEX=1
 VARIANT_INFO="Predecessor model version to p2. p2 contains a range of bug fixes and an updated model tuning: it comprises the introduction of a brine plume parameterization in MPI-OM and a modification of sea-ice lead closure parameters, improvement of the automated topography generation scripts (e.g. baroclinic exchange in shallow fjords), and a replacement of the direct use of topographic roughness in ECHAM6 by a parameterized treatment like the one in the standard ECHAM6 model. Tuning parameters are taken from Mauritsen and Roeckner (2020) and comprise the surface value of the critical humidity profile for cloud formation (crs), the threshold between cloud water and ice (csecfrl) and the fall speed of cloud ice (cvtfall). Additionally, a mean volcanic aerosol forcing was taken into account. For details see Supporting Information of Kapsch et al., 2022."
 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>
@@ -27,7 +27,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
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 45c83a3f66298cb120ef676ae85d94c497151f4b..ec23b3e9f56003026268a7a48e99d9cbb78c7731 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp
@@ -10,8 +10,8 @@
 #SBATCH --mem=100G
 #SBATCH --time=01:30:00
 #SBATCH --mail-type=none
-#SBATCH --account=bm0021
-#SBATCH --qos=esgf
+#SBATCH --account=bk1192
+##SBATCH --qos=esgf
 #SBATCH --exclusive
 #SBATCH --output=transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp_%j.log
 
@@ -46,6 +46,9 @@ alias ncatted=/sw/spack-levante/nco-5.0.6-3xkdth/bin/ncatted
 EXP_ID=transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR
 EXP_DIR_TRUNK=/work/kd1292/k204212/palmod_errataBatch/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR
 
+#initial files for aggregation
+INPUT_DIR_TRUNK=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments
+
 # RAW Data - path and sub-simulations
 RAW_EXP_DIR_TRUNK=/work/bk1192/from_Mistral/bk1192/WG1/WP1.1/MPIM_transient_deglaciation_prescribed/
 # As usually the Experiment is composed of several sub-simulations, define them as array
@@ -78,9 +81,7 @@ SCRIPT_DIR=/work/bm0021/PalMod2/cmor/mpiesm/scripts
 ############################
 
 #Where to find cdo (incl. CMOR operator), eg.:
-#cdo="/home/k/k204210/cdo-git/bin/cdo -v"  # PalMod cdo cmor version
-#cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-05-11_cmor3.6.0_gcc/bin/cdo -v" # hopefully the correct version
-cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v"
+cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v" # latest version
 
 # Base directory for DataRequest related headers/scripts
 SCRIPT_ROOT=${SCRIPT_DIR}
@@ -102,6 +103,11 @@ chunks+=([echam6]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([jsbach]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([mpiom]=" 0 ")
 chunks+=([jsbachdec]=" 0 ")
+atmmod_dec_chunk="???0123124"
+srfmod_dec_chunk="???0123124"
+srfmod_input_dec_chunk="???1123124"
+ocemod_dec_chunk="???1123124"
+icemod_dec_chunk="???1123124"
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -351,10 +357,10 @@ function run_agg
         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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
         fi
 
         # Perform aggregation for each submodel
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 f2980afad47a3bc28e98d18b8e5607357f670b60..0322221a3fadaedab88dcbd90a817827c61572f6 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR_cdocmorinfo
@@ -18,7 +18,7 @@ PHYSICS_INDEX=2
 FORCING_INDEX=2
 VARIANT_INFO="For details see Supporting Information of Kapsch et al., 2022."
 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>
@@ -27,7 +27,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
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 ae316514cad5ec2a0a1a7024fc46e7b655e4f552..2d9e694f3688e2550f91287312db3534190b8939 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp
@@ -10,8 +10,8 @@
 #SBATCH --mem=100G
 #SBATCH --time=01:30:00
 #SBATCH --mail-type=none
-#SBATCH --account=bm0021
-#SBATCH --qos=esgf
+#SBATCH --account=bk1192
+##SBATCH --qos=esgf
 #SBATCH --exclusive
 #SBATCH --output=transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp_%j.log
 
@@ -46,6 +46,9 @@ alias ncatted=/sw/spack-levante/nco-5.0.6-3xkdth/bin/ncatted
 EXP_ID=transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR
 EXP_DIR_TRUNK=/work/kd1292/k204212/palmod_errataBatch/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR
 
+#initial files for aggregation
+INPUT_DIR_TRUNK=/work/mh0110/m211003/mpiesm-1.2.00p1/experiments
+
 # RAW Data - path and sub-simulations
 RAW_EXP_DIR_TRUNK=/work/bk1192/from_Mistral/bk1192/WG1/WP1.1/MPIM_transient_deglaciation_prescribed/
 # As usually the Experiment is composed of several sub-simulations, define them as array
@@ -78,9 +81,7 @@ SCRIPT_DIR=/work/bm0021/PalMod2/cmor/mpiesm/scripts
 ############################
 
 #Where to find cdo (incl. CMOR operator), eg.:
-#cdo="/home/k/k204210/cdo-git/bin/cdo -v"  # PalMod cdo cmor version
-#cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-05-11_cmor3.6.0_gcc/bin/cdo -v" # hopefully the correct version
-cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v"
+cdo="/work/bm0021/cdo_incl_cmor/cdo-2022-09-20_cmor3.6.0_gcc/bin/cdo -v" # latest version
 
 # Base directory for DataRequest related headers/scripts
 SCRIPT_ROOT=${SCRIPT_DIR}
@@ -102,6 +103,11 @@ chunks+=([echam6]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([jsbach]=" 0 1 2 3 4 5 6 7 8 9 ")
 chunks+=([mpiom]=" 0 ")
 chunks+=([jsbachdec]=" 0 ")
+atmmod_dec_chunk="???0123124"
+srfmod_dec_chunk="???0123124"
+srfmod_input_dec_chunk="???1123124"
+ocemod_dec_chunk="???1123124"
+icemod_dec_chunk="???1123124"
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -351,10 +357,10 @@ function run_agg
         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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${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
+            input_file_jsbach=$INPUT_DIR_TRUNK/${RAW_EXP_IDS[$y]}/restart/topo/jsbach_T31GR30_11tiles_5layers_natural-veg_${inv_chunk_year}k.nc
         fi
 
         # Perform aggregation for each submodel
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 05cdf3e2e49c2b33bf18706634bc87e2724eaa88..a91d93bbdb8c0a8ead6e39f08df7890c1a2bd506 100644
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR_cdocmorinfo
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR_cdocmorinfo
@@ -18,7 +18,7 @@ PHYSICS_INDEX=3
 FORCING_INDEX=2
 VARIANT_INFO="Equivalent to p2 but with a different model tuning. For details see Supporting Information of Kapsch et al., 2022."
 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>
@@ -27,7 +27,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