<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> </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. 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.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </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> </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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<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></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> </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> </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 =
2 ;for three days<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>ndata = 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 = "1month"<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>PRM = "totDRE_surface"
;SOBS_RAD";"T_2M";<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>UNIT = " [delta W/m^2]"<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 ;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> </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 = "12" ; MONTH DEC OR JAN
(01)<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>filename =
("winter."+IYR+"."+PRM+name+"")<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></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> lvl = 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> </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></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>; 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> </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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>dum1 =
dim_avg_n(vardumbs,0)<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>plotvar = (dum1) ; in ppb<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>aveplotvar =
dim_avg_n(plotvar,0) ;<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> </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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US> lat =
f1[0]->lat({nbeglat:nendlat},{nbeglon:nendlon}) ;to retrieve lat <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> lon =
f1[0]->lon({nbeglat:nendlat},{nbeglon:nendlon}) ;to retrieve lon<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> nlat =
dimsizes(lat)
;this gives nlat=330<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> nlon =
dimsizes(lon)
;this gives nlon=384<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></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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US> wks = gsn_open_wks
("png",filename)<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></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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>cmap =
(/"white","black","white","lemonchiffon1",\<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
"lightblue1",
"lightblue2","deepskyblue1","deepskyblue3","orchid3",\<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
"darkorchid","blueviolet","darkorchid4","deeppink2"/)<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>gsn_define_colormap(wks,"temp_19lev")
;meteoro<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></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 = ""+PRM+"
"+UNIT+""<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>tdate = "20.12.2013 - 21.1.2014 "<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></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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res
= True<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@tiMainString
= " "<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnLevelSelectionMode =
"ManualLevels"
; set manual contour levels<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@lbLabelAngleF
=
45
; angle labels<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnMinLevelValF
= SCL1<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnMaxLevelValF
= SCL2<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnLevelSpacingF
= SCL3<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnLabelBarEndStyle
= "ExcludeOuterBoxes"<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> res@cnFillPalette =
"temp_19lev" ;meteor<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> res@gsnLeftString
= " "<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnRightString
= " "<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> </o:p></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)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnFrame
= False
; Do not advance frame<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> res@cnFillOn
=
True
; color Fill<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnFillMode
= "AreaFill" ; raster mode<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnRasterSmoothingOn
= True<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnRasterMinCellSizeF
= 0.0005<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnLinesOn
=
True
; Turn off contour lines<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnLineLabelsOn
= False
; Turn off label of contour lines<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnMaxLevelCount
= 100<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> res@cnInfoLabelOn
= False
; do not plot info label<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnAddCyclic
= False
; data is not cyclic<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> res@mpDataBaseVersion
= "MediumRes" ; map resolution<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpProjection
= "CylindricalEquidistant"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpOutlineBoundarySets =
"National"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpFillOn
= True<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpPerimOn
= False<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpMaskAreaSpecifiers
= (/"land"/)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpFillDrawOrder
= "PreDraw"<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> res@mpGridAndLimbOn
= True
; en- or disable lon lat lines<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpGeophysicalLineThicknessF = 2.5<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> res@pmTickMarkDisplayMode =
"Always"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@tmXTMinorOn
= False<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@tmYRMinorOn
= False<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></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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
res@tfDoNDCOverlay = True<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
res@mpDataBaseVersion =
"HighRes"
; map resolution<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
res@mpProjection =
"CylindricalEquidistant" ; map projection<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpOutlineBoundarySets =
"National"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpCenterLonF =
lon({rlon|0},{rlat|0})<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpCenterLatF =
lat({rlon|0},{rlat|0})<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpLimitMode =
"Corners"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpLeftCornerLatF =
lat(0,0)
; latitude left lower corner from model output<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@mpLeftCornerLonF =
lon(0,0)
; longitude left lower corner from model output<o:p></o:p></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<o:p></o:p></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<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> res@sfXCStartV =
nbeglon
; minimal length (rotated): embedd data into map<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@sfXCEndV =
nendlon
; maximal length
(rotated): embedd data into map<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@sfYCStartV =
nbeglat
; minimal length (rotated): embedd data into map<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@sfYCEndV =
nendlat
; maximal length (rotated): embedd data into map<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></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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@vpWidthF =
0.6
; change the aspect ratio<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@vpHeightF = 0.6<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
res@vpXF =
.1
; location of where plot starts<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ;res@vpYF
= 1<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> bot_plot =
gsn_csm_contour_map(wks,plotvar(:,:),res)<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>;-----------add
polymarker-----------------;<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> light_gray =
NhlNewColor(wks,0.85,0.85,0.85)
; Add light gray<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>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 ; Filled circle<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> mkres@tfPolyDrawOrder =
"PreDraw"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 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> <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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>;*****************************************************
<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>;*****************************************************
<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></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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
res_text
= True<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res_text@txFontHeightF =
0.03
; change font size<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> txid_top = gsn_create_text(wks,
title, res_text)<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>
amres
= True<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
amres@amJust =
"BottomCenter"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> amres@amParallelPosF
= 0.0 ; This is the center of the plot.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> amres@amOrthogonalPosF =
-0.72 ; This is above the top edge of the plot.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res_text@txFontHeightF =
0.02
; change font size<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> txid_mid = gsn_create_text(wks,
tdate,res_text)<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> amres@amOrthogonalPosF = -0.62
; This is just below the previous title.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US> pres = True<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> maximize_output(wks,pres)<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>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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </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> </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> </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> </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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US><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"'>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> </o:p></span></p>
</div>
</body>
</html>