E3SM-Project / E3SM

Energy Exascale Earth System Model source code. NOTE: use "maint" branches for your work. Head of master is not validated.
https://docs.e3sm.org/E3SM
Other
346 stars 353 forks source link

land model calling incorrect surface datasets for trigrid under SSP8.5 #3652

Open stephenprice opened 4 years ago

stephenprice commented 4 years ago

When running the A_WCYCLSSP585_CMIP6 compset with the new trigrid configuration, ELM attempts to load files incompatible with the new grids resulting in an error:

Reading in urban input data from fsurdat file ... (GETFIL): attempting to find local file surfdata_ne30np4_simyr1850_2015_c171018.nc (GETFIL): using /lustre/scratch3/turquoise/sprice/E3SM/input_data/lnd/clm2/surfdata_map/surfdat a_ne30np4_simyr1850_2015_c171018.nc Opened existing file /lustre/scratch3/turquoise/sprice/E3SM/input_data/lnd/clm2/surfdata_map/surfdat a_ne30np4_simyr1850_2015_c171018.nc 81 UrbanInput /lustre/scratch3/turquoise/sprice/E3SM/input_data/lnd/clm2/surfdata_map/surfdat a_ne30np4_simyr1850_2015_c171018.nc UrbanInputldomain and input file do not match dims UrbanInputldomain%ni,ni,= 720 48602 UrbanInputldomain%nj,nj,= 360 1 UrbanInputldomain%ns,ns,= 259200 48602 ENDRUN: ERROR in UrbanParamsType.F90 at line 472

The workaround is the add these changes to the user_clm_nl file before running:

flanduse_timeseries = /pathToFile/landuse.timeseries_0.5x0.5_rcp8.5_simyr2015-2100_c200615.nc fsurdat = /pathToFile/surfdata_0.5x0.5_rcp8.5_simyr2015-2100_c200615.nc

Changes are needed to the land build namelist that populates flanduse_timeseries and fsurdat.

stephenprice commented 4 years ago

@jonbob, @bishtgautam -- please let me know if this description looks adequate for the problem.

bishtgautam commented 4 years ago

@stephenprice, What was the value for --res that you use?

stephenprice commented 4 years ago

This was using ne30pg2_r05_oARRM60to10_ICG

mt5555 commented 4 years ago

not sure if this is related, but I ran into a similar issue before - a confusion between 'res' and 'hgrid':

https://github.com/E3SM-Project/E3SM/issues/3464

bishtgautam commented 4 years ago

@mt5555 The fix implemented here was a great one. In ded99d1440ec88a474a92fe77fc655a0fe5bf891, I have updated the remaining composets in components/clm/bld/namelist_files/use_cases to include hgrid. But, I haven't tested if my fix resolves this issue.

Btw, I found out that res is used when a new case is created, while hgrid is used by mksurfdata.

wlin7 commented 4 years ago

A related one: in PR #3449, it chose to move all fsurdat specifications to namelist_defaults_clm4_5.xml. This only affects compset F2010SC5-CMIP6. It works for grids that have corresponding entries in namelist_defaults_clm4_5.xml. But one for r05 and sim_year 2010 is still needed, which is now added via #3655.