<div dir="ltr">Hi Mary & Guilherme,<div><br></div><div>Thanks for the replies. When I run 'make' in my cnvgrid-3.1.0 directory I get the following error:</div><div><br></div><div><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">In file pds2pdt.f:50</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);min-height:13px"><span style="font-variant-ligatures:no-common-ligatures"></span><br></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> use params</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"> 1</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">Fatal Error: Can't open module file 'params.mod' at (1) for reading: No such file or directory</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small">I think the ultimate problem lies with my not being able to install g2lib; when I attempt to do so, I get the error:</span><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">m2c -o g2grids.o g2grids.mod</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">
</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">make: m2c: No such file or directory</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small">and inserting the line %.o : %.mod results in a different error:</span><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">make: *** No rule to make target `libg2.a(g2grids.o)', needed by `libg2.a'. Stop.</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small">The problem is probably simple, but I've yet to make any headway on these errors.</span><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small">Regarding wgrib2, I found examples and documentation of using it to convert grib1 to grib2, but not vice versa. If you could kindly provide guidance on how the conversion to grib1 could be done using wgrib2 (if possible), I'd greatly appreciate it.</span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small">Best regards,</span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small">Colin</span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small"><br></span></span></p></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Mar 16, 2018 at 5:58 AM, <span dir="ltr"><<a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send ncl-talk mailing list submissions to<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:ncl-talk-request@ucar.edu">ncl-talk-request@ucar.edu</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:ncl-talk-owner@ucar.edu">ncl-talk-owner@ucar.edu</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of ncl-talk digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: Missing values after ESMF regridding (Mary Haley)<br>
2. Re: Plotting netCDF (Dave Allured - NOAA Affiliate)<br>
3. Unable to set explicit levels on contour plot<br>
(Nicholas Thomas Luchetti)<br>
4. Re: Unable to set explicit levels on contour plot (Adam Phillips)<br>
5. Re: Possible to convert grib2 or netcdf files to grib1?<br>
(Mary Haley)<br>
6. Re: Possible to convert grib2 or netcdf files to grib1?<br>
(Guilherme Martins)<br>
7. Regridding Moving Nest WRF outputs (Lara Quiti?n Hern?ndez)<br>
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------<br>
<br>
Message: 1<br>
Date: Thu, 15 Mar 2018 13:54:35 -0600<br>
From: Mary Haley <<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>><br>
To: Bian Qingyun <<a href="mailto:bianqy@tea.ac.cn">bianqy@tea.ac.cn</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Missing values after ESMF regridding<br>
Message-ID:<br>
<CACNN_CK6jVgR4S=Z=<wbr>d7fwr4MGdw5QUPhbQ_Ho3-yi=Nc6=<a href="mailto:WJEQ@mail.gmail.com">W<wbr>JEQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Bian,<br>
<br>
I'm looking into this.<br>
<br>
Meanwhile, I see that your original data goes from -180 to 180, but the<br>
lower left corner and upper right corner of the destination grid was set to:<br>
<br>
Opt@DstLLCorner = (/ 0.25d, 0.25d/)<br>
Opt@DstURCorner = (/89.75d, 179.75d/)<br>
<br>
This is restricting your data in the longitude direction to go from 0.25<br>
and 179.75. Perhaps you meant for this to be:<br>
<br>
Opt@DstLLCorner = (/ 0.25d, -180.00d/)<br>
Opt@DstURCorner = (/89.75d, 179.75d/)<br>
<br>
However, this still doesn't fix the issue, so I'm looking into it further.<br>
I'm wondering if this is a missing value issue, because I noticed that your<br>
lat/lon has some missing values in it.<br>
<br>
--Mary<br>
<br>
<br>
<br>
On Thu, Mar 15, 2018 at 2:08 AM, Bian Qingyun <<a href="mailto:bianqy@tea.ac.cn">bianqy@tea.ac.cn</a>> wrote:<br>
<br>
><br>
><br>
> Hello,<br>
><br>
> I'm regridding AMSR-E L3 daily snow following the examples on NCL website<br>
> <a href="http://www.ncl.ucar.edu/Applications/Scripts/ESMF_regrid_13.ncl" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/Scripts/ESMF_<wbr>regrid_13.ncl</a> and hdf<br>
> website <a href="http://hdfeos.org/zoo/index_openNSIDC_Examples.php#AMSR_E" rel="noreferrer" target="_blank">http://hdfeos.org/zoo/index_<wbr>openNSIDC_Examples.php#AMSR_E</a>.<br>
><br>
> All is fine,but there are missing values along 90?E and 90?W after<br>
> regridding while there are no missing values in original data along and<br>
> near 90?E and 90?W ?<br>
><br>
> Below is my script. I also attached the AMSR-E data file in case.<br>
><br>
> ----------------------<br>
><br>
> ----------------------<br>
><br>
> load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/esmf/ESMF_<wbr>regridding.ncl"<br>
><br>
> begin<br>
> ;-----read data<br>
> srcFileName = "AMSR_E_L3_DailySnow_V09_<wbr>20100311.hdf"<br>
> ; Source file<br>
> eos_file = addfile(srcFileName + ".he2","r")<br>
> lon2d = eos_file->GridLon_Northern_<wbr>Hemisphere<br>
> lat2d = eos_file->GridLat_Northern_<wbr>Hemisphere<br>
> ; Process NAN (1e51) fillValue for LAMAZ geolocation.<br>
> lon2d=where(lon2d .gt. 1000, -999.0, lon2d)<br>
> lat2d=where(lat2d .gt. 1000, -999.0, lat2d)<br>
> lon2d@_FillValue = -999.0<br>
> lat2d@_FillValue = -999.0<br>
><br>
><br>
> hdf_file = addfile(srcFileName, "r")<br>
> ; Read the dataset.<br>
> hdf_data = hdf_file->SWE_NorthernDaily<br>
><br>
> ; Filter out invalid range values.<br>
> ; See "Table 2. Pixel Values ofr the SWE Feids" from [2].<br>
> hdf_data = where(hdf_data .gt. 240, hdf_data@_FillValue, hdf_data)<br>
><br>
> ; Prepare data for plotting by converting type.<br>
> swe2d = tofloat(hdf_data)<br>
> swe2d@_FillValue = tofloat(hdf_data@_FillValue)<br>
><br>
> ; Multiply by two according to data spec [2].<br>
> swe2d = 2 * swe2d<br>
><br>
> ; You can get the description of data set from the data spec [2].<br>
> swe2d@long_name = "Northern Hemisphere 5-day Snow Water Equivalent<br>
> ~C~ (" + hdf_data@hdf_name + ")"<br>
> swe2d@units = "mm"<br>
> ; Associate data with lat/lon.<br>
> swe2d@lon2d = lon2d<br>
> swe2d@lat2d = lat2d<br>
><br>
> ;-----ESMF regridding<br>
> Opt = True ; Options for regridding<br>
><br>
> Opt@SrcFileName = "EASE_SCRIP.nc" ; Output files<br>
> Opt@DstFileName = "NH_SCRIP.nc"<br>
> Opt@WgtFileName = "EASE_2_NH_patch.nc"<br>
> Opt@ForceOverwrite = True<br>
> Opt@SrcMask2D = where(ismissing(lat2d),0,1)<br>
> Opt@SrcGridLat = lat2d<br>
> Opt@SrcGridLon = lon2d<br>
><br>
> Opt@DstGridType = "0.25deg" ; Destination grid<br>
> Opt@DstTitle = "Northern Hemisphere 0.25 resolution"<br>
> Opt@DstLLCorner = (/ 0.25d, 0.25d/)<br>
> Opt@DstURCorner = (/89.75d, 179.75d/)<br>
><br>
> Opt@InterpMethod = "bilinear" ; Careful! Patch method takes a<br>
> long time<br>
><br>
> Opt@Debug = True<br>
><br>
> swe_regrid = ESMF_regrid(swe2d,Opt) ; Regrid swe<br>
> printVarSummary(swe_regrid)<br>
><br>
> ;-----------------------------<wbr>------------------------------<wbr>-----------<br>
> ; Plotting section<br>
> ;-----------------------------<wbr>------------------------------<wbr>-----------<br>
> wks = gsn_open_wks("x11","ESMF_<wbr>regrid") ; send graphics to PNG<br>
> file<br>
><br>
> res = True ; Plot mods desired.<br>
> res@gsnMaximize = True ; Maximize plot<br>
><br>
> res@gsnDraw = False<br>
> res@gsnFrame = False<br>
><br>
> res@cnFillOn = True ; color plot desired<br>
> res@cnFillPalette = "amwg" ; set color map<br>
> res@cnLinesOn = False ; turn off contour lines<br>
> res@cnLineLabelsOn = False ; turn off contour labels<br>
> res@cnFillMode = "RasterFill" ; turn raster on<br>
><br>
> res@cnLevelSelectionMode= "ExplicitLevels" ; set explicit contour<br>
> levels<br>
> res@cnLevels = (/-300,-250,-200,-150,-100, \<br>
> 0,1,5,10,25,100,200,300,400/)<br>
><br>
> res@lbLabelBarOn = False ; turn on in panel<br>
><br>
> res@trGridType = "TriangularMesh" ; allow missing<br>
> coordinates<br>
><br>
><br>
> res@gsnPolar = "NH" ; specify the hemisphere<br>
> res@mpMinLatF = 35<br>
><br>
> ;---Plot original data.<br>
> res@gsnAddCyclic = False<br>
> res@sfXArray = lon2d<br>
> res@sfYArray = lat2d<br>
> res@tiMainString = "Original EASE grid (" +<br>
> str_join(dimsizes(lat2d),",") + ")"<br>
> res@trGridType = "TriangularMesh"<br>
> res@cnFillMode = "RasterFill"<br>
><br>
> plot_orig = gsn_csm_contour_map_polar(wks,<wbr>swe2d,res)<br>
><br>
> delete(res@sfXArray)<br>
> delete(res@sfYArray)<br>
><br>
> ;---Plot regridded data.<br>
> res@gsnAddCyclic = True<br>
><br>
> dims = tostring(dimsizes(swe_regrid))<br>
> res@tiMainString = "Regridded to 0.25 degree grid (" + \<br>
> str_join(dims," x ") + ")"<br>
><br>
> plot_regrid = gsn_csm_contour_map_polar(wks,<wbr>swe_regrid,res)<br>
><br>
> ;---Compare the plots in a panel<br>
> pres = True<br>
> pres@gsnMaximize = True<br>
> pres@gsnPanelLabelBar = True<br>
> pres@lbLabelFontHeightF = 0.01<br>
> pres@pmLabelBarWidthF = 0.8<br>
><br>
> gsn_panel(wks,(/plot_orig,<wbr>plot_regrid/),(/1,2/),pres)<br>
> end<br>
><br>
> ---------------------------<br>
><br>
> ---------------------------<br>
><br>
> Any suggestions are appreciated! Thank you!<br>
><br>
> Best,<br>
><br>
> Qingyun<br>
><br>
><br>
><br>
> ______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180315/4e73d5fc/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180315/4e73d5fc/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 15 Mar 2018 14:09:27 -0600<br>
From: Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov">dave.allured@noaa.gov</a>><br>
To: "Smith, Stuart" <<a href="mailto:smit1770@purdue.edu">smit1770@purdue.edu</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Plotting netCDF<br>
Message-ID:<br>
<<a href="mailto:CALqwTFNC1uoYD7t1g71g-LLqCBWadFAot5H6NkSLWe3dUNjz3A@mail.gmail.com">CALqwTFNC1uoYD7t1g71g-<wbr>LLqCBWadFAot5H6NkSLWe3dUNjz3A@<wbr>mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Stuart,<br>
<br>
When reading the whole array, subscripts are not needed. Just use whole<br>
array syntax:<br>
<br>
x = f->ppt<br>
<br>
You did not fix the mismatched _FillVallue problem that I mentioned<br>
before. You will need to either fix the file, or else add this kludge to<br>
make NCL properly recognize the missing values after they are read in:<br>
<br>
x@_FillValue := todouble (x@missing_value)<br>
<br>
With that fix, it looks like your <a href="http://ppt.nc" rel="noreferrer" target="_blank">ppt.nc</a> data file is pre-masked for a<br>
predetermined spatial area. To calculate a simple spatial average time<br>
series for all grid points within that predetermined area, just use a<br>
simple unweighted array averaging function:<br>
<br>
dim_select = (/ 1, 2 /) ; select lat and lon dims to average over<br>
spatial_aves = dim_avg_n_Wrap (x, dim_select) ; makes a 1-D time<br>
series<br>
<br>
Please read the function documentation for dim_avg_n and dim_avg_n_Wrap,<br>
for more details.<br>
<br>
If you need latitude weighting, or you need to select a smaller spatial<br>
area within the existing domain, then more code will be needed.<br>
<br>
--Dave<br>
<br>
<br>
On Wed, Mar 14, 2018 at 8:06 AM, Smith, Stuart <<a href="mailto:smit1770@purdue.edu">smit1770@purdue.edu</a>> wrote:<br>
<br>
> Dave,<br>
><br>
><br>
><br>
> Thanks for the feedback. I was able to get more data. The runoff data<br>
> didn?t look useful, so I have attached another variable (precipitation)<br>
> which had more data over the spatial area, but still missing data (-9999).<br>
> The variables are:<br>
><br>
><br>
><br>
> 0) PPT=Precipitation<br>
><br>
> 1) T = Time<br>
><br>
> 2) X= Longitude<br>
><br>
> 3) Y = Latitude<br>
><br>
><br>
><br>
> In order to calculate spatial averages over time would I read in all<br>
> dimensions, similar to what you have described below and calculate the<br>
> average using a function? Thank you for your time.<br>
><br>
><br>
><br>
> x = f->runoff(:,:,:) ; x will be read as a 3-D array<br>
><br>
><br>
><br>
> Regards,<br>
><br>
><br>
><br>
> -Stuart<br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
> *From:* Dave Allured - NOAA Affiliate [mailto:<a href="mailto:dave.allured@noaa.gov">dave.allured@noaa.gov</a>]<br>
> *Sent:* Tuesday, March 13, 2018 5:12 PM<br>
> *To:* Smith, Stuart <<a href="mailto:smit1770@purdue.edu">smit1770@purdue.edu</a>><br>
> *Cc:* Ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
> *Subject:* Re: [ncl-talk] Plotting netCDF<br>
><br>
><br>
><br>
> Stuart,<br>
><br>
><br>
><br>
> That is a strange Netcdf file. It is a time series of 45 x 29 grids, but<br>
> all grid points are missing values, except for one corner point at<br>
> (*,0,0). You can extract the corner point as a 1-D time series, and make a<br>
> line plot and other 1-D analysis.<br>
><br>
><br>
><br>
> x = f->runoff(:,0,0) ; x will be read as a 1-D array<br>
><br>
><br>
><br>
> However, you requested how to calculate spatial averages. That does not<br>
> make sense if you do not have valid data over multiple grid points. So<br>
> either you can just do 1-D plots and analysis, or you need to go back to<br>
> your data source and get a file with valid data over your desired spatial<br>
> area.<br>
><br>
><br>
><br>
> If you get another data file, also please ask the provider to make the<br>
> _FillValue match the array data type, which is "double" in this case.<br>
><br>
><br>
><br>
> --Dave<br>
><br>
><br>
><br>
><br>
><br>
> On Tue, Mar 13, 2018 at 2:20 PM, Smith, Stuart <<a href="mailto:smit1770@purdue.edu">smit1770@purdue.edu</a>><br>
> wrote:<br>
><br>
> Good afternoon,<br>
><br>
><br>
><br>
> I received a netCDF file that has a format I am unfamiliar with and wanted<br>
> guidance on calculating spatial averages over time. Attached is a copy of<br>
> the .nc file. The variables are:<br>
><br>
> 0) Runoff<br>
><br>
> 1) T=Time<br>
><br>
> 2) X= Longitude<br>
><br>
> 3) Y = Latitude<br>
><br>
> When defining variables for runoff (ex. Var = fi->runoff ), the output<br>
> format is the following:<br>
><br>
> (0,0,0) 0<br>
><br>
> (0,0,1) -9999<br>
><br>
> (0,0,2) -9999<br>
><br>
> (0,0,3) -9999<br>
><br>
> (0,0,4) -9999<br>
><br>
> ??..<br>
><br>
> (1,0,0) 0<br>
><br>
> (1,0,1) -9999<br>
><br>
> (1,0,2) -9999<br>
><br>
><br>
><br>
> The next value will not appear until (1,0,0). I wanted to ask how I would<br>
> plot such data, and calculate averages over the area with the given<br>
> format? Thank you for your time.<br>
><br>
><br>
><br>
> Regards,<br>
><br>
><br>
><br>
> -Stuart<br>
><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180315/625a2408/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180315/625a2408/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 15 Mar 2018 15:16:49 -0600<br>
From: Nicholas Thomas Luchetti <<a href="mailto:Nicholas.Luchetti@colorado.edu">Nicholas.Luchetti@colorado.<wbr>edu</a>><br>
To: ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Subject: [ncl-talk] Unable to set explicit levels on contour plot<br>
Message-ID:<br>
<<a href="mailto:CANDmd5LYvsF1Vf6e_wSDNGfNbvG967SkwBKpheYhWLbqeRBivA@mail.gmail.com">CANDmd5LYvsF1Vf6e_<wbr>wSDNGfNbvG967SkwBKpheYhWLbqeRB<wbr>ivA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi all,<br>
<br>
I am trying to set explicit contour levels on a vertical cross-section plot<br>
however with no luck....<br>
<br>
I've tried both "ExplicitLevels" and "ManualLevels" settings but neither<br>
seem to work.<br>
<br>
I've looked at all sorts of examples online but can't seem to see why mine<br>
won't work<br>
<br>
I've attached my script. As well as two plots from 2 different WRF runs...<br>
I want to be able to compare these, but if I can't set the scale the same<br>
then it's difficult to compare....<br>
<br>
Does anyone have any idea as to why I can not set these levels?<br>
<br>
Thanks!<br>
<br>
Nick Luchetti<br>
<br>
--<br>
Nicholas Luchetti<br>
Graduate Student<br>
Department of Atmospheric and Oceanic Sciences (ATOC)<br>
CU Boulder<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180315/b5a5b41b/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180315/b5a5b41b/<wbr>attachment-0001.html</a>><br>
-------------- next part --------------<br>
<br>
;-----------------------------<wbr>------------------------------<wbr>-----------<br>
; WRF_cn_4.ncl<br>
;<br>
; Concepts illustrated:<br>
; - Plotting WRF data<br>
; - Plotting a cross section<br>
;-----------------------------<wbr>------------------------------<wbr>-----------<br>
; WRF: time-z cross section.<br>
;-----------------------------<wbr>------------------------------<wbr>-----------<br>
; These files are loaded by default in NCL V6.2.0 and newer<br>
load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"<br>
load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"<br>
load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/wrf/WRFUserARW.ncl"<br>
<br>
load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/wrf/WRF_<wbr>contributed.ncl"<br>
<br>
begin<br>
;---Open file; substitute your own WRF output file here<br>
f = addfile ("/Users/nicholasluchetti/<wbr>Build_WRF/WRFV3/test/em_<wbr>quarter_ss/<a href="http://wrfout_d01_control.nc" rel="noreferrer" target="_blank">wrfout_d01_control.<wbr>nc</a>","r")<br>
; f = addfile("/Users/<wbr>nicholasluchetti/Build_WRF/<wbr>WRFV3/test/em_quarter_ss/<a href="http://wrfout_d01_simulation_1.nc" rel="noreferrer" target="_blank">wrfou<wbr>t_d01_simulation_1.nc</a>","r")<br>
<br>
; f = addfile("/Users/<wbr>nicholasluchetti/Build_WRF/<wbr>WRFV3/test/em_quarter_ss/<a href="http://wrfout_d01_simulation_4.nc" rel="noreferrer" target="_blank">wrfou<wbr>t_d01_simulation_4.nc</a>","r")<br>
; f = addfile("/Users/<wbr>nicholasluchetti/Build_WRF/<wbr>WRFV3/test/em_quarter_ss/<a href="http://wrfout_d01_simulation_3.nc" rel="noreferrer" target="_blank">wrfou<wbr>t_d01_simulation_3.nc</a>","r")<br>
; f = addfile("/Users/<wbr>nicholasluchetti/Build_WRF/<wbr>WRFV3/test/em_quarter_ss/<a href="http://wrfout_d01_simulation_2.nc" rel="noreferrer" target="_blank">wrfou<wbr>t_d01_simulation_2.nc</a>","r")<br>
<br>
;<br>
; Read character variable Times<br>
; Convert to units of "hours since" for plotting purposes<br>
;<br>
times = f->Times<br>
Time = wrf_times_c(f->Times, 0) ; convert to "hours since"<br>
print(Time)<br>
;<br>
; Read W(bottom_top_stag,west_east) at lat<br>
; index 15, lon index nlon.<br>
; Read associated levels and longitudes<br>
;<br>
nlat = 85<br>
nlon = 120<br>
w = f->W(:,:,nlat,nlon) ; W(Time,bottom_top_stag)<br>
znw = f->ZNW(0,:) ; znw(bottom_top_stag)<br>
<br>
print("test3")<br>
<br>
printVarSummary(w)<br>
w&Time := Time<br>
printVarSummary(w)<br>
print(w&Time)<br>
<br>
print("test4")<br>
<br>
;---Use simple array syntax [like f90] to change units<br>
w = w ; for demo change units<br>
w@units = "m/s"<br>
<br>
printVarSummary(w)<br>
printMinMax(w,0)<br>
<br>
;---Add/change meta data to conform to netCDF convention standards<br>
w!0 = "Time"<br>
w&Time = Time<br>
w!1 = "lev" ; name dimensions<br>
w&lev = znw ; assign values to named dimensions<br>
<br>
<br>
;---For plot purposes only, read the specific lat/lon point<br>
lat = f->XLAT(0,nlat,nlon)<br>
lon = f->XLONG(0,nlat,nlon)<br>
;<br>
; create plots<br>
; (1) A "BlWhRe" is often selected when plus/minus are of interest<br>
; (2) The "symMinMaxPlt" procedure determines contour limits<br>
; that are symmetric.<br>
; (3) Use the "sprintf" function to format the title<br>
; (4) Because the rightmost dimension will become the "x" axis<br>
; use NCL's "dimension reordering" to reshape<br>
;<br>
; wks = gsn_open_wks("pdf","plt_vv_<wbr>vertical_sim2.pdf")<br>
wks = gsn_open_wks("pdf","plt_vv_<wbr>vertical_control.pdf")<br>
; wks = gsn_open_wks("pdf","plt_vv_<wbr>vertical_sim3.pdf")<br>
; wks = gsn_open_wks("pdf","plt_vv_<wbr>vertical_sim1.pdf")<br>
; wks = gsn_open_wks("pdf","plt_vv_<wbr>vertical_sim4.pdf")<br>
gsn_define_colormap(wks,"MPL_<wbr>PRGn")<br>
<br>
<br>
res = True ; plot mods desired<br>
;res@gsnMaximize = True ; maximize plot size<br>
res@cnLevelSelectionMode = "ManualLevels" ; manually set the contour levels with the following 3 resources<br>
res@cnMinLevelValF = -10. ; set the minimum contour level<br>
res@cnMaxLevelValF = 10. ; set the maximum contour level<br>
res@cnLevelSpacingF = 1. ; set the interval between contours<br>
<br>
;res@cnLevelSelectionMode = "ExplicitLevels" ; set explicit contour levels<br>
;res@cnLevels = (/ -12, -10.,-8.,-6.,-4.,-2.,-1.,1.,2.<wbr>,4.,6.,8.,10.,12./) ; set levels<br>
<br>
res@cnFillOn =True<br>
res@cnLinesOn = False ; turn off contour lines<br>
res@lbOrientation = "vertical" ; vertical label bar<br>
<br>
res@trYReverse = True ; reverse y axis<br>
res@tiXAxisString = Time@units ; label bottom axis with units attribute<br>
symMinMaxPlt(w, 14, False, res) ; contributed.ncl<br>
<br>
res@tiMainString = sprintf("%4.2f", nlat)+ " " \<br>
+ sprintf("%4.2f", fabs(nlon))<br>
print("test")<br>
plot = gsn_csm_contour(wks,w(lev|:,<wbr>Time|:),res)<br>
end<br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: plt_vv_vertical_control.pdf<br>
Type: application/pdf<br>
Size: 63930 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180315/b5a5b41b/attachment-0002.pdf" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180315/b5a5b41b/<wbr>attachment-0002.pdf</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: plt_vv_vertical_sim1.pdf<br>
Type: application/pdf<br>
Size: 62514 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180315/b5a5b41b/attachment-0003.pdf" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180315/b5a5b41b/<wbr>attachment-0003.pdf</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Thu, 15 Mar 2018 15:52:14 -0600<br>
From: Adam Phillips <<a href="mailto:asphilli@ucar.edu">asphilli@ucar.edu</a>><br>
To: <a href="mailto:Nicholas.Luchetti@colorado.edu">Nicholas.Luchetti@colorado.edu</a><br>
Cc: ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Unable to set explicit levels on contour plot<br>
Message-ID:<br>
<<a href="mailto:CACvWuqCSoOyJLvkCHCHHkGfa6yPQdVeJfQYH9bMmQjXbF2xohg@mail.gmail.com">CACvWuqCSoOyJLvkCHCHHkGfa6yPQ<wbr>dVeJfQYH9bMmQjXbF2xohg@mail.<wbr>gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Nick,<br>
I think what is happening is that symMinMaxPlt is overruling your previous<br>
settings of cnLevels. Try uncommenting out these lines:<br>
res@cnLevelSelectionMode = "ExplicitLevels" ; set explicit contour levels<br>
res@cnLevels = (/ -12,<br>
-10.,-8.,-6.,-4.,-2.,-1.,1.,2.<wbr>,4.,6.,8.,10.,12./) ; set levels<br>
and then comment out the symMinMaxPlt line:<br>
;symMinMaxPlt(w, 14, False, res)<br>
<br>
If that doesn't fix it let ncl-talk know..<br>
Adam<br>
<br>
On Thu, Mar 15, 2018 at 3:16 PM, Nicholas Thomas Luchetti <<br>
<a href="mailto:Nicholas.Luchetti@colorado.edu">Nicholas.Luchetti@colorado.edu</a><wbr>> wrote:<br>
<br>
> Hi all,<br>
><br>
> I am trying to set explicit contour levels on a vertical cross-section<br>
> plot however with no luck....<br>
><br>
> I've tried both "ExplicitLevels" and "ManualLevels" settings but neither<br>
> seem to work.<br>
><br>
> I've looked at all sorts of examples online but can't seem to see why mine<br>
> won't work<br>
><br>
> I've attached my script. As well as two plots from 2 different WRF<br>
> runs... I want to be able to compare these, but if I can't set the scale<br>
> the same then it's difficult to compare....<br>
><br>
> Does anyone have any idea as to why I can not set these levels?<br>
><br>
> Thanks!<br>
><br>
> Nick Luchetti<br>
><br>
> --<br>
> Nicholas Luchetti<br>
> Graduate Student<br>
> Department of Atmospheric and Oceanic Sciences (ATOC)<br>
> CU Boulder<br>
><br>
><br>
> ______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
><br>
><br>
<br>
<br>
--<br>
Adam Phillips<br>
Associate Scientist, Climate and Global Dynamics Laboratory, NCAR<br>
<a href="http://www.cgd.ucar.edu/staff/asphilli/" rel="noreferrer" target="_blank">www.cgd.ucar.edu/staff/<wbr>asphilli/</a> <a href="tel:303-497-1726" value="+13034971726">303-497-1726</a><br>
<br>
<<a href="http://www.cgd.ucar.edu/staff/asphilli" rel="noreferrer" target="_blank">http://www.cgd.ucar.edu/<wbr>staff/asphilli</a>><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180315/bbc880fa/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180315/bbc880fa/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Thu, 15 Mar 2018 18:07:33 -0600<br>
From: Mary Haley <<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>><br>
To: Colin Raymond <<a href="mailto:csr859@gmail.com">csr859@gmail.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Possible to convert grib2 or netcdf files to<br>
grib1?<br>
Message-ID:<br>
<CACNN_C+t40fAzx-<wbr>2dimV0QJBFpi2CVfLxLcmQ-<wbr>eKLrcEHa=<a href="mailto:4kQ@mail.gmail.com">4kQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Colin,<br>
<br>
NCL cannot write GRIB data, so I'm afraid we can't help you there.<br>
<br>
But we might be able to help with the cnvgrib compilation issues. Usually<br>
these smaller programs aren't too difficult to compile. Are you getting<br>
some kind of error from it?<br>
<br>
--Mary<br>
<br>
<br>
On Tue, Mar 13, 2018 at 4:54 PM, Colin Raymond <<a href="mailto:csr859@gmail.com">csr859@gmail.com</a>> wrote:<br>
<br>
> Hi all,<br>
><br>
> Does anyone know of a procedure or combination of procedures that can<br>
> convert geospatial data *to* grib1 format using NCL 6.4.0? I'm working on a<br>
> Mac (OS 10.12.6, Darwin Kernel 16.7.0) and consequently I have been unable<br>
> to compile cnvgrib, which would be the only alternative that I know of for<br>
> this conversion.<br>
><br>
><br>
> Thanks,<br>
><br>
> Colin Raymond<br>
><br>
><br>
> ______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180315/f16f2f23/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180315/f16f2f23/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Fri, 16 Mar 2018 00:30:24 +0000<br>
From: Guilherme Martins <<a href="mailto:jgmsantos@gmail.com">jgmsantos@gmail.com</a>><br>
To: Mary Haley <<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>>, Colin Raymond <<a href="mailto:csr859@gmail.com">csr859@gmail.com</a>><br>
Subject: Re: [ncl-talk] Possible to convert grib2 or netcdf files to<br>
grib1?<br>
Message-ID:<br>
<CAH+t_hPttX8uR0xYm5=<a href="mailto:fNCrLG6kLHFsT4OeXAPZToSQayYdj9w@mail.gmail.com">fNCrLG6kL<wbr>HFsT4OeXAPZToSQayYdj9w@mail.<wbr>gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
You can use wgrib2.<br>
<br>
<a href="http://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/" rel="noreferrer" target="_blank">http://www.cpc.ncep.noaa.gov/<wbr>products/wesley/wgrib2/</a><br>
<br>
<br>
<br>
Em qui, 15 de mar de 2018 ?s 21:07, Mary Haley <<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>> escreveu:<br>
<br>
> Hi Colin,<br>
><br>
> NCL cannot write GRIB data, so I'm afraid we can't help you there.<br>
><br>
> But we might be able to help with the cnvgrib compilation issues. Usually<br>
> these smaller programs aren't too difficult to compile. Are you getting<br>
> some kind of error from it?<br>
><br>
> --Mary<br>
><br>
><br>
> On Tue, Mar 13, 2018 at 4:54 PM, Colin Raymond <<a href="mailto:csr859@gmail.com">csr859@gmail.com</a>> wrote:<br>
><br>
>> Hi all,<br>
>><br>
>> Does anyone know of a procedure or combination of procedures that can<br>
>> convert geospatial data *to* grib1 format using NCL 6.4.0? I'm working on a<br>
>> Mac (OS 10.12.6, Darwin Kernel 16.7.0) and consequently I have been unable<br>
>> to compile cnvgrib, which would be the only alternative that I know of for<br>
>> this conversion.<br>
>><br>
>><br>
>> Thanks,<br>
>><br>
>> Colin Raymond<br>
>><br>
>><br>
>> ______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
>><br>
>><br>
> ______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
><br>
--<br>
<br>
<<a href="http://www.inpe.br/queimadas/" rel="noreferrer" target="_blank">http://www.inpe.br/queimadas/</a><wbr>><br>
Guilherme Martins<br>
Meteorologista<br>
<br>
Rodovia Presidente Dutra, km <a href="tel:40%20-%2012630000" value="+14012630000">40 - 12630000</a><br>
Cachoeira Paulista, SP - Brasil<br>
<a href="tel:%2B55%2012%203186-9205" value="+551231869205">+55 12 3186-9205</a> <a href="tel:%2B55%2012%2098111-4292" value="+5512981114292">+55 12 98111-4292</a><br>
<a href="mailto:jgmsantos@gmail.com">jgmsantos@gmail.com</a><br>
<a href="http://www.inpe.br/queimadas" rel="noreferrer" target="_blank">http://www.inpe.br/queimadas</a><br>
[image: skype] guilherme.martins.<br>
<a href="https://www.researchgate.net/profile/Guilherme_Martins6" rel="noreferrer" target="_blank">https://www.researchgate.net/<wbr>profile/Guilherme_Martins6</a><br>
[image: Curr?culo Lattes] <a href="http://lattes.cnpq.br/5997657584785803" rel="noreferrer" target="_blank">http://lattes.cnpq.br/<wbr>5997657584785803</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180316/af705455/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180316/af705455/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 7<br>
Date: Fri, 16 Mar 2018 10:58:04 +0100<br>
From: Lara Quiti?n Hern?ndez <<a href="mailto:laraquitianhernandez@gmail.com">laraquitianhernandez@gmail.<wbr>com</a>><br>
To: <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
Subject: [ncl-talk] Regridding Moving Nest WRF outputs<br>
Message-ID:<br>
<<a href="mailto:CADK-_pk3useTuU96-xPmzL-u5bccCPtBnU5BuQpGAUULciNtww@mail.gmail.com">CADK-_pk3useTuU96-xPmzL-<wbr>u5bccCPtBnU5BuQpGAUULciNtww@<wbr>mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Good morning,<br>
<br>
I'm trying to regrid from a WRF grid (curvilinear) to ECMWF grid<br>
(rectilinear). The issue is that my WRF outputs are now obtained using the<br>
Moving Nest configuration. Therefore, I have a new high-resolution domain<br>
(D02) for each time step (I attach an example of what I'm saying).<br>
When D02 was fixed at specific lat/lon coordinates (i.e., without using the<br>
Moving Nest configuration) I used the rcm2rgrid function. But I do not know<br>
if I can use the same function now.<br>
<br>
Thanks in advance!<br>
<br>
Regards,<br>
Lara.<br>
<br>
--<br>
<br>
<br>
*Lara Quiti?n Hern?ndez, PhD Student*<br>
Departamento de F?sica de la Tierra, Astronom?a y Astrof?sica<br>
Facultad de Ciencias F?sicas<br>
Ciudad Universitaria, Plaza Ciencias, 1<br>
28040 Madrid. Spain<br>
Email: <a href="mailto:lquitian@ucm.es">lquitian@ucm.es</a> <<a href="mailto:mlunar@aemet.es">mlunar@aemet.es</a>><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180316/101a4013/attachment.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180316/101a4013/<wbr>attachment.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: D02_MovingNest.gif<br>
Type: image/gif<br>
Size: 353984 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180316/101a4013/attachment.gif" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20180316/101a4013/<wbr>attachment.gif</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br>
<br>
------------------------------<br>
<br>
End of ncl-talk Digest, Vol 172, Issue 30<br>
******************************<wbr>***********<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><font size="2"><i>Colin Raymond</i></font><div><div>Ph.D. Student, Ocean & Climate Physics</div><div>Columbia University</div></div></div></div>
</div></div>