[ncl-talk] Plot warnings
Anahita Amiri Farahani
aamir003 at ucr.edu
Wed Jan 18 17:19:33 MST 2017
Hello all,
I tried to plot a simple map but I got warnings. Here is my code and
warnings:
;*********** From Observation **********
ana=addfile("NOaerosol_allsky_3DQ.nc","r")
lat2=ana->lat
lon2=ana->lon
lev2=ana->lev/10
Absp3D=ana->Absp3D ; Atmospheric heating (K/sec)
Absp=ana->Absp ; Solar radiation absorption in the atmosphere (W/m^2)
Sfc=ana->Sfc ; Surface net flux (W/m^2)
delete(ana)
;********* T's data *************
;z1 = asciiread("Absp.dat",(/12,128,64,26/),"float")
z1 = asciiread("Absp.dat",(/128,64,26,12/),"float")
z1!1="lat"
z1!0="lon"
;Coordinate
z1&lat=lat2
z1&lon=lon2
wks = gsn_open_wks("eps","T_Atmospheric_heating_26")
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
res = True ; plot mods desired
res at cnFillOn = True ; turn on color
res at gsnSpreadColors = True ; use full color map
res at cnLinesOn = False ; no contour lines
res at cnLineLabelsOn = False ; no line labels
res at mpMaxLatF = 87.8638 ; specify the plot domain
res at mpMinLatF = -87.8638 ;
res at mpMinLonF = 0 ;
res at mpMaxLonF = 357.1875
res at gsnAddCyclic = False
res at cnLevelSelectionMode = "AutomaticLevels" ; AutomaticLevels
;res at cnLevelSelectionMode = "ManualLevels" ; AutomaticLevels
;res at cnMinLevelValF = -15. ; min level
;res at cnMaxLevelValF = 15. ; max level
;res at cnLevelSpacingF = 1. ; interval
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
res at tiMainString = "Atmospheric Heating in Jan at 26th level"
plot= gsn_csm_contour_map(wks,z1(:,:,25,0),res)
delete(wks)
0) check_for_y_lat_coord: Warning: Data either does not contain a valid
latitude coordinate array or doesn't contain one at all.
(0) A valid latitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0) check_for_lon_coord: Warning: Data either does not contain a valid
longitude coordinate array or doesn't contain one at all.
(0) A valid longitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
(0) gsn_csm_map_ce: Warning: you set mpMaxLonF to a value > 180, but
(0) didn't set mpCenterLonF. Setting mpCenterLonF to
178.594
Thanks,
Ana
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170118/1f2aeec1/attachment.html
More information about the ncl-talk
mailing list