<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 id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hello, </p>
<p><br>
</p>
<p>I'm trying to re-grid velocity data from 1km resolution to 5km resolution (2801x1501 to 560x300 resolution) using linint2. </p>
<p><br>
</p>
<p>Right now the script runs and the linint2 function appears to work at first- I attached a screenshot of the script output, including the printVarSummary for the re-gridded velocity (Velmeas), which shows the dimensions are 560x300. However, the plot shows
that instead of decreasing the resolution and showing velocity over all of Greenland, Velmeas appears to have stayed at the same resolution and zoomed into 1:500 x 1:300 of the original 2801x1501 resolution.
<br>
</p>
<p><br>
</p>
<p>While this data is a map of velocity, there's no lat lon data associated with it, so I was trying to simply treat it as an array.
<br>
</p>
<p><br>
</p>
<p>I'm running NCL version 6.3.0 on Yellowstone </p>
<p><br>
</p>
<p>Here's the script: <br>
</p>
<p> g = addfile("Velocity_data.nc","r")<br>
vx = g->vx(0,:,:) ; read in velocity data<br>
vy = g->vy(0,:,:) ; read in velocity data<br>
printVarSummary(vx) <br>
printVarSummary(vy)<br>
Velm = ((vx^2) + (vy^2))^0.5 ; calculate absolute velocity<br>
printVarSummary(Velm)<br>
<br>
xi = ispan(1,1501,1) ; original x resolution<br>
yi = ispan(1,2801,1) ; original y resolution<br>
xo = ispan(1,300,1) ; new x resolution<br>
yo = ispan(1,560,1) ; new y resolution <br>
Velmeas = linint2(xi,yi,Velm,False,xo,yo,0) ; Re-grid<br>
printVarSummary(Velmeas)<br>
<br>
wks = gsn_open_wks("x11","test_Vel")<br>
plot = gsn_contour(wks,Velmeas,False) </p>
<p><br>
</p>
<p><br>
</p>
<p>Thanks for any help, </p>
<p>Lauren <br>
</p>
<p><br>
</p>
<p><br>
</p>
</div>
</body>
</html>