<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Dear helpdesk,
<div class=""><br class="">
</div>
<div class="">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 </div>
<div class=""><br class="">
</div>
<div class="">
<div class=""><font color="#3a88fe" class="">begin </font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">in = "./data/"</font></div>
<div class=""><font color="#3a88fe" class="">out = "./results/"</font></div>
<div class=""><font color="#3a88fe" class="">H_filename = "ice_thickness.nc"</font></div>
<div class=""><font color="#3a88fe" class="">H_file = addfile(in+H_filename,"r")</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">H = H_file->thk</font></div>
<div class=""><font color="#3a88fe" class="">time_sim = ispan(0,16100,100)</font></div>
<div class=""><font color="#3a88fe" class="">ny = dimsizes(H&y1)</font></div>
<div class=""><font color="#3a88fe" class="">nx = dimsizes(H&x1)</font></div>
<div class=""><font color="#3a88fe" class="">ntime = dimsizes(time_sim)</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">H@_FillValue = -999.</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">H = mask(H,H.lt..1,False)</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">lat1 = fspan(58.25,84.61,ny)</font></div>
<div class=""><font color="#3a88fe" class="">lat1@units = "degrees_north"</font></div>
<div class=""><font color="#3a88fe" class="">lat1!0 = "lat"</font></div>
<div class=""><font color="#3a88fe" class="">lat1@long_name = "latitude"</font></div>
<div class=""><font color="#3a88fe" class="">lat1@standard_name = "latitude"</font></div>
<div class=""><font color="#3a88fe" class="">lat1@axis = "Y"</font></div>
<div class=""><font color="#3a88fe" class="">print(lat1)</font></div>
<div class=""><font color="#3a88fe" class="">lon1 = fspan(-95.68,13.33,nx)</font></div>
<div class=""><font color="#3a88fe" class="">lon1@units = "degrees_east"</font></div>
<div class=""><font color="#3a88fe" class="">lon1!0 = "lon"</font></div>
<div class=""><font color="#3a88fe" class="">lon1@long_name = "longitude"</font></div>
<div class=""><font color="#3a88fe" class="">lon1@standard_name = "longitude"</font></div>
<div class=""><font color="#3a88fe" class="">lon1@axis = "X"</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">H!1 = "lat"</font></div>
<div class=""><font color="#3a88fe" class="">H!2 = "lon"</font></div>
<div class=""><font color="#3a88fe" class="">H&lat = lat1</font></div>
<div class=""><font color="#3a88fe" class="">H&lon = lon1</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">minlat = 58.25</font></div>
<div class=""><font color="#3a88fe" class="">maxlat = 84.61</font></div>
<div class=""><font color="#3a88fe" class="">minlon = -95.68</font></div>
<div class=""><font color="#3a88fe" class="">maxlon = 13.33</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class=""> res = True ; plot mods desired</font></div>
<div class=""><font color="#3a88fe" class=""> ;res@gsnDraw = False</font></div>
<div class=""><font color="#3a88fe" class=""> ;res@gsnFrame = False </font></div>
<div class=""><font color="#3a88fe" class=""> res@cnFillOn = True ; turn on color</font></div>
<div class=""><font color="#3a88fe" class=""> res@cnFillMode = "RasterFill"</font></div>
<div class=""><font color="#3a88fe" class=""> res@cnLinesOn = False</font></div>
<div class=""><font color="#3a88fe" class=""> res@cnLineLabelsOn = False</font></div>
<div class=""><font color="#3a88fe" class=""> res@lbLabelBarOn = False ; turn off individual lb's</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class=""> res@trGridType = "TriangularMesh"</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class=""> res@mpProjection = "Stereographic"</font></div>
<div class=""><font color="#3a88fe" class=""> res@mpDataBaseVersion = "mediumres"</font></div>
<div class=""><font color="#3a88fe" class=""> res@mpFillOn = False ; turn off default land map fill</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class=""> res@mpLimitMode = "Corners" </font></div>
<div class=""><font color="#3a88fe" class=""> res@mpLeftCornerLatF = minlat</font></div>
<div class=""><font color="#3a88fe" class=""> res@mpLeftCornerLonF = minlon</font></div>
<div class=""><font color="#3a88fe" class=""> res@mpRightCornerLatF = maxlat</font></div>
<div class=""><font color="#3a88fe" class=""> res@mpRightCornerLonF = maxlon</font></div>
<div class=""><font color="#3a88fe" class=""> res@mpCenterLonF = 0</font></div>
<div class=""><font color="#3a88fe" class=""> res@mpCenterLatF = 90</font></div>
<div class=""><font color="#3a88fe" class=""> res@gsnLeftString = ""</font></div>
<div class=""><font color="#3a88fe" class=""> res@gsnRightString = ""</font></div>
<div class=""><font color="#3a88fe" class=""> res@cnFillPalette = "testcmap"</font></div>
<div class=""><font color="#3a88fe" class=""> res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels</font></div>
<div class=""><font color="#3a88fe" class=""> res@cnMinLevelValF = 0.</font></div>
<div class=""><font color="#3a88fe" class=""> res@cnMaxLevelValF = 1.</font></div>
<div class=""><font color="#3a88fe" class=""> res@cnLevelSpacingF = 0.2 </font></div>
<div class=""><font color="#3a88fe" class=""> res@tfDoNDCOverlay = "DataTransform"</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">do i=0,ntime-1</font></div>
<div class=""><font color="#3a88fe" class=""><span class="Apple-tab-span" style="white-space:pre"></span>wks = gsn_open_wks("pdf","./results/image-"+sprinti("%03i",i))</font></div>
<div class=""><font color="#3a88fe" class=""> print("simulation time(" + i + ") = " + time_sim(i))</font></div>
<div class=""><font color="#3a88fe" class=""> res@gsnRightString = "Simulation time: " + time_sim(i) + " (yr)"</font></div>
<div class=""><font color="#3a88fe" class=""><span class="Apple-tab-span" style="white-space:pre"></span>ice = new((/ny,nx/),double)</font></div>
<div class=""><font color="#3a88fe" class=""><span class="Apple-tab-span" style="white-space:pre"></span>ice(:,:) = H(i,:,:)</font></div>
<div class=""><font color="#3a88fe" class=""><span class="Apple-tab-span" style="white-space:pre"></span>plot = new(1,graphic)</font></div>
<div class=""><font color="#3a88fe" class=""><span class="Apple-tab-span" style="white-space:pre"></span>plot(0) = gsn_csm_contour_map(wks,ice,res)</font></div>
<div class=""><font color="#3a88fe" class=""><span class="Apple-tab-span" style="white-space:pre"></span>delete(wks)</font></div>
<div class=""><font color="#3a88fe" class="">end do</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">cmd = "cd ./results/ ; pdflatex movie.tex"</font></div>
<div class=""><font color="#3a88fe" class="">system(cmd)</font></div>
<div class=""><font color="#3a88fe" class=""><br class="">
</font></div>
<div class=""><font color="#3a88fe" class="">end </font></div>
</div>
<div class=""><br class="">
</div>
<div class="">whereas the error I get is the following </div>
<div class=""><br class="">
</div>
<div class="">
<div style="margin: 0px; font-stretch: normal; font-size: 13px; line-height: normal; font-family: Menlo; color: rgb(248, 195, 36); background-color: rgb(146, 146, 146);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">fatal:Or: operation not supported on type (string)</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 13px; line-height: normal; font-family: Menlo; color: rgb(248, 195, 36); background-color: rgb(146, 146, 146);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">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</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 13px; line-height: normal; font-family: Menlo; color: rgb(248, 195, 36); background-color: rgb(146, 146, 146); min-height: 15px;" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class="">
</div>
<div style="margin: 0px; font-stretch: normal; font-size: 13px; line-height: normal; font-family: Menlo; color: rgb(248, 195, 36); background-color: rgb(146, 146, 146);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">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</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 13px; line-height: normal; font-family: Menlo; color: rgb(248, 195, 36); background-color: rgb(146, 146, 146); min-height: 15px;" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class="">
</div>
<div style="margin: 0px; font-stretch: normal; font-size: 13px; line-height: normal; font-family: Menlo; color: rgb(248, 195, 36); background-color: rgb(146, 146, 146);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">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</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 13px; line-height: normal; font-family: Menlo; color: rgb(248, 195, 36); background-color: rgb(146, 146, 146); min-height: 15px;" class="">
<span style="font-variant-ligatures: no-common-ligatures" class=""></span><br class="">
</div>
<div style="margin: 0px; font-stretch: normal; font-size: 13px; line-height: normal; font-family: Menlo; color: rgb(248, 195, 36); background-color: rgb(146, 146, 146);" class="">
<span style="font-variant-ligatures: no-common-ligatures" class="">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 80 in file contour_movie.ncl</span></div>
</div>
<div class=""><span style="font-variant-ligatures: no-common-ligatures" class=""><br class="">
</span></div>
<div class="">I can’t find any suggestion by looking at previous questions and I really don’t understand what the problem can be. </div>
<div class="">thanks!</div>
<div class=""><br class="">
</div>
<div class="">best,</div>
<div class="">Michele</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;">
******************************************************************<br class="">
Michele Petrini<br class="">
<br class="">
PhD in Earth Science and Fluid Mechanics<br class="">
Department of Geoscience and Remote Sensing, TUDelft<br class="">
Delft, The Netherlands<br class="">
<br class="">
<a href="mailto:M.Petrini@tudelft.nl" class="">M.Petrini@tudelft.nl</a> <br class="">
+39 3398367372<br class="">
michele.petrins (Skype alias)<br class="">
<br class="">
<br class="">
</div>
</div>
<br class="">
</div>
</body>
</html>