Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

How to read and configure external CO2 data?

Ambition_hhh

Ambition_hhh
New Member
./create_newcase --case Myco2case --res f09_g16 --compset IHistClm50BgcCrop --mach myintel --compiler intel --run-unsupported

I created this case, I want the model to use my own file download of co2 (cmip6 download rcp2.6 data):

/home/hlh/cesm/inputdata/atm/datm7/CO2/co2_AERmon_CESM2_historical_r10i1p1f1_gn_200001-201412.nc

1. I changed it in the run_run.xml file:
<entry id="DATM_CO2_TSERIES" value="none">
<type>char</type>
<valid_values>none,20tr,20tr.latbnd,rcp2.6,rcp4.5,rcp6.0,rcp8.5</valid_values>
<desc>DATM CO2 time series</desc>
</entry>

2. I changed it in the run_run.xml file:
<entry id="DATM_CLMNCEP_YR_START" value="2010">
<type>integer</type>
<desc>starting year to loop data over</desc>
</entry>
<entry id="DATM_CLMNCEP_YR_END" value="2010">
<type>integer</type>
<desc>ending year to loop data over</desc>
</entry>

3. I changed it in the run_run.xml file:
<entry id="CLM_CO2_TYPE" value="diagnostic">
<type>char</type>
<valid_values>constant,diagnostic,prognostic</valid_values>

4. I added to the user_nl_clm:
&datm
finidat = '/home/hlh/cesm/inputdata/atm/datm7/CO2/co2_AERmon_CESM2_historical_r10i1p1f1_gn_200001-201412.nc',
/

The lnd_in file is displayed as follows:

&clm_inparm
albice = 0.50,0.30
co2_type = 'diagnostic'
create_crop_landunit = .true.
dtime = 1800
fatmlndfrc = '/home/hlh/cesm/inputdata/share/domains/domain.lnd.fv0.9x1.25_gx1v6.090309.nc'
finidat = '/home/hlh/cesm/inputdata/atm/datm7/CO2/co2_AERmon_CESM2_historical_r10i1p1f1_gn_200001-201412.nc'
fsnowaging = '/home/hlh/cesm/inputdata/lnd/clm2/snicardata/snicar_drdt_bst_fit_60_c070416.nc'
fsnowoptics = '/home/hlh/cesm/inputdata/lnd/clm2/snicardata/snicar_optics_5bnd_c090915.nc'
fsurdat = '/home/hlh/cesm/inputdata/lnd/clm2/surfdata_map/surfdata_0.9x1.25_78pfts_CMIP6_simyr1850_c170824.nc'
glc_do_dynglacier = .true.
glc_snow_persistence_max_days = 0
h2osno_max = 10000.0
int_snow_max = 2000.
irrigate = .true.
maxpatch_glcmec = 10
maxpatch_pft = 79
n_melt_glcmec = 10.0d00
nlevsno = 12
nsegspc = 35
paramfile = '/home/hlh/cesm/inputdata/lnd/clm2/paramdata/clm5_params.c171117.nc'
......


Is this the right way to handle?
How can I see if the model uses the CO2 data I downloaded?
Thanks for your answer!!
 
Top