Skip to content

Commit eb9efce

Browse files
bishtgautambishtgautam
committed
Merge branch 'bishtgautam/lnd/interpinic-notes' (PR #6582)
Adds notes on using the following ELM tools: - `mksurfdata_map`: For creating a surface dataset - `interpinic`: For creating an initial condition. [BFB]
2 parents 4d85553 + 7ba6c85 commit eb9efce

File tree

5 files changed

+185
-42
lines changed

5 files changed

+185
-42
lines changed

‎components/elm/bld/namelist_files/namelist_defaults.xml‎

Lines changed: 41 additions & 41 deletions
Original file line numberDiff line numberDiff line change
@@ -334,7 +334,7 @@ lnd/clm2/surfdata_map/surfdata_64x128_simyr2000_c170111.nc</fsurdat>
334334
<fsurdat hgrid="NLDAS" sim_year="2000" irrig=".false." use_crop=".false." >
335335
lnd/clm2/surfdata_map/surfdata_nldas2_simyr2000_c181207.nc</fsurdat>
336336

337-
<fsurdat hgrid="r0125xr0125" sim_year="2000" use_crop=".false." >
337+
<fsurdat hgrid="r0125" sim_year="2000" use_crop=".false." >
338338
lnd/clm2/surfdata_map/surfdata_0.125x0.125_simyr2000_c190730.nc</fsurdat>
339339
<fsurdat hgrid="0.9x1.25" sim_year="2000" use_crop=".false." >
340340
lnd/clm2/surfdata_map/surfdata_0.9x1.25_simyr2000_c180404.nc</fsurdat>
@@ -405,7 +405,7 @@ lnd/clm2/surfdata_map/surfdata_360x720cru_simyr1850_c180216.nc</fsurdat>
405405
<fsurdat hgrid="48x96" sim_year="1850" use_crop=".false." >
406406
lnd/clm2/surfdata_map/surfdata_48x96_simyr1850_c130927.nc</fsurdat>
407407

408-
<fsurdat hgrid="r0125xr0125" sim_year="1850" use_crop=".false." >
408+
<fsurdat hgrid="r0125" sim_year="1850" use_crop=".false." >
409409
lnd/clm2/surfdata_map/surfdata_0.125x0.125_simyr1850_c190730.nc</fsurdat>
410410
<fsurdat hgrid="0.9x1.25" sim_year="1850" use_crop=".false." >
411411
lnd/clm2/surfdata_map/surfdata_0.9x1.25_simyr1850_c180306.nc</fsurdat>
@@ -1695,13 +1695,6 @@ this mask will have smb calculated over the entire global land surface
16951695
<map frm_hgrid="1km-merge-10min" frm_lmask="HYDRO1K-merge-nomask" to_hgrid="2.5x3.33" to_lmask="nomask"
16961696
>lnd/clm2/mappingdata/maps/2.5x3.33/map_1km-merge-10min_HYDRO1K-merge-nomask_to_2.5x3.33_nomask_aave_da_c130405.nc</map>
16971697

1698-
<!-- mapping files for 0.25x0.25 START added on Wed May 1 11:10:45 2024-->
1699-
1700-
<map frm_hgrid="0.01x0.01" frm_lmask="nomask" to_hgrid="0.25x0.25" to_lmask="nomask"
1701-
>lnd/clm2/mappingdata/maps/0.25x0.25/map_0.01x0.01_nomask_to_0.25x0.25_nomask_aave_da_c240501.nc</map>
1702-
1703-
<!-- mapping files for 0.25x0.25 END -->
1704-
17051698
<!-- mapping files for 0.5x0.5 START added on Wed Sep 19 09:05:45 2012-->
17061699
<!-- Created by lnd/clm/bld/namelist_files/createMapEntry.pl-->
17071700

@@ -2065,76 +2058,83 @@ this mask will have smb calculated over the entire global land surface
20652058
<map frm_hgrid="0.5x0.5" frm_lmask="GSDTG2000" to_hgrid="ne0np4_northamericax4v1.pg2" to_lmask="nomask"
20662059
>lnd/clm2/mappingdata/maps/northamericax4v1pg2/map_0.5x0.5_GSDTG2000_to_northamericax4v1pg2_nomask_aave_da_c210112.nc</map>
20672060

2068-
<!-- mapping files for r0125xr0125 START added on Fri Aug 1 15:44:32 2014-->
2061+
<!-- mapping files for r0125 START added on Fri Aug 1 15:44:32 2014-->
20692062
<!-- Created by lnd/clm/bld/namelist_files/createMapEntry.pl-->
20702063

2071-
<map frm_hgrid="0.5x0.5" frm_lmask="AVHRR" to_hgrid="r0125xr0125" to_lmask="nomask"
2064+
<map frm_hgrid="0.5x0.5" frm_lmask="AVHRR" to_hgrid="r0125" to_lmask="nomask"
20722065
>lnd/clm2/mappingdata/maps/0.125x0.125/map_0.5x0.5_AVHRR_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2073-
<map frm_hgrid="0.5x0.5" frm_lmask="MODIS" to_hgrid="r0125xr0125" to_lmask="nomask"
2066+
<map frm_hgrid="0.5x0.5" frm_lmask="MODIS" to_hgrid="r0125" to_lmask="nomask"
20742067
>lnd/clm2/mappingdata/maps/0.125x0.125/map_0.5x0.5_MODIS_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2075-
<map frm_hgrid="0.9x1.25" frm_lmask="GRDC" to_hgrid="r0125xr0125" to_lmask="nomask"
2068+
<map frm_hgrid="0.9x1.25" frm_lmask="GRDC" to_hgrid="r0125" to_lmask="nomask"
20762069
>lnd/clm2/mappingdata/maps/0.125x0.125/map_0.9x1.25_GRDC_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2077-
<map frm_hgrid="10x10min" frm_lmask="IGBPmergeICESatGIS" to_hgrid="r0125xr0125" to_lmask="nomask"
2070+
<map frm_hgrid="10x10min" frm_lmask="IGBPmergeICESatGIS" to_hgrid="r0125" to_lmask="nomask"
20782071
>lnd/clm2/mappingdata/maps/0.125x0.125/map_10x10min_IGBPmergeICESatGIS_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2079-
<map frm_hgrid="10x10min" frm_lmask="nomask" to_hgrid="r0125xr0125" to_lmask="nomask"
2072+
<map frm_hgrid="10x10min" frm_lmask="nomask" to_hgrid="r0125" to_lmask="nomask"
20802073
>lnd/clm2/mappingdata/maps/0.125x0.125/map_10x10min_nomask_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2081-
<map frm_hgrid="1km-merge-10min" frm_lmask="HYDRO1K-merge-nomask" to_hgrid="r0125xr0125" to_lmask="nomask"
2074+
<map frm_hgrid="1km-merge-10min" frm_lmask="HYDRO1K-merge-nomask" to_hgrid="r0125" to_lmask="nomask"
20822075
>lnd/clm2/mappingdata/maps/0.125x0.125/map_1km-merge-10min_HYDRO1K-merge-nomask_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2083-
<map frm_hgrid="360x720cru" frm_lmask="cruncep" to_hgrid="r0125xr0125" to_lmask="nomask"
2076+
<map frm_hgrid="360x720cru" frm_lmask="cruncep" to_hgrid="r0125" to_lmask="nomask"
20842077
>lnd/clm2/mappingdata/maps/0.125x0.125/map_360x720cru_cruncep_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2085-
<map frm_hgrid="3x3min" frm_lmask="GLOBE-Gardner-mergeGIS" to_hgrid="r0125xr0125" to_lmask="nomask"
2078+
<map frm_hgrid="3x3min" frm_lmask="GLOBE-Gardner-mergeGIS" to_hgrid="r0125" to_lmask="nomask"
20862079
>lnd/clm2/mappingdata/maps/0.125x0.125/map_3x3min_GLOBE-Gardner-mergeGIS_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2087-
<map frm_hgrid="3x3min" frm_lmask="GLOBE-Gardner" to_hgrid="r0125xr0125" to_lmask="nomask"
2080+
<map frm_hgrid="3x3min" frm_lmask="GLOBE-Gardner" to_hgrid="r0125" to_lmask="nomask"
20882081
>lnd/clm2/mappingdata/maps/0.125x0.125/map_3x3min_GLOBE-Gardner_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2089-
<map frm_hgrid="3x3min" frm_lmask="LandScan2004" to_hgrid="r0125xr0125" to_lmask="nomask"
2082+
<map frm_hgrid="3x3min" frm_lmask="LandScan2004" to_hgrid="r0125" to_lmask="nomask"
20902083
>lnd/clm2/mappingdata/maps/0.125x0.125/map_3x3min_LandScan2004_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2091-
<map frm_hgrid="3x3min" frm_lmask="MODIS" to_hgrid="r0125xr0125" to_lmask="nomask"
2084+
<map frm_hgrid="3x3min" frm_lmask="MODIS" to_hgrid="r0125" to_lmask="nomask"
20922085
>lnd/clm2/mappingdata/maps/0.125x0.125/map_3x3min_MODIS_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2093-
<map frm_hgrid="3x3min" frm_lmask="USGS" to_hgrid="r0125xr0125" to_lmask="nomask"
2086+
<map frm_hgrid="3x3min" frm_lmask="USGS" to_hgrid="r0125" to_lmask="nomask"
20942087
>lnd/clm2/mappingdata/maps/0.125x0.125/map_3x3min_USGS_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2095-
<map frm_hgrid="5x5min" frm_lmask="IGBP-GSDP" to_hgrid="r0125xr0125" to_lmask="nomask"
2088+
<map frm_hgrid="5x5min" frm_lmask="IGBP-GSDP" to_hgrid="r0125" to_lmask="nomask"
20962089
>lnd/clm2/mappingdata/maps/0.125x0.125/map_5x5min_IGBP-GSDP_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2097-
<map frm_hgrid="5x5min" frm_lmask="ISRIC-WISE" to_hgrid="r0125xr0125" to_lmask="nomask"
2090+
<map frm_hgrid="5x5min" frm_lmask="ISRIC-WISE" to_hgrid="r0125" to_lmask="nomask"
20982091
>lnd/clm2/mappingdata/maps/0.125x0.125/map_5x5min_ISRIC-WISE_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2099-
<map frm_hgrid="5x5min" frm_lmask="nomask" to_hgrid="r0125xr0125" to_lmask="nomask"
2092+
<map frm_hgrid="5x5min" frm_lmask="nomask" to_hgrid="r0125" to_lmask="nomask"
21002093
>lnd/clm2/mappingdata/maps/0.125x0.125/map_5x5min_nomask_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
2101-
<map frm_hgrid="0.5x0.5" frm_lmask="GSDTG2000" to_hgrid="r0125xr0125" to_lmask="nomask"
2094+
<map frm_hgrid="0.5x0.5" frm_lmask="GSDTG2000" to_hgrid="r0125" to_lmask="nomask"
21022095
>lnd/clm2/mappingdata/maps/0.125x0.125/map_0.5x0.5_GSDTG2000_to_0.125x0.125_nomask_aave_da_c190725.nc</map>
21032096

2104-
<!-- mapping files for r0125xr0125 END -->
2097+
<!-- mapping files for r0125 END -->
21052098

2106-
<map frm_hgrid="0.5x0.5" frm_lmask="AVHRR" to_hgrid="r025xr025" to_lmask="nomask"
2099+
<map frm_hgrid="0.5x0.5" frm_lmask="AVHRR" to_hgrid="r025" to_lmask="nomask"
21072100
>lnd/clm2/mappingdata/maps/0.25x0.25/map_0.5x0.5_AVHRR_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2108-
<map frm_hgrid="0.5x0.5" frm_lmask="MODIS" to_hgrid="r025xr025" to_lmask="nomask"
2101+
<map frm_hgrid="0.5x0.5" frm_lmask="MODIS" to_hgrid="r025" to_lmask="nomask"
21092102
>lnd/clm2/mappingdata/maps/0.25x0.25/map_0.5x0.5_MODIS_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2110-
<map frm_hgrid="0.9x1.25" frm_lmask="GRDC" to_hgrid="r025xr025" to_lmask="nomask"
2103+
<map frm_hgrid="0.9x1.25" frm_lmask="GRDC" to_hgrid="r025" to_lmask="nomask"
21112104
>lnd/clm2/mappingdata/maps/0.25x0.25/map_0.9x1.25_GRDC_to_0.25x0.25_nomask_aave_da_c240124.nc</map>
21122105
<map frm_hgrid="10x10min" frm_lmask="IGBPmergeICESatGIS" to_hgrid="r025xr025" to_lmask="nomask"
21132106
>lnd/clm2/mappingdata/maps/0.25x0.25/map_10x10min_IGBPmergeICESatGIS_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2114-
<map frm_hgrid="10x10min" frm_lmask="nomask" to_hgrid="r025xr025" to_lmask="nomask"
2107+
<map frm_hgrid="10x10min" frm_lmask="nomask" to_hgrid="r025" to_lmask="nomask"
21152108
>lnd/clm2/mappingdata/maps/0.25x0.25/map_10x10min_nomask_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2116-
<map frm_hgrid="1km-merge-10min" frm_lmask="HYDRO1K-merge-nomask" to_hgrid="r025xr025" to_lmask="nomask"
2109+
<map frm_hgrid="1km-merge-10min" frm_lmask="HYDRO1K-merge-nomask" to_hgrid="r025" to_lmask="nomask"
21172110
>lnd/clm2/mappingdata/maps/0.25x0.25/map_1km-merge-10min_HYDRO1K-merge-nomask_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
21182111
<map frm_hgrid="360x720cru" frm_lmask="cruncep" to_hgrid="r025xr025" to_lmask="nomask"
21192112
>lnd/clm2/mappingdata/maps/0.25x0.25/map_360x720cru_cruncep_to_0.25x0.25_nomask_aave_da_c240124.nc</map>
2120-
<map frm_hgrid="3x3min" frm_lmask="GLOBE-Gardner-mergeGIS" to_hgrid="r025xr025" to_lmask="nomask"
2113+
<map frm_hgrid="3x3min" frm_lmask="GLOBE-Gardner-mergeGIS" to_hgrid="r025" to_lmask="nomask"
21212114
>lnd/clm2/mappingdata/maps/0.25x0.25/map_3x3min_GLOBE-Gardner-mergeGIS_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2122-
<map frm_hgrid="3x3min" frm_lmask="GLOBE-Gardner" to_hgrid="r025xr025" to_lmask="nomask"
2115+
<map frm_hgrid="3x3min" frm_lmask="GLOBE-Gardner" to_hgrid="r025" to_lmask="nomask"
21232116
>lnd/clm2/mappingdata/maps/0.25x0.25/map_3x3min_GLOBE-Gardner_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2124-
<map frm_hgrid="3x3min" frm_lmask="LandScan2004" to_hgrid="r025xr025" to_lmask="nomask"
2117+
<map frm_hgrid="3x3min" frm_lmask="LandScan2004" to_hgrid="r025" to_lmask="nomask"
21252118
>lnd/clm2/mappingdata/maps/0.25x0.25/map_3x3min_LandScan2004_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2126-
<map frm_hgrid="3x3min" frm_lmask="MODIS" to_hgrid="r025xr025" to_lmask="nomask"
2119+
<map frm_hgrid="3x3min" frm_lmask="MODIS" to_hgrid="r025" to_lmask="nomask"
21272120
>lnd/clm2/mappingdata/maps/0.25x0.25/map_3x3min_MODIS_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2128-
<map frm_hgrid="3x3min" frm_lmask="USGS" to_hgrid="r025xr025" to_lmask="nomask"
2121+
<map frm_hgrid="3x3min" frm_lmask="USGS" to_hgrid="r025" to_lmask="nomask"
21292122
>lnd/clm2/mappingdata/maps/0.25x0.25/map_3x3min_USGS_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2130-
<map frm_hgrid="5x5min" frm_lmask="IGBP-GSDP" to_hgrid="r025xr025" to_lmask="nomask"
2123+
<map frm_hgrid="5x5min" frm_lmask="IGBP-GSDP" to_hgrid="r025" to_lmask="nomask"
21312124
>lnd/clm2/mappingdata/maps/0.25x0.25/map_5x5min_IGBP-GSDP_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2132-
<map frm_hgrid="5x5min" frm_lmask="ISRIC-WISE" to_hgrid="r025xr025" to_lmask="nomask"
2125+
<map frm_hgrid="5x5min" frm_lmask="ISRIC-WISE" to_hgrid="r025" to_lmask="nomask"
21332126
>lnd/clm2/mappingdata/maps/0.25x0.25/map_5x5min_ISRIC-WISE_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2134-
<map frm_hgrid="5x5min" frm_lmask="nomask" to_hgrid="r025xr025" to_lmask="nomask"
2127+
<map frm_hgrid="360x720cru" frm_lmask="cruncep" to_hgrid="r025" to_lmask="nomask"
2128+
>lnd/clm2/mappingdata/maps/0.25x0.25/map_360x720_cruncep_to_0.25x0.25_nomask_aave_da_c240124.nc</map>
2129+
<map frm_hgrid="5x5min" frm_lmask="nomask" to_hgrid="r025" to_lmask="nomask"
21352130
>lnd/clm2/mappingdata/maps/0.25x0.25/map_5x5min_nomask_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2136-
<map frm_hgrid="0.5x0.5" frm_lmask="GSDTG2000" to_hgrid="r025xr025" to_lmask="nomask"
2131+
<map frm_hgrid="0.5x0.5" frm_lmask="GSDTG2000" to_hgrid="r025" to_lmask="nomask"
21372132
>lnd/clm2/mappingdata/maps/0.25x0.25/map_0.5x0.5_GSDTG2000_to_0.25x0.25_nomask_aave_da_c240123.nc</map>
2133+
<map frm_hgrid="0.01x0.01" frm_lmask="nomask" to_hgrid="r025" to_lmask="nomask"
2134+
>lnd/clm2/mappingdata/maps/0.25x0.25/map_0.01x0.01_nomask_to_0.25x0.25_nomask_aave_da_c240501.nc</map>
2135+
<map frm_hgrid="0.1x0.1" frm_lmask="nomask" to_hgrid="r025" to_lmask="nomask"
2136+
>lnd/clm2/mappingdata/maps/0.25x0.25/map_0.1x0.1_nomask_to_0.25x0.25_nomask_aave_da_c240308.nc</map>
2137+
21382138

21392139
<map frm_hgrid="0.5x0.5" frm_lmask="AVHRR" to_hgrid="ne0np4_antarcticax4v1" to_lmask="nomask"
21402140
>lnd/clm2/mappingdata/maps/antarcticax4v1/map_0.5x0.5_AVHRR_to_antarcticax4v1_nomask_aave_da_c210130.nc</map>

‎components/elm/bld/namelist_files/namelist_definition.xml‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1508,7 +1508,7 @@ If TRUE, irrigation will be active (find surface datasets with active irrigation
15081508

15091509
<entry id="sim_year" type="integer" category="default_settings"
15101510
group="default_settings" valid_values=
1511-
"1000,850,1100,1350,1600,1850,1855,1865,1875,1885,1895,1905,1915,1925,1935,1945,1950,1955,1965,1975,1985,1995,2000,2005,2010,2015,2025,2035,2045,2055,2065,2075,2085,2095,2100,2105">
1511+
"1000,850,1100,1350,1600,1850,1855,1865,1875,1885,1895,1905,1915,1925,1935,1945,1950,1955,1965,1975,1980,1985,1995,2000,2005,2010,2015,2025,2035,2045,2055,2065,2075,2085,2095,2100,2105">
15121512
Year to simulate and to provide datasets for (such as surface datasets, initial conditions, aerosol-deposition, Nitrogen deposition rates etc.)
15131513
A sim_year of 1000 corresponds to data used for testing only, NOT corresponding to any real datasets.
15141514
A sim_year greater than 2005 corresponds to rcp scenario data

‎components/elm/docs/user-guide/index.md‎

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -76,3 +76,11 @@ Using the above-mentioned settings:
7676
[FATES](fates.md) can be run in various modes with ELM through the use of namelist settings.
7777
The [FATES User's Guide section on namelist options](https://fates-users-guide.readthedocs.io/en/latest/user/Namelist-Options-and-Run-Time-Modes.html)
7878
provides guidance on enabling these different FATES run modes.
79+
80+
## Create land surface dataset
81+
82+
A new surface dataset for ELM is generated using `mksurfdata_map` and the notes about it are available [here](surface_dataset.md)
83+
84+
## Generate land initial condition
85+
86+
Initial ELM condition can be generated using `interpinic` and the notes about it are available [here](interpinic.md).

‎components/elm/docs/user-guide/interpinic.md‎

Lines changed: 79 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,79 @@
1+
# Creating an ELM initial condition file
2+
3+
An ELM initial condition (IC) file can be created by remapping an existing IC file from
4+
one resolution to another using the `interpinic`, located at
5+
`components/elm/tools/interpinic`. An ELM IC file is in the same format as an ELM restart file.
6+
The composet of the remapped IC file will be the same as that of the input IC file.
7+
So, for a new ELM SP-mode IC file, use an ELM input file corresponding to the SP-mode.
8+
9+
The steps involved in creating a new IC files are as follows:
10+
11+
1. Identifying an input ELM IC or restart file that will be remapped.
12+
2. Obtaining an ELM restart file at the new resolution.
13+
3. Compiling `interpinic` on the machine of interest.
14+
4. Running `interpinic` to perform the interpolation.
15+
16+
The notes below provide an example of creating 1850 ELM IC file for the NARRM grid using E3SM v3 LR piControl from year = 0101. These notes are provided for Chrysalis.
17+
18+
## 1. Identification of the input ELM IC file
19+
20+
The identified input land condition file for this case is the following:
21+
22+
```bash
23+
/lcrc/group/e3sm2/ac.golaz/E3SMv3/v3.LR.piControl/archive/rest/0101-01-01-00000/v3.LR.piControl.elm.r.0101-01-01-00000.nc
24+
```
25+
26+
## 2. Obtaining an ELM restart file
27+
28+
Using an existing NARRM land IC and making a copy of it
29+
30+
```bash
31+
cd components/elm/tools/interpinic
32+
33+
cp /lcrc/group/e3sm/data/inputdata/lnd/clm2/initdata_map/elmi.v3-NARRM.northamericax4v1pg2_r025_IcoswISC30E3r5.1870-01-01-00000.c20240704.nc \
34+
elmi.v3-NARRM.northamericax4v1pg2_r025_IcoswISC30E3r5.1850-01-01-00000.c`date "+%Y%m%d"`.nc
35+
```
36+
37+
## 3. Compiling `interpinic`
38+
39+
```bash
40+
# Load relevant modules
41+
cd <e3sm-dir>
42+
eval $(./cime/CIME/Tools/get_case_env)
43+
44+
# change directory
45+
cd components/elm/tools/interpinic/src
46+
47+
export USER_LDFLAGS="-L$NETCDF_C_DIR/lib -lnetcdf -L$NETCDF_F_DIR/lib -lnetcdff -L$HDF5_DIR/lib -lhdf5"
48+
49+
USER_FC=ifort LIB_NETCDF="`nc-config --flibs`" INC_NETCDF="`nf-config --includedir`" make VERBOSE=1
50+
51+
cd ../
52+
```
53+
54+
## 4. Run `interpinic`
55+
56+
The `interpinic` can then be run via the following batch job (e.g., `remap.r025_RRSwISC6to18E3r4.1850.batch`) to generate the initial condition.
57+
58+
```bash
59+
>cat remap.r025_RRSwISC6to18E3r4.1850.batch
60+
61+
#!/bin/sh
62+
#SBATCH --job-name=remap
63+
#SBATCH --nodes=1
64+
#SBATCH --exclusive
65+
#SBATCH --time 24:00:00
66+
#SBATCH -p slurm
67+
#SBATCH --account esmd
68+
69+
# Load relevant modules.
70+
cd <e3sm-dir>
71+
eval $(./cime/CIME/Tools/get_case_env)
72+
73+
# Change dir to `interpinic`
74+
cd components/elm/tools/interpinic/
75+
76+
srun -n 1 ./interpinic \
77+
-i /lcrc/group/e3sm2/ac.golaz/E3SMv3/v3.LR.piControl/archive/rest/0101-01-01-00000/v3.LR.piControl.elm.r.0101-01-01-00000.nc \
78+
-o elmi.v3-NARRM.northamericax4v1pg2_r025_IcoswISC30E3r5.1850-01-01-00000.c20240903.nc
79+
```

‎components/elm/docs/user-guide/surface_dataset.md‎

Lines changed: 56 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,56 @@
1+
# Creating an ELM surface dataset
2+
3+
The notes describe the steps in creating an ELM surface dataset at 0.5x0.5 resolution for 1950 on Perlmutter.
4+
5+
## 1. Load the appropriate modules
6+
7+
```bash
8+
cd <e3ms-dir>
9+
eval $(./cime/CIME/Tools/get_case_env)
10+
```
11+
12+
## 2. Compile `mksurfdata_map`
13+
14+
```bash
15+
cd components/elm/tools/mksurfdata_map/src/
16+
17+
make clean
18+
export USER_LDFLAGS="-L$NETCDF_DIR/lib -lnetcdf -lnetcdff -lnetcdf_intel"
19+
export USER_LDFLAGS=$USER_LDFLAGS" -L$HDF5_DIR/lib -lhdf5 -lhdf5_fortran -lhdf5_hl_intel -lhdf5hl_fortran_intel"
20+
21+
USER_FC=ifort LIB_NETCDF="`nc-config --flibs`" INC_NETCDF="`nf-config --includedir`" make VERBOSE=1
22+
```
23+
24+
## Build the namelist
25+
26+
This step assumes that the resolution for which the new surface dataset is being created is a supported resolution.
27+
If the surface dataset is being created for an unsupported resolution, 16 mapping files will have to be created to map the raw datasets
28+
onto this unsupported resolution. The `namelist` file with default number of glaciers (equal to zero) can be generated as:
29+
30+
```bash
31+
cd ../
32+
33+
RES=0.5x0.5
34+
YR=1950
35+
DIN_LOC_ROOT=/global/cfs/cdirs/e3sm/inputdata
36+
37+
./mksurfdata.pl -res $RES -years $YR -d -dinlc $DIN_LOC_ROOT
38+
```
39+
40+
An example of generating the namelist for 0.25 deg (`r025`) resolution for 1980 with 10 glacier layers is as follows:
41+
42+
```bash
43+
RES=r025
44+
YR=1980
45+
DIN_LOC_ROOT=/global/cfs/cdirs/e3sm/inputdata
46+
47+
./mksurfdata.pl -res $RES -years $YR -d -dinlc $DIN_LOC_ROOT -glc_nec 10
48+
```
49+
50+
## Run `mksurfdata_map` via an interactive job
51+
52+
```bash
53+
salloc --nodes 1 --qos interactive --time 01:00:00 --constraint cpu --account e3sm
54+
55+
srun -n 1 ./mksurfdata_map < namelist
56+
```

0 commit comments

Comments
 (0)