...
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
...