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

Questions Regarding SSP Simulations

L

Li Ming

Li Ming
Member
Hi,
I'm using version release-clm5.0.35, and I want to simulate NPP changes from 2025 to 2100. So, I used the following command to create a case: `./create_newcase --case ice-4 --res f19_g16 --compset SSP126_DATM%GSWP3v1_CLM50%BGC_SICE_SOCN_SROF_SGLC_SWAV --compiler intel --mach myintel --run-unsupported`. I made the following main modifications:
<entry id="DATM_CLMNCEP_YR_START" value="1977">
<entry id="DATM_CLMNCEP_YR_END" value="2014">
<entry id="RUN_STARTDATE" value="2025-01-01">
<entry id="STOP_OPTION" value="nyears">
<entry id="STOP_N" value="76">
<entry id="HIST_OPTION" value="nyears">
<entry id="HIST_N" value="1">
And I added anomaly forcing information in user_nl_datm.
anomaly_forcing = 'Anomaly.Forcing.Precip','Anomaly.Forcing.Temperature','Anomaly.Forcing.Pressure','Anomaly.Forcing.Humidity','Anomaly.Forcing.Uwind','Anomaly.Forcing.Vwind','Anomaly.Forcing.Shortwave','Anomaly.Forcing.Longwave'
However, when I run `./check_input_data`, the model throws an error indicating that the relevant files are missing.
Model datm missing file file1 = '/gpfs/gpfs_fs/home/xmx/cesm/inputdata/atm/datm7/anomaly_forcing/af.uas.ccsm4.rcp45.2006-2300.nc'
Model datm missing file file1 = '/gpfs/gpfs_fs/home/xmx/cesm/inputdata/atm/datm7/anomaly_forcing/af.vas.ccsm4.rcp45.2006-2300.nc'
Model datm missing file file1 = '/gpfs/gpfs_fs/home/xmx/cesm/inputdata/atm/datm7/anomaly_forcing/af.rlds.ccsm4.rcp45.2006-2300.nc'
I have the following questions:

1. Is my current setup for future simulations appropriate, or are there areas that need improvement?
2. I downloaded all the files from the "svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/atm/datm7/anomaly_forcing/" folder but couldn't find the missing files. Could you please advise on where to locate these files? I noticed a folder named "CMIP6-SSP1-2.6" within the "anomaly_forcing" directory. Do I need to manually configure the use of the NC files from this folder?
3. If manual configuration is necessary to utilize the files from the "CMIP6-SSP1-2.6" folder, how should I call these files when creating a case for a low emissions scenario, such as "ISSP119Clm50BgcCrop"?
4. I've observed that the files in the "anomaly_forcing" directory have a resolution of 0.9°×1.25°. Is it possible to combine these files with my own atmospheric forcing data at a 0.1° resolution for simulating future scenarios?

I eagerly await your response and am sincerely grateful for any assistance you can provide.
 

Attachments

  • clm5.0.35.jpg
    clm5.0.35.jpg
    130.8 KB · Views: 22
  • check_input_data.jpg
    check_input_data.jpg
    103.2 KB · Views: 23
slevis

slevis

Moderator
On cheyenne (/glade/p/cesmdata/cseg/inputdata/atm/datm7/anomaly_forcing) I see the rcp45 files that the warning describes as missing. Since you wish to use them, I am importing these three files to the repository. Now you should be able to download these three files.
 
L

Li Ming

