ESCOMP
Don't allow users to try to subset datm for a region with subset_data: doesn't work yet
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 2000 --datm-eyr 2001 --create-user-mods --crop --outdir /glade/scratch/oleson/my_subset_data_Bo1_region2 --verbose
doesn't produce any subsetted datm files.
General bug information
CTSM version you are using: ctsm5.1.dev133
Does this bug cause significantly incorrect results in the model's science? No
Details of bug
The command above produces:
INFO: User = oleson INFO: Current directory = /glade/work/oleson/ctsm5.1.dev133/tools/site_and_regional INFO: dir_input_datm : /glade/p/cesmdata/inputdata/atm/datm7/atm_forcing.datm7.GSWP3.0.5d.v1.c170516 INFO: dir_output_datm: /glade/scratch/oleson/my_subset_data_Bo1_region2/datmdata INFO: ---------------------------------------------------------------------------- INFO: This script extracts a region from the global CTSM datasets. INFO: Creating surface dataset file at region: CO_region INFO: fsurf_in: /glade/p/cesmdata/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_0.9x1.25_hist_78pfts_CMIP6_simyr2000_c190214.nc INFO: fsurf_out: /glade/scratch/oleson/my_subset_data_Bo1_region2/surfdata_CO_region_hist_78pfts_CMIP6_simyr2000_c230815.nc INFO: created file (fsurf_out) /glade/scratch/oleson/my_subset_data_Bo1_region2/surfdata_CO_region_hist_78pfts_CMIP6_simyr2000_c230815.nc INFO: Successfully ran script for a regional case.
and no datm files.
From @adrifoster - this is not working yet. One reason is that you need the mesh file.
We should remove this option and/or give an informative error message if you try to do this.
Solution (from @adrifoster ) - remove argument forcing from the parser for regional, just keeping it in the single point arg parsing.
@billsacks @adrifoster I'm looking into this because I'd like to make some regional-subsetted datm files for testing purposes. Could you explain why the mesh file is needed for this? Is it because the climate forcings are at half-degree?
I'm not sure... I have a feeling that I was just recording comments from @adrifoster ....