<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN">
<html><body>
<p>Dear Denis,</p>
<p>Thank you for the answer.</p>
<p>I tried to make in another way. I interpolated from hybrid-sigma levels to pressure using the NCL script vinth2p_ecmwf.ncl, that is for CESM hybrid coordinates to pressure coordinates but uses an ECMWF formulation to extrapolate values below ground (https://bb.cgd.ucar.edu/converting-hybrid-levels-pressures). This interpolation was good, so I tried to do the interpolation from pressure to height, like your suggestions in this post: https://www.ncl.ucar.edu/Support/talk_archives/2004/0394.html<br /><br />And the interpolation was not good. Is there any NCL script to interpolate wind from pressure levels to 10m height? I need the wind in 10m just over the ocean.</p>
<p><span id="result_box" class="short_text"><span>My script is attached.</span></span></p>
<p><span class="short_text"><span>Best regards.<br /></span></span></p>
<p><span class="short_text"><span>Leilane</span></span></p>
<p>&nbsp;</p>
<p>Em 22-05-2017 16:49, Dennis Shea escreveu:</p>
<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px; width:100%"><!-- html ignored --><!-- head ignored --><!-- meta ignored -->
<div dir="ltr">
<div>The wrf_* functions were developed for the WRF model(s) with the WRF 'eta' vertical coordinate system.<br />[ <a href="https://www.researchgate.net/post/Does_WRF_simulation_use_sigma_or_eta_vertical_coordinate_system">https://www.researchgate.net/post/Does_WRF_simulation_use_sigma_or_eta_vertical_coordinate_system</a> ]<br /><br /></div>
<div>Hence, these functions are not appropriate for CAM data.</div>
<div>---</div>
<div>The calculation of CAM U10 is non-trivial. Many factors are involved: surface roughness, atmospheric stability, ...<br /><br /></div>
<div>Over the oceans, for neutral stability something like the following is used:</div>
<div><br />&nbsp;&nbsp;&nbsp; U10 = UL / [ 1 + (sqrt(CN10)/k)*log(ZL/10-PSI)]&nbsp;&nbsp; ; &nbsp; [1] <br /><br /></div>
<div>The following variables are readily available at each grid point</div>
<div>&nbsp;</div>
<div>UL&nbsp; - zonal component of the wind at the Lowest model level<br />ZL&nbsp; - geopotential height at the Lowest model level</div>
<div>PSI - Sfc. geopotential<br /><br /></div>
<div>&nbsp;The CN10 must be derived:<br /><br /></div>
<div>&nbsp;&nbsp; CN10 - drag Coefficient at 10-m height and Neutral stability
<div style="font-size: 19.9253px; font-family: sans-serif;">&nbsp;</div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;">where:</span></div>
<div style="font-family: serif;"><span style="font-size: small;">&nbsp;&nbsp;&nbsp; CN10 = (C4/U10) + C5 + C6*U10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; [2]<br /><br /></span></div>
<div style="font-family: serif;"><span style="font-size: small;">and, over the ocean,<br /></span></div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;">&nbsp;&nbsp;&nbsp; C4=0.0027 m/s , C5=0.000142, C6=0.0000764 m/s<br /><br /></span></div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;">Note that U10 appears in both [1] and [2]. Hence the solution requires iteration.<br />----<br /></span></div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;">Over land, ????? <br />---<br /></span></div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;">If you need 'numbers' as a filler then maybe using a simple&nbsp; u10 = U*log(..) could be used. <br /><br /></span></div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;">Attached is a VERY SIMPLE function. The function requires 'surface roughness length' which can be a function of time (season).&nbsp; This information may be available online.<br /><br /></span></div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;">Also, as noted within the function, the approximation is most valid in the lowest 100m of the atmosphere.<br /></span></div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;">&nbsp;</span></div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;">Good Luck<br /></span></div>
<div style="font-size: 19.9253px; font-family: sans-serif;"><span style="font-size: small;"><br /><br /></span></div>
</div>
</div>
<div class="gmail_extra"><br />
<div class="gmail_quote">On Mon, May 15, 2017 at 1:09 PM, <span>&lt;<a href="mailto:leilane.passos@inpe.br">leilane.passos@inpe.br</a>&gt;</span> wrote:<br />
<blockquote class="gmail_quote" style="margin: 0  0  0  .8ex; border-left: 1px  #ccc  solid; padding-left: 1ex;"><span style="text-decoration: underline;"></span>
<div>
<p>Hi,</p>
<p>I need to interpolate the wind results from CCSM4(CAM) in hybrid sigma level coordinate to height of 10 m. Then, I tried to use the function wrf_user_intrp3d with the u, v and the geopotential height (Z3), but the interpolated field has undefined values and I don't know the reason. I have the following questions:</p>
<p>1) I think that I should find the interpolated values in the entire domain. But why did I find undefined values in certain areas like Andes Chain?</p>
<p>2) If I use the command vinth2p to interpolate from hybrid to pressure coordinates, what is a similar function that I could use to convert from pressure to height of 10 m?</p>
<p>&nbsp;</p>
<p>; --------------&nbsp; LOAD FUNCTIONS AND PROCEDURES ----------------<br /><br />&nbsp;&nbsp;&nbsp; load "/usr/share/ncarg/nclscripts/csm/gsn_code.ncl"<br />&nbsp;&nbsp;&nbsp; load "/usr/share/ncarg/nclscripts/csm/gsn_csm.ncl"<br />&nbsp;&nbsp;&nbsp; load "/usr/share/ncarg/nclscripts/csm/contributed.ncl"<br />&nbsp;&nbsp;&nbsp; load "/usr/share/ncarg/nclscripts/csm/shea_util.ncl"<br />&nbsp;&nbsp;&nbsp; load "/usr/share/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br />&nbsp;&nbsp;&nbsp; load "/usr/share/ncarg/nclscripts/acentos.ncl"<br /><br />&nbsp;&nbsp;&nbsp; ; --------------&nbsp; BEGINING OF NCL SCRIPT ----------------<br /><br />&nbsp;&nbsp;&nbsp; begin<br /><br />&nbsp;&nbsp;&nbsp; outfile = "<a href="http://test_U10mV10m.1985.20thC.nc">test_U10mV10m.1985.20thC.nc</a>"<br /><br />&nbsp;&nbsp;&nbsp; if (isfilepresent(outfile)) then<br />&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; system("rm -rf "+outfile)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- make sure that file does not exist<br />&nbsp;&nbsp;&nbsp; end if<br /><br />&nbsp;&nbsp;&nbsp; ;-- open data file<br /><br />&nbsp;&nbsp;&nbsp; in_u = addfile("<a href="http://224331.U.1985.20thC.U.nc">224331.U.1985.20thC.U.nc</a>","r")<br />&nbsp;&nbsp;&nbsp; in_v = addfile("<a href="http://224331.V.1985.20thC.V.nc">224331.V.1985.20thC.V.nc</a>","r")<br />&nbsp;&nbsp;&nbsp; in_z = addfile("<a href="http://239928.Z3.1985.20thC.Z3.nc">239928.Z3.1985.20thC.Z3.nc</a>","r")<br /><br />&nbsp;&nbsp;&nbsp; ;-- get variable t and its dimensions and dimension sizes<br /><br />&nbsp;&nbsp;&nbsp; u&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; in_u-&gt;U ;(time, lev, lat, lon)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get variable u<br />&nbsp;&nbsp;&nbsp; v&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; in_v-&gt;V ;(time, lev, lat, lon)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get variable v<br />&nbsp;&nbsp;&nbsp; z&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; in_z-&gt;Z3 ;(time, lev, lat, lon)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get variable u<br /><br />&nbsp;&nbsp;&nbsp; time&nbsp; =&nbsp; in_u-&gt;time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get dimension time<br />&nbsp;&nbsp;&nbsp; lev&nbsp;&nbsp; =&nbsp; in_u-&gt;lev&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get dimension lev<br />&nbsp;&nbsp;&nbsp; lat&nbsp;&nbsp; =&nbsp; in_u-&gt;lat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get dimension lat<br />&nbsp;&nbsp;&nbsp; lon&nbsp;&nbsp; =&nbsp; in_u-&gt;lon&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get dimension lon<br /><br />&nbsp;&nbsp;&nbsp; ntim&nbsp; =&nbsp; dimsizes(time)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get dimension sizes of time<br />&nbsp;&nbsp;&nbsp; nlev&nbsp; =&nbsp; dimsizes(lev)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get dimension sizes of lev<br />&nbsp;&nbsp;&nbsp; nlat&nbsp; =&nbsp; dimsizes(lat)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get dimension sizes of lat<br />&nbsp;&nbsp;&nbsp; nlon&nbsp; =&nbsp; dimsizes(lon)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- get dimension sizes of lon<br /><br />&nbsp;&nbsp;&nbsp; ;--&nbsp; Interpolate u and v to 10m:<br />&nbsp;&nbsp;&nbsp; ;u10m&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = u&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- copy variable and its dimensions and attributes<br />&nbsp;&nbsp;&nbsp; u10m&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = wrf_user_intrp3d(u,z,"h",10.0,0.,False)&nbsp;&nbsp; ;-- interpolate u<br />&nbsp;&nbsp;&nbsp; u10m@units&nbsp; = "m/s"&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- define new units<br />&nbsp;&nbsp;&nbsp; ;v10m&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = v&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- copy variable and its dimensions and attributes<br />&nbsp;&nbsp;&nbsp; v10m&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = wrf_user_intrp3d(v,z,"h",10.0,0.,False) ;-- interpolate v<br />&nbsp;&nbsp;&nbsp; v10m@units&nbsp; = "m/s"&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- define new units<br />&nbsp;&nbsp;&nbsp; ;print(u10m)<br /><br />&nbsp;&nbsp;&nbsp; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br />&nbsp;&nbsp;&nbsp; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br />&nbsp;&nbsp;&nbsp; ;-- create new netCDF file<br />&nbsp;&nbsp;&nbsp; fout = addfile(outfile,"c")<br /><br />&nbsp;&nbsp;&nbsp; ;-- begin output file settings<br />&nbsp;&nbsp;&nbsp; setfileoption(fout,"DefineMode",True) ;-- explicitly declare file definition mode<br /><br />&nbsp;&nbsp;&nbsp; ;-- create global attributes of the file<br />&nbsp;&nbsp;&nbsp; fAtt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- assign file attributes<br />&nbsp;&nbsp;&nbsp; fAtt@title&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "NCL Efficient Approach to netCDF Creation"<br />&nbsp;&nbsp;&nbsp; fAtt@source_file&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "CCSM4.nc"<br />&nbsp;&nbsp;&nbsp; fAtt@Conventions&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "CF"<br />&nbsp;&nbsp;&nbsp; fAtt@creation_date&nbsp;&nbsp;&nbsp; =&nbsp; systemfunc ("date")<br />&nbsp;&nbsp;&nbsp; fAtt@history&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; "NCL script: ex_netcdf_2.ncl"<br />&nbsp;&nbsp;&nbsp; fAtt@comment&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "Convert variable t: Kelvin to Celsius"<br />&nbsp;&nbsp;&nbsp; fileattdef(fout,fAtt)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- copy file attributes<br /><br />&nbsp;&nbsp;&nbsp; ;-- predefine the coordinate variables and their dimensionality<br />&nbsp;&nbsp;&nbsp; dimNames = (/"time", "lat", "lon"/)<br />&nbsp;&nbsp;&nbsp; dimSizes = (/ -1&nbsp;&nbsp; , nlat,&nbsp; nlon /)<br />&nbsp;&nbsp;&nbsp; dimUnlim = (/ True , False, False/)<br />&nbsp;&nbsp;&nbsp; filedimdef(fout,dimNames,dimSizes,dimUnlim)<br /><br />&nbsp;&nbsp;&nbsp; ;-- predefine the the dimensionality of the variables to be written out<br />&nbsp;&nbsp;&nbsp; filevardef(fout, "time" ,typeof(time),getvardims(time))<br />&nbsp;&nbsp;&nbsp; filevardef(fout, "lat"&nbsp; ,typeof(lat), getvardims(lat))<br />&nbsp;&nbsp;&nbsp; filevardef(fout, "lon"&nbsp; ,typeof(lon), getvardims(lon))<br />&nbsp;&nbsp;&nbsp; filevardef(fout, "u10m" ,typeof(u10m),&nbsp; getvardims(u10m))<br />&nbsp;&nbsp;&nbsp; filevardef(fout, "v10m" ,typeof(v10m),&nbsp; getvardims(v10m))<br /><br />&nbsp;&nbsp;&nbsp; ;-- copy attributes associated with each variable to the file<br />&nbsp;&nbsp;&nbsp; filevarattdef(fout,"time" ,time)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- copy time attributes<br />&nbsp;&nbsp;&nbsp; filevarattdef(fout,"lat"&nbsp; ,lat)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- copy lat attributes<br />&nbsp;&nbsp;&nbsp; filevarattdef(fout,"lon"&nbsp; ,lon)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- copy lon attributes<br />&nbsp;&nbsp;&nbsp; filevarattdef(fout,"u10m" ,u10m)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- copy u10m attributes<br />&nbsp;&nbsp;&nbsp; filevarattdef(fout,"v10m" ,v10m)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- copy v10m attributes<br /><br />&nbsp;&nbsp;&nbsp; ;-- explicitly exit file definition mode (not required)<br />&nbsp;&nbsp;&nbsp; setfileoption(fout,"DefineMode",False)<br /><br />&nbsp;&nbsp;&nbsp; ;-- output only the data values since the dimensionality and such have been predefined.<br />&nbsp;&nbsp;&nbsp; ;-- The "(/", "/)" syntax tells NCL to only output the data values to the predefined<br />&nbsp;&nbsp;&nbsp; ;-- locations on the file.<br />&nbsp;&nbsp;&nbsp; fout-&gt;time&nbsp;&nbsp; =&nbsp; (/time/)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- write time to new netCDF file<br /><br />&nbsp;&nbsp;&nbsp; fout-&gt;lat&nbsp;&nbsp;&nbsp; =&nbsp; (/lat/)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- write lat to new netCDF file<br />&nbsp;&nbsp;&nbsp; fout-&gt;lon&nbsp;&nbsp;&nbsp; =&nbsp; (/lon/)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- write lon to new netCDF file<br />&nbsp;&nbsp;&nbsp; fout-&gt;u10m&nbsp;&nbsp; =&nbsp; (/u10m/)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- write variable to new netCDF file<br />&nbsp;&nbsp;&nbsp; fout-&gt;v10m&nbsp;&nbsp; =&nbsp; (/v10m/)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;-- write variable to new netCDF file<br /><br />end<br /><br /></p>
<p>&nbsp;</p>
<p>&nbsp;</p>
<p>Best regards,</p>
<p>Leilane</p>
<div>&nbsp;</div>
</div>
<br />_______________________________________________<br /> ncl-talk mailing list<br /><a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br /> List instructions, subscriber options, unsubscribe:<br /><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br /><br /></blockquote>
</div>
</div>
</blockquote>
<p>&nbsp;</p>
<div>&nbsp;</div>
</body></html>