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

Question about subsetingDATM files for a regional case?

Status
Not open for further replies.
M

mengqi

mj
Member
Hi,

I am working on a regional simulation using CLM5. Concretely, I tried to subset DATM files but I failed. Here is my script below:

:~/my_ctsm_sandbox/tools/site_and_regional> qcmd -- ./subset_data region --lat1 -40.0 --lat2 15.0 --lon1 275.0 --lon2 330.0 --reg CO_region --create-surface --create-datm --create-user-mods --datm-syr 2000 --datm-eyr 2001 --crop --outdir /glade/scratch/mengqij/my_subset_data_Bo1_region2

I can successfully subset surface data; however, I cannot subset the DATM data (i.e., there are no data files in the datmdata file), although I used "--create-datm", "--datm-syr" and "--datm-eyr".

I am not sure why DATM data does not work. Could anyone please give me some advice?

Thanks!

M.
 
M

mengqi

mj
Member
oleson said:
It didn't work for me either. I've submitted an issue here:

github.com

Don't allow users to try to subset datm for a region with subset_data: doesn't work yet · Issue #2110 · ESCOMP/CTSM

Brief summary of bug As reported by a user and confirmed by myself: ./subset_data region --lat1 -40.0 --lat2 15.0 --lon1 275.0 --lon2 330.0 --reg CO_region --create-surface --create-datm --datm-syr...
github.com github.com
Click to expand...
Thank you, Oleson! To solve the issue, I have an idea below. Please let me know if it does work. Thanks!

If I run a single-point simulation, then I can successfully subset the DATM data. So I was wondering if I can copy the DATM data (single-point case) to the datmdata file (regional case). In other words, I am not sure if the DATM data is the same between the single-point case and the regional case.
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
I don't think that will work. The single-point datm data is just for one grid cell. You could just run with the global GSWP3 data for now, subsetting the data just makes it run faster.
 
M

mengqi

mj
Member
oleson said:
I don't think that will work. The single-point datm data is just for one grid cell. You could just run with the global GSWP3 data for now, subsetting the data just makes it run faster.
Click to expand...
Thank you, Oleson. I was wondering if you refer to running a global simulation using CLM5. If so, could you let me know or offer a few tutorials on how to subset the global GSWP3. I have ZERO experience with global simulations.

Actually, I found a material about the global simulation (CTSM-Tutorial/notebooks/Day1a_GlobalCase.ipynb at main · NCAR/CTSM-Tutorial). However, I did not notice how to subset datm data. That is, in this tutorial, there is no subsetting DATM files for a global simulation.

Thank you in advance!

M.
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
What I meant was that you can run a regional simulation using the global forcing data, you don't necessarily need to subset the forcing data.
 
M

mengqi

mj
Member
Hi Oleson,

Sorry for bothering again. I tried to run a regional simulation using the global forcing data but there is an error. Please see below:

1. Create a new regional case
~/my_ctsm_sandbox/cime/scripts> ./create_newcase --case ~/A_corn_simulation_region/I2000_CTSM_singlept_clm_av_corn_test_1 --res CLM_USRDAT --compset I2000Clm50BgcDvCrop --project UIUC0025 --run-unsupported

2 Set up the case and build the executable
./case.setup
./xmlchange MPILIB=mpi-serial
qcmd -- ./case.build

However, there is an error. I do not understand the use_cndv =T. Actually, I checked the clm5 code. The use_cndv = 'false'. Please give me some input here. Thanks!

ERROR: Command /glade/u/home/mengqij/my_ctsm_sandbox/bld/build-namelist failed rc=255
out=
err=Warning : CLM build-namelist::CLMBuildNamelist::setup_cmdl_dynamic_vegetation() : The use_cndv=T option is deprecated. We do NOT recommend using it. It's known to have issues and it's not calibrated.
-- Add -ignore_warnings option to CLM_BLDNML_OPTS to ignore this warning
 
M

mengqi

mj
Member
oleson said:
What I meant was that you can run a regional simulation using the global forcing data, you don't necessarily need to subset the forcing data.
Click to expand...
One more naive question: how about the surface data? In the single-point simulation, I think that both surface and DATM data can be subsetted together.
So, I was wondering how to subset the surface data. Thanks!
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
The use_cndv is being set by the compset you chose (I2000Clm50BgcDvCrop). Maybe you want I2000Clm50BgcCrop?
 
O

oleson

Keith Oleson
CSEG and Liaisons
Staff member
You've already subset the surface data above right?
/glade/scratch/mengqij/my_subset_data_Bo1_region2/surfdata_0.9x1.25_hist_78pfts_CMIP6_simyr2000_CO_region_c230815.nc
which looks like it is subsetted:
lsmlon = 45 ;
lsmlat = 58 ;
 
M

mengqi

mj
Member
oleson said:
It didn't work for me either. I've submitted an issue here:

github.com

Don't allow users to try to subset datm for a region with subset_data: doesn't work yet · Issue #2110 · ESCOMP/CTSM

Brief summary of bug As reported by a user and confirmed by myself: ./subset_data region --lat1 -40.0 --lat2 15.0 --lon1 275.0 --lon2 330.0 --reg CO_region --create-surface --create-datm --datm-syr...
github.com github.com
Click to expand...
Hi Oleson,

I want to follow up on the issue.

I looked through Bill's advice on the issue. Here is my understanding below (please correct me if I am wrong):

(a) We cannot subset the DATM data for a regional case because there is no mesh file.
(b) In other words, if we use "./subset_data region", then we only subset the surface data.
(c) When I ran a single-point simulation, I did not notice a mesh file. I was wondering if you could give me some advice on how to get a mesh file for a regional case.

Thanks,

M.
 
Status
Not open for further replies.
Top