<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> </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> </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> </o:p></span></p>
<div>
<p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p> </o:p></span></p>
</div>
<p class=MsoNormal><span style='color:#1F497D'><o:p> </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> </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> </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> </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> </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> </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> </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>;-------------------------------<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ; 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> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>PATHb =
"/work/pr001/eathana/orestis/KIT/KIT/COSMO-ART/COSMO_output/Out_dust.storm_base_indOFF2015013012_05Dec16_1709"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>lfile =
addfile(PATHb+"/"+"lffd2015020112.nc","r")<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>cfile =
addfile(PATHb+"/"+"lffd2015013012c.nc","r")<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 temperature<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>jval = 0 ; rotated latitude index for
cross-section (M_Thission)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>;t =
jmb_getvar(lfile,"T",False) a cosmolib function, currently out
of order 9-12-2016<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>temp =
lfile->T(:,:,jval,:) ;to reduce
the dimensons by 1<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>t = dim_sum_n_Wrap(temp,
0) ;to reduce the dimensions by 1<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>;t := t-273.15 ; degrees Celsius<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>hhl = cfile->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> </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> nlev =
dimsizes(t(:,0)) ;this gives
nlev=40 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> nlon =
dimsizes(t(0,:)) ;this gives nlev=130<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>; close files<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> delete(cfile)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> delete(lfile)<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>; convert units<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>; t = t-273.15 ; degrees Celsius<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> hhl = 0.001 * hhl ; km<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>; compute data positions<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>x2d = conform_dims((/nlev,nlon/),
t&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> </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> ptype = "png"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> wks =
gsn_open_wks(ptype,"cosmo")<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>; SETUP IRREGULAR MESH<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@trGridType
= "TriangularMesh" ; used for irregular mesh triangulation<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@sfXArray
= x2d<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@sfYArray
= y2d<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@tiXAxisString
= "Rotated Longitude [deg]"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@tiYAxisString
= "Height [km]"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ;;res@trXMinF
= <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ;;res@trXMaxF
= <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@trYMinF
= 0.0<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@trYMaxF
= 23.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><o:p> </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> </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> res@vpWidthF
= 0.85<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@vpHeightF
= 0.5<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ;;res@cnFillMode
= "RasterFill"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnMaximize
= True
; maxmize plot in frame<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnFillOn
= True
; turn on color <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnLinesOn
= False ; no
contour lines<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnLineLabelsOn
= False ; no
contour labels<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnLevelSelectionMode
= "ManualLevels" ; manual level selection<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnMinLevelValF
=
180.0
;<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@cnMaxLevelValF
= 300.0 ; <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnSpreadColors
= True
; use full color map<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnSpreadColorEnd
=
2
; skipt black&white<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnSpreadColorEnd
=
-1
; end of color table<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@pmTickMarkDisplayMode =
"conditional"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnAddCyclic
= False<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@lbOrientation
= "vertical"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>;;
res@lbLabelFontHeightF = 0.015<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>;;
res@lbLabelStride = 2<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>; postpone drawing<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnDraw
= False<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsnFrame
= 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>; make contour + map plot<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> pl = gsn_csm_contour(wks, t, res)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> delete(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><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ; add topography<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@gsLineColor
= "black"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> res@gsLineThicknessF
= 1.0<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> pl@topopoly = gsn_add_polyline(wks,
pl, t&rlon, hhl(nlev,:), res)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> delete(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>; draw and frame<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> draw(pl)<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><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US> ; cleanup<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> delete(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>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 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> </o:p></span></p>
</div>
</body>
</html>