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

RCP atmospheric forcing for CLM5 projections

Status
Not open for further replies.
G

Gabriel.Hes

Hes
Member
oleson said:
Ok, thanks. Syntax seems correct.
Are the contents of these files as expected, i.e., do they point to files for 1990-2014:

datm.streams.txt.CLMGSWP3v1.Precip
datm.streams.txt.CLMGSWP3v1.Solar
datm.streams.txt.CLMGSWP3v1.TPQW
Click to expand...
No, this is actually my concern. For example the beginning of file
datm.streams.txt.CLMGSWP3v1.Precip is:
<?xml version="1.0"?>
<file id="stream" version="1.0">
<dataSource>
GENERIC
</dataSource>
<domainInfo>
<variableNames>
time time
xc lon
yc lat
area area
mask mask
</variableNames>
<filePath>
/scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516
</filePath>
<fileNames>
domain.lnd.360x720_gswp3.0v1.c170606.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
PRECTmms precn
</variableNames>
<filePath>
/scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/Precip
</filePath>
<fileNames>
clmforc.GSWP3.c2011.0.5x0.5.Prec.1901-01.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1901-02.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1901-03.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1901-04.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1901-05.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1901-06.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1901-07.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1901-08.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1901-09.nc
Click to expand...
and ends like this:
clmforc.GSWP3.c2011.0.5x0.5.Prec.1920-10.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1920-11.nc
clmforc.GSWP3.c2011.0.5x0.5.Prec.1920-12.nc
</fileNames>
<offset>
0
</offset>
</fieldInfo>
</file>
Click to expand...
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
Do you have any user_* versions of these in your case directory that could be setting the years incorrectly, e.g.,

user_datm.streams.txt.CLMGSWP3v1.Precip
 
G

Gabriel.Hes

Hes
Member
No. But my script is in the same directory as the script for the Anomaly forcing simulation. Could this be a problem?
oleson said:
Do you have any user_* versions of these in your case directory that could be setting the years incorrectly, e.g.,

user_datm.streams.txt.CLMGSWP3v1.Precip
Click to expand...
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
I don't know what script you mean.
But you might try resetting your case (./case.setup --reset) and then just run ./preview_namelists and see if things are correct.
Without access to your case, I'm not going to being able to help troubleshoot.
 
G

Gabriel.Hes

Hes
Member
oleson said:
I don't know what script you mean.
But you might try resetting your case (./case.setup --reset) and then just run ./preview_namelists and see if things are correct.
Without access to your case, I'm not going to being able to help troubleshoot.
Click to expand...
I meant my scripts which creates, setups and submits the case. However I have tried to move my script in another folder where the model code is aswell and this doesn't improve anything.
 
G

Gabriel.Hes

Hes
Member
oleson said:
I don't know what script you mean.
But you might try resetting your case (./case.setup --reset) and then just run ./preview_namelists and see if things are correct.
Without access to your case, I'm not going to being able to help troubleshoot.
Click to expand...
When I try the ./case.setup I get the following:
Successfully cleaned .case.run
Successfully cleaned env_mach_specific.xml
Successfully cleaned Macros.make
Successfully cleaned Macros.cmake
ERROR: module command /opt/modules/default/bin/modulecmd python load craype cray-mpich cray-netcdf-hdf5parallel cray-parallel-netcdf daint-mc PrgEnv-intel craype cray-mpich cray-netcdf-hdf5parallel cray-parallel-netcdf daint-mc perftools-base perftools-preload failed with message:
cray-netcdf-hdf5parallel/4.7.4.4(6):ERROR:150: Module 'cray-netcdf-hdf5parallel/4.7.4.4' conflicts with the currently loaded module(s) 'cray-netcdf/4.7.4.4'
cray-netcdf-hdf5parallel/4.7.4.4(6):ERROR:102: Tcl command execution failed: conflict cray-netcdf

daint-mc/21.09(56):ERROR:150: Module 'daint-mc/21.09' conflicts with the currently loaded module(s) 'daint-gpu/21.09'
daint-mc/21.09(56):ERROR:102: Tcl command execution failed: conflict daint-gpu

