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

Something wrong with my own forcing data

Status
Not open for further replies.
Z

zzwflood

zzwflood
New Member
I'm using the CMFD dataset,I made it be suitable for CESM,主要是将原始时间步长为3小时的数据进行重采样为6小时的数据，然后在


Code: 
./create_newcase --case 5SHJ0.1X0.1 --res CLM_USRDAT --compset HIST_DATM%CRUv7_CLM50%SP_SICE_SOCN_SROF_SGLC_SWAV --mach tc4600 --run-unsupported
./xmlchange ATM_DOMAIN_PATH=/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files,LND_DOMAIN_PATH=/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files
./xmlchange LND_DOMAIN_PATH=/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files,LND_DOMAIN_PATH=/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files
./xmlchange ATM_DOMAIN_FILE=domain.lnd.SHJ0.1X0.1_noocean.230327.nc
./xmlchange LND_DOMAIN_FILE=domain.lnd.SHJ0.1X0.1_noocean.230327.nc
./xmlchange CLM_USRDAT_NAME=CMFD
./xmlchange DATM_MODE='CLM1PT'
./xmlchange CPL_NTASKS=15,ATM_NTASKS=25,LND_NTASKS=45,ICE_NTASKS=15,OCN_NTASKS=15,ROF_NTASKS=10,GLC_NTASKS=20,WAV_NTASKS=15,ESP_NTASKS=1
./xmlchange DATM_CLMNCEP_YR_ALIGN=1996
./xmlchange DATM_CLMNCEP_YR_START=1996
./xmlchange DATM_CLMNCEP_YR_END=1999
./xmlchange DIN_LOC_ROOT_CLMFORC=/home/zzw/software/inputdata/atm/datm7
./xmlchange RUN_STARTDATE=1996-01-01
./xmlchange RESUBMIT=1
./xmlchange STOP_OPTION=nyears
./xmlchange STOP_N=2
./xmlchange CLM_ACCELERATED_SPINUP=on
change env_build.xml---line155----<entry id="CALENDAR" value="GREGORIAN">
./case.setup
#vim user_nl_clm
hist_empty_htapes = .true.
hist_fincl1 = 'QFLX_EVAP_TOT:A', 'QRUNOFF:A', 'QH2OSFC:A', 'QOVER:A','RAIN:A', 'QRUNOFF_R:A', 'QRUNOFF_U:A'
hist_fincl2 = 'QFLX_EVAP_TOT:A', 'QRUNOFF:A', 'QH2OSFC:A', 'QOVER:A','RAIN:A', 'QRUNOFF_R:A', 'QRUNOFF_U:A'
hist_mfilt = 1461,12
hist_nhtfrq = -24,0
fsurdat = '/home/zzw/software/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_SHJ0.1X0.1_hist_16pfts_Irrig_CMIP6_simyr2000_c230328.nc'
flanduse_timeseries = '/home/zzw/software/inputdata/lnd/clm2/surfdata_map/landuse.timeseries_SHJ0.1X0.1_hist_16pfts_Irrig_CMIP6_simyr1850-2015_c230328.nc'
finidat='/home/zzw/software/inputdata/lnd/clm2/initdata_map/clmi.I2000Clm50BgcCrop.2011-01-01.1.9x2.5_gx1v7_gl4_simyr2000_c190312.nc'
use_init_interp = .true.
check_dynpft_consistency = .false.
#vim user_nl_datm
&datm_nml
  decomp = "1d"
  factorfn = "null"
  force_prognostic_true = .false.
  iradsw = 1
  presaero = .true.
  restfilm = "undefined"
  restfils = "undefined"
  wiso_datm = .false.
