Setting up a coupled climate run for pre-industrial times

The page follows the workflow of HowTo/CoupledDofoco and HowTo/SetupZoomEurope. For the commands and actions consult those pages. On this page the settings are specified.

The aim is to run an equilibrium simulation for the year 1750. Orchidee off-line line was used to simulate the surface area in 1750 and now we want to know how the climate looked liked for such a surface. The coupled run is zoomed over Europe, nudged at the edges of Europe uses the CWRR hydrology and all the DOFOCO features.

The boundary conditions

LMDZOR was installed twice to avoid any confusion between the boundary files of the pre-industrial and the present day simulations. Hence, one installation will be used for the pre-industrial runs (and will thus use the pre-industrial boundary conditions), the other installation is to be used for present day simulations. For both sets of simulations an AMIP set-up was chosen. This set-up comes with acceptable BC between 1870 and 2100 which allows us to focus on the climate and reduce issues with sea surface temperature, see ice and aerosol concentrations. So copy the an amip config.card to define the job that will create the BC

cd config/LMDZOR_v5.2/
cp EXPERIMENTS/LMDZ/CREATE_amip/config.card ./

The following changes to were made to the config.card:

(1) Add a suffix, i.e., -1750 to the job name so that the pre-industrial and present day boundary files will not interfere with each other

JobName= ELI-128x118x39-1750

(2) Set the begin date to the earliest year available. In this case we are limited by the SST and sea ice files to 1870. So we assume that the SST in 1870 represents pre-industrial conditions. Further we are planning 20 year long equilibrium simulations but we are not sure whether it will take 20 or 30 years so we simulate the BC for 30 years to be on the safe side.

DateBegin=1870-01-01
DateEnd=1899-12-31

There is no need to manually create a directory for this experiment. When you use ./util/ins_job the script will create a directory with a copy of the config.card as well as the COMP, DRIVER and PARAM folder. Note that many of the data files that are used to describe the boundary conditions can be found at

/ccc/work/cont003/dsm/p86ipsl/IGCM/INIT/ATM/LMDZOR/

Following, more settings need to be changed and/or confirmed before launching the job. The following changes to were made to the COMP/lmdz.card:

(1) Specify in the list of parameter files that the BC should be based on the predefined pre-industrial settings. These settings contain the TOA solar incoming radiation, atmospheric CO2, CH4 and N2O among other settings. Have a loop in PARAM/config.def_preind to get a better idea of these settings.

(${SUBMIT_DIR}/PARAM/config.def_preind     ,  config.def),    \

(2) In the DRIVER/lmdz.driver it is specified that in the absence of the Oasis coupler, which is the case for LMDZOR, the ocean is forced. Hence the need of SST and sea ice files for all years in the experiment.

(3) Most likely this is redundant at this stage but better safe than sorry so in PARAM/run.def set

VEGET=y

Now the ELI-128x118x39-1750 can be launched. It takes about 70 seconds per year so in this example it will take over 30 minutes to complete this 30-year job

ccc_msub Job_ELI-128x118x39-1750

The pre-industrial experiment

Where for the BC run, an LMD experiment was copied, we now have to copy a LMDZOR experiment. We haven't tried whether LMDZ/clim BC files can be combined with LMDZOR/amip experiments but to keep it simple we propose to be consistent and thus combine LMDZ/amip BC files with an LMDZOR/amip experiment

cp ./EXPERIMENTS/LMDZOR/clim/config.card .

Edit the config card. Change the JobName?, dates

../../util/ins_job 

This creates a new directory and copies some from elsewhere... our account is gen6328

cd JOBNAME

change PARAM/orchidee.def_CWRR to include the info from the offline orchidee runs (i.e., the run.def used offline) change PARAM/orchidee.def_CWRR: set river routing to "RIVER_ROUTING = n" set'DefSuffix?=CWRR' in COMP/orchidee.card

We want to use the new physics. Therefore, we need to change a few things in the input files.

# in PARAM/orchidee.def_CWRR
---> add 'evapot_corr=y'

# in COMP/lmdz.card
---> LMDZ_Physics=NPv3.1


# in PARAM/physiq.def_L39_NPv3.1
--->
iflag_thermals=17
iflag_pbl=11

# Make sure that stomate is on by looking in PARAM/orchidee.def_CWRR
STOMATE_OK_STOMATE = y
# and in config.card (if this line is missing, there will be no STOMATE)
SBG= ("", "")

If you are interested in a zoomed set-up follow the further instructions on HowTo/SetupZoomEurope. For the regular grid you are done and can try to launch the job:

ccc_msub Job_JOBNAME

Notice that I had to ask the hotline to be included in this group, since when I tried to launch the job it was refused