[ncl-talk] Not getting the original lat-lon label in contour map

Atul Saini atulsainimail at gmail.com
Thu Dec 5 07:38:43 MST 2019


*lat=a->latprintVarSummary(lat)*

The above given code in ncl itself results into following detail, which
contain the units and other things-













*Variable: latType: floatTotal Size: 200 bytes            50 valuesNumber
of Dimensions: 1Dimensions and sizes: [lat | 50]Coordinates:
lat: [ 8..20.25]Number Of Attributes: 4  standard_name : latitude
long_name : latitude  units : degrees_north  axis : Y*

Regards,





On Thu, Dec 5, 2019 at 7:59 PM Atul Saini <atulsainimail at gmail.com> wrote:

> Hi dear Karin,
>               I followed the same as suggested by you but it could not
> result into expected output. On the other hand I did 'ncdump -h' on the
> data file and it shows that data file contain the detail of unit of
> coordinate variable. May be I am still doing mistake somewhere, please
> suggest what should be done now, when setting lat lon unit do not work?
>
>
> a=addfile("WesternAgroClimRegion_MonthlySum_1901_2017Jan.nc","r")
> b=a->rf
> lat at units = "degrees_north"
> lon at units = "degrees_east"
> aveg=dim_avg_n(b,0)
> stdev=dim_stddev_n(b,0)
> ;stdev = where(stdev.ne.0,stdev,stdev at _FillValue)
> cv=(stdev/aveg)*100
> ;*****************************************
> ;create plot
> ;*****************************************
>   wks  = gsn_open_wks("png","CV_map")   ; send graphics to PNG file
>   res                     = True          ; plot mods desired
>   res at cnFillOn            = True          ; turn on color
>   res at cnLinesOn           = False
>   res at cnFillMode          = "RasterFill"       ; Raster Mode
>   res at cnFillPalette       = "gui_default" ; set color map
>   ;res at vpWidthF           = 0.5          ; change aspect ratio of plot
>   ;res at vpHeightF          = 0.3
>   res at gsnMaximize         = False          ; maximize plot size
>   ;res at gsnDraw            = False
>   ;res at gsnFrame           = False
>   plot = gsn_csm_contour(wks,cv,res)       ; contour the variable
>   frame(wks)
>   end
>
> Regards,
>
>
>
>
>
> On Thu, Dec 5, 2019 at 7:16 PM Karin Meier-Fleischer via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
>> The lat and lon variables don't have a units attribute.
>>
>> lat at units = "degrees_north"
>> lon at units = "degrees_east"
>>
>> -Karin
>>
>> Am 05.12.19 um 13:50 schrieb Atul Saini via ncl-talk:
>>
>> Hi NCL,
>>               In spite of getting the real lat-lon label, I am getting
>> the serial number of lat-lon. Please suggest the solution!
>>
>> Please find the script and data attached with this mail.
>> Thanks in advance.
>>
>> Regards,
>>
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing listncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191205/b0c7678f/attachment.html>


More information about the ncl-talk mailing list