/
&shr_strdata_nml
  datamode = "CLMNCEP"
  domainfile = "/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files/domain.lnd.SHJ0.1X0.1_noocean.230327.nc"
  dtlimit = 1.5, 1.5, 1.5, 1.5, 1.5, 1.5
  fillalgo = "nn", "nn", "nn", "nn", "nn", "nn"
  fillmask = "nomask", "nomask", "nomask", "nomask", "nomask", "nomask"
  fillread = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
  fillwrite = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
  mapalgo = "bilinear", "bilinear", "bilinear", "bilinear", "bilinear", "nn"
  mapmask = "nomask", "nomask", "nomask", "nomask", "nomask", "nomask"
  mapread = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
  mapwrite = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
  readmode = "single", "single", "single", "single", "single", "single"
  streams = "user_datm.streams.txt.CMFD.Solar 1996 1996 1997",
      "user_datm.streams.txt.CMFD.Precip 1996 1996 1997",
      "user_datm.streams.txt.CMFD.TPQW 1996 1996 1997",
      "user_datm.streams.txt.presaero.trans_1850-2000 1849 1849 1997",
      "user_datm.streams.txt.topo.observed 1 1 1",
      "user_datm.streams.txt.co2tseries.20tr 1850 1850 1997"
  taxmode = "cycle", "cycle", "cycle", "cycle", "cycle", "extend"
  tintalgo = "coszen", "nearest", "linear", "linear", "lower", "linear"
  vectors = "null"
/
./check_input_data
./case.build --skip-provenance-check
./case.submit
Then I get error,I put logs in the attach files.
Could you give some advices?Thanks very much!
 

Attachments

  • atm.log.230420-214507.txt
    237.4 KB · Views: 8
  • cesm.log.230420-214507.txt
    19.6 KB · Views: 7
  • cpl.log.230420-214507.txt
    70.8 KB · Views: 0
  • lnd.log.230420-214507.txt
    200 KB · Views: 3
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
It looks like there might be something wrong with the time stamps on your forcing files. I would check those with ncdump -v time $file and make sure the time stamps are monotonically increasing and are consistent with what appears to be 6-hourly data.
I also see you are using Gregorian calendar. Do you have a calendar = "gregorian" attribute on your time variable?
 
L

LiwenWu

LiwenWu
Member
zzwflood said:
I'm using the CMFD dataset,I made it be suitable for CESM,主要是将原始时间步长为3小时的数据进行重采样为6小时的数据，然后在


Code: 
./create_newcase --case 5SHJ0.1X0.1 --res CLM_USRDAT --compset HIST_DATM%CRUv7_CLM50%SP_SICE_SOCN_SROF_SGLC_SWAV --mach tc4600 --run-unsupported
./xmlchange ATM_DOMAIN_PATH=/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files,LND_DOMAIN_PATH=/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files
./xmlchange LND_DOMAIN_PATH=/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files,LND_DOMAIN_PATH=/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files
./xmlchange ATM_DOMAIN_FILE=domain.lnd.SHJ0.1X0.1_noocean.230327.nc
./xmlchange LND_DOMAIN_FILE=domain.lnd.SHJ0.1X0.1_noocean.230327.nc
./xmlchange CLM_USRDAT_NAME=CMFD
./xmlchange DATM_MODE='CLM1PT'
./xmlchange CPL_NTASKS=15,ATM_NTASKS=25,LND_NTASKS=45,ICE_NTASKS=15,OCN_NTASKS=15,ROF_NTASKS=10,GLC_NTASKS=20,WAV_NTASKS=15,ESP_NTASKS=1
./xmlchange DATM_CLMNCEP_YR_ALIGN=1996
./xmlchange DATM_CLMNCEP_YR_START=1996
./xmlchange DATM_CLMNCEP_YR_END=1999
./xmlchange DIN_LOC_ROOT_CLMFORC=/home/zzw/software/inputdata/atm/datm7
./xmlchange RUN_STARTDATE=1996-01-01
./xmlchange RESUBMIT=1
./xmlchange STOP_OPTION=nyears
./xmlchange STOP_N=2
./xmlchange CLM_ACCELERATED_SPINUP=on
change env_build.xml---line155----<entry id="CALENDAR" value="GREGORIAN">
./case.setup
#vim user_nl_clm
hist_empty_htapes = .true.
hist_fincl1 = 'QFLX_EVAP_TOT:A', 'QRUNOFF:A', 'QH2OSFC:A', 'QOVER:A','RAIN:A', 'QRUNOFF_R:A', 'QRUNOFF_U:A'
hist_fincl2 = 'QFLX_EVAP_TOT:A', 'QRUNOFF:A', 'QH2OSFC:A', 'QOVER:A','RAIN:A', 'QRUNOFF_R:A', 'QRUNOFF_U:A'
hist_mfilt = 1461,12
hist_nhtfrq = -24,0
fsurdat = '/home/zzw/software/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_SHJ0.1X0.1_hist_16pfts_Irrig_CMIP6_simyr2000_c230328.nc'
flanduse_timeseries = '/home/zzw/software/inputdata/lnd/clm2/surfdata_map/landuse.timeseries_SHJ0.1X0.1_hist_16pfts_Irrig_CMIP6_simyr1850-2015_c230328.nc'
finidat='/home/zzw/software/inputdata/lnd/clm2/initdata_map/clmi.I2000Clm50BgcCrop.2011-01-01.1.9x2.5_gx1v7_gl4_simyr2000_c190312.nc'
use_init_interp = .true.
check_dynpft_consistency = .false.
#vim user_nl_datm
&datm_nml
  decomp = "1d"
  factorfn = "null"
  force_prognostic_true = .false.
  iradsw = 1
  presaero = .true.
  restfilm = "undefined"
  restfils = "undefined"
  wiso_datm = .false.
