Commit 4a2aecc0 authored by Nadine Wieters's avatar Nadine Wieters

Merge branch 'ndepo_fix' into 'develop'

Ndepo fix

See merge request esm-tools/esm-runscripts!197
parents 029a7e70 47d57c6b
......@@ -34,6 +34,13 @@ awicm_set_defaults()
;;
esac
CMIP6_namelists=0
if [ "x$CMIP6_namelists" == "x1" ]; then
CMIP6_nml_tag="_CMIP6"
else
CMIP6_nml_tag=""
fi
omp_num_threads_compute_awicm=1
omp_num_threads_post_awicm=1
......
......@@ -598,7 +598,7 @@ echam_prepare_config()
if [[ "x${RUN_NUMBER_echam}" = "x1" ]]; then
echo " Taking echam namelist from NAMELIST_DIR_echam=${NAMELIST_DIR_echam}"
fi
add_to ${NAMELIST_DIR_echam}/namelist.echam namelist.echam
add_to ${NAMELIST_DIR_echam}/namelist.echam${CMIP6_nml_tag} namelist.echam
else
mecho "Unknown scenario or scenario not set. Please choose out of HIST, SCEN, PI-CTRL, RCP*, SSP* and PALEO"
mecho "and change values according to your needs."
......
......@@ -523,15 +523,13 @@ fesom_cleanup()
if files=$(ls *_fesom_${START_DATE_fesom}.nc 2>/dev/null); then
for file in $files; do
# Christoper Danek
#
# fix filenames that get moved from work/ to outdata/fesom/
#
# fesoms output filenames are constructed in fesom_cpl/g_meandata.F90:
# CD: fix filenames that get moved from work/ to outdata/fesom/
# fesom 1.4 filenames are constructed in fesom_cpl/g_meandata.F90:
# filename=trim(ResultPath)//this%name//'_'//runid//'_'//firstdate//'.nc'
# where 'this%name' is the variable name according to the yaml list
# substream=`echo $file | cut -d_ -f1` # here, e.g. virtual_salt is not considered
substream=${file%$"_fesom"*} # actually, 'fesom' is 'runid' from namelist.config
# (actually, 'fesom' is 'runid' from namelist.config but always set to 'fesom')
# substream=`echo $file | cut -d_ -f1` does not consider 'virtual_salt' due to the '_'
substream=${file%$"_fesom"*}
if output_files=$(ls ${substream}_fesom_*.nc 2>/dev/null); then
for ofile in $output_files; do
add_to ${WORK_DIR}/${ofile} ${EXP_ID}_${ofile} output
......
......@@ -18,6 +18,9 @@ jsbach_set_defaults()
if [[ "x${DYNVEG_jsbach}" = "xdynveg" ]]; then
jsbach_nc="_dynveg"
jsbach_message="\t\tRunning WITH dynamic vegetation"
# CD: workaround for wrong filenames in pool
jsbach_nc=""
jsbach_message="${jsbach_message}; CD: had to set \$jsbach_nc manually to \"\""
else
jsbach_nc="_no-dynveg"
jsbach_message="\t\tRunning WITHOUT dynamic vegetation"
......@@ -103,6 +106,12 @@ jsbach_prepare_init()
;;
esac
;;
mpiesm )
case $echam_VERSION in
6.3.02p4 )
;;
esac
;;
mpiesm )
case $echam_VERSION in
6.3.02p4 )
......@@ -175,148 +184,163 @@ jsbach_prepare_init()
fi
}
jsbach_modify_forcing()
{
jsbach_modify_forcing() {
if [[ ${jsbach_VERSION} = "3.20" ]]; then
case ${RES_jsbach} in
echo -e "\t\tRES_jsbach=${RES_jsbach}"
echo -e "\t\tSCENARIO_jsbach=${SCENARIO_jsbach}"
echo -e "\t\tFORCING_DIR_jsbach=${FORCING_DIR_jsbach}"
case ${RES_jsbach} in
T63 | T31 )
case $SCENARIO_jsbach in
PALEO | 1850 | PI-CTRL* | 4CO2 | 1percCO2 | SCEN | HIST)
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr))
do
echo " jsbach: Setting year to $yr in forcing file Ndepo.${yr}.nc"
cdo setyear,${yr} ${FORCING_DIR_jsbach}/Ndepo.${yr}.nc ${FORCING_DIR_jsbach}/foo${yr}.nc >> $dumpfile 2>&1
mv ${FORCING_DIR_jsbach}/foo${yr}.nc ${FORCING_DIR_jsbach}/Ndepo.${yr}.nc >> $dumpfile
done
case $SCENARIO_jsbach in
PALEO | 1850 | PI-CTRL* | 4CO2 | 1percCO2 | SCEN)
for ((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr)) do
echo " jsbach: Setting year to $yr in forcing file Ndepo.${yr}.nc"
cdo setyear,${yr} ${FORCING_DIR_jsbach}/Ndepo.${yr}.nc ${FORCING_DIR_jsbach}/foo${yr}.nc >> $dumpfile 2>&1
mv ${FORCING_DIR_jsbach}/foo${yr}.nc ${FORCING_DIR_jsbach}/Ndepo.${yr}.nc >> $dumpfile
done
;;
HIST)
for ((yr = YR0_jsbach; yr <= YRN_jsbach + 0; ++yr)) do
echo -e "\t\tjsbach: selecting year $yr from ${FORCING_DIR_jsbach}/Ndepo.${yr}.nc"
cdo selyear,${yr} ${FORCING_DIR_jsbach}/Ndepo.${yr}.nc ${FORCING_DIR_jsbach}/foo${yr}.nc >> $dumpfile 2>&1
mv ${FORCING_DIR_jsbach}/foo${yr}.nc ${FORCING_DIR_jsbach}/Ndepo.${yr}.nc >> $dumpfile
done
;;
1950)
echo "Nothing to do in jsbach_modify_forcing for scenario $SCENARIO_jsbach"
echo "Nothing to do in jsbach_modify_forcing for scenario $SCENARIO_jsbach"
;;
*)
echo "Unknown scenario $SCENARIO_jsbach"
exit 42
;;
esac
esac
*)
echo "Unknown scenario $SCENARIO_jsbach"
exit 42
;;
esac
;;
esac
fi
}
jsbach_prepare_forcing()
{
if [[ "x${iterative_coupling}" == "x1" ]]; then
if [ -f "${CONFIG_DIR_jsbach}/jsbach_init_override.dat" ]; then
echo "PG: Sourcing ${CONFIG_DIR_jsbach}/jsbach_init_override.dat"
cat ${CONFIG_DIR_jsbach}/jsbach_init_override.dat
source ${CONFIG_DIR_jsbach}/jsbach_init_override.dat
fi
if [[ "x${iterative_coupling}" == "x1" ]]; then
if [ -f "${CONFIG_DIR_jsbach}/jsbach_init_override.dat" ]; then
echo "PG: Sourcing ${CONFIG_DIR_jsbach}/jsbach_init_override.dat"
cat ${CONFIG_DIR_jsbach}/jsbach_init_override.dat
source ${CONFIG_DIR_jsbach}/jsbach_init_override.dat
fi
if [ ! -z $LAND_BOUNDARY_CONDITIONS_jsbach ]; then
add_to ${LAND_BOUNDARY_CONDITIONS_jsbach} jsbach.nc
else
case $SCENARIO_jsbach in
fi
if [ ! -z $LAND_BOUNDARY_CONDITIONS_jsbach ]; then
add_to ${LAND_BOUNDARY_CONDITIONS_jsbach} jsbach.nc
else
case $SCENARIO_jsbach in
1850 | RCP* | PI-CTRL* | PALEO | 4CO2 | 1percCO2)
add_to ${ADJ_INPUT_DIR_jsbach}/jsbach_${RES_jsbach}${OCERES_jsbach}_11tiles_5layers_1850${jsbach_nc}.nc \
jsbach.nc
jsbach.nc
;;
1950 | SCEN)
add_to ${ADJ_INPUT_DIR_jsbach}/jsbach_${RES_jsbach}${OCERES_jsbach}_11tiles_5layers_1976${jsbach_nc}.nc \
jsbach.nc
jsbach.nc
;;
HIST)
add_to ${ADJ_INPUT_DIR_jsbach}/jsbach_${RES_jsbach}${OCERES_jsbach}_11tiles_5layers_1850${jsbach_nc}.nc \
jsbach.nc
for ((yr = YR0_jsbach + -1; yr <= YRN_jsbach + 2; ++yr)); do
if [ $yr -lt 1850 ]; then
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_1850.nc cover_fract.${yr}.nc
else
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
fi
done
;;
SSP*)
add_to ${ADJ_INPUT_DIR_jsbach}/jsbach_${RES_jsbach}${OCERES_jsbach}_11tiles_5layers_1850${jsbach_nc}.nc \
jsbach.nc
for ((yr = YR0_jsbach + -1; yr <= YRN_jsbach + 2; ++yr)); do
if [ $yr -lt 2015 ]; then
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
else
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/ssp${SCEN_NB_jsbach}/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
fi
done
;;
*)
echo "Unknown scenario $SCENARIO_jsbach in jsbach"
exit 42
;;
HIST)
add_to ${ADJ_INPUT_DIR_jsbach}/jsbach_${RES_jsbach}${OCERES_jsbach}_11tiles_5layers_1850${jsbach_nc}.nc \
jsbach.nc
for ((yr = YR0_jsbach + -1; yr <= YRN_jsbach + 2; ++yr)); do
if [ $yr -lt 1850 ]; then
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_1850.nc cover_fract.${yr}.nc
else
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
fi
done
;;
SSP*)
add_to ${ADJ_INPUT_DIR_jsbach}/jsbach_${RES_jsbach}${OCERES_jsbach}_11tiles_5layers_1850${jsbach_nc}.nc \
jsbach.nc
for ((yr = YR0_jsbach + -1; yr <= YRN_jsbach + 2; ++yr)); do
if [ $yr -lt 2015 ]; then
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
else
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/ssp${SCEN_NB_jsbach}/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
fi
done
;;
*)
echo "Unknown scenario $SCENARIO_jsbach in jsbach"
exit 42
;;
esac
fi
case ${RES_jsbach} in
T63 | T31 )
if [[ ${jsbach_VERSION} = "3.20" ]]; then
case ${RES_jsbach} in
T63 | T31 )
if [[ ${jsbach_VERSION} = "3.20" ]]; then
# PG: New input files needed for "spitfire" module
# These appear to be independent of the scenario
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/spitfire/a_nd_${RES_jsbach}.nc \
a_nd_file.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/spitfire/LISOTD_HRMC_V2.2_${RES_jsbach}.nc \
lightning.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/spitfire/population_density_HYDE_${RES_jsbach}.nc \
population_density.nc
fi
case $SCENARIO_jsbach in
PALEO | 1850 | PI-CTRL* | 4CO2 | 1percCO2)
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr))
do
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_1850.nc \
landuseHarvest.${yr}.nc
if [[ ${jsbach_VERSION} = "3.20" ]]; then
# PG: New for mpiesm-1.2.01, nitrogren input files
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/${RES_jsbach}_ndepo_CMIP_NCAR_CCMI-1-0_gr_185001-185012-clim.nc \
Ndepo.${yr}.nc
fi
done
;;
1950 )
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr))
do
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_${SCENARIO_jsbach}.nc \
landuseHarvest.${yr}.nc
done
;;
SCEN | RCP* | HIST)
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr))
do
if [ $yr -le 1849 ] ; then
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_1850.nc \
landuseHarvest.$yr.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_1850.nc cover_fract.${yr}.nc
elif [ $yr -le 2015 ] ; then
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_$yr.nc \
landuseHarvest.$yr.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
else
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/rcp${SCEN_NB_jsbach}/LUH_harvest_${RES_jsbach}_rcp${SCEN_NB_jsbach}_$yr.nc \
landuseHarvest.$yr.nc
fi
done
;;
SSP*)
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr))
do
if [ $yr -lt 2015 ] ; then
# add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_${yr}.nc \
# landuseHarvest.$yr.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
else
# add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/ssp${SCEN_NB_jsbach}_${DYNVEG_jsbach}/LUH_harvest_${RES_jsbach}_rcp${SCEN_NB_jsbach}_$yr.nc \
# landuseHarvest.$yr.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/ssp${SCEN_NB_jsbach}/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
fi
done
;;
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/spitfire/a_nd_${RES_jsbach}.nc \
a_nd_file.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/spitfire/LISOTD_HRMC_V2.2_${RES_jsbach}.nc \
lightning.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/spitfire/population_density_HYDE_${RES_jsbach}.nc \
population_density.nc
fi
case $SCENARIO_jsbach in
PALEO | 1850 | PI-CTRL* | 4CO2 | 1percCO2)
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr)) do
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_1850.nc \
landuseHarvest.${yr}.nc
if [[ ${jsbach_VERSION} = "3.20" ]]; then
# PG: New for mpiesm-1.2.01, nitrogren input files
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/${RES_jsbach}_ndepo_CMIP_NCAR_CCMI-1-0_gr_185001-185012-clim.nc \
Ndepo.${yr}.nc
fi
done
;;
1950)
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr))
do
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_${SCENARIO_jsbach}.nc \
landuseHarvest.${yr}.nc
done
;;
SCEN | RCP* | HIST)
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr))
do
if [ $yr -le 1849 ] ; then
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_1850.nc \
landuseHarvest.$yr.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_1850.nc cover_fract.${yr}.nc
elif [ $yr -le 2015 ] ; then
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_$yr.nc \
landuseHarvest.$yr.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
else
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/rcp${SCEN_NB_jsbach}/LUH_harvest_${RES_jsbach}_rcp${SCEN_NB_jsbach}_$yr.nc \
landuseHarvest.$yr.nc
fi
done
if [[ ${jsbach_VERSION} = "3.20" ]]; then
# PG: New for mpiesm-1.2.01, nitrogren input files
for ((yr = YR0_jsbach; yr <= YRN_jsbach + 0; ++yr)) do
echo -e "\t\tjsbach: selecting year $yr from ${jsbach_INPUT_DIR}/${RES_jsbach}/${RES_jsbach}_ndepo_CMIP_NCAR_CCMI-1-0_gr_185001-201412.nc"
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/${RES_jsbach}_ndepo_CMIP_NCAR_CCMI-1-0_gr_185001-201412.nc \
Ndepo.${yr}.nc
done
fi
;;
SSP*)
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr))
do
if [ $yr -lt 2015 ] ; then
# add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/${hist_fold}/LUH_harvest_${RES_jsbach}_${yr}.nc \
# landuseHarvest.$yr.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
else
# add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/ssp${SCEN_NB_jsbach}_${DYNVEG_jsbach}/LUH_harvest_${RES_jsbach}_rcp${SCEN_NB_jsbach}_$yr.nc \
# landuseHarvest.$yr.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/land_use_maps/ssp${SCEN_NB_jsbach}/cover_fract_${RES_jsbach}_11tiles_${yr}.nc cover_fract.${yr}.nc
fi
done
;;
*)
echo "Unknown scenario $SCENARIO_jsbach"
exit 42
......@@ -326,8 +350,8 @@ jsbach_prepare_forcing()
esac
for((yr = YR0_jsbach + -2; yr <= YRN_jsbach + 0; ++yr))
do
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/no_LUH_transitions_${RES_jsbach}.nc \
landuseTransitions.${yr}.nc
add_to ${jsbach_INPUT_DIR}/${RES_jsbach}/New_Hampshire_LCC/no_LUH_transitions_${RES_jsbach}.nc \
landuseTransitions.${yr}.nc
done
echo " Taking jsbach forcing data from jsbach_INPUT_DIR=${jsbach_INPUT_DIR}"
......@@ -376,7 +400,7 @@ jsbach_prepare_config()
if [[ "x${RUN_NUMBER_jsbach}" = "x1" ]]; then
echo " Taking jsbach namelist from jsbach_NAMELIST_DIR=${jsbach_NAMELIST_DIR}"
fi
add_to ${jsbach_NAMELIST_DIR}/namelist.jsbach namelist.jsbach
add_to ${jsbach_NAMELIST_DIR}/namelist.jsbach${CMIP6_nml_tag} namelist.jsbach
else
mecho "Unknown scenario or scenario not set. Please choose out of HIST, SCEN, PI-CTRL, RCP*, SSP* and PALEO"
mecho "and change values according to your needs."
......@@ -544,14 +568,7 @@ jsbach_post_process()
tar --append $tar_compress_program ${TAR_DIR}/${EXP_ID}_jsbach_raw_restarts_${END_DATE_jsbach}.tar.gz restart_${EXP_ID}_*_${END_DATE_jsbach}.nc
fi
cd ${DATA_DIR_jsbach}
if ! [[ -e ${EXP_ID}_jsbach_raw_output_${END_DATE_jsbach}.tar.gz ]]; then
# --use-compress-program=pigz -cvf
tar --create $tar_compress_program ${TAR_DIR}/${EXP_ID}_jsbach_raw_output_${END_DATE_jsbach}.tar.gz ${EXP_ID}_jsbach_*_${END_DATE_jsbach}.nc
else
tar --append $tar_compress_program ${TAR_DIR}/${EXP_ID}_jsbach_raw_output_${END_DATE_jsbach}.tar.gz ${EXP_ID}_jsbach_*_${END_DATE_jsbach}.nc
fi
fi
fi
ARCHIVE_TAR_OUTPUT_jsbach=${ARCHIVE_TAR_OUTPUT_jsbach:-0}
if [[ "x${ARCHIVE_TAR_OUTPUT_jsbach}" = "x1" ]]; then
......
&parctl
NPROCA = 24
NPROCB = 48
/
&runctl
NPROMA = 72
OUT_EXPNAME = "TST"
OUT_DATAPATH = "/work/bk0988/awicm/a270062/fesom_echam/glob//cpl_work_test/"
OUT_FILETYPE = 2
OUT_ZTYPE = 0
RERUN_FILETYPE = 4
LRESUME = T
LTIMER = .false.
LMIDATM = .true.
LFRACTIONAL_MASK=.false.
! LHD = .true.
LCOUPLE = .true.
LCOUPLE_CO2 = .true.
! LIPCC = .true.
LAMIP = .false.
DT_START = 1850,01,01,0,0,0
DT_STOP = 2402,1,1,0,0,0
DELTA_TIME = 200
PUTOCEAN = 1,'days','last',0
GETOCEAN = 1,'days','last',0
PUTRERUN = 1,'years','first',0
PUTDATA = 3,'hours','last',0
default_output = .true.
TRIGFILES = 12,'months','exact',0
LDEBUGEV = .false.
! lmeltpond = .false.
ltdiag = .true.
/
&cfdiagctl
locfdiag = .true.
/
&submdiagctl
vphysc_lpost = .false.
/
&dynctl
ENSTDIF = 1.
/
&submodelctl
lmethox = .true.
lco2 = .false.
/
&radctl
iaero = 8
io3 = 4
isolrad = 1
ico2 = 4
ighg = 1
icfc = 4
ich4 = 4
in2o = 4
ldiag_aeropt = .true.
/
&co2ctl
LCO2_SCENARIO = .true.
/
! CD: the rest of this namelist was added by hand from
! /pf/a/a270061/mpiesm-1.2.01p1_CMIP6/mpiesm-1.2.01p1/experiments/LIG_127k_CMIP6_output/scripts/LIG_127k_CMIP6_output.run
! /work/ba1066/a270061/mpiesm-1.2.01p1_CMIP6/mpiesm-1.2.01p1/experiments/LIG_127k_CMIP6_output/outdata/echam6
&mvstreamctl ! 'default spm glm g3bm'
filetag = 'echam'
source = 'sp', 'gl', 'g3b'
/
&mvstreamctl ! 'default co2m'
filetag = 'co2'
source = 'co2'
/
&mvstreamctl ! 'default surfm'
filetag = 'surf'
source = 'surf'
/
&mvstreamctl ! 'default tracerm'
filetag = 'tracer'
source = 'tracer'
/
&mvstreamctl ! 'default jsbachm'
filetag = 'jsbach'
source = 'jsbach'
/
&mvstreamctl ! 'default landm'
filetag = 'land'
source = 'land'
/
&mvstreamctl ! 'default vegm'
filetag = 'veg'
source = 'veg'
/
&mvstreamctl ! 'default yassom'
filetag = 'yasso'
source = 'yasso'
/
&set_stream ! 'default jsbach'
stream = 'jsbach'
lpost = 0
/
&set_stream ! 'default land'
stream = 'land'
lpost = 0
/
&set_stream ! 'default veg'
stream = 'veg'
lpost = 0
/
&set_stream ! 'default yasso'
stream = 'yasso'
lpost = 0
/
&set_stream ! 'scm'
stream = 'scm'
lpost = 0
/
&mvstreamctl ! 'cmip6 accw3hr'
filetag = 'echam3hr'
source = 'accw'
target = '*3hr'
interval = 3, 'hours', 'last', 0
variables = 'drainage', 'glacier_melt', 'runoff', 'snow_melt'
/
&set_stream ! 'cmip6 accw3hr'
stream = 'accw3hr'
lrerun = 0
/
&mvstreamctl ! 'cmip6 aeroptday'
filetag = 'echamday'
source = 'aeropt'
target = '*day'
interval = 1, 'days', 'last', 0
variables = 'tau_aero_550'
/
&set_stream ! 'cmip6 aeroptday'
stream = 'aeroptday'
lrerun = 0
/
&mvstreamctl ! 'cmip6 aeroptmon'
filetag = 'echammon'
source = 'aeropt'
target = '*mon'
interval = 1, 'months', 'last', 0
variables = 'omega_aero_550', 'tau_aero_550'
/
&set_stream ! 'cmip6 aeroptmon'
stream = 'aeroptmon'
lrerun = 0
/
&mvstreamctl ! 'cmip6 cfdiag3hr'
source = 'cfdiag'
target = '*3hr'
interval = 3, 'hours', 'last', 0
variables = 'irsucs:inst>=6'
/
&set_stream ! 'cmip6 cfdiag3hr'
stream = 'cfdiag3hr'
lrerun = 0
/
&set_stream_element ! 'cmip6 cfdiag3hr irsucs'
stream = 'cfdiag3hr'
name = 'irsucs'
code = 6
lpost = 1
/
&mvstreamctl ! 'cmip6 cfdiagday'
source = 'cfdiag'
target = '*day'
interval = 1, 'days', 'last', 0
variables = 'mc'
/
&set_stream ! 'cmip6 cfdiagday'
stream = 'cfdiagday'
lrerun = 0
/
&mvstreamctl ! 'cmip6 cfdiagmon'
source = 'cfdiag'
target = '*mon'
interval = 1, 'months', 'last', 0
variables = 'dmc', 'irld>rld=3', 'irldcs>rldcs=7', 'irlu>rlu=1',
'irlucs>rlucs=5', 'irsd>rsd=4', 'irsu>rsu=2', 'irsucs>rsucs=6',
'mc', 'mcd', 'mcu', 'smc'
/
&set_stream ! 'cmip6 cfdiagmon'
stream = 'cfdiagmon'
lrerun = 0
/
&mvstreamctl ! 'cmip6 co2mon'
filetag = 'echammon'
source = 'co2'
target = '*mon'
interval = 1, 'months', 'last', 0
variables = 'co2_burden', 'co2_flx_anthro', 'co2_flx_harvest', 'co2_flux',
'co2_flx_land', 'co2_flx_lcc', 'co2_flx_ocean'
/
&set_stream ! 'cmip6 co2mon'
stream = 'co2mon'
lrerun = 0
/
&mvstreamctl ! 'cmip6 echam3hr'
filetag = 'echam3hr'
source = 'sp', 'gl', 'g3b'
target = '*3hr'
interval = 3, 'hours', 'last', 0
variables = 'aclcov_na>aclcov=164', 'ahfl_na>ahfl=147', 'ahfs_na>ahfs=146',
'aprc_na>aprc=143', 'aprl_na>aprl=142', 'aprs_na>aprs=144',
'aps:inst', 'geosp', 'lsp', 'lsp:inst', 'q2m:inst', 'q:inst',
'rsdscs_na>rsdscs=50', 'rsuscs_na:inst', 'rsuscs_na>rsuscs=51',
'sd:inst', 'srad0d_na>srad0d=184', 'srad0u_na>srad0u=203',
'srads_na>srads=176', 'sradsu_na>sradsu=204',
'sraf0_na>sraf0=187', 'st', 'st:inst', 'svo:inst', 'temp2:inst',
'trad0_na>trad0=179', 'trads_na>trads=177',
'tradsu_na>tradsu=205', 'traf0_na>traf0=188',
'trafs_na>trafs=186', 'u10', 'u10:inst', 'v10', 'v10:inst',
'wind10_na>wind10=171'
/
&set_stream ! 'cmip6 sp3hr'
stream = 'sp3hr'
lrerun = 0
/
&set_stream ! 'cmip6 gl3hr'