<html><head></head><body><div class="yahoo-style-wrap" style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div><div><div><span><span>Hello,</span></span></div><div><span><span>This is a repeat of the post I made yesterday 06 April 2019. To make my problem clear, I get the error below. From output of "printVarSummary(p)" given after the error message below, I have never seen coordinates given in the way they are given in this file.</span></span></div><div><span><span><br></span></span></div><div><span><span>check_for_y_lat_coord: Warning: Data either does not contain</span></span></div><span><span>(0) a valid latitude coordinate array or doesn't contain one at all.<br>(0) A valid latitude coordinate array should have a 'units'<br>(0) attribute equal to one of the following values: <br>(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'<br><br></span><span>Variable: p<br>Type: float<br>Total Size: 155976 bytes<br> 38994 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [y | 201] x [x | 194]<br><br></span></span><span></span><div>I would appreciate if someone can assist me modify my script below so that I can plot the data for the first time step. Part of output from ncdump is also given after the script.</div><div><br></div><div>I will appreciate any help.<br></div><div><br></div><div><br></div></div><div><span>begin</span></div><span> <br> f = addfile("~/CORDEX/DATA/pr_AFR-44_MPI-M-MPI-ESM-LR_historical_r1i1p1_MPI-CSC-REMO2009_v1_day_1971-2000.nc","r")<br><br>; p=f->pr(y|:,x|:,time|:) ;With this command the system hangs<br> p=f->pr(0,:,:) ; This command gives the error indicated above<br></span><div><span><br></span></div><div><span> printVarSummary(p)</span></div><div><span></span><br></div><div><span> wks = gsn_open_wks("png", "Simple_Contour_Map")<br><span><br> res = True<br><br> res@gsnAddCyclic = False<br><br> res@cnInfoLabelOn = False<br><br> res@gsnRightString = ""<br> res@gsnLeftString = ""<br><br> res@tiMainString= "Sea Level Pressure" ; draw a title<br><br> res@tiMainString = "Regional Map Plot"<br><br> res@mpDataBaseVersion = "MediumRes" ; use finer database<br> res@mpDataSetName = "Earth..4"<br> res@mpOutlineBoundarySets = "AllBoundaries"<br></span></span><div><span><span> <br></span></span></div><div><span><span><span> res@mpOutlineOn = True</span></span></span></div><div><span><span><span></span></span></span></div><span><span><span> res@mpNationalLineThicknessF = 1.5 ; turn on country boundaries<br><br> plot= gsn_csm_contour_map(wks, p, res) ; create the plot<br><br></span></span></span><div><span><span><span> end</span></span></span></div><div><span><span><span></span></span></span><br></div><div>Part Output of ncdump -c filename</div><div><br></div><span><span><span><span>dimensions:<br> time = UNLIMITED ; // (10958 currently)<br> bnds = 2 ;<br> x = 194 ;<br> y = 201 ;<br> vertices = 4 ;<br>variables:<br> double time(time) ;<br> time:standard_name = "time" ;<br> time:long_name = "time" ;<br> time:bounds = "time_bnds" ;<br> time:units = "days since 1949-12-01 00:00:00" ;<br> time:calendar = "proleptic_gregorian" ;<br> time:axis = "T" ;<br> double time_bnds(time, bnds) ;<br> float lon(y, x) ;<br><span> lon:standard_name = "longitude" ;<br> lon:long_name = "longitude coordinate" ;<br> lon:units = "degrees_east" ;<br> lon:_CoordinateAxisType = "Lon" ;<br> lon:bounds = "lon_bnds" ;<br> float lon_bnds(y, x, vertices) ;<br> float lat(y, x) ;<br> lat:standard_name = "latitude" ;<br> lat:long_name = "latitude coordinate" ;<br> lat:units = "degrees_north" ;<br> lat:_CoordinateAxisType = "Lat" ;<br> lat:bounds = "lat_bnds" ;<br> float lat_bnds(y, x, vertices) ;<br> float pr(time, y, x) ;<br><br></span><br></span><br><br></span><br></span><br></span><br></div><span><br></span><br></div></div></body></html>