cray-netcdf-hdf5parallel/4.7.4.4(6):ERROR:150: Module 'cray-netcdf-hdf5parallel/4.7.4.4' conflicts with the currently loaded module(s) 'cray-netcdf/4.7.4.4'
cray-netcdf-hdf5parallel/4.7.4.4(6):ERROR:102: Tcl command execution failed: conflict cray-netcdf

daint-mc/21.09(56):ERROR:150: Module 'daint-mc/21.09' conflicts with the currently loaded module(s) 'daint-gpu/21.09'
daint-mc/21.09(56):ERROR:102: Tcl command execution failed: conflict daint-gpu
Click to expand...

Yes, I understand that without the case its hard to help. Thank you for having guided me through this.
 
G

Gabriel.Hes

Hes
Member
oleson said:
Ah, I was wondering where people were getting instructions to modify user_nl_cpl. That part is not necessary any more, I'm not sure when that changed.
Since you are using the mct coupler, please disregard my instructions for getting the anomaly forcing data from the inputdata repository. The data there is the same as that at /glade/p/cgd/tss/people/oleson/anomaly_forcing but the individual variable files have been converted into a single file, which will not work with the mct coupler version that you are using.
I've put the anomaly forcing data on our ftp site:

ftp://ftp.cgd.ucar.edu/pub/oleson/anomaly_forcing

To run an anomaly-forced SSP simulation (e.g, SSP370), you could start with something like this:

./create_newcase --compset SSP370_DATM%GSWP3v1_CLM51%SP_SICE_SOCN_MOSART_SGLC_SWAV --res hcru_hcru --case ctsm51sp_ctsm51d029_hcru_GSWP3V1_SSP3-7_Anomaly --run-unsupported

The base forcing (GSWP3V1) that the anomaly forcing is applied to comes out of the box with:

<entry id="DATM_CLMNCEP_YR_ALIGN" value="$DATM_CLMNCEP_YR_START">
<entry id="DATM_CLMNCEP_YR_START" value="2001">
<entry id="DATM_CLMNCEP_YR_END" value="2013">

Since we are starting the model here at 2015, do this:

./xmlchange DATM_CLMNCEP_YR_ALIGN=2015

That will align model year 2015 with the 2001 GSWP3V1 forcing.

