<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:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 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;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
span.EmailStyle18
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@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 style='color:#1F497D'>Kind reminder,<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US style='color:#1F497D'>Thank you all<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US style='color:#1F497D'>Orestis<o:p></o:p></span></p>

<p class=MsoNormal><span style='color:#1F497D'><o:p>&nbsp;</o:p></span></p>

<div>

<p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p>&nbsp;</o:p></span></p>

<p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p>&nbsp;</o:p></span></p>

</div>

<p class=MsoNormal><span style='color:#1F497D'><o:p>&nbsp;</o:p></span></p>

<div>

<div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'>

<p class=MsoNormal><b><span lang=EN-US style='font-size:10.0pt;font-family:
"Tahoma","sans-serif"'>From:</span></b><span lang=EN-US style='font-size:10.0pt;
font-family:"Tahoma","sans-serif"'> ncl-talk-bounces@ucar.edu
[mailto:ncl-talk-bounces@ucar.edu] <b>On Behalf Of </b>Orestis<br>
<b>Sent:</b> Friday, December 09, 2016 2:53 PM<br>
<b>To:</b> ncl-talk@ucar.edu<br>
<b>Subject:</b> [ncl-talk] Slice from original rotated coordinates to
geographical coordinates<o:p></o:p></span></p>

</div>

</div>

<p class=MsoNormal><o:p>&nbsp;</o:p></p>

<p class=MsoNormal><span lang=EN-US>Hi everyone,<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 running the COSMO model and I'm trying
to do a slice. This model outputs results with rotated coordinates e.g. T
(time, level, rlat, rlon).<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>I closely followed the example on <a
href="http://www.ncl.ucar.edu/Applications/Scripts/cosmo_2.ncl">http://www.ncl.ucar.edu/Applications/Scripts/cosmo_2.ncl</a>
and I do get a slice of the temperature (image attached).<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>The problem is the x-axis is in rotated
coordinates and therefore difficult to decipher.<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>Could you please help me in putting the
x-axis in geographical coordinates?<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>Thank you all,<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>Orestis<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>The script is as follows:<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></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>begin<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>&nbsp;; Read Data<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>PATHb =
&quot;/work/pr001/eathana/orestis/KIT/KIT/COSMO-ART/COSMO_output/Out_dust.storm_base_indOFF2015013012_05Dec16_1709&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>lfile =
addfile(PATHb+&quot;/&quot;+&quot;lffd2015020112.nc&quot;,&quot;r&quot;)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>cfile =
addfile(PATHb+&quot;/&quot;+&quot;lffd2015013012c.nc&quot;,&quot;r&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>; read temperature<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>jval&nbsp; = 0 ; rotated latitude index for
cross-section (M_Thission)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>;t&nbsp; =
jmb_getvar(lfile,&quot;T&quot;,False)&nbsp; a cosmolib function, currently out
of order 9-12-2016<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>temp&nbsp; =
lfile-&gt;T(:,:,jval,:)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;to reduce
the dimensons by 1<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>t&nbsp; = dim_sum_n_Wrap(temp,
0)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;;to reduce the dimensions by 1<o:p></o:p></span></p>

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

<p class=MsoNormal><span lang=EN-US>hhl = cfile-&gt;HHL(0,:,jval,:)<o:p></o:p></span></p>

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

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

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

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

<p class=MsoNormal><span lang=EN-US>;print (hhl(40,:))<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>; get dimensions<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;nlev =
dimsizes(t(:,0))&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;this gives
nlev=40&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;nlon =
dimsizes(t(0,:))&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;this gives nlev=130<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>; close files<o:p></o:p></span></p>

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

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

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

<p class=MsoNormal><span lang=EN-US>&nbsp;hhl = 0.001 * hhl ; km<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>; compute data positions<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>x2d = conform_dims((/nlev,nlon/),
t&amp;rlon, 1) ;conform_dims (dims, r, ndim)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>y2d = 0.5*(hhl(0:nlev-1,:)+hhl(1:nlev,:))<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>; open graphic port<o:p></o:p></span></p>

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

<p class=MsoNormal><span lang=EN-US>&nbsp;wks =
gsn_open_wks(ptype,&quot;cosmo&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>; SETUP IRREGULAR MESH<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>res&nbsp;&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@trGridType&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= &quot;TriangularMesh&quot; ; used for irregular mesh triangulation<o:p></o:p></span></p>

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

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

<p class=MsoNormal><span lang=EN-US>&nbsp;res@tiXAxisString&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= &quot;Rotated Longitude [deg]&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@tiYAxisString&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= &quot;Height [km]&quot;<o:p></o:p></span></p>

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

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

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

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

<p class=MsoNormal><span lang=EN-US>; COUNTOUR PLOT RESOURCES<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>; setup contour plot resources<o:p></o:p></span></p>

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

<p class=MsoNormal><span lang=EN-US>&nbsp;res@vpHeightF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= 0.5<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;RasterFill&quot;<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@gsnMaximize&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; maxmize plot in frame<o:p></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;
; turn on color <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;
= False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; no
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; ; no
contour labels<o:p></o:p></span></p>

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

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

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

<p class=MsoNormal><span lang=EN-US>&nbsp;res@gsnSpreadColors&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
= True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; use full color map<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@gsnSpreadColorEnd&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
=
2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; skipt black&amp;white<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@gsnSpreadColorEnd&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
=
-1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
; end of color table<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;res@pmTickMarkDisplayMode&nbsp; =
&quot;conditional&quot;<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<o:p></o:p></span></p>

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

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

<p class=MsoNormal><span lang=EN-US>;;
res@lbLabelStride&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2<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>; postpone drawing<o:p></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<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<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>; make contour + map plot<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;pl = gsn_csm_contour(wks, t, res)<o:p></o:p></span></p>

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

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

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

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

<p class=MsoNormal><span lang=EN-US>&nbsp;pl@topopoly = gsn_add_polyline(wks,
pl, t&amp;rlon, hhl(nlev,:), res)<o:p></o:p></span></p>

<p class=MsoNormal><span lang=EN-US>&nbsp;delete(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>; draw and frame<o:p></o:p></span></p>

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

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

<p class=MsoNormal><span lang=EN-US>&nbsp;delete(wks)<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>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 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><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<br>
Phone: +30 210 8109170</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>