Li Ming
Member
Thank you for your response. After downloading the missing components, the model is running smoothly. However, I have some new questions that I need your assistance with.
I used the following command to create a case for simulating the changes in NPP (Net Primary Productivity) under future scenarios: "./create_newcase --case test10 --res CLM_USRDAT --compset SSP585_DATM%GSWP3v1_CLM50%BGC_SICE_SOCN_SROF_SGLC_SWAV --compiler gnu --mach mygnu --run-unsupported." I also referenced a forum post for configuring the model (CMIP6 Anomaly Forcing Details | DiscussCESM Forums (ucar.edu)), and the model is running without any issues. But I've noticed some differences between my settings and those in this forum post.. Here is important information from my user_datm.streams.txt.Anomaly.Forcing.Humidity and other user_datm.streams.txt.Anomaly.Forcing.* files:
<filePath>
/gpfs/gpfs_fs/home/xmx/cesm/inputdata/atm/datm7/anomaly_forcing
</filePath>
<fileNames>
domain.permafrostRCN_P2.c2013.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
huss shum_af
</variableNames>
<filePath>
/gpfs/gpfs_fs/home/xmx/cesm/inputdata/atm/datm7/anomaly_forcing
</filePath>
<fileNames>
af.huss.ccsm4.rcp45.2006-2300.nc
Click to expand...
However, in the forum post, the example file information provided by @oleson is as follows:
<filePath>
/glade/p/cgd/tss/people/oleson/anomaly_forcing
</filePath>
<fileNames>
domain.b.e21.BHIST.f09_g17.CMIP6-historical.010.c200203.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
huss shum_af
</variableNames>
<filePath>
/glade/p/cgd/tss/people/oleson/anomaly_forcing
</filePath>
<fileNames>
af.huss.CESM.SSP370.2015-2100.nc
Click to expand...
The file 'af.huss.CESM.SSP370.2015-2100.nc' seems more like a part of 'af.allvars.CESM.SSP1-2.6.2015-2100_c20220628.nc' located in the '/trunk/inputdata/atm/datm7/anomaly_forcing/CMIP6-SSP3-7.0' directory. Therefore, I have the following questions:
Q1: I couldn't find the 'domain.b.e21.BHIST.f09_g17.CMIP6-historical.010.c200203.nc' file. By default, it is using 'domain.permafrostRCN_P2.c2013.nc.' Is this reasonable?
Q2: I want to simulate multiple SSP scenarios, but I noticed that all SSP scenarios are using the 'af.huss.ccsm4.rcp45.2006-2300.nc' by default. Is this reasonable? If not, do I need to use 'af.allvars.CESM.SSP1-2.6.2015-2100_c20220628.nc' from CMIP6-SSP3-7.0 to simulate SSP3.70?
Q3: In the post titled How to incorporate CMIP6 climate anomaly in ISSP-based land model under CESM2.2? | DiscussCESM Forums (ucar.edu), it was mentioned that the model can run successfully by directly using 'af.allvars.CESM.SSP1-2.6.2015-2100_c20220628.nc' in 'user_datm.streams.txt.Anomaly.*' files. Is this safe? Do I need to manually create different NC component files?
Q4: I intend to simulate changes in regional NPP under future scenarios, so I created the compset 'SSP585_DATM%GSWP3v1_CLM50%BGC_SICE_SOCN_SROF_SGLC_SWAV.' However, all SSP scenario compsets supported by clm5.0.35 are of type 'BGC-CROP,' such as 'ISSP585Clm50BgcCrop.' Can I create a case without the Crop component? Do you have any other recommendations regarding this?
Q5: There is no 'CMIP6-SSP1-1.9' folder under '/trunk/inputdata/atm/datm7/anomaly_forcing.' If I want to simulate the SSP1.19 scenario, what files should I use?

My questions are quite numerous, and I'm eagerly looking forward to your response!
 

Attachments

  • user_datm.streams.txt.Anomaly.Forcing.Humidity.txt
    710 bytes · Views: 5
  • user_nl_datm.txt
    3.9 KB · Views: 7
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
To answer Q1-Q3, I'm going to point you to the set of data that I've put on our anonymous ftp at ftp://ftp.cgd.ucar.edu/pub/oleson/anomaly_forcing. This is the set of data you should use since you are using the CLM5 release version of the code. There were various updates to the anomaly forcing because of some reruns of the simulations that the anomaly forcing was based on. The files I've provided have the correct data and are the same data as those on the "allvars" files. The allvars files were created for running with the development version of CESM and I'm not sure the CLM5 release code will handle that.

Q4. The compset you described should work. It's up to you whether you want to run with the crop model or not.
Q5. To my knowledge, anomaly forcing for SSP1-1.9 was not created because CESM was not run for that scenario (officially at least).
 
J

jinping F

冯进平
New Member
oleson said:
To answer Q1-Q3, I'm going to point you to the set of data that I've put on our anonymous ftp at ftp://ftp.cgd.ucar.edu/pub/oleson/anomaly_forcing. This is the set of data you should use since you are using the CLM5 release version of the code. There were various updates to the anomaly forcing because of some reruns of the simulations that the anomaly forcing was based on. The files I've provided have the correct data and are the same data as those on the "allvars" files. The allvars files were created for running with the development version of CESM and I'm not sure the CLM5 release code will handle that.

