<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">

<head>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
 /* Font Definitions */
 @font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:"Trebuchet MS";
        panose-1:2 11 6 3 2 2 2 2 2 4;}
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;}
@page Section1
        {size:612.0pt 792.0pt;
        margin:72.0pt 90.0pt 72.0pt 90.0pt;}
div.Section1
        {page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
 <o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
 <o:shapelayout v:ext="edit">
  <o:idmap v:ext="edit" data="1" />
 </o:shapelayout></xml><![endif]-->
</head>

<body lang=EL link=blue vlink=purple>

<div class=Section1>

<p class=MsoNormal><span lang=EN-US>Dear all,<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></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. &nbsp;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 &quot;dynamically&quot;). The frustrating/interesting part is that
with overlay (ovelay countour lines to a basic countour map plot in another
script) the dot appears.<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></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>.<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; horizontal contour plot with station<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>load
&quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>load
&quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;
================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>begin<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; choose wks <o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>nbeglon = (/-0.20/)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>nbeglat = (/-0.2/) ; chooses begin of
lon/lat(rotated coordinates of the data)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>nendlon = (/0.55/)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>nendlat = (/0.58/)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>starthour = 0<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>endhour&nbsp;&nbsp; =
2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;for three days<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>ndata&nbsp;&nbsp;&nbsp;&nbsp; = 72<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>aa = starthour+1;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>bb = endhour;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>name = &quot;1month&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>PRM = &quot;totDRE_surface&quot;
;SOBS_RAD&quot;;&quot;T_2M&quot;;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>UNIT = &quot;&nbsp; [delta W/m^2]&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>SCL1 = -2.6;2.5;.4;0.2 ;-0.5;-1<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>SCL2 = 2.6;2.5;.4;0.2&nbsp; ;0.5;1<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>SCL3 = 0.2;0.05;0.25 ;0.1<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>IYR = 2013;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>MON = &quot;12&quot; ; MONTH DEC OR JAN
(01)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>filename =
(&quot;winter.&quot;+IYR+&quot;.&quot;+PRM+name+&quot;&quot;)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; choose level<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; lvl&nbsp; = 39<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; open file and read in data, I cut this
out as It does not matter<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; read in data&nbsp;&nbsp;&nbsp; I cut this
out as It does not matter<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></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)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>dum1&nbsp;&nbsp;&nbsp; =
dim_avg_n(vardumbs,0)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>plotvar = (dum1)&nbsp;&nbsp; ; in ppb<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>aveplotvar =
dim_avg_n(plotvar,0)&nbsp;&nbsp; ;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>aveplotvar2 = dim_avg_n(aveplotvar,0)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>printVarSummary (plotvar)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>printMinMax (plotvar, True)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>delete(dum1)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; lat =
f1[0]-&gt;lat({nbeglat:nendlat},{nbeglon:nendlon}) ;to retrieve lat <o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; lon =
f1[0]-&gt;lon({nbeglat:nendlat},{nbeglon:nendlon}) ;to retrieve lon<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; nlat =
dimsizes(lat)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
;this gives nlat=330<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; nlon =
dimsizes(lon)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
;this gives nlon=384<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; Create wks according to definition<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;wks&nbsp;&nbsp; = gsn_open_wks
(&quot;png&quot;,filename)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; Define Colormap<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>cmap =
(/&quot;white&quot;,&quot;black&quot;,&quot;white&quot;,&quot;lemonchiffon1&quot;,\<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&quot;lightblue1&quot;,
&quot;lightblue2&quot;,&quot;deepskyblue1&quot;,&quot;deepskyblue3&quot;,&quot;orchid3&quot;,\<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&quot;darkorchid&quot;,&quot;blueviolet&quot;,&quot;darkorchid4&quot;,&quot;deeppink2&quot;/)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>gsn_define_colormap(wks,&quot;temp_19lev&quot;)&nbsp;&nbsp;
;meteoro<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;===================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; Manage ressources and title<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;===================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>title = &quot;&quot;+PRM+&quot;
&quot;+UNIT+&quot;&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>tdate = &quot;20.12.2013 - 21.1.2014 &quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; Set some resources that will apply to the
base<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; contour/map plot <o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= True<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@tiMainString&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= &quot; &quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnLevelSelectionMode&nbsp; =
&quot;ManualLevels&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; set manual contour levels<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@lbLabelAngleF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
=
45&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; angle labels<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnMinLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
=&nbsp; SCL1<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnMaxLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
=&nbsp; SCL2<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnLevelSpacingF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
=&nbsp; SCL3<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnLabelBarEndStyle&nbsp;&nbsp;&nbsp;
=&nbsp; &quot;ExcludeOuterBoxes&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; <o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnFillPalette =
&quot;temp_19lev&quot;&nbsp;&nbsp; ;meteor<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@gsnLeftString&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= &quot; &quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@gsnRightString&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= &quot; &quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;===============================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; standard stuff<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;===============================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@gsnDraw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; Do not draw plot now (helps adding stuff to the plot or producing labels etc)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@gsnFrame&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; Do not advance frame<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
=
True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; color Fill<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnFillMode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= &quot;AreaFill&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; raster mode<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnRasterSmoothingOn&nbsp;&nbsp;&nbsp;
= True<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnRasterMinCellSizeF&nbsp;&nbsp;
= 0.0005<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnLinesOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
=
True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; Turn off contour lines<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnLineLabelsOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&nbsp;&nbsp;&nbsp;&nbsp;; Turn off label of contour lines<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnMaxLevelCount&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= 100<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@cnInfoLabelOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; do not plot info label<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@gsnAddCyclic&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; data is not cyclic<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@mpDataBaseVersion&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= &quot;MediumRes&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; map resolution<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@mpProjection&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= &quot;CylindricalEquidistant&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@mpOutlineBoundarySets&nbsp; =
&quot;National&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@mpFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= True<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@mpPerimOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= False<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@mpMaskAreaSpecifiers&nbsp;&nbsp;
= (/&quot;land&quot;/)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@mpFillDrawOrder&nbsp;&nbsp;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;= &quot;PreDraw&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@mpGridAndLimbOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; en- or disable lon lat lines<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@mpGeophysicalLineThicknessF = 2.5<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@pmTickMarkDisplayMode&nbsp; =
&quot;Always&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@tmXTMinorOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= False<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@tmYRMinorOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= False<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; map projection<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;
res@tfDoNDCOverlay&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;
res@mpDataBaseVersion&nbsp;&nbsp;&nbsp;&nbsp; =
&quot;HighRes&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; map resolution<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;
res@mpProjection&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
&quot;CylindricalEquidistant&quot;&nbsp; ; map projection<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@mpOutlineBoundarySets =
&quot;National&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@mpCenterLonF =
lon({rlon|0},{rlat|0})<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@mpCenterLatF =
lat({rlon|0},{rlat|0})<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@mpLimitMode =
&quot;Corners&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@mpLeftCornerLatF =
lat(0,0)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&nbsp;&nbsp;; latitude left lower corner from model output<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@mpLeftCornerLonF =
lon(0,0)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; longitude left lower corner from model output<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@mpRightCornerLatF =
lat(nlat(0)-1,nlat(1)-1)&nbsp;&nbsp;&nbsp; ; latitude right upper corner from
model output<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@mpRightCornerLonF =
lon(nlon(0)-1,nlat(1)-1)&nbsp;&nbsp;&nbsp; ; longitutde right upper corner from
model output<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@sfXCStartV =
nbeglon&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; minimal length (rotated): embedd data into map<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@sfXCEndV =&nbsp;
nendlon&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;; maximal length
(rotated): embedd data into map<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@sfYCStartV =
nbeglat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; minimal length (rotated): embedd data into map<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@sfYCEndV =
nendlat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; maximal length (rotated): embedd data into map<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; Scaling<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;=================================================;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@vpWidthF&nbsp; =
0.6&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; change the aspect ratio<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res@vpHeightF = 0.6<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;
res@vpXF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
.1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; location of where plot starts<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; ;res@vpYF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= 1<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;bot_plot =
gsn_csm_contour_map(wks,plotvar(:,:),res)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;-----------add
polymarker-----------------;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; light_gray =
NhlNewColor(wks,0.85,0.85,0.85)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; Add light gray<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>mkres =True<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>thissio_lat = 0.2732<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>thissio_lon =0.1720<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>mkres@gsMarkerIndex =
17&nbsp;&nbsp;&nbsp;&nbsp; ; Filled circle<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; mkres@tfPolyDrawOrder =
&quot;PreDraw&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; mkres@gsMarkerSizeF = 0.03 <o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>Station =
gsn_add_polymarker(wks,bot_plot,thissio_lon,thissio_lat,mkres)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>draw (bot_plot)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;*****************************************************&nbsp;&nbsp;
<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; Manually create and attach titles<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;*****************************************************&nbsp;&nbsp;
<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>; Attach some titles at the top.<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;
res_text&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= True<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res_text@txFontHeightF =
0.03&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; change font size<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; txid_top = gsn_create_text(wks,
title, res_text)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;
amres&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= True<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;
amres@amJust&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
&quot;BottomCenter&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; amres@amParallelPosF&nbsp;&nbsp;
=&nbsp; 0.0&nbsp;&nbsp;&nbsp; ; This is the center of the plot.<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; amres@amOrthogonalPosF =
-0.72&nbsp;&nbsp; ; This is above the top edge of the plot.<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; annoid_top = gsn_add_annotation(bot_plot,
txid_top, amres)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; res_text@txFontHeightF =
0.02&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; change font size<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; txid_mid = gsn_create_text(wks,
tdate,res_text)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; amres@amOrthogonalPosF = -0.62&nbsp;
; This is just below the previous title.<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; annoid_mid =
gsn_add_annotation(bot_plot, txid_mid, amres)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; pres = True<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp; maximize_output(wks,pres)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>frame (wks)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>delete(plotvar)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>end<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US>Once again thank you,<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></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"'><o:p></o:p></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><o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US><o:p>&nbsp;</o:p></span></p>

</div>

</body>

</html>