"lnd_import ERROR: Longwave down sent from the atmosphere model is negative or zero"

[Han Lubin]

Hello Mr.Oleson,

Many thanks for your reply! I have successfully generated the surface file surfdata_0.47x0.63_hist_78pfts_CMIP6_simyr1850_c230128.nc. I try to get the results for surface crop yields corresponding to 0.47x0.63 resolution by CESM2.1.3. However, the following error message appears.

ERROR:
ERROR in /public1/home/sc93599/CESM2.1.3/my_cesm_sandbox/components/clm/src/mai
n/ncdio_pio.F90.in at line 2388
ncd_getiodesc ERROR in vsize 829440 221184 3

the code snippet with line 2388 in ncdio_pio.F90.in are as follows:

vsize = fullsize / gsize
if (mod(fullsize,gsize) /= 0) then
write(iulog,*) subname,' ERROR in vsize ',fullsize,gsize,vsize
call shr_sys_abort(errMsg(sourcefile, __LINE__))
endif

fullsize is "size of entire array on cdf"
gsize is "global size of clmlevel"

I have browsed the forum for similar error messages and I guess the error may be caused by a mismatch between flanduse_timeseries and fatmlndfrc. The flanduse_timeseries I set is landuse.timeseries_0.9x1.25_hist_78pfts_CMIP6_simyr1850-2015_c170824.nc, and fatmlndfrc is domain.lnd.fv0.47 x0.63_gx1v7.180521.nc. Fullsize (829440) is a multiple of 0.9x1.25° resolution corresponding to latitude and longitude, 192x288x15, while gsize is a product of the number of resolutions corresponding to 0.47x0.63°, 384x576. But I don't find the flanduse_timeseries file for 0.47x0.63° resolution with 78 pfts. How can I deal with this error message?

The specific steps and documents are attached.

./create_newcase --case maize_f05_test --compset HIST_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_SGLC_SWAV --res f05_g17 --run-unsupported --compiler intel --mach bscca3
cd maize_f05_test

vi user_nl_clm
hist_nhtfrq=0,0
hist_mfilt=1,12
hist_fincl2='GRAINC_TO_FOOD'
hist_dov2xy=.true.,.false.
irrigate=.false.
fsurdat='/public1/home/sc93599/CESM2.1.3/my_cesm_sandbox/components/clm/tools/mksurfdata_map/surfdata_0.47x0.63_hist_78pfts_CMIP6_simyr1850_c230128.nc'
flanduse_timeseries='/public1/home/sc93599/CESM/my_cesm_sandbox/inputdata/lnd/clm2/surfdata_map/landuse.timeseries_0.9x1.25_hist_78pfts_CMIP6_simyr1850-2015_c170824.nc'
use_init_interp = .true.

./xmlchange STOP_OPTION=nyears
./xmlchange STOP_N=1
./xmlchange RESUBMIT=0
./xmlchange RUN_STARTDATE=1979-01-01
./xmlchange DATM_CLMNCEP_YR_START=1979
./xmlchange DATM_CLMNCEP_YR_END=1979
./xmlchange DATM_CLMNCEP_YR_ALIGN=1979
./xmlchange DOUT_S=TRUE

./case.setup
./case.build --skip-provenance-check
./check_case
./check_input_data
./case.submit

 

[oleson]

You will also need to generate the landuse file in addition to the 1850 surface dataset.

./mksurfdata.pl -res 0.47x0.63 -years 1850-2000

Despite the end year set here to 2000, this should generate a landuse dataset 1850-2015.

 

[Han Lubin]

Hello Mr.Oleson,

Sorry to bother you again. Firstly, many thanks for your advice. I have generated landuse dataset 1850-2105 with 0.47x0.63 resolution and got the output of this case (compset: HIST_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_SGLC_SWAV, res: f05_g17).

