[ncl-talk] Error with lat and lon corrdinates

[Rick Brownrigg]

Hi,

What do you mean by ".... lat lon dimensions don't show up"?  What do you
get instead?

I think the reason for the warnings is that you are not plotting on top of
a map, and thus those resources don't apply.  Perhaps you want to plot
using gsn_csm_contour_map()  ?

HTH...
Rick


On Sun, Nov 20, 2016 at 5:17 PM, Ipshita Majhi <ipmajhi at alaska.edu> wrote:

> Dear NCL,
>
> I am trying to plot SWE data and the lat lon dimensions don't show up.
> I would be grateful if you could guide me on this
> *********************************************************
> ;********************************************
> ;This is reading in Snow water equivalent which goes from 0 to 250 mm
> ; from 1982 to 2012
> ;********************************************
>
> dira = "/Users/ipshita/Desktop/GlobsnowV2.0/version02/"
> dirb = dira + "2011/"
> ;dira = "./"
> ;dirb = "./"
> a = addfile(dira+"monthly_82_13.nc","r")
> b = addfile(dirb+"GlobSnow_SWE_L3A_20110101_v2.0.nc","r")
>
> swe  = a->swe
> time = a->time
> lat  = b->lat
> lon  = b->lon
> swe at lat2d = lat
> swe at lon2d = lon
>
> print(dimsizes(lat))
> print(dimsizes(lon))
>
> swe_FillValue=-9.99e+08
>
> utc_date = cd_calendar(time, 0)
> year=utc_date(:,0)
> month=utc_date(:,1)
>
> ;********************************************
> ;This is to extract 1982-2012 of swe data
> ;********************************************
> swe_82_12=swe(0:371,:,:)
>
> ;********************************************
>
> ;********************************************
> ;This is to read in All India Rainfall
> ;For March , March , June
> ;********************************************
>
> MJ_1982_2014=asciiread("~/Documents/NCL_files/SST/
> Monthly/MJ.txt",(/34,1/),"float")
> MJ_1982_2012=MJ_1982_2014(0:30,0)
>
> print(dimsizes(MJ_1982_2012))
>
> ;*******************************************
> ;Extracting monthly data for sce
> ;*******************************************
>
> March_swe = new((/31,721,721/),float)
>
> March_swe at _FillValue = -9.99e+08
>
> do nyr=0,371,12
>
> March_swe(nyr/12,:,:) =swe_82_12(nyr+2,:,:)
>
> end do
> ;=================================================
> March_swe at _FillValue = -9.99e+08
> ;**********************************************
>
> March_dt = dtrend_msg(March_swe&time,March_swe(lat|:,lon|:,time|:),
> True,False)
>
> printMinMax(March_dt,False)
>
> ;****************************************
> ;Calculating Correlation
> ;***************************************
> x1d=ndtooned(MJ_1982_2012)
>
> corr_March=escorc(March_dt,x1d)
> corr_March!0="lat"
> corr_March!1="lon"
> corr_March&lat=lat
> corr_March&lon=lon
>
> copy_VarCoords(swe,corr_March)
>
>
> ;**********************************************
>
> ;************************************************
> ; create plot
> ;************************************************
>   wks = gsn_open_wks("x11","SWE_March_MJ_correlation_significance")
>         ;
>
>     res                     = True              ; Plot modes desired.
>     res at gsnMaximize         = True              ; Maximize plot
>     res at cnFillOn            = True              ; color plot desired
>     res at cnFillPalette       = "BrownBlue12"
>     res at cnLinesOn           = False             ; turn off contour lines
>     res at pmLabelBarWidthF    = 0.9               ; make wider
>     res at pmLabelBarHeightF   = 0.1               ; default is taller
>     res at lbLabelFontHeightF  = .018              ; default is HUGE
>
>
> ;*******************************************
> ; georeferencing: plot on polar projection
> ;*******************************************
>     res at trGridType   = "TriangularMesh"        ; allow missing coordinates
>     res at gsnAddCyclic = False
>     res at cnFillMode   = "RasterFill"
>
>     res at gsnPolar   = "NH"                          ; specify the
> hemisphere
>     res at mpMinLatF  = 35
>
>     res at tiMainString = "March SWE and MJ correlation"
>
>   plot = gsn_csm_contour(wks,corr_March,res)
>
> ;*******************************************************
>
> warning:gsnPolar is not a valid resource in SWE_March_MJ_correlation_significance_contour
> at this time
> warning:mpMinLatF is not a valid resource in SWE_March_MJ_correlation_significance_contour
> at this time
>
>
>
>
> _______________________________________________
> 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/20161121/e56b4b5d/attachment.html