/
&shr_strdata_nml
  datamode = "CLMNCEP"
  domainfile = "/home/zzw/cesmwork/cesm/cime/tools/mapping/gen_domain_files/domain.lnd.SHJ0.1X0.1_noocean.230327.nc"
  dtlimit = 1.5, 1.5, 1.5, 1.5, 1.5, 1.5
  fillalgo = "nn", "nn", "nn", "nn", "nn", "nn"
  fillmask = "nomask", "nomask", "nomask", "nomask", "nomask", "nomask"
  fillread = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
  fillwrite = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
  mapalgo = "bilinear", "bilinear", "bilinear", "bilinear", "bilinear", "nn"
  mapmask = "nomask", "nomask", "nomask", "nomask", "nomask", "nomask"
  mapread = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
  mapwrite = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
  readmode = "single", "single", "single", "single", "single", "single"
  streams = "user_datm.streams.txt.CMFD.Solar 1996 1996 1997",
      "user_datm.streams.txt.CMFD.Precip 1996 1996 1997",
      "user_datm.streams.txt.CMFD.TPQW 1996 1996 1997",
      "user_datm.streams.txt.presaero.trans_1850-2000 1849 1849 1997",
      "user_datm.streams.txt.topo.observed 1 1 1",
      "user_datm.streams.txt.co2tseries.20tr 1850 1850 1997"
  taxmode = "cycle", "cycle", "cycle", "cycle", "cycle", "extend"
  tintalgo = "coszen", "nearest", "linear", "linear", "lower", "linear"
  vectors = "null"
/
./check_input_data
./case.build --skip-provenance-check
./case.submit
Then I get error,I put logs in the attach files.
Could you give some advices?Thanks very much!
Click to expand...
Hi zzw, seeing your post above, I found that you are using CMFD as your forcing data. Recently I have been searching for relevant code to process CMFD into the forcing data required by CLM for a long time since I am not familiar with matlab. Could you please show me a piece of your code as a reference?My email is wuliwen@stumail.nwu.edu.cn. Of cource I know my request may be offensive. Begging for your understanding. Any suggestions would be greatly appreciated! Thank you!
 
H

Houhhu

