<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I'm attempting to generate weights for regridding a 4 km grid to a 14 km grid using ESMF_regrid_gen_weights.
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Source grid:</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Univ. Arizon Snow Depth and SWE grids<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
4 km rectilinear, [lat | 621] x [lon |1405]</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Destination grid:</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
VR-CESM CONUS (already regridded using ESMF)</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
14 km rectilinear, [lat | 361] x [lon |560]</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
In the below NCL script, I use rectilinear_to_SCRIP on source and destination grids, then ESMF_regrid_gen_weights to create the weight file. No errors occur, however, when I try to regrid using ncremap, I get a bunch of zeros. I'd like to use the weight file
on hundreds of files, hence using ncremap. <br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
When I test the weight file on one source variable using NCL's ESMF_regrid_with_weights, I get an argument type mismatch error, but I can't figure out why. Any ideas how to fix this? Files are on cheyenne.
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Script to generate weights:</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>;======================================================================<br>
</span>
<div>; ESMF Regrid: UA-SNOW-4km (rectilinear) to VRCESM-CONUS (rectilinear)<br>
</div>
<div>;======================================================================<br>
</div>
<div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
</div>
<div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
</div>
<div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>
</div>
<div>load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"<br>
</div>
<div><br>
</div>
<div>;======================================================================<br>
</div>
<div>; The main code<br>
</div>
<div>;======================================================================<br>
</div>
<div><br>
</div>
<div>begin<br>
</div>
<div>;---Specify interpolation to be used<br>
</div>
<div> method = "patch"<br>
</div>
<div><br>
</div>
<div>;---Specify weight file name<br>
</div>
<div> wgtFileName = "UASnow4km2VRCESMCONUS14km_"+method+".nc"<br>
</div>
<div><br>
</div>
<div>;---Input file (Source)<br>
</div>
<div> srcFileDir = "/glade/work/burakows/plot/validation/UA_Snow_4km/"<br>
</div>
<div> srcFileName = "4km_SWE_Depth_WY2014_v01.nc" ; any UA-Snow file<br>
</div>
<div> srcGridName = "UASnow_4km_SCRIP.nc"<br>
</div>
<div><br>
</div>
<div>;---Retrieve a sample variable (used in plot)<br>
</div>
<div> sfile = addfile(srcFileDir+srcFileName,"r")<br>
</div>
<div> srcVarName = "SWE" ; [lat | 621] x [x | 1405]<br>
</div>
<div> var_in = sfile->$srcVarName$<br>
</div>
<div> printVarSummary(var_in)<br>
</div>
<div><br>
</div>
<div> Opt = True<br>
</div>
<div> Opt@ForceOverwrite = True<br>
</div>
<div> Opt@Title = "UA Snow 4km Grid"<br>
</div>
<div><br>
</div>
<div>;--- Convert UA Snow 4 km grid to SCRIP convention file<br>
</div>
<div> rectilinear_to_SCRIP(srcGridName,var_in&lat,var_in&lon,Opt)<br>
</div>
<div> delete(Opt)<br>
</div>
<div><br>
</div>
<div>;================================================================<br>
</div>
<div><br>
</div>
<div>;---File that contains the destination file coordinates<br>
</div>
<div> dstFileDir = "/glade/scratch/burakows/archive/VR-CESM_FHISTSTUB.1984-2015.CONUS30x4_prod/lnd/remap/"<br>
</div>
<div> dstFileName = "VR-CESM_FHISTSTUB.1984-2015.CONUS30x4_prod.clm2.h0.2014-01-01-00000_regrid.nc" ; any VR-CESM file<br>
</div>
<div> dstGridName = "VRCESMCONUS_14km_SCRIP.nc"<br>
</div>
<span></span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>;---open VR-CESM file that contains the VR-CESM grid coordinates<br>
</span>
<div> dfile = addfile(dstFileDir+dstFileName,"r")<br>
</div>
<div> dstVarName = "H2OSNO"<br>
</div>
<div> dvar_in = dfile->$dstVarName$ ; [lat | 361 ] x [lon | 560]<br>
</div>
<div> printVarSummary(dvar_in)<br>
</div>
<div><br>
</div>
<div> Opt = True<br>
</div>
<div> Opt@ForceOverwrite = True<br>
</div>
<div> Opt@Title = "VR-CESM CONUS (14 km) grid"<br>
</div>
<div><br>
</div>
<div>;--- Convert VR-CESM CONUS 14 km to SCRIP convention file<br>
</div>
<div> rectilinear_to_SCRIP(dstGridName,dvar_in&lat,dvar_in&lon,Opt)<br>
</div>
<div> delete(Opt)<br>
</div>
<div><br>
</div>
<div><span>;--- Generate the weights to convert UA Snow 4 km grid to VR-CESM CONUS 14 km grid<br>
</span>
<div><br>
</div>
<div> Opt = True<br>
</div>
<div> Opt@InterpMethod = "patch"<br>
</div>
<div> Opt@ForceOverwrite = "True"<br>
</div>
<div><br>
</div>
<div> ESMF_regrid_gen_weights(srcGridName,dstGridName,wgtFileName,Opt)<br>
</div>
<div> delete(Opt)<br>
</div>
<span></span>end</div>
<div><br>
</div>
<span>;===============================================================</span><br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="Signature">
<div></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font face="Corbel, sans-serif">Dr. Elizabeth A Burakowski</font></div>
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<font face="Corbel, sans-serif">She/Her/Hers<br>
</font>
<div><font face="Corbel, sans-serif">Research Assistant Professor</font></div>
<div><font face="Corbel, sans-serif">Earth Systems Research Center & </font></div>
<div><font face="Corbel, sans-serif">Institute for the Study of Earth, Oceans, and Space</font></div>
<div><font face="Corbel, sans-serif">University of New Hampshire</font></div>
<div><font face="Corbel, sans-serif">8 College Road Durham, NH 03824</font></div>
<div><font face="Corbel, sans-serif">603-862-1796</font></div>
<div><font face="Corbel, sans-serif">elizabeth.burakowski@unh.edu</font></div>
</div>
</div>
</body>
</html>