[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