<div dir="ltr"><div><div><div><div><div>[1]<br> PLEASE ... only send a minimal script. ncl-talk does not want to hack thru many commented out lines. *Our time is valuable too!"<br><br>[2]<br></div>Look at the dimensionality of the array<br><br> surf_pptn = f->surfacePrecipitation ; (:,:) => (nlon,nlat) ==> (720,160)<br> printVarSummary(surf_pptn)<br><br></div>Incidently, it is more efficient to *not* use the (:,:) when reading an array.<br><br>[3]<br></div>NCL's high level graphics routines *require* (lat,lon) ordering. Hence, you must reorder. From your code (I do not have nor do I want sent a 3A12 file)<br><br>[4]<br> Y_axis = fspan(-40,40,nlat)<br> X_axis = fspan(-180.0,180.0,nlon)<br> Y_axis@units = "degrees-north"<br> X_axis@units = "degrees-east"<br> Y_axis!0 = "latitude"<br> X_axis!0 = "longitude"<br><br> surf_pptn = f->surfacePrecipitation ; (:,:) => (nlon,nlat)<br> printVarSummary(surf_pptn)<br> printMinMax(surf_pptn,True)<br><br> surf_pptn@_FillValue = -9999.9<br> surf_pptn!0 = "longitude" ;Name the DIMS <br> surf_pptn!1 = "latitude" ;"latitude"<br> surf_pptn&latitude = Y_axis<br> surf_pptn&longitude = X_axis<br> printVarSummary(surf_pptn)<br> printMinMax(surf_pptn,True)<br> ;;;asciiwrite("surf.dat",surf_pptn)<br><br>; NCL high-level graphics *require* (lat,lon) order <br><br> work = surf_pptn<br> surf_pptn := work(latitude|:,longitude|:)<br> printVarSummary(surf_pptn)<br><br>[5]<br></div>I have no idea what you are trying to do with your graphics.<br><br></div>Good luck<br><div><div><br><div><div><br></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 18, 2016 at 7:07 AM, Geeta Geeta <span dir="ltr"><<a href="mailto:geetag54@yahoo.com" target="_blank">geetag54@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,Sans-Serif;font-size:16px"><div dir="ltr"><span>Hi all. </span></div><div dir="ltr"><span><br></span></div><div dir="ltr"><span>I have to plot the TRMM 3A12 data which is gridded. I am getting following error even after defining the units of the </span></div><div dir="ltr"><span>cordinate variables.</span></div><div dir="ltr"><span><br></span></div><div dir="ltr">Variable: surf_pptn</div><div dir="ltr">Type: float</div><div dir="ltr">Total Size: 460800 bytes</div><div dir="ltr"> 115200 values</div><div dir="ltr">Number of Dimensions: 2</div><div dir="ltr">Dimensions and sizes:<span style="white-space:pre-wrap"> </span>[longitude | 720] x [latitude | 160]</div><div dir="ltr">Coordinates: </div><div dir="ltr"> longitude: [-180..180]</div><div dir="ltr"> latitude: [-40..40]</div><div dir="ltr">Number Of Attributes: 5</div><div dir="ltr"> longitude :<span style="white-space:pre-wrap"> </span>degrees_east</div><div dir="ltr"> latitude :<span style="white-space:pre-wrap"> </span>degrees_north</div><div dir="ltr"> _FillValue :<span style="white-space:pre-wrap"> </span>-9999.9</div><div dir="ltr"> units :<span style="white-space:pre-wrap"> </span>mm/hr</div><div dir="ltr"> hdf_name :<span style="white-space:pre-wrap"> </span>surfacePrecipitation</div><div dir="ltr">(0)<span style="white-space:pre-wrap"> </span> </div><div dir="ltr">(0)<span style="white-space:pre-wrap"> </span>min=0 max=2.99362</div><div dir="ltr">(0)<span style="white-space:pre-wrap"> </span>check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all.</div><div dir="ltr">(0)<span style="white-space:pre-wrap"> </span>A valid latitude coordinate array should have a 'units' attribute equal to one of the following values: </div><div dir="ltr">(0)<span style="white-space:pre-wrap"> </span> 'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'</div><div dir="ltr">(0)<span style="white-space:pre-wrap"> </span>check_for_lon_coord: Warning: Data either does not contain a valid longitude coordinate array or doesn't contain one at all.</div><div dir="ltr">(0)<span style="white-space:pre-wrap"> </span>A valid longitude coordinate array should have a 'units' attribute equal to one of the following values:(0)<span style="white-space:pre-wrap"> </span> 'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'</div><div dir="ltr"><br></div><div dir="ltr">I am attaching the figure which is distorted. </div><div dir="ltr"><span><br></span></div><div dir="ltr">Secondly the documentation says that the range of the data is from SurfaceRain Float SDS 4 nlat*nlon -‐ 0 – 3000 mm/hr </div><div dir="ltr"><a href="http://pps.gsfc.nasa.gov/Documents/filespec.TRMM.V7.pdf" target="_blank">http://pps.gsfc.nasa.gov/Documents/filespec.TRMM.V7.pdf</a>) Page 203. <br></div><div dir="ltr">when I plot the MinMAx values of this variable, </div><div dir="ltr">hdf_name :<span style="white-space:pre-wrap"> </span>surfacePrecipitation</div><div dir="ltr">(0)<span style="white-space:pre-wrap"> </span> </div><div dir="ltr">(0)<span style="white-space:pre-wrap"> </span>min=0 max=2.99362</div><span class="HOEnZb"><font color="#888888"><div dir="ltr"><br></div><div>Geeta.</div></font></span></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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>