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

Calculate carbon dioxide removed by forest in CLM4

Newt

DZ
New Member
Hi all, I want to calculate the amount of carbon dioxide removed by forest in a given area. I wonder if CLM4 or 4.5 can achieve this?
And why is the variable not found in the normal output file? How should it be output?
Thank you!
 

erik

Erik Kluzek
CSEG and Liaisons
Staff member
FCO2 is the net CO2 flux from the land to the atmosphere. It's normally inactive, but you can add it your history files using history namelist options (for example hist_fincl1 to add it to the first history series). It's a gridcell quantity so if you just want the the amount removed from forest you'd need to make some changes to to the code to only output over tree PFT's.
 

Newt

DZ
New Member
FCO2 is the net CO2 flux from the land to the atmosphere. It's normally inactive, but you can add it your history files using history namelist options (for example hist_fincl1 to add it to the first history series). It's a gridcell quantity so if you just want the the amount removed from forest you'd need to make some changes to to the code to only output over tree PFT's.
Thank you for the clear explanation!!
 

Debanjana

Dr Debanjana Das
New Member
Hi, How I can change the forest fraction in the model and check the impact of forest cover in future climates which compset would be the best fit?
Also how what model use for forest information? Could you please navigate the path and how modify it?
 

Debanjana

Dr Debanjana Das
New Member
I want to change the land input files with FINN files... Where I can change or any example of tutorial practical notes could be a great help.
 

oleson

Keith Oleson
CSEG and Liaisons
Staff member
I don't know what FINN files are, but you can change land cover and land use on the surface and landuse datasets.
Those are the files pointed to by "fsurdat" and "flanduse_timeseries" in your CaseDocs/lnd_in file in your case directory.
See Chapters 2 and 27 of the technical note for a description of the default land cover/use

 

Debanjana

Dr Debanjana Das
New Member
I don't know what FINN files are, but you can change land cover and land use on the surface and landuse datasets.
Those are the files pointed to by "fsurdat" and "flanduse_timeseries" in your CaseDocs/lnd_in file in your case directory.
See Chapters 2 and 27 of the technical note for a description of the default land cover/use

Thanks, Keith. Do you have any hands-on tutorial examples? Also what compset would be appropriate to study the climate change signal by altering land use input?
 

oleson

Keith Oleson
CSEG and Liaisons
Staff member
I don't think there are any tutorial examples for that specifically (changing landcover). The SSP clm compsets are appropriate for studying landcover change:

ISSP585Clm50BgcCrop : SSP585_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP126Clm50BgcCrop : SSP126_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP119Clm50BgcCrop : SSP119_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP245Clm50BgcCrop : SSP245_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP370Clm50BgcCrop : SSP370_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP434Clm50BgcCrop : SSP434_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP460Clm50BgcCrop : SSP460_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP534Clm50BgcCrop : SSP534_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV

The most used ones are SSP585, SSP245, SSP126, SSP370. Each of those will have a different landuse_timeseries file. There are some plots here that show how landcover changes in these compsets:


If you want to study how landcover change interacts with climate change, then you can look at fully-coupled (active land, ocean, atmosphere, sea-ice, etc) versions of these simulations (ScenarioMIP). You can find those results on the Earth System Grid. Just use ScenarioMIP and CESM2 as search criteria.

 

Debanjana

Dr Debanjana Das
New Member
I don't think there are any tutorial examples for that specifically (changing landcover). The SSP clm compsets are appropriate for studying landcover change:

ISSP585Clm50BgcCrop : SSP585_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP126Clm50BgcCrop : SSP126_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP119Clm50BgcCrop : SSP119_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP245Clm50BgcCrop : SSP245_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP370Clm50BgcCrop : SSP370_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP434Clm50BgcCrop : SSP434_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP460Clm50BgcCrop : SSP460_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV
ISSP534Clm50BgcCrop : SSP534_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_CISM2%NOEVOLVE_SWAV

The most used ones are SSP585, SSP245, SSP126, SSP370. Each of those will have a different landuse_timeseries file. There are some plots here that show how landcover changes in these compsets:


If you want to study how landcover change interacts with climate change, then you can look at fully-coupled (active land, ocean, atmosphere, sea-ice, etc) versions of these simulations (ScenarioMIP). You can find those results on the Earth System Grid. Just use ScenarioMIP and CESM2 as search criteria.

Thanks a lot! I was just wondering ScenarioMIP outputdata is are available in Cheyenne or Casper? If so, Please share
 