Nash
Member
@oleson , I try to run the model with my own forcing data, surface data and domain. But in the ./case.build, I meet some error about the namelist.
The domain file and surfdata longitude range from 70°E to 140°E, and latitude range from 15°N to 55°N
my own forcing data CMFD is of the same range as domain file, and in the same format as GSWP3v1.
---------------------------------------------- after ./case.build
Creating component namelists
Calling /public3/home/a2s001011/code/clm5new/cime/src/components/data_comps/datm/cime_config/buildnml
Calling /public3/home/a2s001011/code/clm5new/cime_config/buildnml
ERROR: Command /public3/home/a2s001011/code/clm5new/bld/build-namelist failed rc=255
out=
err=ERROR : CLM build-namelist::CLMBuildNamelist::process_namelist_commandline_infile() : Invalid namelist variable in '-infile' /public3/home/a2s001011/code/clm5new/cime/scripts/my_userdata/Buildconf/clmconf/namelist.
ERROR: in _validate_pair (package Build::NamelistDefinition): Variable name fsurfdat not found in /public3/home/a2s001011/code/clm5new/bld/namelist_files/namelist_definition_drv.xml, /public3/home/a2s001011/code/clm5new/bld/namelist_files/namelist_definition_drv_flds.xml, /public3/home/a2s001011/code/clm5new/bld/namelist_files/namelist_definition_clm4_5.xml
------------------------------------------
Here are some setting:
./create_newcase --case my_userdata --res CLM_USRDAT --compset 2000_DATM%GSWP3v1_CLM50%BGC_SICE_SOCN_SROF_SGLC_SWAV --mach bscca2 --run-unsupported
./xmlchange ATM_DOMAIN_PATH=/code/clm5new/inputdata/share/domains
./xmlchange LND_DOMAIN_PATH=/code/clm5new/inputdata/share/domains
./xmlchange ATM_DOMAIN_FILE=domain.lnd.0.1_0.1_china_0.1_0.1_china.230626.nc
./xmlchange LND_DOMAIN_FILE=domain.lnd.0.1_0.1_china_0.1_0.1_china.230626.nc
./xmlchange CLM_USRDAT_NAME=CMFD
./xmlchange DATM_CLMNCEP_YR_ALIGN=1979
./xmlchange DATM_CLMNCEP_YR_START=1979
./xmlchange DATM_CLMNCEP_YR_END=1979
./xmlchange DIN_LOC_ROOT_CLMFORC=/public3/home/a2s001011/code/clm5new/inputdata/lmwg/atm_forcing.datm7.GSWP3.0.5d.v1.c170516/CMFD
./xmlchange RUN_STARTDATE=1979-01-01
./xmlchange STOP_N=30
./xmlchange STOP_OPTION=ndays
./xmlchange NTASKS=64
./case.setup
----------------------------
vim user_nl_clm
fsurfdat="/public3/home/a2s001011/code/clm5new/inputdata/lnd/clm2/surfdata_map/surfdata_0.1_0.1_china_hist_16pfts_Irrig_CMIP6_simyr2000_c230627.nc"
finidat="/public3/home/a2s001011/code/clm5new/inputdata/lnd/clm2/initdata_map/clmi.I2000Clm50BgcCrop.2011-01-01.1.9x2.5_gx1v7_gl4_simyr2000_c190312.nc"
use_init_interp = .true.
hist_fincl1 = 'TLAI', 'TSA'
hist_nhtfrq = -24
hist_mfilt = 10
vim user_nl_datm
&datm_nml
decomp = "1d"
factorfn = "null"
force_prognostic_true = .false.
iradsw = 1
presaero = .true.
restfilm = "undefined"
restfils = "undefined"
wiso_datm = .false.
/
&shr_strdata_nml
datamode = "CLMNCEP"
domainfile = "/public3/home/a2s001011/code/clm5new/inputdata/share/domains/domain.lnd.0.1_0.1_china_0.1_0.1_china.230626.nc"
dtlimit = 1.5, 1.5, 1.5, 1.5, 1.5,
fillalgo = "nn", "nn", "nn", "nn", "nn"
fillmask = "nomask", "nomask", "nomask", "nomask", "nomask"
fillread = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET",
fillwrite = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET",
mapalgo = "bilinear", "bilinear", "bilinear", "bilinear", "bilinear"
mapmask = "nomask", "nomask", "nomask", "nomask", "nomask"
mapread = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
mapwrite = "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET", "NOT_SET"
readmode = "single", "single", "single", "single", "single"
streams = "user_datm.streams.txt.CMFD.Solar 1979 1979 1979",
"user_datm.streams.txt.CMFD.Precip 1979 1979 1979",
"user_datm.streams.txt.CMFD.TPQW 1979 1979 1979",
"user_datm.streams.txt.presaero.clim_2000 1 2000 2000",
"user_datm.streams.txt.topo.observed 1 1 1"
taxmode = "cycle", "cycle", "cycle", "cycle", "cycle"
tintalgo = "coszen", "nearest", "linear", "linear", "lower"
vectors = "null"
/
--------------------
vim env_build.xml

