<div dir="ltr"><div>It is the user's responsibility to look at the file being used.</div><div>---<br></div><div>WRF files are (generally) netCDF. <br></div><div>NARR files are GRIB files. <br></div><div><br></div><div>Different information is available for each:<br></div><div>---<br></div><div>%> <b>ncl_filedump</b> merged_AWIP32.1979010100.3D.NARR.<b>grb</b></div><div><b><br></b></div><div>provides the following information for the GRIB file variables<br></div><div><b><br></b></div><div>[SNIP]<br></div><div> float gridlat_221 ( gridx_221, gridy_221 )<br> corners : ( 1, 0.8972797, 46.35283, 46.63458 )<br> long_name : latitude<br> grid_description : AWIPS - Regional - NOAMHI - High Resolution North American Master Grid (Lambert Conformal)<br> units : degrees_north<br><b> mpLambertMeridianF : 253<br> mpLambertParallel2F : 50<br> mpLambertParallel1F : 50<br> mpProjection : LAMBERTCONFORMAL</b><br> Dy : 32463.41<br> Dx : 32463.41<br> Lov : 253<br> Lo1 : 214.5<br> La1 : 1</div><div>=====================================</div><div>This information is *NOT* available on the WRF file.</div><div>Please examine the WRF file's global attributes:</div><div><br></div><div>%> ncdump -h wrfout...</div><div>or</div><div>%> ncl_filedump wrfout ...</div><div><br></div><div>Maybe something like:</div><div><br></div><div>[SNIP]</div><div> :CEN_LAT = 43.99802f ;<br> :CEN_LON = -98.f ;<br> :TRUELAT1 = 30.f ;<br> :TRUELAT2 = 60.f ;[SNIP]<br></div><div>[SNIP]</div><div><br></div><div><div>wrf = addfile("wrf...nc","r")</div><div> ...</div><div>Something like:<br></div><div> <br></div>res@mpLambertParallel1F = fwrf->TRUELAT1 <br><div>res@mpLambertParallel2F = fwrf->TRUELAT2 <br>res@mpLambertMeridianF = fwrf->CEN_LON <br></div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 28, 2020 at 9:13 AM Zhifeng Yang <<a href="mailto:yangzf01@gmail.com">yangzf01@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Hi Dennis,<div>Thank you. I tried to add</div><div><div> res@mpLimitMode = "Corners" ; choose range of map</div><div> res@mpLeftCornerLatF = lat2d(0,0)</div><div> res@mpLeftCornerLonF = lon2d(0,0)</div><div> res@mpRightCornerLatF = lat2d(nlat-1,mlon-1)</div><div> res@mpRightCornerLonF = lon2d(nlat-1,mlon-1)</div></div><div>It can select map range. However, it seems that contour does not move along with map. So there is a shift between contour and map. Then I want to add</div><div><div> res@mpProjection = "LambertConformal" </div><div> res@mpLambertParallel1F = lat2d@mpLambertParallel1F</div><div> res@mpLambertParallel2F = lat2d@mpLambertParallel2F</div><div> res@mpLambertMeridianF = lat2d@mpLambertMeridianF</div></div><div><br></div><div>The problem is there is no latitude attributes (mpLambertParallel1F, mpLambertParallel2F, mpLambertMeridianF) in wrfout files.</div><div><br></div><div>Best,</div><div>Zhifeng</div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jan 27, 2020 at 6:06 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>NARR are Lambert Conformat:</div><div><br></div><div><a href="http://www.ncl.ucar.edu/Applications/narr.shtml" target="_blank"><b>http://www.ncl.ucar.edu/Applications/narr.shtml</b></a></div><div>ncl_3.ncl</div><div><br></div><div>See plot section</div><div>===</div><div>WRF projection information is contained within the file attributes</div><div><br></div><div> fwrf = addfile("wrf...nc","r")</div><div> ...</div><div>Something like:<br></div><div> <br></div>res@mpLambertParallel1F = fwrf->Latin1 <br><div>res@mpLambertParallel2F = fwrf->Latin2 <br>res@mpLambertMeridianF = fwrf->LoV </div><div><br></div><div>========================================</div><div><pre> res@gsnAddCyclic = False ; regional data
res@mpLimitMode = "Corners" ; choose range of map
res@mpLeftCornerLatF = lat2d(0,0)
res@mpLeftCornerLonF = lon2d(0,0)
res@mpRightCornerLatF = lat2d(nlat-1,mlon-1)
res@mpRightCornerLonF = lon2d(nlat-1,mlon-1)
res@tfDoNDCOverlay = True
res@mpProjection = "LambertConformal"
res@mpLambertParallel1F = lat2d@mpLambertParallel1F
res@mpLambertParallel2F = lat2d@mpLambertParallel2F
res@mpLambertMeridianF = lat2d@mpLambertMeridianF
</pre></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jan 27, 2020 at 2:02 PM Zhifeng Yang via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">Hi all,<div>I have WRF output files with Lambert Conformal projection. When I plot images with contour on maps, I tried to limit the plot over a specific region, not the whole domain. I tried the following</div><div><div> res@mpLimitMode = "latlon"</div><div> res@mpMinLatF = minlat</div><div> res@mpMaxLatF = maxlat</div><div> res@mpMinLonF = minlon</div><div> res@mpMaxLonF = maxlon</div></div><div><br></div><div>However, it seems it does not work. Do you have any ideas on it? Thank you!</div><div><br></div><div>Best,</div><div>Zhifeng</div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div>
</blockquote></div>