Q4. The compset you described should work. It's up to you whether you want to run with the crop model or not.
Q5. To my knowledge, anomaly forcing for SSP1-1.9 was not created because CESM was not run for that scenario (officially at least).
Click to expand...
Hello everyone! Could someone please re-upload the following four .nc files from the /cesm/inputdata/atm/datm7/anomaly_forcing/ directory? The files are:

af.allvars.CESM.SSP1-2.6.2015-2100_c20220628.nc
af.allvars.CESM.SSP2-4.5.2015-2100_c20220628.nc
af.allvars.CESM.SSP3-7.0.2015-2100_c20220628.nc
af.allvars.CESM.SSP5-8.5.2015-2100_c20220628.nc
The reason is that I found these four .nc files are showing a NetCDF error. Thank you in advance for uploading them!

[fengjinping@swarm02 anomaly_forcing]$ ncdump -h af.allvars.CESM.SSP1-2.6.2015-2100_c20220628.nc
ncdump: af.allvars.CESM.SSP1-2.6.2015-2100_c20220628.nc: NetCDF: Unknown file format
[fengjinping@swarm02 anomaly_forcing]$ ncdump -h af.allvars.CESM.SSP2-4.5.2015-2100_c20220628.nc
ncdump: af.allvars.CESM.SSP2-4.5.2015-2100_c20220628.nc: NetCDF: Unknown file format
[fengjinping@swarm02 anomaly_forcing]$ ncdump -h af.allvars.CESM.SSP3-7.0.2015-2100_c20220628.nc
ncdump: af.allvars.CESM.SSP3-7.0.2015-2100_c20220628.nc: NetCDF: Unknown file format
[fengjinping@swarm02 anomaly_forcing]$ ncdump -h af.allvars.CESM.SSP5-8.5.2015-2100_c20220628.nc
ncdump: af.allvars.CESM.SSP5-8.5.2015-2100_c20220628.nc: NetCDF: Unknown file format
 
J

jinping F

冯进平
New Member
Could you please share the following 8 anomaly forcing data files with me? If possible, please also share the 8 anomaly forcing data files for the Shared Socioeconomic Pathways (SSP) scenarios SSP126, SSP245, and SSP585. Thank you in advance.
  • Anomaly.Forcing.Temperature: datafiles = /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.tas.CESM.SSP370.2015-2100.nc
  • Anomaly.Forcing.Precip: datafiles = /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.pr.CESM.SSP370.2015-2100.nc
  • Anomaly.Forcing.Humidity: datafiles = /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.huss.CESM.SSP370.2015-2100.nc
  • Anomaly.Forcing.Pressure: datafiles = /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.ps.CESM.SSP370.2015-2100.nc
  • Anomaly.Forcing.Longwave: datafiles = /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.rlds.CESM.SSP370.2015-2100.nc
  • Anomaly.Forcing.Shortwave: datafiles = /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.rsds.CESM.SSP370.2015-2100.nc
  • Anomaly.Forcing.Uwind: datafiles = /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.uas.CESM.SSP370.2015-2100.nc
  • Anomaly.Forcing.Vwind: datafiles = /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.vas.CESM.SSP370.2015-2100.nc."
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
ftp://ftp.cgd.ucar.edu/pub/oleson/anomaly_forcing_mct_for_ftp.tar
 
J

jinping F

冯进平
New Member
I have the following questions regarding the configuration of user_nl_datm:
1.Why is the time setting for datm.streams.txt.presaero.SSP3-7.0 set to 2015-2101, while my file aerodep_clm_SSP126_b.e21.BSSP126cmip6.f09_g17.CMIP6-SSP1-2.6.001_2014-2101_monthly_0.9x1.25_c190523.nc starts from 2014? I would like to verify if the files I am using are reasonable and accurate. Also, is the domain still the domain.b.e21.BHIST.f09_g17.CMIP6-historical.010.c200203.nc you shared with me?
2.I haven't seen where the forcing data "datm.streams.txt.topo.observed 1 1 1" is set. What settings are required?
3.The time setting for "datm.streams.txt.co2tseries.SSP3-7.0" is 2015-2500, while my file (fco2_datm_globalSSP1-2.6__simyr_2014-2501_CMIP6_c190506.nc) is set to 2014-2501. Is my file applied correctly? Also, is the domain still the domain.b.e21.BHIST.f09_g17.CMIP6-historical.010.c200203.nc you shared with me?
3.I have uploaded my settings as an attachment. Could you share the settings for the three files you have:

  • "datm.streams.txt.presaero.SSP3-7.0 2015 2015 2101",
  • "datm.streams.txt.topo.observed 1 1 1",
  • "datm.streams.txt.co2tseries.SSP3-7.0 2015 2015 2500"?