oleson

Keith Oleson
CSEG and Liaisons
Staff member
Output for selected variables is here:

/glade/collections/cdg/data/CMIP6/ScenarioMIP/NCAR/CESM2

The timeseries output is here:

/glade/campaign/collections/cmip/CMIP6/timeseries-cmip6

There might be other holdings elsewhere on glade. I'd post on the Community Projects Forum.
 

Debanjana

Dr Debanjana Das
New Member
Output for selected variables is here:

/glade/collections/cdg/data/CMIP6/ScenarioMIP/NCAR/CESM2

The timeseries output is here:

/glade/campaign/collections/cmip/CMIP6/timeseries-cmip6

There might be other holdings elsewhere on glade. I'd post on the Community Projects Forum.
Thank you ! Can I have a look at the submission script? For my experiment with 1 degree resolution or 0.5 if possible, I want to put present-day land cover change with present-day aerosol and run the coupled model.
 

oleson

Keith Oleson
CSEG and Liaisons
Staff member
The case directories for the ScenarioMIP simulations are here:

/glade/work/cmip6/cases/ScenarioMIP
 

Debanjana

Dr Debanjana Das
New Member
Thank you so much! It was very helpful. Next, can we run the model with a .5 degree resolution? Then what should I use in the command line?
 

oleson

Keith Oleson
CSEG and Liaisons
Staff member
You can run ./query_config --grids from cime/scripts to see what grid/resolutions are available.
I see: f05_g17 (which is 0.47x0.63) and hcru_hcru (which is 0.5deg x 0.5deg).
You should look at the CESM tutorial to see how to specify resolutions when creating a case:

 

Debanjana

Dr Debanjana Das
New Member
Thank you so much! In the scenariomip what kind of input data is used- model data or observation based?
Where I can get the details of the input files and path
 

Debanjana

Dr Debanjana Das
New Member
To be more specific - I am looking at this" fsurdat " -More details of the data- (obs-based or model etc?)
What kind of data / data description

ncdump -h /glade/p/cesmdata/cseg/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_0.9x1.25_hist_78pfts_CMIP6_simyr1850_c190214.nc| more

