Hi, I was using ~/CTSM/tools/site_and_regional to get surface data/domain file/datm files, and trying to run a single point simulations.
I can successfully run it while creating cases from CTSM script root, but met errors while creating cases from CESM script root (it said CLM_USRDAT_DIR, PTS_LON, and PTS_LAT do not exist). I asked this questions in CTSM tutorial, seems like CESM was not updated to be able to use this tool for single point simulations.
Then I followed the user guide (1.6.3. Running Single Point Configurations — ctsm CTSM master documentation), but using the domain file and surface data obtained from the subset_data tool.
Here are my code:
# use subset_data tool
cd ~/CTSM/tools/site_and_regional
module load python
ncar_pylib
./subset_data point --lat 55.5 --lon 11.6 --site DK_SOR --dompft 8 --pctpft 100 --create-surface --create-domain --create-datm --datm-syr 2000 --datm-eyr 2000 --create-user-mods --outdir /glade/scratch/jiamengl/subset_data
deactivate
# create case
cd ~/my_cesm_sandbox/cime/scripts
./create_newcase --case ~/case/I2000Clm50BgcCropGs.v3.gm.ocs.c13.point --compset 2000_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_SGLC_SWAV --res CLM_USRDAT --project UCOR0047 --user-mods-dir /glade/scratch/jiamengl/subset_data/user_mods --run-unsupported
cd ~/case/I2000Clm50BgcCropGs.v3.gm.ocs.c13.point
./xmlchange ATM_DOMAIN_PATH=/glade/scratch/jiamengl/subset_data,LND_DOMAIN_PATH=/glade/scratch/jiamengl/subset_data
./xmlchange ATM_DOMAIN_FILE=domain.lnd.fv0.9x1.25_gx1v7_DK_SOR_c220527.nc,LND_DOMAIN_FILE=domain.lnd.fv0.9x1.25_gx1v7_DK_SOR_c220527.nc
./xmlchange CLM_USRDAT_NAME=DK_SOR
./case.setup
./case.build
./xmlchange DIN_LOC_ROOT=/glade/scratch/jiamengl/subset_data
./xmlchange STOP_OPTION=nyears,STOP_N=1
./xmlchange RUN_STARTDATE=2000-01-01
./xmlchange DOUT_S_ROOT=~/case/output/I2000Clm50BgcCropGs.v1.ocs.Ogee.c13.point
./xmlchange DATM_CLMNCEP_YR_START=2000
./xmlchange DATM_CLMNCEP_YR_END=2000
./xmlchange DATM_CLMNCEP_YR_ALIGN=2000
./case.submit
I can successfully submit, but error occurs in reading the stream file for topo_forcing. Please see below for error message in atm.log (also attached)
(shr_strdata_init) stream 5
(shr_strdata_init) filePath = 5
/glade/scratch/jiamengl/subset_data/atm/datm7/topo_forcing/
(shr_strdata_init) fileName = 5
/glade/scratch/jiamengl/subset_data/atm/datm7/topo_forcing/topodata_0.9x1.25_US
GS_070110_stream_c151201.nc
(shr_strdata_init) timeName = 5 time
(shr_strdata_init) lonName = 5 LONGXY
(shr_strdata_init) latName = 5 LATIXY
(shr_strdata_init) hgtName = 5 unknownname
(shr_strdata_init) maskName = 5 mask
(shr_strdata_init) areaName = 5 area
(shr_ncread_domain) ERROR xc ndim = 0
ERROR: (shr_ncread_domain) ERROR ndim must be 1 or 2 for xc
I guess something goes wrong with my datm setting. Also it is strange to me that my user_nl_datm_streams has been changed to point to my subset datm files (I manually substitute the $CLM_USRDAT_DIR with my datm file directory); but seems the datm.streams in Buildconf/datmconf still remain unchanged?
Can anyone help me with this?
I can successfully run it while creating cases from CTSM script root, but met errors while creating cases from CESM script root (it said CLM_USRDAT_DIR, PTS_LON, and PTS_LAT do not exist). I asked this questions in CTSM tutorial, seems like CESM was not updated to be able to use this tool for single point simulations.
Then I followed the user guide (1.6.3. Running Single Point Configurations — ctsm CTSM master documentation), but using the domain file and surface data obtained from the subset_data tool.
Here are my code:
# use subset_data tool
cd ~/CTSM/tools/site_and_regional
module load python
ncar_pylib
./subset_data point --lat 55.5 --lon 11.6 --site DK_SOR --dompft 8 --pctpft 100 --create-surface --create-domain --create-datm --datm-syr 2000 --datm-eyr 2000 --create-user-mods --outdir /glade/scratch/jiamengl/subset_data
deactivate
# create case
cd ~/my_cesm_sandbox/cime/scripts
./create_newcase --case ~/case/I2000Clm50BgcCropGs.v3.gm.ocs.c13.point --compset 2000_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_SGLC_SWAV --res CLM_USRDAT --project UCOR0047 --user-mods-dir /glade/scratch/jiamengl/subset_data/user_mods --run-unsupported
cd ~/case/I2000Clm50BgcCropGs.v3.gm.ocs.c13.point
./xmlchange ATM_DOMAIN_PATH=/glade/scratch/jiamengl/subset_data,LND_DOMAIN_PATH=/glade/scratch/jiamengl/subset_data
./xmlchange ATM_DOMAIN_FILE=domain.lnd.fv0.9x1.25_gx1v7_DK_SOR_c220527.nc,LND_DOMAIN_FILE=domain.lnd.fv0.9x1.25_gx1v7_DK_SOR_c220527.nc
./xmlchange CLM_USRDAT_NAME=DK_SOR
./case.setup
./case.build
./xmlchange DIN_LOC_ROOT=/glade/scratch/jiamengl/subset_data
./xmlchange STOP_OPTION=nyears,STOP_N=1
./xmlchange RUN_STARTDATE=2000-01-01
./xmlchange DOUT_S_ROOT=~/case/output/I2000Clm50BgcCropGs.v1.ocs.Ogee.c13.point
./xmlchange DATM_CLMNCEP_YR_START=2000
./xmlchange DATM_CLMNCEP_YR_END=2000
./xmlchange DATM_CLMNCEP_YR_ALIGN=2000
./case.submit
I can successfully submit, but error occurs in reading the stream file for topo_forcing. Please see below for error message in atm.log (also attached)
(shr_strdata_init) stream 5
(shr_strdata_init) filePath = 5
/glade/scratch/jiamengl/subset_data/atm/datm7/topo_forcing/
(shr_strdata_init) fileName = 5
/glade/scratch/jiamengl/subset_data/atm/datm7/topo_forcing/topodata_0.9x1.25_US
GS_070110_stream_c151201.nc
(shr_strdata_init) timeName = 5 time
(shr_strdata_init) lonName = 5 LONGXY
(shr_strdata_init) latName = 5 LATIXY
(shr_strdata_init) hgtName = 5 unknownname
(shr_strdata_init) maskName = 5 mask
(shr_strdata_init) areaName = 5 area
(shr_ncread_domain) ERROR xc ndim = 0
ERROR: (shr_ncread_domain) ERROR ndim must be 1 or 2 for xc
I guess something goes wrong with my datm setting. Also it is strange to me that my user_nl_datm_streams has been changed to point to my subset datm files (I manually substitute the $CLM_USRDAT_DIR with my datm file directory); but seems the datm.streams in Buildconf/datmconf still remain unchanged?
Can anyone help me with this?