Secondly, In order to compare the difference between CAM and DATM results, I try to get further output with FHIST_BGC and f05_g17 in CESM2.1.3, but when executing ./case.build with an error message ERROR: Need to provide valid mapping file between atm and ice in xml variable atm2ice_smapname. I use map_fv0.47x0.63_TO_gx1v7_blin.230203.nc, generated by ./interpic -t fv_0.47x0.63.nc map_fv0.9x1.25_TO_gx1v7_blin.151008.nc map_fv0.47x0.63_TO_gx1v7_blin.230203.nc, as the atm2ice_smapname. But when I execute ./case.build again, it doesn't work and report an error message
ERROR: in _validate_pair (package Build::NamelistDefinition): Variable name atm2ice_smapname not found in /public1/home/sc93599/CESM2.1.3/my_cesm_sandbox/components/cam/bld/namelist_files/namelist_definition.xml

The steps and documents are attached. Could you give me some advice on how to fix this error? Thanks again.

./create_newcase --case maize_CAM_f05_irr --compset HIST_CAM60_CLM50%BGC-CROP_CICE%PRES_DOCN%DOM_MOSART_CISM2%NOEVOLVE_SWAV --res f05_g17 --run-unsupported --compiler intel --mach bscca3
cd maize_CAM_f05_irr
./case.setup

vi user_nl_clm
hist_nhtfrq=0,0
hist_mfilt=1,12
hist_fincl2='GRAINC_TO_FOOD'
hist_dov2xy=.true.,.false.
irrigate=.true.
fsurdat='/public1/home/sc93599/CESM2.1.3/my_cesm_sandbox/components/clm/tools/mksurfdata_map/surfdata_0.47x0.63_hist_78pfts_CMIP6_simyr1850_c230128.nc'
flanduse_timeseries='/public1/home/sc93599/CESM2.1.3/my_cesm_sandbox/components/clm/tools/mksurfdata_map/landuse.timeseries_0.47x0.63_hist_78pfts_CMIP6_simyr1850-2015_c230131.nc'

vi user_nl_cam
ncdata="public1/home/sc93599/CESM/my_cesm_sandbox/inputdata/atm/cam/inic/fv/cami-mam3_0000-01-01_0.47x0.63_L30_c100929.nc"

./xmlchange STOP_OPTION=nyears
./xmlchange STOP_N=1
./xmlchange RESUBMIT=0
./xmlchange RUN_STARTDATE=1979-01-01
./xmlchange DOUT_S=TRUE

./case.build --skip-provenance-check

 

[oleson]

I'm not sure, you should try posting on the CAM forum. I got the same error as you when I tried it. It's possible that resolution/compset combination is not supported in CAM/CLM.
You could try f05_f05_mg17, but I got mapping errors with that also, just different.

 

[Han Lubin]

I'm not sure, you should try posting on the CAM forum. I got the same error as you when I tried it. It's possible that resolution/compset combination is not supported in CAM/CLM.
You could try f05_f05_mg17, but I got mapping errors with that also, just different.

Thank you so much for your patience and advice! I will test further and post the issue to CAM forum.

 

[Ceil]

I've placed all of the mapping files we have for that resolution on our ftp site under the subdirectory mapping_files:

ftp://ftp.cgd.ucar.edu/pub/oleson/mapping_files

Let me know if anything is missing.

Hi Oleson,

I am also using --res f05_g17. But the above files disappear from the ftp. Could you please upload them again?
What's more, could you upload the two files: landuse.timeseries_0.47x0.63_78pfts*nc, surfdata_0.47x0.63_78pfts*nc? Because I only found the following files for 16pfts from isvn-ccsm-nputdata websites:inputdata/lnd/clm2/surfdata_map/landuse.timeseries_0.47x0.63_hist_16pfts_Irrig_CMIP6_simyr1850-2015_c171025.nc, inputdata/lnd/clm2/surfdata_map/surfdata_0.47x0.63_16pfts_Irrig_CMIP6_simyr1850_c170919.nc

Many thanks in advance!

 

[oleson]

The ftp space is aggressively scrubbed, about every 2 weeks or so.
I've put everything we have for that resolution on the ftp site.