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

Errors with yd outside when running a monthly forcing dataset

Status
Not open for further replies.
S

syang

Shawn Yang
Member
Dear all,

I am currently running a single-point case (I2000CLM50BgcCruGs) for consecutive years with my own forcing data (modified CMIP6 dataset). I used the "GWSP3v1" for "DATM_MODE".

However, the temporal resolution of my own dataset is monthly, which means there is only one average value for a single month. Thus, I integrated all monthly values from year 2005 to year 2014 into a single file as a input file, named "clmfc_cmiphist_heihe_Prec_2005-2014.nc". Time dimension=120 (12x10).

The domain files are successfully made by CESM tools. The ATM domain file has the same grids with forcing grids. lon: 96 to 103, lat: 37 to 43. nx=70, ny=60
The LND domain file is a single-point file.

Further, I modified the user_nl_datm and user_nl_clm to have a yearly output.

The streamfiles are also modified as shown below

However, once I submit the case, an error showed up with "ERROR: yd outside bounds 38.8555100000000". Does anyone can help me with this problem? and I am not sure whether my setups are correct. THANKS!!!
 

Attachments

  • atm.log.3865240.220629-152336.txt
    4.7 KB · Views: 4
  • user_nl_datm.txt
    894 bytes · Views: 3
  • user_nl_clm.txt
    1.5 KB · Views: 4
  • user_datm.streams.txt.CLMGSWP3v1.TPQW.txt
    815 bytes · Views: 2
  • user_datm.streams.txt.CLMGSWP3v1.Solar.txt
    723 bytes · Views: 2
  • user_datm.streams.txt.CLMGSWP3v1.Precip.txt
    728 bytes · Views: 3
  • own forcing data.png
    own forcing data.png
    25.2 KB · Views: 13
  • env_run.png
    env_run.png
    28 KB · Views: 10
  • cesm.log.3865240.220629-152336.txt
    9.1 KB · Views: 3
S

syang

Shawn Yang
Member
oleson said:
Maybe your error is related to this:

bb.cgd.ucar.edu

My coordinate system was viewed as Polar coordinates, but it is now

Dear NCAR scientist, I recently created my own forcing and domain data for running CLM45 which I got the following error: 'ERROR: ('shr_map_getWts') ERROR yd outside 90 degree bounds'. in shr_abour_mod.f90. I found my coordinate units were treated as radiance but they are not, they are...
bb.cgd.ucar.edu bb.cgd.ucar.edu
Click to expand...
Thank you. It looks like the CLM must have at least 6hr forcing to run. I only have anomaly monthly dataset.
 
S

syang

Shawn Yang
Member
oleson said:
Maybe your error is related to this:

bb.cgd.ucar.edu

My coordinate system was viewed as Polar coordinates, but it is now

Dear NCAR scientist, I recently created my own forcing and domain data for running CLM45 which I got the following error: 'ERROR: ('shr_map_getWts') ERROR yd outside 90 degree bounds'. in shr_abour_mod.f90. I found my coordinate units were treated as radiance but they are not, they are...
bb.cgd.ucar.edu bb.cgd.ucar.edu
Click to expand...
And I would like to ask how can I setup anomaly forcing (monthly dataset) to run the case? From the guide book, it seems that I should change the "user_nl_cpl", "user_nl_datm", "user_nl_clm", and modify stream files. Is there an example I can follow? THANKS！
 
S

syang

Shawn Yang
Member
oleson said:
The only example/instructions I'm aware of are here:


I don't believe you need any changes to user_nl_cpl or user_nl_clm for anomaly forcing specifically.
Click to expand...
Thanks for the link of the documentation. I followed the instruction and modified user_nl_datm with "anomaly forcing stuff". The case can successfully run this time! However, I still have some questions. I set up the start time of 2020. However, in the atm.log, it looks like it only opens the present reference dataset (GSWP3 dataset) without anomaly dataset as a recycling data. My goal is to simulate a regional carbon/nitrogen budgets between 2000 to 2100 under different scenarios (e.g., SSP126, SSP585). Thus, I would like to know how can I check if the anomaly forcing is used and added into the present dataset. Is that automatic? I attached the relevant files below.
 

Attachments

  • env_run.png
    env_run.png
    29 KB · Views: 11
  • user_nl_datm.png
    user_nl_datm.png
    22 KB · Views: 11
  • user_nl_clm.png
    user_nl_clm.png
    4.9 KB · Views: 11
  • datm.streams.txt.Anomaly.Forcing.Temperature.png
    datm.streams.txt.Anomaly.Forcing.Temperature.png
    21.2 KB · Views: 11
  • atm.log.png
    atm.log.png
    64.8 KB · Views: 14
S

syang

Shawn Yang
Member
syang said:
Thanks for the link of the documentation. I followed the instruction and modified user_nl_datm with "anomaly forcing stuff". The case can successfully run this time! However, I still have some questions. I set up the start time of 2020. However, in the atm.log, it looks like it only opens the present reference dataset (GSWP3 dataset) without anomaly dataset as a recycling data. My goal is to simulate a regional carbon/nitrogen budgets between 2000 to 2100 under different scenarios (e.g., SSP126, SSP585). Thus, I would like to know how can I check if the anomaly forcing is used and added into the present dataset. Is that automatic? I attached the relevant files below.
Click to expand...
 

Attachments

  • datm_in.png
    datm_in.png
    24.7 KB · Views: 13
  • lnd_in.png
    lnd_in.png
    37.5 KB · Views: 12
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
There should be some indication in the atm log that the anomaly forcing is being read in.
Do you have something like this in your datm_in (the screen shot only shows part of the file)?

streams = "datm.streams.txt.CLMGSWP3v1.Solar 2005 2005 2014",
"datm.streams.txt.CLMGSWP3v1.Precip 2005 2005 2014",
"datm.streams.txt.CLMGSWP3v1.TPQW 2005 2005 2014",
"datm.streams.txt.presaero.clim_2000 1 2000 2000",
"datm.streams.txt.topo.observed 1 1 1",
"datm.streams.txt.Anomaly.Forcing.Precip 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Temperature 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Pressure 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Humidity 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Uwind 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Vwind 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Shortwave 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Longwave 2005 2090 2090"
 
S

syang

Shawn Yang
Member
oleson said:
There should be some indication in the atm log that the anomaly forcing is being read in.
Do you have something like this in your datm_in (the screen shot only shows part of the file)?

streams = "datm.streams.txt.CLMGSWP3v1.Solar 2005 2005 2014",
"datm.streams.txt.CLMGSWP3v1.Precip 2005 2005 2014",
"datm.streams.txt.CLMGSWP3v1.TPQW 2005 2005 2014",
"datm.streams.txt.presaero.clim_2000 1 2000 2000",
"datm.streams.txt.topo.observed 1 1 1",
"datm.streams.txt.Anomaly.Forcing.Precip 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Temperature 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Pressure 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Humidity 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Uwind 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Vwind 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Shortwave 2005 2090 2090",
"datm.streams.txt.Anomaly.Forcing.Longwave 2005 2090 2090"
Click to expand...
Yes! there are streams like that and I found the indication in the atm_log. Thank you so much!
 
Status
Not open for further replies.
Top