176 | | Within the 20-year coupled run the land surface will stay constant to that described by the restart file that will be used, i.e., 1750. Forest management and litter raking are still activated (else we will affect the C-pools). We will keep these settings constant to their values for the year 1750. We simply hard code the years for these maps as detailed below. Note that because we start from a restart file the [!InitialStateFiles] are not used. We have to define the litter and FM maps as [!BoundaryFiles] which are read only once during the year. Note that 1:12: means that the file will be copied to the working directory at the first integration step and then every 12 iterations until the simulation is finished. More information on the difference between !InitialStateFiles, !BoundaryFilesn and !SmoothFiles can be found at http://forge.ipsl.jussieu.fr/igcmg_doc/wiki/DocEsetup#COMPdirectory |
| 178 | Within the 20-year coupled run the land surface will stay constant to that described by the restart file that will be used, i.e., 2010. Forest management and litter raking are still activated (else we will affect the C-pools). We will keep these settings constant to their values for the year 2010. We simply hard code the years for these maps as detailed below. Note that because we start from a restart file the [!InitialStateFiles] are not used. We have to define the litter and FM maps as [!BoundaryFiles] which are read only once during the year. Note that 1:12: means that the file will be copied to the working directory at the first integration step and then every 12 iterations until the simulation is finished. More information on the difference between !InitialStateFiles, !BoundaryFilesn and !SmoothFiles can be found at http://forge.ipsl.jussieu.fr/igcmg_doc/wiki/DocEsetup#COMPdirectory |
| 179 | CAREFUL: the example here is for BCLF, when BC, BCF or BCL are to be set different fm_start maps need to be used. |
190 | | === PARAM/config.def_preind === |
| 193 | === Ground work for adjusting dates === |
| 194 | Because we have BC files for a present day experiment but also for pre-industrial experiment we changed the name of the BC folder by adding the suffic '-2010'. We have to make sure that LMDZOR finds the correct BC by adding this suffix. The ELI files go from 1981 to 2010 but we often don't want to use the 1981 conditions so in those cases we will hard code the date to 2010. |
| 195 | |
| 196 | ?????????????????????? |
| 197 | The nudging files |
| 198 | and fixing the year to 2010 rather than the date in the config.card (which was in turn adjusted to ensure that the ELI simulation could be run for the length of our experiment, i.e., 30 years): |
| 199 | Following, more settings need to be changed and/or confirmed before launching the job: |
| 200 | (1) We now have to set-up the job to create files for the boundary conditions (BC) for the years 1751 to 1780 but we have no driver files for these years so we use the earliest years available. That varies from data stream to data stream so we used a creative solution. We are limited, for example by the SST and sea ice files that only go back until 1870. So we assume that the SST in 1870 represents pre-industrial conditions. |
| 201 | |
| 202 | Open '''.../libIGCM/libIGCM_config/libIGCM_config.tsh''' and search for the phrase '# BEGIN: SHOULD GO IN A FUNCTION FROM libIGCM_date.ksh'. Then add: |
| 203 | {{{ |
| 204 | # BEGIN: SHOULD GO IN A FUNCTION FROM libIGCM_date.ksh |
| 205 | # Compute year_m1 and year_p1 (year minus 1Y and year plus 1Y) |
| 206 | year_m1=$(( year - 1 )) |
| 207 | year_p1=$(( year + 1 )) |
| 208 | # First year for SST 1870 = 1751 + 119 |
| 209 | year_sst=$((year + 119)) |
| 210 | # First year for O3 is 1855 = 1751 + 104 |
| 211 | year_o3=$((year + 104)) |
| 212 | # First year for nudged wind fields 1980 = 1751 + 229 |
| 213 | year_nudge=$((year + 229)) |
| 214 | # First year of pre-industrial forcing 1765 = 1751 + 14 |
| 215 | year_preind=$((year + 14)) |
| 216 | }}} |
| 217 | (2) These newly defined years now have to be used in COMP/lmdz.card. |
| 218 | ????????????????? |
| 219 | |
| 220 | === PARAM/config.def_actuel === |
237 | | Ideally wind fields for the year 1750 should be used for the pre-industrial run. However, ERAI data only start in 1958 so such files do not exist. Consequently, the 1980 to 2011 wind fields are used to nudge the pre-industrial simulation. the assumption we make here is that the winter patterns outside the nudged area have not changed dramatically between 1980 and 1750. For a strongly zoomed simulation over Europe that is probably an acceptable assumption. The years 1980 to 2011 have been prepared for the 128x118 zoomed grid and are stored at /ccc/work/cont003/dsm/p86ipsl/IGCM/BC/ATM/LMDZ/LMD128118/NUDGE_FILES/ERAI. If you are not using a zoomed grid that has been used previously, or you want to use the 128*118 grid for different years, you will have to prepare your own nudging files. ERAI nudging files for wind (u and v) can be interpolated to the zoomed grid using the method described here : [https://forge.ipsl.jussieu.fr/igcmg/wiki/platform/en/documentation/I_modeles/B_LMDZ#Interpolationofnudgefiles]. Note that you must be in the group subipsl to have the permission to access these files. Contact Anne Cozic for TGCC to be added to subipsl. |
238 | | |
239 | | We still need to tell the model which nudging files it should use. This requires a change to COMP/lmdz.card so that libIGCM knows to copy the files before the run. Note that here we substitute the argument ${year} by ${year+130} so we will run the simulation from 1750 to 1780 (as specified in the config.card) but use the wind fields from 1980 to 2000 during the 20 years of the equilibrium climate experiment. |
| 267 | The 1980 to 2011 wind fields are used to nudge the present day simulation. The years 1980 to 2011 have been prepared for the 128x118 zoomed grid and are stored at /ccc/work/cont003/dsm/p86ipsl/IGCM/BC/ATM/LMDZ/LMD128118/NUDGE_FILES/ERAI. If you are not using a zoomed grid that has been used previously, or you want to use the 128*118 grid for different years, you will have to prepare your own nudging files. ERAI nudging files for wind (u and v) can be interpolated to the zoomed grid using the method described here : [https://forge.ipsl.jussieu.fr/igcmg/wiki/platform/en/documentation/I_modeles/B_LMDZ#Interpolationofnudgefiles]. Note that you must be in the group subipsl to have the permission to access these files. Contact Anne Cozic for TGCC to be added to subipsl. |
| 268 | |
| 269 | We need to tell the model which nudging files it should use. This requires a change to COMP/lmdz.card so that libIGCM knows to copy the files before the run. |
255 | | ListNonDel= (${R_BC}/ATM/${config_UserChoices_TagName}/IPCC_AR5/HISTORIQUE/SOLARANDVOLCANOES_1610_2008.txt, SOLARANDVOLCANOES.txt),\ |
256 | | (${R_BC}/ATM/${config_UserChoices_TagName}/IPCC_AR5/HISTORIQUE/CO2_1765_2005.txt, CO2.txt),\ |
257 | | (${R_BC}/ATM/${config_UserChoices_TagName}/IPCC_AR5/HISTORIQUE/CH4_1765_2005.txt, CH4.txt),\ |
258 | | (${R_BC}/ATM/${config_UserChoices_TagName}/IPCC_AR5/HISTORIQUE/N2O_1765_2005.txt, N2O.txt),\ |
259 | | (${R_BC}/ATM/${config_UserChoices_TagName}/IPCC_AR5/HISTORIQUE/CFC11_1765_2005.txt, CFC11.txt),\ |
260 | | (${R_BC}/ATM/${config_UserChoices_TagName}/IPCC_AR5/HISTORIQUE/CFC12_1765_2005.txt, CFC12.txt),\ |
| 285 | ListNonDel= (/ccc/work/cont003/dofoco/dofoco/INPUT_FILES/ATM_COMP/SOLAR_ACTUAL.txt, SOLARANDVOLCANOES.txt),\ |
| 286 | (/ccc/work/cont003/dofoco/dofoco/INPUT_FILES/ATM_COMP/CO2_ACTUAL.txt, CO2.txt),\ |
| 287 | (/ccc/work/cont003/dofoco/dofoco/INPUT_FILES/ATM_COMP/CH4_ACTUAL.txt, CH4.txt),\ |
| 288 | (/ccc/work/cont003/dofoco/dofoco/INPUT_FILES/ATM_COMP/N2O_ACTUAL.txt, N2O.txt),\ |
| 289 | (/ccc/work/cont003/dofoco/dofoco/INPUT_FILES/ATM_COMP/CFC11_ACTUAL.txt, CFC11.txt),\ |
| 290 | (/ccc/work/cont003/dofoco/dofoco/INPUT_FILES/ATM_COMP/CFC12_ACTUAL.txt, CFC12.txt),\ |