# # Setup for pre-industrial control run at T63L??/ORCA05. # using pre-initialized carbon pools for JSBACH, and ECHAM6 restartt # files to avoid an initial shock to the system. The ocean is # initialized from Levitus, and the sea ice uses an initalization # file from a present day, forced ocean run from Klaus Getzlaff, TM # the jsbach carbon pools (Cpools.nc) are from a 8500 year holocene # run done at MPI # # Model information VERSION_ = Hist_scn_CMIP6-foci-mops_daily_output_restart_from_CC009.config foci1.20.0 2020-09-29 Chia-Te Chien MODEL_ID = FOCI1.20.0 EXP_TYPE = foci #EXP_OPTIONS = echam_daily_output_with_tdiag # non-ESM model EXP_OPTIONS = echam_mops_daily_output_with_tdiag # ESM model SCENARIO = ssp585 # if year < 2015, historical is used NITERATION = 1 #nb of year done in one job, important to use for HLRN4 #JVD # save veg stream at it's original resolution #LAND_KEEP_ORIG = veg # CPU Layout, if commented a default layout for the current machine (see # standard_environments/.config will be used) ATM_DOMAINS_X = 24 # 24 is also the default for nprocb ATM_DOMAINS_Y = 24 # 12 is also the default for nproca #JVD ATM_NPROCIO = 12 #for HLRN4, 24 for HLRN3 #JVD OCE_DOMAINS_X = 28 # 3 is the default #JVD OCE_DOMAINS_Y = 24 #for HLRN4, 24 for HLRN3 # 24 is currently the default # Atmosphere setup ATM_RES = T63 ATM_LEVELS = 95 ATM_GRID_SIZE_X = 192 ATM_GRID_SIZE_Y = 96 ATM_TIME_STEP = 450 ATM_BLOCK_SIZE = 32 #for HLRN4, 72 for HLRN3 # nproma #JVD ATM_THREADS = 1 # uncomment below two lines if run in ESM mode ATMO_RESTART_TAGS = echam co2 accw tracer ATMO_FILE_TAGS = echam co2 accw tracer tdiag # Ocean setup OCE_RES = ORCA05 OCE_LEVELS = 46 OCE_GRID_SIZE_X = 722 OCE_GRID_SIZE_Y = 511 OCE_TIME_STEP = 2160 #JVD OCE_THREADS = 1 OCE_USE_TRACER = true OCE_CONF = ORCA05_LIM2_FOCI_MOPS OCEAN_EXE = nemo_ORCA05_LIM2_FOCI_MOPS.x OCEAN_FILE_TAGS = "grid_T grid_U grid_V grid_W icemod ptrc_T diad_T scalar" # Coupling setup OCE2ATM_INTERVAL = 10800 ATM2OCE_INTERVAL = $OCE2ATM_INTERVAL # Initial and final date INITIAL_DATE = 2015-01-01 FINAL_DATE = 2015-12-31 # example settings for restarting from another experiment, they are used if you start the model # using .run_start # to use your own restart data, replace PARENT_DIR with the correct path, data is expected # in subdirectories nemo, jsbach, echam, oasis3mct # in this case use .run_start to start the model PARENT_EXP_ID = FOCI1.20.0-CC105_RCP_ESM_spinup2099 #PARENT_EXP_ID = FOCI1.20.0-CC009_piControl_ESM PARENT_DATE = 20141231 PARENT_ORIG_DATE = $PARENT_DATE OCEAN_PARENT_NEMOTIMESTEP = 02410600 #CSW no need to set PARENT_DIR, it defaults to that location #PARENT_DIR = /scratch/usr/shkifmsw/foci_input2/FOCI_RESTART/FOCI1.20.0-JD100_spinup #PARENT_DIR = /scratch/usr/shkchien/foci_input2/FOCI_RESTART/FOCI1.20.0-CC009_piControl_ESM PARENT_DIR = /scratch/usr/shkchien/foci_input2/FOCI_RESTART/FOCI1.20.0-CC105_RCP_ESM_spinup2099 #PARENT_DIR = /scratch/usr/shkchien/models/foci1.20.0/experiments/FOCI1.20.0-CC104_RCP_ESM_spinup2089/restart #CSW: uncomment below if run in ESM mode ATMO_PARENT_DIR = $PARENT_DIR/echam6 #PARENT_DIR = $INPUT_ROOT/FOCI_RESTART/FOCI1.3-SW038_MOPSRESTART # on HLRN this resolves to $WORK2/foci_input/FOCI/restart/$PARENT_EXP_ID # data is expected in subdir nemo, jsbach, echam, ... #xios mode, if commented the default setting for the current machine (see # standard_environments/.config will be used) #XIOS_DETACHED = false #extra XIOS server (true) or not (false), false under HLRN4, true under HLRN3 # Model failure work-around - enable this if echam hangs and returns "cuadjtq (1): lookup table overflow" #JVD ATMO_DISTURBED_YEARS = 1872 1879 #space separated year(s) where the model fails #JVD #ATMO_DISTURBANCE = 1.00001 #JVD [namelists] [[namcouple_foci]] nlogprt = 0 # possible values 0(silent),1,2,5,10,12,15,20,30(verbose) xios = eval( 'xios.x' if '${XIOS_DETACHED}'=='true' else '') nmod = eval( 3 if '${XIOS_DETACHED}'=='true' else 2) [[namelist.echam]] [[[runctl]]] lcouple_co2 = .true. default_output = .false. ltdiag = .true. #default_output and ltdiag are new for CO2 emission files [[[radctl]]] iaero = 8 ighg = 1 io3 = 4 isolrad = 1 ico2 = 1 #CSW: use ico2 = 1 if in ESM mode, = 4 in non-ESM model ich4 = 4 in2o = 4 icfc = 4 # uncomment if you run in ESM mode [[[submodelctl]]] lco2 = .true. lmethox = .true. # lmethox is new for ESM mode scenarios [[[co2ctl]]] lco2_emis = .true. lco2_2perc = .true. # There two are new for ESM mode scenarios [[namelist.jsbach]] [[[jsbach_ctl]]] use_dynveg = true use_disturbance = true with_nitrogen = true lcc_forcing_type = transitions lcc_scheme = 2 test_cconservation = true test_nconservation = true [[[dynveg_ctl]]] dynveg_feedback = .true. [[[hydrology_ctl]]] ldebughd = false diag_water_budget = true [[[cbalance_ctl]]] read_ndepo = true [[[disturbance_ctl]]] fire_name = thonicke fuel_name = thonicke lburn_pasture = true [[[input_ctl lightning_frq]]] var_name = lightning_frq var_file_name = HRMC_COM_FR file_name = lightning.nc action_cycle = all [[[input_ctl population_density]]] # only read first record of population_density file (i.e. year 1850) var_name = population_density # dt_update = 0 # init_rec = 1 [[namelist_top_cfg]] ln_trcrad = .false. # artificially correct negative concentrations (T) or not (F) !DE change 25.04.2018 [files] [[echam]] [[[multi_year]]] .offsets = 0 .set label = _$scn .set hist_or_scn = $scn greenhouse_gases.nc = greenhouse_$${hist_or_scn}.nc MAC-SP.nc = MAC-SP$${label}.nc [[[mapped]]] [[[[aero]]]] .sub_dir = aero # echam.radctl.iaero = 8 needs 1850 aerosols aero_fine_1849.nc = $${ATM_RES}_aeropt_kinne_sw_b14_fin_1865.nc aero_coarse_1849.nc = $${ATM_RES}_aeropt_kinne_sw_b14_coa.nc aero_farir_1849.nc = $${ATM_RES}_aeropt_kinne_lw_b16_coa.nc aero_fine_1850.nc = $${ATM_RES}_aeropt_kinne_sw_b14_fin_1865.nc aero_coarse_1850.nc = $${ATM_RES}_aeropt_kinne_sw_b14_coa.nc aero_farir_1850.nc = $${ATM_RES}_aeropt_kinne_lw_b16_coa.nc aero_fine_1851.nc = $${ATM_RES}_aeropt_kinne_sw_b14_fin_1865.nc aero_coarse_1851.nc = $${ATM_RES}_aeropt_kinne_sw_b14_coa.nc aero_farir_1851.nc = $${ATM_RES}_aeropt_kinne_lw_b16_coa.nc # CSW: carbon emissions for historical simulations with carbon cycle [[[[carbonemissions]]]] .sub_dir = carbon_emissions carbon_emissions.nc = $${ATM_RES}_carbon_emissions_${scn}.nc [[[varying]]] [[[[solar_irradiance]]]] .sub_dir = solar_irradiance .set label = 1850, 1850:%year swflux_%year.nc = swflux_14band_$${label}.nc [[[[mapped]]]] [[[[[aero]]]]] aero_fine_%year.nc = aero_coarse_%year.nc = aero_farir_%year.nc = [[[[[ozone]]]]] .set historical_or_scn = historical, 2015:$scn .set label = 1850, 1850:%year, 2100:2099 ozon%year = $${ATM_RES}_ozone_$${historical_or_scn}_$${label}.nc [[[[[volcano_aerosols]]]]] .sub_dir = volcano_aerosols .set label = 1850, 1850:%year, 2025:2024 strat_aerosol_ir_%year.nc = strat_aerosol_ir_$${ATM_RES}_$${label}.nc strat_aerosol_sw_%year.nc = strat_aerosol_sw_$${ATM_RES}_$${label}.nc [[jsbach]] [[[mapped]]] jsbach.nc = jsbach_$${ATM_RES}$${OCE_RES}${.fractional}_${LAND_TILES}tiles${.layers}_${LAND_REF_YEAR}_dynveg.nc [[[[landuse]]]] .offsets = -1, 0 .sub_dir = New_Hampshire_LCC .set hist_or_scn = hist, 2015:$scn landuseTransitions.%year.nc = $${hist_or_scn}_dynveg/LUH_scaled_transitions_$${ATM_RES}_%year.nc landuseHarvest.%year.nc = $${hist_or_scn}_harvest/LUH_harvest_$${ATM_RES}_%year.nc [[[[nitrogen]]]] .offsets = -1, 0 .set label = 1850, 1850:%year .set hist_or_scn = _CMIP_NCAR_CCMI-1-0_gr_185001-201412, 2015:.input4MIPs_surfaceFluxes_ScenarioMIP-NCAR-CCMI-$scn .method = 'cdo selyear,$${label}' Ndepo.%year.nc = $${ATM_RES}_ndepo$${hist_or_scn}.nc [[[[spitfire]]]] .sub_dir = spitfire a_nd_file.nc = a_nd_$${ATM_RES}.nc lightning.nc = LISOTD_HRMC_V2.2_$${ATM_RES}.nc population_density.nc = population_density_HYDE_$${ATM_RES}.nc [jobs] project = shk00056 #relpace shk??? with your project id and uncomment #JVD # [[run, run_init, run_start]] #USE on HLRN3 # time_limit = 02:00:00 #USE on HLRN3 # [[run]] #USE on HLRN4 # time_limit = 12:00:00 #USE on HLRN4 [[run]] time_limit = 8:00:00 #USE on HLRN4 [[[files]]] [[[[echam]]]] [[[[[restart]]]]] .base_dir = $ATMO_RESTART_DIR .method = cp -pv restart_${EXP_ID}_tracer.nc = restart_$${EXP_ID}_tracer_$${prevdate}.nc [[run_start, run_init]] #USE on HLRN4 time_limit = 01:00:00 #USE on HLRN4 #JVD [[run_start]] [[[namelists]]] [[[[namelist.jsbach]]]] [[[[[cbalance_ctl]]]]] read_npools = false init_npools = false [[[files]]] [[[[echam]]]] [[[[[restart]]]]] .base_dir = $ATMO_PARENT_DIR .method = cp -pv restart_${EXP_ID}_tracer.nc = restart_$${ATMO_PARENT_EXP_ID}_tracer_$${ATMO_PARENT_DATE}.nc [[run_init]] [[[namelists]]] # [[[[namelist_cfg]]]] # ln_rstart = .false. # ln_tsd_init = .true. # Initialisation of ocean T & S with T &S input data (T) or not (F) # ln_limini = .true. # don't use ice initialization in restart case [[[[namelist.jsbach]]]] [[[[[cbalance_ctl]]]]] read_cpools = .true. # read yasso carbon pools from Cpools.nc read_ndepo = .true. # read initial nitrogen dep rates from Ndepo.nc [[[[[dynveg_ctl]]]]] read_fpc = .true. [[[files]]] [[[[jsbach]]]] [[[[[mapped]]]]] Cpools.nc = cpools_$${ATM_RES}$${OCE_RES}_pdv0002_51251231.nc Ndepo.nc = cpools_$${ATM_RES}$${OCE_RES}_pdv0002_51251231.nc fpc.nc = cpools_$${ATM_RES}$${OCE_RES}_pdv0002_51251231.nc [[[[oasis]]]] [[[[[init]]]]] #.base_dir = $COUPLER_INPUT_DIR/$${ATM_RES}_$${OCE_RES} #.method = cp -pv flxatmos = flxatmos_FOCI1.4.1-KJH001_18501231 sstocean = sstocean_FOCI1.4.1-KJH001_18501231_CO2OCEAN_CO2TRA_FCO2_SOSALINE_CO2FLX_FF_CO2