Re: Error with drawing a polygon in a polar sterographic map

Hi Ipshita,

Could you send the complete script. The polygon part is missing hence we =
can’t say anything about it.

The warning about mpMinLonF and MaxLonF tells you what it is, you can’t =
set it. Just mpMinLatF and
mpMaxLatF can be used. See example polar_5.ncl on =
http://www.ncl.ucar.edu/Applications/polar.shtml

Bye,
Karin


Am 23.07.2014 um 21:44 schrieb Ipshita Majhi <ipmajhi@alaska.edu>:

> Hi,
>
> I am trying put a box around the region of interest. But I get two =
errors and I am not sure how to get a handle around it.
>
> One is
> Warning: You should not set the mpMinLonf and /or MarLonF resources.
>
> fatal: NhlCvStringToEnum:Unable to convert string "missing" to =
requested type
>
> Warning : Error retrieving resource om Just from args-Ignoring Args.
>
> Here is the code:-
>
> I want to put a box around a particular region and the extract the =
data and take average of that region. Do it for all the time series
>
> ;Date:07/23/2014
> ;Ipshita Majhi
>
> ;************************************************
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> ;************************************************
> begin
> ;************************************************
> ; read in netCDF file
> ;************************************************
> a =
=addfile("~/Documents/Snow_cover_IMS/nhsce_v01r01_19661004_20140602.nc",=
"r")
>
> sce=a->snow_cover_extent(0,:,:)
> lat=a->latitude
> lon=a->longitude
>
>
> ;Attaching latitude and longitude to sce
> sce@lat2d=lat
> sce@lon2d=lon
>
>
> ;************************************************
> ; create plot
> ;************************************************
> wks = gsn_open_wks("ps","polyg") ; open a ncgm file
>
> res = True ; plot mods =
desired
> res@gsnDraw = False ; don't draw yet
> res@gsnFrame = False ; don't advance =
yet
>
> res@mpMinLatF = 40 ; zoom in on a =
subregion
> res@mpMaxLatF = 60
> res@mpMinLonF = 10
> res@mpMaxLonF = 30
>
> res@gsnPolar = "NH" ; specify the =
hemisphere
> plot = gsn_csm_contour_map_polar(wks,sce,res) ; create the plot
>
> ;************************************************
> ; create points for box
> ;************************************************
> ypts = (/ 30.0, 30.0, 0.0, 0.0, 30.0/)
> xpts = (/-90.0, -45.0,-45.0, -90.0,-90.0/)
> ;************************************************
> ; add the box
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk


Received on Wed Jul 23 08:46:57 2014