Then place the following files that are in the ftp directory, into your case directory (you can follow the instructions at 1.5.9. Running with anomaly forcing — ctsm CTSM master documentation to generate these files, but I'm providing them here for completeness):

user_datm.streams.txt.Anomaly.Forcing.Humidity
user_datm.streams.txt.Anomaly.Forcing.Longwave
user_datm.streams.txt.Anomaly.Forcing.Precip
user_datm.streams.txt.Anomaly.Forcing.Pressure
user_datm.streams.txt.Anomaly.Forcing.Shortwave
user_datm.streams.txt.Anomaly.Forcing.Temperature
user_datm.streams.txt.Anomaly.Forcing.Uwind
user_datm.streams.txt.Anomaly.Forcing.Vwind
user_nl_datm

You'll need to change the paths for the domain and anomaly forcing files to the path where your data is located.
Click to expand...
Dear Keith Oleson,
I am still trying to run my future anomaly forcing simulations. I am not sure the model is currently taking the anomaly forcing into account: is there a way to check this? Maybe in the datm_in file?
In this post (copy above), you say that I have to change the domain paths but I am not sure where to do this.
I tried adding the domainfile at the beginning of the user_nl_datm file as below:

cat > user_nl_datm << EOF
domainfile = "/users/ghes/anomaly_forcing/domain.b.e21.BHIST.f09_g17.CMIP6-historical.010.c200203_cdf5.nc"
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'
streams = "datm.streams.txt.CLMGSWP3v1.Solar 2015 2001 2013",
"datm.streams.txt.CLMGSWP3v1.Precip 2015 2001 2013",
"datm.streams.txt.CLMGSWP3v1.TPQW 2015 2001 2013",
"datm.streams.txt.presaero.SSP5-8.5 2015 2015 2101",
"datm.streams.txt.topo.observed 1 1 1",
"datm.streams.txt.co2tseries.SSP5-8.5 2015 2015 2500",
"datm.streams.txt.Anomaly.Forcing.Precip 2015 2015 2100",
"datm.streams.txt.Anomaly.Forcing.Temperature 2015 2015 2100",
"datm.streams.txt.Anomaly.Forcing.Pressure 2015 2015 2100",
"datm.streams.txt.Anomaly.Forcing.Humidity 2015 2015 2100",
"datm.streams.txt.Anomaly.Forcing.Uwind 2015 2015 2100",
"datm.streams.txt.Anomaly.Forcing.Vwind 2015 2015 2100",
"datm.streams.txt.Anomaly.Forcing.Shortwave 2015 2015 2100",
"datm.streams.txt.Anomaly.Forcing.Longwave 2015 2015 2100"


EOF
Click to expand...

Doing so I get the following error in the atm.log file:
(shr_strdata_readnml) reading input namelist file: datm_in
(shr_stream_init) Reading file datm.streams.txt.CLMGSWP3v1.Solar
(shr_stream_init) Reading file datm.streams.txt.CLMGSWP3v1.Precip
(shr_stream_init) Reading file datm.streams.txt.CLMGSWP3v1.TPQW
(shr_stream_init) Reading file datm.streams.txt.presaero.SSP5-8.5
(shr_stream_init) Reading file datm.streams.txt.topo.observed
(shr_stream_init) Reading file datm.streams.txt.co2tseries.SSP5-8.5
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Precip
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Temperature
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Pressure
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Humidity
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Uwind
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Vwind
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Shortwave
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Longwave
(datm_comp_init) datm datamode = CLMNCEP
ERROR: (shr_ncread_domain) ERROR var does not exist frac
Click to expand...

Do you know what this can come from? Should the specification of the domainfile be done differently?
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
The specification of the domain file is done in the user_datm.streams.txt.Anomaly.Forcing.* files (e.g., user_datm.streams.txt.Anomaly.Forcing.Humidity).
You will need to change the path to the file in those user_datm.streams.txt.Anomaly.Forcing.* files. For the examples I've provided, the path to the domain file is specified as:

<filePath>
/glade/p/cgd/tss/people/oleson/anomaly_forcing
</filePath>

Once the model is running you should see something like this in the atm log to indicate that the anomaly forcing is being used:

(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Precip
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Temperature
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Pressure
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Humidity
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Uwind
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Vwind
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Shortwave
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Longwave

and then later in the log file:

(shr_dmodel_readstrm) open : /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.pr.CESM.SSP370.2015-2100.nc
(shr_dmodel_readstrm) file lb: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.pr.CESM.SSP370.2015-2100.nc 1032
(shr_dmodel_readstrm) file ub: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.pr.CESM.SSP370.2015-2100.nc 1
(shr_dmodel_readstrm) open : /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.tas.CESM.SSP370.2015-2100.nc
(shr_dmodel_readstrm) file lb: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.tas.CESM.SSP370.2015-2100.nc 1032
(shr_dmodel_readstrm) file ub: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.tas.CESM.SSP370.2015-2100.nc 1
(shr_dmodel_readstrm) open : /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.ps.CESM.SSP370.2015-2100.nc
(shr_dmodel_readstrm) file lb: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.ps.CESM.SSP370.2015-2100.nc 1032
(shr_dmodel_readstrm) file ub: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.ps.CESM.SSP370.2015-2100.nc 1
(shr_dmodel_readstrm) open : /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.huss.CESM.SSP370.2015-2100.nc
(shr_dmodel_readstrm) file lb: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.huss.CESM.SSP370.2015-2100.nc 1032
(shr_dmodel_readstrm) file ub: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.huss.CESM.SSP370.2015-2100.nc 1
(shr_dmodel_readstrm) open : /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.uas.CESM.SSP370.2015-2100.nc
(shr_dmodel_readstrm) file lb: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.uas.CESM.SSP370.2015-2100.nc 1032
(shr_dmodel_readstrm) file ub: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.uas.CESM.SSP370.2015-2100.nc 1
(shr_dmodel_readstrm) open : /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.vas.CESM.SSP370.2015-2100.nc
(shr_dmodel_readstrm) file lb: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.vas.CESM.SSP370.2015-2100.nc 1032
(shr_dmodel_readstrm) file ub: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.vas.CESM.SSP370.2015-2100.nc 1
(shr_dmodel_readstrm) open : /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.rsds.CESM.SSP370.2015-2100.nc
(shr_dmodel_readstrm) file lb: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.rsds.CESM.SSP370.2015-2100.nc 1032
(shr_dmodel_readstrm) file ub: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.rsds.CESM.SSP370.2015-2100.nc 1
(shr_dmodel_readstrm) open : /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.rlds.CESM.SSP370.2015-2100.nc
(shr_dmodel_readstrm) file lb: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.rlds.CESM.SSP370.2015-2100.nc 1032
(shr_dmodel_readstrm) file ub: /glade/p/cgd/tss/people/oleson/anomaly_forcing/af.rlds.CESM.SSP370.2015-2100.nc 1
 
G

Gabriel.Hes

Hes
Member
Thank you!
I have checked the atm.log file and I do have this:

(shr_strdata_readnml) reading input namelist file: datm_in
(shr_stream_init) Reading file datm.streams.txt.CLMGSWP3v1.Solar
(shr_stream_init) Reading file datm.streams.txt.CLMGSWP3v1.Precip
(shr_stream_init) Reading file datm.streams.txt.CLMGSWP3v1.TPQW
(shr_stream_init) Reading file datm.streams.txt.presaero.SSP5-8.5
(shr_stream_init) Reading file datm.streams.txt.topo.observed
(shr_stream_init) Reading file datm.streams.txt.co2tseries.SSP5-8.5
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Precip
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Temperature
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Pressure
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Humidity
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Uwind
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Vwind
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Shortwave
(shr_stream_init) Reading file datm.streams.txt.Anomaly.Forcing.Longwave
(datm_comp_init) datm datamode = CLMNCEP
Click to expand...

and this:

(datm_comp_init) flds_strm = strm_swdn:strm_precn:strm_tbot:strm_wind:strm_pbot:strm_shum:strm_lwdn:strm_co2diag:strm_prec_af:strm_tbot_af:strm_pbot_af:strm_shum_af:strm_u_af:strm_v_af:strm_swdn_af:strm_lwdn_af
(shr_dmodel_readstrm) open : /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/Solar/clmforc.GSWP3.c2011.0.5x0.5.Solr.2001-01.nc
(shr_dmodel_readstrm) file lb: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/Solar/clmforc.GSWP3.c2011.0.5x0.5.Solr.2001-01.nc 1
(shr_dmodel_readstrm) file ub: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/Solar/clmforc.GSWP3.c2011.0.5x0.5.Solr.2001-01.nc 2
(shr_dmodel_readstrm) open : /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/Precip/clmforc.GSWP3.c2011.0.5x0.5.Prec.2013-12.nc
(shr_dmodel_readstrm) file lb: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/Precip/clmforc.GSWP3.c2011.0.5x0.5.Prec.2013-12.nc 248
(shr_dmodel_readstrm) close : /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/Precip/clmforc.GSWP3.c2011.0.5x0.5.Prec.2013-12.nc
(shr_dmodel_readstrm) open : /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/Precip/clmforc.GSWP3.c2011.0.5x0.5.Prec.2001-01.nc
(shr_dmodel_readstrm) file ub: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/Precip/clmforc.GSWP3.c2011.0.5x0.5.Prec.2001-01.nc 1
(shr_dmodel_readstrm) open : /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/TPHWL/clmforc.GSWP3.c2011.0.5x0.5.TPQWL.2013-12.nc
(shr_dmodel_readstrm) file lb: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/TPHWL/clmforc.GSWP3.c2011.0.5x0.5.TPQWL.2013-12.nc 248
(shr_dmodel_readstrm) close : /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/TPHWL/clmforc.GSWP3.c2011.0.5x0.5.TPQWL.2013-12.nc
(shr_dmodel_readstrm) open : /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/TPHWL/clmforc.GSWP3.c2011.0.5x0.5.TPQWL.2001-01.nc
(shr_dmodel_readstrm) file ub: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/TPHWL/clmforc.GSWP3.c2011.0.5x0.5.TPQWL.2001-01.nc 1
(shr_dmodel_readstrm) open : /scratch/snx3000/ghes/cesm_inputdata/atm/cam/chem/trop_mozart_aero/aero/aerodep_clm_SSP585_b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.001_2014-2101_monthly_0.9x1.25_c190419.nc
(shr_dmodel_readstrm) file lb: /scratch/snx3000/ghes/cesm_inputdata/atm/cam/chem/trop_mozart_aero/aero/aerodep_clm_SSP585_b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.001_2014-2101_monthly_0.9x1.25_c190419.nc 1056
(shr_dmodel_readstrm) file ub: /scratch/snx3000/ghes/cesm_inputdata/atm/cam/chem/trop_mozart_aero/aero/aerodep_clm_SSP585_b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.001_2014-2101_monthly_0.9x1.25_c190419.nc 13
(shr_dmodel_readstrm) open : /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/topo_forcing/topodata_0.9x1.25_USGS_070110_stream_c151201.nc
(shr_dmodel_readstrm) file lb: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/topo_forcing/topodata_0.9x1.25_USGS_070110_stream_c151201.nc 1
(shr_dmodel_readstrm) file ub: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/topo_forcing/topodata_0.9x1.25_USGS_070110_stream_c151201.nc 1
(shr_dmodel_readstrm) open : /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/CO2/fco2_datm_globalSSP5-8.5__simyr_2014-2501_CMIP6_c190506.nc
(shr_dmodel_readstrm) file lb: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/CO2/fco2_datm_globalSSP5-8.5__simyr_2014-2501_CMIP6_c190506.nc 2
(shr_dmodel_readstrm) file ub: /scratch/snx3000/ghes/cesm_inputdata/atm/datm7/CO2/fco2_datm_globalSSP5-8.5__simyr_2014-2501_CMIP6_c190506.nc 2
(shr_dmodel_readstrm) open : /users/ghes/anomaly_forcing/af.pr.CESM.SSP585.2015-2100_cdf5.nc
(shr_dmodel_readstrm) file lb: /users/ghes/anomaly_forcing/af.pr.CESM.SSP585.2015-2100_cdf5.nc 1032
(shr_dmodel_readstrm) file ub: /users/ghes/anomaly_forcing/af.pr.CESM.SSP585.2015-2100_cdf5.nc 1
(shr_dmodel_readstrm) open : /users/ghes/anomaly_forcing/af.tas.CESM.SSP585.2015-2100_cdf5.nc
(shr_dmodel_readstrm) file lb: /users/ghes/anomaly_forcing/af.tas.CESM.SSP585.2015-2100_cdf5.nc 1032
(shr_dmodel_readstrm) file ub: /users/ghes/anomaly_forcing/af.tas.CESM.SSP585.2015-2100_cdf5.nc 1
(shr_dmodel_readstrm) open : /users/ghes/anomaly_forcing/af.ps.CESM.SSP585.2015-2100_cdf5.nc
(shr_dmodel_readstrm) file lb: /users/ghes/anomaly_forcing/af.ps.CESM.SSP585.2015-2100_cdf5.nc 1032
(shr_dmodel_readstrm) file ub: /users/ghes/anomaly_forcing/af.ps.CESM.SSP585.2015-2100_cdf5.nc 1
(shr_dmodel_readstrm) open : /users/ghes/anomaly_forcing/af.huss.CESM.SSP585.2015-2100_cdf5.nc
(shr_dmodel_readstrm) file lb: /users/ghes/anomaly_forcing/af.huss.CESM.SSP585.2015-2100_cdf5.nc 1032
(shr_dmodel_readstrm) file ub: /users/ghes/anomaly_forcing/af.huss.CESM.SSP585.2015-2100_cdf5.nc 1
(shr_dmodel_readstrm) open : /users/ghes/anomaly_forcing/af.uas.CESM.SSP585.2015-2100_cdf5.nc
(shr_dmodel_readstrm) file lb: /users/ghes/anomaly_forcing/af.uas.CESM.SSP585.2015-2100_cdf5.nc 1032
(shr_dmodel_readstrm) file ub: /users/ghes/anomaly_forcing/af.uas.CESM.SSP585.2015-2100_cdf5.nc 1
(shr_dmodel_readstrm) open : /users/ghes/anomaly_forcing/af.vas.CESM.SSP585.2015-2100_cdf5.nc
(shr_dmodel_readstrm) file lb: /users/ghes/anomaly_forcing/af.vas.CESM.SSP585.2015-2100_cdf5.nc 1032
(shr_dmodel_readstrm) file ub: /users/ghes/anomaly_forcing/af.vas.CESM.SSP585.2015-2100_cdf5.nc 1
(shr_dmodel_readstrm) open : /users/ghes/anomaly_forcing/af.rsds.CESM.SSP585.2015-2100_cdf5.nc
(shr_dmodel_readstrm) file lb: /users/ghes/anomaly_forcing/af.rsds.CESM.SSP585.2015-2100_cdf5.nc 1032
(shr_dmodel_readstrm) file ub: /users/ghes/anomaly_forcing/af.rsds.CESM.SSP585.2015-2100_cdf5.nc 1
(shr_dmodel_readstrm) open : /users/ghes/anomaly_forcing/af.rlds.CESM.SSP585.2015-2100_cdf5.nc
(shr_dmodel_readstrm) file lb: /users/ghes/anomaly_forcing/af.rlds.CESM.SSP585.2015-2100_cdf5.nc 1032
(shr_dmodel_readstrm) file ub: /users/ghes/anomaly_forcing/af.rlds.CESM.SSP585.2015-2100_cdf5.nc 1
(datm_comp_run) max values = 317.767639160156 0.000000000000000E+000
Click to expand...

I guess I was disturbed by the fact that there were also GSWP3 climate forcing present, but now I realise that this is the "present day" climate (probably looping on the 2001-2013 period) on which the anomaly forcing is applied to, so it is normal to have it here too. Am I correct ? Is there a preferable period to choose for the "observed climate" to loop on ? Would a ten year period like 1996-2005 work?

Related to this, I tried to modify the env_run.xml file by setting the following parameters (defining the loop):
./xmlchange DATM_CLMNCEP_YR_START=1996, DATM_CLMNCEP_YR_END=2005
However, these modifications were never taken into account. I suspect this is due to the fact that I am using "DATM_MODE" value="CLMGSWP3v1". For instance it is stated that the DATM_CLMNCEP_YR_START and DATM_CLMNCEP_YR_END parameters are used for for CLM_QIAN or CLMCRUNCEP modes (see 1.2.2. Customizing CLM’s Configuration — ctsm release-clm5.0 documentation). Do you think this could be the issue? Were I to use the CRUNCEP data could I set these looping parameters properly?
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
Yes, it is the present day climate that the anomaly forcing is applied to. I think I would stick with the suggested present day climate instead of 1996-2005 since the anomaly forcings are based on the more recent data.
Changing DATM_CLMNCEP_YR_START and DATM_CLMNCEP_YR_END should work for GSWP3v1, you must have something else overriding those settings.
 
G

Gabriel.Hes

Hes
Member
Gabriel.Hes said:
No, this is actually my concern. For example the beginning of file
datm.streams.txt.CLMGSWP3v1.Precip is:

and ends like this:
Click to expand...
Dear Keith,
I believe that I have solved the problem by setting the following commands after the ./case.build:
./xmlchange DATM_CLMNCEP_YR_ALIGN=1990
./xmlchange DATM_CLMNCEP_YR_START=1990
./xmlchange DATM_CLMNCEP_YR_END=2014
Previously, I had put them after the ./creat_newcase but this was apparently not the proper order.
Now, the forcings dates are consistent with what I ask.

Thank you for your help!
 
H

Hemraj

Hemraj Bhattarai
Member
Hi,
upon running

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'

I can see datm.streams.txt.Anomaly.Forcing.* files in my run directly and also in CaseDocs. However, I can't see user_nl_datm.streams.txt.Anomaly.Forcing.* in my case directory? What's the way of getting them in our case directory?

Many thanks!
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
Can you create a new post that describes exactly what you are trying to do? This thread is long with several different subjects/projects and it's time to end it.
 
Y

Yuan Sun

Yuan Sun
Member
Hemraj said:
Hi,
upon running

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'

I can see datm.streams.txt.Anomaly.Forcing.* files in my run directly and also in CaseDocs. However, I can't see user_nl_datm.streams.txt.Anomaly.Forcing.* in my case directory? What's the way of getting them in our case directory?

Many thanks!
Click to expand...
You can run the command ./preview_namelists afer adding anomaly_forcing description to generate the datm.streams.txt.Anomaly.Forcing.* files

It works in my case.
 
Status
Not open for further replies.
Top