Looking forward to your response, and thank you in advance!"
 

Attachments

  • user_datm.streams.txt.co2tseries.SSP1-2.6.txt
    747 bytes · Views: 0
  • user_datm.streams.txt.presaero.SSP1-2.6.txt
    1.1 KB · Views: 0
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
1. Dec 2014 is used to interpolate to get Jan 1-14 2015 values. Yes, the domain file provided is for the anomaly forcing streams.
2. Should be set by default in the compset of your choice, e.g., SSP126_DATM%GSWP3v1_CLM50%BGC-CROP...... There is no domain file required for co2.
3. I only have a case setup for SSP126 with release-cesm2.1.3. These are the settings:

cat datm.streams.txt.co2tseries.SSP1-2.6

<?xml version="1.0"?>
<file id="stream" version="1.0">
<dataSource>
GENERIC
</dataSource>
<domainInfo>
<variableNames>
time time
lonc lon
latc lat
area area
mask mask
</variableNames>
<filePath>
/glade/p/cesmdata/cseg/inputdata/atm/datm7/CO2
</filePath>
<fileNames>
fco2_datm_globalSSP1-2.6__simyr_2014-2501_CMIP6_c190506.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
CO2 co2diag
</variableNames>
<filePath>
/glade/p/cesmdata/cseg/inputdata/atm/datm7/CO2
</filePath>
<fileNames>
fco2_datm_globalSSP1-2.6__simyr_2014-2501_CMIP6_c190506.nc
</fileNames>
<offset>
0
</offset>
</fieldInfo>
</file>

cat datm.streams.txt.presaero.SSP1-2.6

<?xml version="1.0"?>
<file id="stream" version="1.0">
<dataSource>
GENERIC
</dataSource>
<domainInfo>
<variableNames>
time time
lon lon
lat lat
area area
mask mask
</variableNames>
<filePath>
/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/aero
</filePath>
<fileNames>
aerodep_clm_SSP126_b.e21.BSSP126cmip6.f09_g17.CMIP6-SSP1-2.6.001_2014-2101_monthly_0.9x1.25_c190523.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
BCDEPWET bcphiwet
BCPHODRY bcphodry
BCPHIDRY bcphidry
OCDEPWET ocphiwet
OCPHIDRY ocphidry
OCPHODRY ocphodry
DSTX01WD dstwet1
DSTX01DD dstdry1
DSTX02WD dstwet2
DSTX02DD dstdry2
DSTX03WD dstwet3
DSTX03DD dstdry3
DSTX04WD dstwet4
DSTX04DD dstdry4
</variableNames>
<filePath>
/glade/p/cesmdata/cseg/inputdata/atm/cam/chem/trop_mozart_aero/aero
</filePath>
<fileNames>
aerodep_clm_SSP126_b.e21.BSSP126cmip6.f09_g17.CMIP6-SSP1-2.6.001_2014-2101_monthly_0.9x1.25_c190523.nc
</fileNames>
<offset>
0
</offset>
</fieldInfo>
</file>

cat datm.streams.txt.topo.observed

<?xml version="1.0"?>
<file id="stream" version="1.0">
<dataSource>
GENERIC
</dataSource>
<domainInfo>
<variableNames>
time time
LONGXY lon
LATIXY lat
area area
mask mask
</variableNames>
<filePath>
/glade/p/cesmdata/cseg/inputdata/atm/datm7/topo_forcing
</filePath>
<fileNames>
topodata_0.9x1.25_USGS_070110_stream_c151201.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
TOPO topo
</variableNames>
<filePath>
/glade/p/cesmdata/cseg/inputdata/atm/datm7/topo_forcing
</filePath>
<fileNames>
topodata_0.9x1.25_USGS_070110_stream_c151201.nc
</fileNames>
<offset>
0
</offset>
</fieldInfo>
</file>
 
You must log in or register to reply here.
Top