<div dir="ltr"><div dir="ltr"><div><b><br></b></div><div><b>[1]</b></div><div>I have tried opening the <br></div><div> <b>3A-DAY.</b>GPM.DPR.GPM-SLH.20161231-S000000-E235959.366.<b>V06A.HDF5</b></div><div>file with NCL versions 6.3.0, 6.5.0 and 6.6.2 and have encountered the same errors as Rick.</div><div><br></div><div>Please place the file you are using "<b>3A-DAY.20161231.V06A.</b><b>HDF5</b>" onto to ftp site</div><div><br></div><div>ftp <a href="http://ftp.cgd.ucar.edu">ftp.cgd.ucar.edu</a></div><div>anonymous</div><div>your_email</div><div>cd incoming</div><div>put 3A-DAY.20161231.V06A.HDF5</div><div>quit<br></div><div><br></div><div>===</div><div>I tried to open a monthly '3A' file with no problems:</div><div><br></div><div>%> <a href="https://www.ncl.ucar.edu/Document/Tools/ncl_filedump.shtml" target="_blank"><b>ncl_filedump</b></a> <b>3A-MO</b>.GCOMW1.AMSR2.GRID2014R2.20150101-S000000-E235959.01.<b>V03C.HDF5</b></div><div>===<b><br></b></div><div><b>[2]<br></b></div><div>A netCDF<a href="https://www.ncl.ucar.edu/Document/Language/cv.shtml"><b> coordinate variable</b></a> is a one-dimensional array where the name of the variable is the same as the dimension name: EG: time(time), level(level), lat(lat), longitude(longitude)</div><div>---</div><div> layer=fspan(0.25,20,80)<br> lon=fspan(-180,180,720)<br> lat=fspan(-67,67,268)<br> lat@long_name = "latitude"<br> lon@long_name = "longitude"<br> lat@units = "degrees_north"<br> lon@units = "degrees_east"</div><div> lon!0 = "lon"</div><div> lat!0 = "lat"</div><div> layer!0 = "layer"</div><div><br></div><div><b></b></div></div>QR1 = "/Grid/LHCndMean"<br> QR=f->$QR1$<br> QR!0="layer" ; name dimensions<br> QR!1="lon"<br> QR!2="lat"<br> QR&layer=layer ; assign variables to named dimensions<br> QR&lon=lon<br><div> QR&lat=lat</div><div><br></div><div> printVarSummary(QR)<br></div> printVarSummary(QR&lat)<br><div> printVarSummary(QR&non)</div><div><br></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Mar 8, 2020 at 8:16 AM Rick Brownrigg 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>Unfortunately, I am unable to open this file. My version of NCL complains about mismatched headers/libraries, which would be an issue on my part. I'll have to wait until tomorrow when I have access to a different version of NCL. That said, the various HDF5 command-line utilities I have installed also say they can not open the file (???)</div><div><br></div><div>In the meantime, perhaps someone else can take a crack at it....</div><div><br></div><div>Rick</div><div><br></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Mar 7, 2020 at 10:13 PM Ipsita Putatunda <<a href="mailto:ipsita.putatunda@gmail.com" target="_blank">ipsita.putatunda@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>Hi,</div><div>I think the length of the array "lat" is corresponding the dimension of the variable. But am not sure where exactly I am making mistake to plot the variable.</div><div>For your reference I have sent my file (3A-DAY.GPM.DPR.GPM-SLH.20161231-S000000-E235959.366.V06A.HDF5) through ftp in <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a>.</div><div><br></div><div>Thanks in advance,</div><div>Ipsita<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Mar 8, 2020 at 10:01 AM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@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>Hi,</div><div><br></div><div>Offhand, it looks like you are constructing the coordinate arrays properly. I notice you calculate the lat coordinate array with this statement:</div><div><br></div><div> lat=fspan(-67,67,268)</div><div><br></div><div>Are you certain the length of this array lines up with the corresponding dimension of the variable QR?</div><div><br></div><div>HTH...Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Mar 7, 2020 at 10:58 AM Ipsita Putatunda 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>Dear NCL users,</div><div>Getting warning message while trying to plot the variable.</div><div><br></div><div>"Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all."</div><div>My script is as follows:</div><div>load "/usr/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "/usr/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br> begin<br> netCDF=True<br> plot=True<br> diri="./"<br>;==============================================================<br> fn1="3A-DAY.20161231.V06A.HDF5"<br> printVarSummary(fn1)<br> lev=fspan(0.25,20,80)<br> lon=fspan(-180,180,720)<br> lat=fspan(-67,67,268)<br> lat@long_name = "latitude"<br> lon@long_name = "longitude"<br> lat@units = "degrees_north"<br> lon@units = "degrees_east"<br>;-------------------------------------------------------------------<br> f= addfile(fn1, "r")<br> ;g=f->Grid<br> QR1 = "/Grid/LHCndMean"<br> QR=f->$QR1$<br> QR!0="nlayer"<br> QR!1="nlon"<br> QR!2="nlat"<br> QR&nlayer=lev<br> QR&nlon=lon<br> QR&nlat=lat<br> QR&nlat@units="degrees_north"<br> QR&nlon@units="degrees_east"<br> printVarSummary(QR&nlat)<br> printVarSummary(QR&nlon)<br>;-------------------------------------------------------------------<br> wks = gsn_open_wks ("png","QR")<br> gsn_define_colormap(wks,"rainbow")<br><br>res = True<br>res@gsnMaximize = True<br>res@cnFillOn = True<br>res@cnLinesOn = False<br>res@cnLineLabelsOn = False<br>res@gsnSpreadColors = True<br>res@gsnAddCyclic = False<br> res@mpMinLatF = min(QR&nlat)<br> res@mpMaxLatF = max(QR&nlat)<br> res@mpMinLonF = min(QR&nlon)<br> res@mpMaxLonF = max(QR&nlon)<br>plot = gsn_csm_contour_map(wks,QR(0,:,:),res)<br> end</div><div><br></div><div><br></div><div>Any help in this issue will be appreciated.</div><div>Thanks in advance,</div><div>Ipsita<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></blockquote></div>
</blockquote></div>
</blockquote></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></div>