Abstract: There is a small task out of LLNL LDRD-SI ``Climate Resilience for National Security’’ — we are curious if we can run SCREAM RRM at 100m. However, the existing USGS cubed topo is all in 3 km processed from the USGS GTOPO30 30 arch-second (1 km) source DEM (10.5194/gmd-8-3975-2015). Therefore, we need to generate new cubed topo using ultra-high resolution source DEM data. Some of us are interested in the process of generating new global uniform cubed topo used for whatever new grids with dx < 3 km @MarkT Mark Taylor, @PeterB Peter Bogenschutz, @Jishi Jishi Zhang. Specifically, this page documents the practice on LC of generating 800m cubed topo using 500m / 250m USGS GMTED2010 source DEM dataset.
...
Other processing: a) no polar regions, stitched from 500m GMTED2010 DEM (first remapped to 250m resolution); b) fill Greenland from 500m GMTED2010 DEMl; c) use landfract interpolated from 500m GMTED2010 DEM.
Code Block |
---|
cdo --no_remap_weights remapbil,r172800x86400 S5P_OPER_REF_DEM_15_00000000T000000_99999999T999999_20160111T104226.NCL_24-3.nc S5P_OPER_REF_DEM_15_00000000T000000_99999999T999999_20160111T104226.NCL_24-3.r172800x86400.nc |
grid_WL.topo.DEM.src_glb.stitch.ncl
Code Block |
---|
dir="/pscratch/sd/z/zhang73/DATA/DEM/" base="S5P_OPER_REF_DEM_15_00000000T000000_99999999T999999_20160111T104226.NCL_24-3.r172800x86400.nc" target="be75_grd_1.nc" a=addfile(dir+base,"r") b=addfile(dir+target,"r") htopo=a->htopo htopo_stitch=htopo htopo_stitch(16320:83519,0:86399)=(/b->Band1(:,86400:172799)/) print("1 timer: "+get_cpu_time()) htopo_stitch(16320:83519,86400:172799)=(/b->Band1(:,0:86399)/) print("2 timer: "+get_cpu_time()) htopo_stitch(69436:83238,137136:168188)=\ where(htopo(69436:83238,137136:168188).ne.0.and.htopo_stitch(69436:83238,137136:168188).eq.0, \ htopo(69436:83238,137136:168188), htopo_stitch(69436:83238,137136:168188)) ; ************************************************************************************************ setfileoption("nc", "Format", "LargeFile") diro=dir filo = "GMTED2010_7.5_stitch_S5P_OPER_REF_DEM_15_NCL_24-3.r172800x86400.nc" system("if ! test -f " + diro+filo +" ; then cp "+ diro+base +" "+ diro+filo +" ; fi") fout = addfile (diro + filo, "w") fout->htopo = htopo_stitch |
Step 2. bin_to_cube:
build issues
...
250m source DEM to cubed topo failed due to OOM. It requires ~638G to cube48000 (200m), ~507G to cube40000 (250m), ~316G to cube24000 (400m), not enough for both LC and perlmutter.
...