[ncarg-talk] contour over map (from CISM output)

[Mary Haley]

Hi Michele,

The quick answer is that in NCL V6.4.0 there's some code that still expects
tfDoNDCOverlay to just be set to True or False, so instead of:

res at tfDoNDCOverlay = "DataTransform"

use:

res at tfDoNDCOverlay = False

This resource is False by default, however, so you can also just not set it
at all.

In the future, please post NCL questions to ncl-talk and not ncarg-talk.
The ncarg-talk talk list is for people working with the old NCAR Graphics
libraries, and there's not much activity on this email list.  You will need
to subscribe:

http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Regards,

--Mary


On Tue, May 8, 2018 at 3:21 PM, Michele Petrini - CITG <M.Petrini at tudelft.nl
> wrote:

> Dear helpdesk,
>
> I am trying to plot Greenland ice thickness (from CISM model) contour over
> a map, but I get an error I never had before - so far, I plotted ice
> thickness contours using outputs from different ice sheet models with no
> problems. Below you can find the code
>
> begin
>
> in = "./data/"
> out = "./results/"
> H_filename = "ice_thickness.nc"
> H_file = addfile(in+H_filename,"r")
>
> H            = H_file->thk
> time_sim         = ispan(0,16100,100)
> ny =  dimsizes(H&y1)
> nx =  dimsizes(H&x1)
> ntime = dimsizes(time_sim)
>
> H at _FillValue = -999.
>
> H = mask(H,H.lt..1,False)
>
> lat1 = fspan(58.25,84.61,ny)
> lat1 at units = "degrees_north"
> lat1!0 = "lat"
> lat1 at long_name = "latitude"
> lat1 at standard_name = "latitude"
> lat1 at axis = "Y"
> print(lat1)
> lon1 = fspan(-95.68,13.33,nx)
> lon1 at units = "degrees_east"
> lon1!0 = "lon"
> lon1 at long_name = "longitude"
> lon1 at standard_name = "longitude"
> lon1 at axis = "X"
>
> H!1 = "lat"
> H!2 = "lon"
> H&lat = lat1
> H&lon = lon1
>
> minlat = 58.25
> maxlat = 84.61
> minlon = -95.68
> maxlon = 13.33
>
>   res                       = True              ; plot mods desired
>   ;res at gsnDraw               = False
>   ;res at gsnFrame              = False
>   res at cnFillOn              = True              ; turn on color
>   res at cnFillMode            = "RasterFill"
>   res at cnLinesOn             = False
>   res at cnLineLabelsOn        = False
>   res at lbLabelBarOn          = False             ; turn off individual lb's
>
>   res at trGridType            = "TriangularMesh"
>
>   res at mpProjection          = "Stereographic"
>   res at mpDataBaseVersion     = "mediumres"
>   res at mpFillOn              = False             ; turn off default land
> map fill
>
>   res at mpLimitMode           = "Corners"
>   res at mpLeftCornerLatF      = minlat
>   res at mpLeftCornerLonF      = minlon
>   res at mpRightCornerLatF     = maxlat
>   res at mpRightCornerLonF     = maxlon
>   res at mpCenterLonF          = 0
>   res at mpCenterLatF          = 90
>   res at gsnLeftString         = ""
>   res at gsnRightString        = ""
>   res at cnFillPalette        = "testcmap"
>   res at cnLevelSelectionMode = "ManualLevels"     ; set manual contour
> levels
>   res at cnMinLevelValF       =     0.
>   res at cnMaxLevelValF       =     1.
>   res at cnLevelSpacingF      =     0.2
>   res at tfDoNDCOverlay      = "DataTransform"
>
> do i=0,ntime-1
> wks = gsn_open_wks("pdf","./results/image-"+sprinti("%03i",i))
>         print("simulation time(" + i + ") = " + time_sim(i))
>         res at gsnRightString = "Simulation time: " + time_sim(i) + " (yr)"
> ice = new((/ny,nx/),double)
> ice(:,:) = H(i,:,:)
> plot = new(1,graphic)
> plot(0) = gsn_csm_contour_map(wks,ice,res)
> delete(wks)
> end do
>
> cmd = "cd ./results/ ; pdflatex movie.tex"
> system(cmd)
>
> end
>
> whereas the error I get is the following
>
> fatal:Or: operation not supported on type (string)
> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 3642 in
> file /usr/local/ncl-6.4.0/lib/ncarg/nclscripts/csm/gsn_csm.ncl
>
> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 7985 in
> file /usr/local/ncl-6.4.0/lib/ncarg/nclscripts/csm/gsn_csm.ncl
>
> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 11287 in
> file /usr/local/ncl-6.4.0/lib/ncarg/nclscripts/csm/gsn_csm.ncl
>
> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 80 in
> file contour_movie.ncl
>
> I can’t find any suggestion by looking at previous questions and I really
> don’t understand what the problem can be.
> thanks!
>
> best,
> Michele
>
>
>
> ******************************************************************
> Michele Petrini
>
> PhD in Earth Science and Fluid Mechanics
> Department of Geoscience and Remote Sensing, TUDelft
> Delft, The Netherlands
>
> M.Petrini at tudelft.nl
> +39 3398367372
> michele.petrins (Skype alias)
>
>
>
>
> _______________________________________________
> ncarg-talk mailing list
> ncarg-talk at ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncarg-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncarg-talk/attachments/20180509/739819e7/attachment.html>