<div dir="ltr"><div>Thanks a lot Dave. It's helpful. I tried running without any coordinate script, it is running and I was able to generate the EOF and EOF-time series netcdf files for the entire domain. I am interested in coordinate subscripting as in the EOF example 1. <br><br>When I tried to enter the lon lat ranges, the script does not run. It says lon is not a dimension or variable. It accepts x or y, but x and y are some strange values in the data. If I type ncl_filedump file -v lon, I see the values. <br><br></div><div>Here is the ncl_filedump summary.<br></div><div><br>short air ( time, y, x )<br>..<br>float lat ( y, x )<br>....<br></div><div>float lon ( y, x )<br><br>1) Is there a way to do the subscripting in the EOF code by defining lon and lat instead of x and y.<br><br></div><div>2) Also, the output netcdf file of EOF does not have the lon, lat meta data. Instead it has x and y, which cannot does not overlay when I tried with NARR plot.<br><br></div><div>Here is my script.<br><br>;Read the NARR data (Air temperature at 2m)<br>f = addfile ("<a href="http://air.2m.mon.mean.nc">air.2m.mon.mean.nc</a>", "r")<br> <br>;Define the parameter, convert the format short to float<br> p = short2flt(f->air(::12,:,:))<br><br>;Calculate the weights<br>w = sqrt(cos(0.01745329*p&y)) ; p&y works, but if i enter p&lat, it does not work.<br>wp = p*conform(p, w, 1)<br>copy_VarCoords(p, wp)<br><br>t = wp(y|:,x|:,time|:) ; x and y works, but if i enter lon or lat, it does not work.<br><br>;Number of EOF to be calculated<br>neof = 3<br><br>;NCL EOF functions <br>eof = eofunc_Wrap(t, neof, False) ;EOF (3modesXlonXlat)<br>eof_ts = eofunc_ts_Wrap (t, eof, False) ;EOF (3modesXtime)<br><br>;Export EOF spatial data as netcdfile<br>system ("/bin/rm -f EOFnarr.nc") <br>fout = addfile ("EOFnarr.nc", "c")<br>fout@title = "EOFs of pre 1971-2013"<br>fout->EOF = eof<br><br></div><div>Regards<br></div><div>Jothiganesh<br></div><div><br></div><div><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 16, 2015 at 9:58 PM, Dave Allured - NOAA Affiliate <span dir="ltr"><<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Jothiganesh,<div><br></div><div>Yes you can directly read NARR data and compute EOF's without regridding. The underlying EOF functions eofunc and eofunc_ts do not use any spatial coordinates. They operate purely on numeric arrays.<br></div><div><br></div><div>The first EOF example uses coordinate subscripting and the lonFlip function. If you remove both of these methods from the example, then you can compute EOF's over the entire NARR grid.</div><div><br></div><div>However, to generate valid plots, you will need to retain the NARR coordinate metadata, and modify the plot section to use one of the plotting methods shown on the NARR examples page.</div><div><br></div><div>I recommend one of the projection methods in NARR example 4, because they use the retained 2-D coordinate grids to directly align the data with the underlying map. After you get that working, you can switch back to Lambert or some other projection of your choice.</div><div class="gmail_extra"><br></div><div class="gmail_extra">--Dave<br><br><br><div class="gmail_quote"><div><div class="h5">On Tue, Jun 16, 2015 at 6:39 PM, Jothi Ganesh Shanmuga Sundaram <span dir="ltr"><<a href="mailto:joshanmugasundaram@mix.wvu.edu" target="_blank">joshanmugasundaram@mix.wvu.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div><div><div><div><div><div>Hello,<br><br></div>I use NARR gridded data (349x277 Lambert Conformal Conic grid) in netcdf format in NCL. I am able to plot graphics without any problem using the existing script in NCL website.<br><br></div>I want to do EOF analysis with NARR data. When I tried to use the EOF script available from NCL website to analyze NARR data, it is not working. When I checked the data using ncl_file dump the lon, lat formats are very different. My question is<br><br></div>1) Whether I can do EOF analysis in NCL by directly reading NARR data? If so, I cannot define the LON LAT value, as NARR has a projected coordinates value?<br><br>2) Do I need to regrid the NARR curvilinear to rectlinear data before doing EOF? or Is there a way to convert the format of LON, LAT in grid without regridding inorder to do EOF analysis in NCL<br><br></div>Thanks for your help.<br><br></div>Regards<br></div>Jothiganesh<br></div>
<br></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><br>
<br></blockquote></div><br></div></div>
</blockquote></div><br></div>