<entry id="GRID" value="a%CLM_USRDAT_l%CLM_USRDAT_oi%null_r%null_g%null_w%null_m%null">
<entry id="ATM_GRID" value="CLM_USRDAT">
<entry id="ATM_NX" value="1">
<entry id="ATM_NY" value="1">
<entry id="LND_GRID" value="CLM_USRDAT">
<entry id="LND_NX" value="1">
<entry id="LND_NY" value="1">
------------------------
The domain file longitude range from 70°E to 140°E, and latitude range from 15°N to 55°N
ncdump -h domain.lnd.0.1_0.1_china_0.1_0.1_china.230626.nc
dimensions:
n = 280000 ;
ni = 700 ;
nj = 400 ;
nv = 4 ;
variables:
double xc(nj, ni) ;
xc:long_name = "longitude of grid cell center" ;
xc:units = "degrees_east" ;
xc:bounds = "xv" ;
double yc(nj, ni) ;
yc:long_name = "latitude of grid cell center" ;
yc:units = "degrees_north" ;
yc:bounds = "yv" ;
double xv(nj, ni, nv) ;
xv:long_name = "longitude of grid cell verticies" ;
xv:units = "degrees_east" ;
double yv(nj, ni, nv) ;
yv:long_name = "latitude of grid cell verticies" ;
yv:units = "degrees_north" ;
int mask(nj, ni) ;
mask:long_name = "domain mask" ;
mask:note = "unitless" ;
mask:coordinates = "xc yc" ;
mask:comment = "0 value indicates cell is not active" ;
double area(nj, ni) ;
area:long_name = "area of grid cell in radians squared" ;
area:coordinates = "xc yc" ;
area:units = "radian2" ;
double frac(nj, ni) ;
frac:long_name = "fraction of grid cell that is active" ;
frac:coordinates = "xc yc" ;
frac:note = "unitless" ;
frac:filter1 = "error if frac> 1.0+eps or frac < 0.0-eps; eps = 0.1000000E-11" ;
frac:filter2 = "limit frac to [fminval,fmaxval]; fminval= 0.1000000E-02 fmaxval= 1.000000" ;

ncdump -h surfdata_0.1_0.1_china_hist_16pfts_Irrig_CMIP6_simyr2000_c230627.nc
dimensions:
lsmlon = 700 ;
lsmlat = 400 ;
nglcec = 10 ;
nglcecp1 = 11 ;
numurbl = 3 ;
nlevurb = 5 ;
numrad = 2 ;
nchar = 256 ;
nlevsoi = 10 ;
time = UNLIMITED ; // (12 currently)
lsmpft = 17 ;
natpft = 15 ;
cft = 2 ;
variables:
int mxsoil_color ;
mxsoil_color:long_name = "maximum numbers of soil colors" ;
mxsoil_color:units = "unitless" ;
int SOIL_COLOR(lsmlat, lsmlon) ;
SOIL_COLOR:long_name = "soil color" ;
SOIL_COLOR:units = "unitless" ;
double PCT_SAND(nlevsoi, lsmlat, lsmlon) ;
........
.........
Could you give me some advice on how to solve the problem, Should I modify the setting in /public3/home/a2s001011/code/clm5new/bld/namelist_files/namelist_definition_drv.xml and so on?

Best,
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
Your specification of "fsurfdat" is incorrect, it should be "fsurdat".
This is indicated in the error you got:
ERROR: in _validate_pair (package Build::NamelistDefinition): Variable name fsurfdat not found in
 
H

Houhhu

Nash
Member
oleson said:
Your specification of "fsurfdat" is incorrect, it should be "fsurdat".
This is indicated in the error you got:
ERROR: in _validate_pair (package Build::NamelistDefinition): Variable name fsurfdat not found in
Click to expand...
Oops! I am too carless. Thanks for you help
 
Status
Not open for further replies.
Top