<div dir="ltr"><div class="gmail_default" style="font-size:small">Orestis,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Whenever people try to draw markers and lines on a map plot using lat/lon coordinates and they don't show up, it's usually because the lat/lon values they used were out-of-range of the plot.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You are using these values:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US">thissio_lat = 0.2732<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US">thissio_lon =0.1720</span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><br></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US">Which seems a bit suspicious to me. You must use values of degrees, and not radians or any other units. Maybe what you meant to do was this:</span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><br></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><span style="font-size:12.8px">thissio_lat = </span><span style="font-size:12.8px">f1[0]->lat({</span><span style="font-size:12.8px">0.2732</span><span style="font-size:12.8px">},{</span><span style="font-size:12.8px">0.1720</span><span style="font-size:12.8px">})</span><br></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><span style="font-size:12.8px">thissio_lon = </span><span style="font-size:12.8px">f1[0]->lon({</span><span style="font-size:12.8px">0.2732</span><span style="font-size:12.8px">},{</span><span style="font-size:12.8px">0.1720</span><span style="font-size:12.8px">})</span><br></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><span style="font-size:12.8px"><br></span></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><span style="font-size:12.8px">The second issue is that you are setting this, which tells NCL to draw the marker before drawing other plot components:</span></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><span style="font-size:12.8px"><br></span></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"> mkres@tfPolyDrawOrder = "PreDraw"<u></u><u></u></span></p><div><span lang="EN-US"><br></span></div><div>Try setting this instead:</div><div><br></div><div><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"> mkres@tfPolyDrawOrder = "PostDraw"<u></u><u></u></span></p></div><div><span lang="EN-US"><br></span></div><div><span lang="EN-US"><br></span></div><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US">The third issue may be the size of the marker you chose:</span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><br></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"> mkres@gsMarkerSizeF = 0.03 <u></u><u></u></span></p><div><span lang="EN-US"><br></span></div><div><span lang="EN-US">I would comment this out initially, and if you don't see anything, try a value like:</span></div><div><span lang="EN-US"><br></span></div><div><span lang="EN-US"><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"> mkres@gsMarkerSizeF = 5.0</span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><br></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US">There are a number of other issues with your script that makes me worry that the data is not being plotted correctly.</span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><br></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US">For starters, if you set this:</span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><br></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"> res@tfDoNDCOverlay = True<u></u><u></u></span></p><div><span lang="EN-US"><br></span></div><div><span lang="EN-US">Then you are telling NCL that you know the exact map projection that your data is in, and hence you do NOT want to set these resources:</span></div><div><span lang="EN-US"><br></span></div><div><span lang="EN-US"><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"> res@sfXCStartV = nbeglon <wbr> ; minimal length (rotated): embedd data into map<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"> res@sfXCEndV = nendlon ; maximal length (rotated): embedd data into map<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"> res@sfYCStartV = nbeglat <wbr> ; minimal length (rotated): embedd data into map<u></u><u></u></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"> res@sfYCEndV = nendlat <wbr> ; maximal length (rotated): embedd data into map</span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><br></span></p><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><br></span></p></span></div><p class="MsoNormal" style="font-size:12.8px"><span lang="EN-US"><br></span></p><p class="MsoNormal" style="font-size:12.8px">--Mary</p><p class="MsoNormal" style="font-size:12.8px"><br></p><p class="MsoNormal" style="font-size:12.8px"><br></p></span></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 28, 2016 at 12:37 AM, Orestis Speyer <span dir="ltr"><<a href="mailto:ospeyer@noa.gr" target="_blank">ospeyer@noa.gr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div lang="EL" link="blue" vlink="purple">
<div class="m_8244606710584712229Section1">
<p class="MsoNormal"><span lang="EN-US">Dear all,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">I am trying to plot a dot (a station) in a
map but the dot does not appear. I looked into the archive search but
could not find a solution. I attach the relevant code below. I am not sure if
the problem is a map/transform/projection issue, an issue of Draw priority or
simply NCL does not know where to plot the dot (as the main plot coordinates
are created "dynamically"). The frustrating/interesting part is that
with overlay (ovelay countour lines to a basic countour map plot in another
script) the dot appears.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span lang="EN-US">Here is the code (sorry for the length).
Thank you for your time</span></b><span lang="EN-US">.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; horizontal contour plot with station<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">load
"$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">load
"$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/contributed.<wbr>ncl"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;
==============================<wbr>==================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">begin<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; choose wks <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">nbeglon = (/-0.20/)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">nbeglat = (/-0.2/) ; chooses begin of
lon/lat(rotated coordinates of the data)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">nendlon = (/0.55/)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">nendlat = (/0.58/)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">starthour = 0<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">endhour =
2 ;for three days<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">ndata = 72<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">aa = starthour+1;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">bb = endhour;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">name = "1month"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">PRM = "totDRE_surface"
;SOBS_RAD";"T_2M";<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">UNIT = " [delta W/m^2]"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">SCL1 = -2.6;2.5;.4;0.2 ;-0.5;-1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">SCL2 = 2.6;2.5;.4;0.2 ;0.5;1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">SCL3 = 0.2;0.05;0.25 ;0.1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">IYR = 2013;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">MON = "12" ; MONTH DEC OR JAN
(01)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">filename =
("winter."+IYR+"."+PRM+name+""<wbr>)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; choose level<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> lvl = 39<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; open file and read in data, I cut this
out as It does not matter<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; read in data I cut this
out as It does not matter<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;vardumbs is the variable I'm going to
countour (The basic countour plot to which I want to add a dot)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">dum1 =
dim_avg_n(vardumbs,0)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">plotvar = (dum1) ; in ppb<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">aveplotvar =
dim_avg_n(plotvar,0) ;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">aveplotvar2 = dim_avg_n(aveplotvar,0)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">printVarSummary (plotvar)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">printMinMax (plotvar, True)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">delete(dum1)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> lat =
f1[0]->lat({nbeglat:nendlat},{<wbr>nbeglon:nendlon}) ;to retrieve lat <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> lon =
f1[0]->lon({nbeglat:nendlat},{<wbr>nbeglon:nendlon}) ;to retrieve lon<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> nlat =
dimsizes(lat) <wbr>
;this gives nlat=330<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> nlon =
dimsizes(lon) <wbr>
;this gives nlon=384<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; Create wks according to definition<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> wks = gsn_open_wks
("png",filename)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; Define Colormap<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">cmap =
(/"white","black","white","<wbr>lemonchiffon1",\<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">
"lightblue1",
"lightblue2","deepskyblue1","<wbr>deepskyblue3","orchid3",\<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">
"darkorchid","blueviolet","<wbr>darkorchid4","deeppink2"/)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">gsn_define_colormap(wks,"temp_<wbr>19lev")
;meteoro<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>======;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; Manage ressources and title<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>======;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">title = ""+PRM+"
"+UNIT+""<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">tdate = "20.12.2013 - 21.1.2014 "<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; Set some resources that will apply to the
base<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; contour/map plot <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>===================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res
= True<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@tiMainString
= " "<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnLevelSelectionMode =
"ManualLevels" <wbr>
; set manual contour levels<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@lbLabelAngleF
=
45 <wbr>
; angle labels<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnMinLevelValF
= SCL1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnMaxLevelValF
= SCL2<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnLevelSpacingF
= SCL3<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnLabelBarEndStyle
= "ExcludeOuterBoxes"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnFillPalette =
"temp_19lev" ;meteor<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@gsnLeftString
= " "<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@gsnRightString
= " "<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>==================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; standard stuff<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>==================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@gsnDraw
= False
; Do not draw plot now (helps adding stuff to the plot or producing labels etc)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@gsnFrame
= False
; Do not advance frame<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnFillOn
=
True
; color Fill<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnFillMode
= "AreaFill" ; raster mode<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnRasterSmoothingOn
= True<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnRasterMinCellSizeF
= 0.0005<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnLinesOn
=
True
; Turn off contour lines<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnLineLabelsOn
= False
; Turn off label of contour lines<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnMaxLevelCount
= 100<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@cnInfoLabelOn
= False
; do not plot info label<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@gsnAddCyclic
= False
; data is not cyclic<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpDataBaseVersion
= "MediumRes" ; map resolution<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpProjection
= "CylindricalEquidistant"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpOutlineBoundarySets =
"National"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpFillOn
= True<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpPerimOn
= False<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpMaskAreaSpecifiers
= (/"land"/)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpFillDrawOrder
= "PreDraw"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpGridAndLimbOn
= True
; en- or disable lon lat lines<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@<wbr>mpGeophysicalLineThicknessF = 2.5<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@pmTickMarkDisplayMode =
"Always"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@tmXTMinorOn
= False<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@tmYRMinorOn
= False<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; map projection<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">
res@tfDoNDCOverlay = True<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">
res@mpDataBaseVersion =
"HighRes"
; map resolution<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">
res@mpProjection =
"CylindricalEquidistant" ; map projection<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpOutlineBoundarySets =
"National"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpCenterLonF =
lon({rlon|0},{rlat|0})<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpCenterLatF =
lat({rlon|0},{rlat|0})<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpLimitMode =
"Corners"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpLeftCornerLatF =
lat(0,0)
; latitude left lower corner from model output<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpLeftCornerLonF =
lon(0,0)
; longitude left lower corner from model output<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpRightCornerLatF =
lat(nlat(0)-1,nlat(1)-1) ; latitude right upper corner from
model output<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@mpRightCornerLonF =
lon(nlon(0)-1,nlat(1)-1) ; longitutde right upper corner from
model output<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@sfXCStartV =
nbeglon <wbr>
; minimal length (rotated): embedd data into map<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@sfXCEndV =
nendlon
; maximal length
(rotated): embedd data into map<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@sfYCStartV =
nbeglat <wbr>
; minimal length (rotated): embedd data into map<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@sfYCEndV =
nendlat <wbr>
; maximal length (rotated): embedd data into map<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; Scaling<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;=============================<wbr>====================;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@vpWidthF =
0.6
; change the aspect ratio<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res@vpHeightF = 0.6<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">
res@vpXF =
.1
; location of where plot starts<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> ;res@vpYF
= 1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> bot_plot =
gsn_csm_contour_map(wks,<wbr>plotvar(:,:),res)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;-----------add
polymarker-----------------;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> light_gray =
NhlNewColor(wks,0.85,0.85,0.<wbr>85)
; Add light gray<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">mkres =True<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">thissio_lat = 0.2732<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">thissio_lon =0.1720<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">mkres@gsMarkerIndex =
17 ; Filled circle<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> mkres@tfPolyDrawOrder =
"PreDraw"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> mkres@gsMarkerSizeF = 0.03 <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Station =
gsn_add_polymarker(wks,bot_<wbr>plot,thissio_lon,thissio_lat,<wbr>mkres)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">draw (bot_plot)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;*****************************<wbr>************************
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; Manually create and attach titles<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;*****************************<wbr>************************
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">; Attach some titles at the top.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">
res_text
= True<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res_text@txFontHeightF =
0.03
; change font size<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> txid_top = gsn_create_text(wks,
title, res_text)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">
amres
= True<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">
amres@amJust =
"BottomCenter"<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> amres@amParallelPosF
= 0.0 ; This is the center of the plot.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> amres@amOrthogonalPosF =
-0.72 ; This is above the top edge of the plot.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> annoid_top = gsn_add_annotation(bot_plot,
txid_top, amres)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> res_text@txFontHeightF =
0.02
; change font size<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> txid_mid = gsn_create_text(wks,
tdate,res_text)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> amres@amOrthogonalPosF = -0.62
; This is just below the previous title.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> annoid_mid =
gsn_add_annotation(bot_plot, txid_mid, amres)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> pres = True<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"> maximize_output(wks,pres)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">frame (wks)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">delete(plotvar)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">end<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Once again thank you,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span lang="EN-US" style="font-size:12.0pt;font-family:"Trebuchet MS","sans-serif"">Orestis
Speyer,</span><span lang="EN-US" style="font-size:12.0pt;font-family:"Times New Roman","serif""><u></u><u></u></span></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span lang="EN-US" style="font-size:12.0pt;font-family:"Trebuchet MS","sans-serif"">Research
Fellow, National Observatory of Athens<br>
Institute for Environmental Research and Sustainable Development</span><span lang="EN-US"><u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
</div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>