Updates for vilma and mpism variable processing, update of...

Updates for vilma and mpism variable processing, update of transient-deglaciation-interactive cdocmorinfo

cmor/mpiesm/scripts/cdocmorinfo_mpism_CR

 #NOMINAL_RESOLUTION read from external file
 #For ice (mPISM)
-NOMINAL_RESOLUTION="25 km"
+NOMINAL_RESOLUTION="10 km"
 #Level/Depth Bounds
 grid_info="/work/bm0021/PalMod2/cmor/mpiesm/scripts/gridinfo_mPISM_L121.nc"
 switch_z=y

cmor/mpiesm/scripts/palmod2_aggregation_MPI-ESM1-2_vilma.h

@@ -2,10 +2,12 @@
 # VILMA ########################
 ################################
 
-{ (if_requested $member $slmod fx areacellsl $chunk && {
-  # Create fx variable
-  $cdo setname,glac $input_file_glac_fx $sdir/out_aggr/SLfx_glac_${period_bk}.nc || echo ERROR
-}; )&; }>$err.fx.areacellsl 2>&1
+# Not required - instead the gridarea of the direct VILMA output will be calculated in the
+#  diagnostic script
+#{ (if_requested $member $slmod fx areacellsl $chunk && {
+#  # Create fx variable
+#  $cdo setname,glac $input_file_glac_fx $sdir/out_aggr/SLfx_glac_${period_bk}.nc || echo ERROR
+#}; )&; }>$err.fx.areacellsl 2>&1
 
 
 { (if_requested $member $slmod vilmadecfx input $chunk "$slmod_input_dec_chunk" && {

cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_mpism_auto.h

 #-- 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.
+# 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 '${icemod_dec_chunk}' && {
   find_file -e            "$sdir" "pism_ant_${period}*.nc" ifile
   echo "$cdo -f nc -O \
-    -expr,'libmassbffl=bmelt*917*86400*365;' \
+    -expr,'libmassbffl=bmelt*917/86400/365;' \
        $ifile ${sdir}/out_diag/IdecAnt_libmassbffl_${period}.nc"
   $cdo -f nc -O \
-    -expr,'libmassbffl=bmelt*917*86400*365;' \
+    -expr,'libmassbffl=bmelt*917/86400/365;' \
        $ifile ${sdir}/out_diag/IdecAnt_libmassbffl_${period}.nc || echo ERROR
 }; )&; }>$err.libmassbffl.IdecAnt 2>&1
 
 #-- 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.
+# 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 '${icemod_dec_chunk}' && {
   find_file -e            "$sdir" "pism_nh_${period}*.nc" ifile
   echo "$cdo -f nc -O \
-    -expr,'libmassbffl=bmelt*917*86400*365;' \
+    -expr,'libmassbffl=bmelt*917/86400/365;' \
        $ifile ${sdir}/out_diag/IdecGre_libmassbffl_${period}.nc"
   $cdo -f nc -O \
-    -expr,'libmassbffl=bmelt*917*86400*365;' \
+    -expr,'libmassbffl=bmelt*917/86400/365;' \
        $ifile ${sdir}/out_diag/IdecGre_libmassbffl_${period}.nc || echo ERROR
 }; )&; }>$err.libmassbffl.IdecGre 2>&1
 

cmor/mpiesm/scripts/palmod2_diagnostic_MPI-ESM1-2_vilma_auto.h

@@ -16,9 +16,9 @@
 
 #-- Diagnostic for vilma (ESM: MPI-ESM1-2) variable areacellsl / table fx
 { (if_requested $member $slmod fx areacellsl $chunk && {
-  find_file -e            "$sdir" "SLfx_glac_${period_bk}.nc" ifile
+  find_file -e            "$sdir" "delta_rsl_${period_bk}.nc" ifile
   $cdo -f nc -O \
-    expr,'areacellsl=gridarea(glac);' \
+    expr,'areacellsl=gridarea(rsl);' \
     $ifile ${sdir}/out_diag/fx_areacellsl_${period_bk}.nc || echo ERROR
 }; )&; }>$err.areacellsl.fx 2>&1
 

cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p1f1-CR_cdocmorinfo → cmor/mpiesm/scripts/transient-deglaciation-interactive_r1i1p4f2-CR_cdocmorinfo

@@ -11,12 +11,12 @@ _history_template="%s ; CMOR rewrote data to be consistent with <activity_id>, <
 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=r1i1p4f1
+VARIANT_LABEL=r1i1p4f2
 REALIZATION_INDEX=1
 INITIALIZATION_INDEX=1
 PHYSICS_INDEX=4
-FORCING_INDEX=1
-VARIANT_INFO="Note, forcing includes volcanoes. 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."
+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)."
 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  
@@ -45,7 +45,7 @@ SUB_EXPERIMENT=none
 #
 SOURCE_ID=MPI-ESM1-2-1-CR
 #MODEL_ID=MPI-ESM1-2-1-CR
-REFERENCES="Kapsch, M.-L., Mikolajewicz, U., Ziemen, F., & Schannwell, C. (2022). Ocean response in transient simulations of the last deglaciation dominated by underlying ice-sheet reconstruction and method of meltwater distribution. Geophysical Research Letters, 49, e2021GL096767. https://doi.org/10.1029/2021GL096767\nMPI-ESM: Mauritsen, T. et al. (2019), Developments in the MPI‐M Earth System Model version 1.2 (MPI‐ESM1.2) and Its Response to Increasing CO2, J. Adv. Model. Earth Syst.,11, 998-1038, doi:10.1029/2018MS001400,\nMueller, W.A. et al. (2018): A high‐resolution version of the Max Planck Institute Earth System Model MPI‐ESM1.2‐HR. J. Adv. Model. EarthSyst.,10,1383–1413, doi:10.1029/2017MS001217"
+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"
 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_TYPE="AOGCM ISM SLM"
 #CALENDAR=proleptic_gregorian
@@ -60,7 +60,7 @@ GRID_LABEL=gn
 #NOMINAL_RESOLUTION="500 km"
 #For ECHAM6/JSBACH CR / T31
 #NOMINAL_RESOLUTION="500 km"
-CHAR_AXIS_TYPELICE="ice sheet","ice shelf","ocean","ice-free land"
+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