netcdf surfdata_0.9x1.25_hist_78pfts_CMIP6_simyr1850_c190214 {

dimensions:

lsmlon = 288 ;

lsmlat = 192 ;

nglcec = 10 ;

nglcecp1 = 11 ;

numurbl = 3 ;

nlevurb = 5 ;

numrad = 2 ;

nchar = 256 ;

nlevsoi = 10 ;

time = UNLIMITED ; // (12 currently)

lsmpft = 79 ;

natpft = 15 ;

cft = 64 ;
 

oleson

Keith Oleson
CSEG and Liaisons
Staff member
For ScenarioMIP, these are obviously future projections of landcover, so not based on observations. The file you point to is just an 1850 surface dataset. ScenarioMIP simulations run 2015-2100. The 1850 data matches the 1850 data in the "flanduse_timeseries" file. But the model will use the 2015 data in the flanduse_timeseries file to start and then change landcover every year per that file.
For the source and methods involving that data, please see:

Lawrence, D. M., G. C. Hurtt, A. Arneth, V. Brovkin, K. V. Calvin, A. D. Jones, C. D. Jones, P. J. Lawrence, N. de Noblet-Ducoudre, J. Pongratz, S. I. Seneviratne and E. Shevliakova, 2016. The Land Use Model Intercomparison Project (LUMIP): Rationale and experimental design. Geoscientific Model Development. DOI: 10.5194/gmd-9-2973-2016
 

Debanjana

Dr Debanjana Das
New Member
Greatly appreciate!

Just making sure I understand correctly - so the following are the two land input files used for 4.5 and 8.5 scenario runs?
Correct me if I am wrong.
And if I want to run the model with a new land input file - here I need to change it?


diff b.e21.BSSP245cmip6.f09_g17.CMIP6-SSP2-4.5.101/CaseDocs/lnd_in b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.101/CaseDocs/lnd_in

79a80

> use_init_interp = .true.

93c94

< stream_fldfilename_ndep = '/glade/p/cesmdata/cseg/inputdata/lnd/clm2/ndepdata/fndep_clm_f09_g17.CMIP6-SSP2-4.5-WACCM_1849-2101_monthly_c191007.nc'

---

> stream_fldfilename_ndep = '/glade/p/cesmdata/cseg/inputdata/lnd/clm2/ndepdata/fndep_clm_f09_g17.CMIP6-SSP5-8.5-WACCM_1849-2101_monthly_c191007.nc'

100c101

< stream_fldfilename_popdens = '/glade/p/cesmdata/cseg/inputdata/lnd/clm2/firedata/clmforc.Li_2018_SSP2_CMIP6_hdm_0.5x0.5_AVHRR_simyr1850-2100_c181205.nc'

---

> stream_fldfilename_popdens = '/glade/p/cesmdata/cseg/inputdata/lnd/clm2/firedata/clmforc.Li_2018_SSP5_CMIP6_hdm_0.5x0.5_AVHRR_simyr1850-2100_c181205.nc'

150c151

< flanduse_timeseries = '/glade/p/cesmdata/cseg/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/landuse.timeseries_0.9x1.25_SSP2-4.5_78pfts_CMIP6_simyr1850-2100_c190214.nc'

---

> flanduse_timeseries = '/glade/p/cesmdata/cseg/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/landuse.timeseries_0.9x1.25_SSP5-8.5_78pfts_CMIP6_simyr1850-2100_c190214.nc'

159a161

> init_interp_method = 'general'

312c314

< !# /gpfs/u/home/cmip6/cesm_tags/cesm2.1.4-exp02/components/clm/bld/CLM build-namelist -cimeroot /gpfs/u/home/cmip6/cesm_tags/cesm2.1.4-exp02/cime/scripts/Tools/../.. -infile /glade/work/nanr/cmip6/cases/ScenarioMIP/b.e21.BSSP245cmip6.f09_g17.CMIP6-SSP2-4.5.101/Buildconf/clmconf/namelist -csmdata /glade/p/cesmdata/cseg/inputdata -inputdata /glade/work/nanr/cmip6/cases/ScenarioMIP/b.e21.BSSP245cmip6.f09_g17.CMIP6-SSP2-4.5.101/Buildconf/clm.input_data_list -ignore_ic_year -namelist &clm_inparm start_ymd=20150101 / -no-megan -use_case 1850-2100_SSP2-4.5_transient -res 0.9x1.25 -clm_start_type startup -envxml_dir /glade/work/nanr/cmip6/cases/ScenarioMIP/b.e21.BSSP245cmip6.f09_g17.CMIP6-SSP2-4.5.101 -l_ncpl 48 -lnd_frac /glade/p/cesmdata/cseg/inputdata/share/domains/domain.lnd.fv0.9x1.25_gx1v7.151020.nc -glc_nec 10 -co2_ppmv 284.7 -co2_type diagnostic -config /glade/work/nanr/cmip6/cases/ScenarioMIP/b.e21.BSSP245cmip6.f09_g17.CMIP6-SSP2-4.5.101/Buildconf/clmconf/config_cache.xml -bgc bgc -crop -clm_accelerated_spinup off -lnd_tuning_mode clm5_0_cam6.0 -mask gx1v7

---

> !# /gpfs/u/home/cmip6/cesm_tags/cesm2.1.4-exp02/components/clm/bld/CLM build-namelist -cimeroot /gpfs/u/home/cmip6/cesm_tags/cesm2.1.4-exp02/cime/scripts/Tools/../.. -infile /glade/work/nanr/cmip6/cases/ScenarioMIP/b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.101/Buildconf/clmconf/namelist -csmdata /glade/p/cesmdata/cseg/inputdata -inputdata /glade/work/nanr/cmip6/cases/ScenarioMIP/b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.101/Buildconf/clm.input_data_list -ignore_ic_year -namelist &clm_inparm start_ymd=20150101 / -no-megan -use_case 1850-2100_SSP5-8.5_transient -res 0.9x1.25 -clm_start_type startup -envxml_dir /glade/work/nanr/cmip6/cases/ScenarioMIP/b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.101 -l_ncpl 48 -lnd_frac /glade/p/cesmdata/cseg/inputdata/share/domains/domain.lnd.fv0.9x1.25_gx1v7.151020.nc -glc_nec 10 -co2_ppmv 284.7 -co2_type diagnostic -config /glade/work/nanr/cmip6/cases/ScenarioMIP/b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.101/Buildconf/clmconf/config_cache.xml -bgc bgc -crop -clm_accelerated_spinup off -lnd_tuning_mode clm5_0_cam6.0 -mask gx1v7

-bash-4.2$
 
Top