<div dir="ltr"><div class="gmail_default" style="font-size:small">Michael,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Thanks for including the script and data offline.</div>
<div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It turns out that you *do* have coordinate arrays attached to your original data, but they were getting stripped off when you did the variables. To see the coordinate arrays, do a "printVarSummary" on the "tho" variable, for example:</div>
<div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style><div class="gmail_default" style><font face="courier new, monospace">Variable: tho</font></div><div class="gmail_default" style>
<font face="courier new, monospace">Type: float</font></div><div class="gmail_default" style><font face="courier new, monospace">Total Size: 409920 bytes</font></div><div class="gmail_default" style><font face="courier new, monospace"> 102480 values</font></div>
<div class="gmail_default" style><font face="courier new, monospace">Number of Dimensions: 4</font></div><div class="gmail_default" style><font face="courier new, monospace">Dimensions and sizes:<span class="" style="white-space:pre"> </span>[time | 1] x [depth | 20] x [lat | 122] x [lon | 42]</font></div>
<div class="gmail_default" style><font face="courier new, monospace">Coordinates: </font></div><div class="gmail_default" style><font face="courier new, monospace"> time: [8760..8760]</font></div><div class="gmail_default" style>
<font face="courier new, monospace"> depth: [ 7..2435]</font></div><div class="gmail_default" style><font face="courier new, monospace"> lat: [45.375..-45.375]</font></div><div class="gmail_default" style>
<font face="courier new, monospace"> lon: [-51.375..-20.625]</font></div><div class="gmail_default" style><font face="courier new, monospace">Number Of Attributes: 6</font></div><div class="gmail_default" style>
<font face="courier new, monospace"> standard_name :<span class="" style="white-space:pre"> </span>sea_water_potential_temperature</font></div><div class="gmail_default" style><font face="courier new, monospace"> long_name :<span class="" style="white-space:pre"> </span>Sea water potential temperature</font></div>
<div class="gmail_default" style><font face="courier new, monospace"> units :<span class="" style="white-space:pre"> </span>C</font></div><div class="gmail_default" style><font face="courier new, monospace"> code :<span class="" style="white-space:pre"> </span>2</font></div>
<div class="gmail_default" style><font face="courier new, monospace"> _FillValue :<span class="" style="white-space:pre"> </span>-9e+33</font></div><div class="gmail_default" style><font face="courier new, monospace"> missing_value :<span class="" style="white-space:pre"> </span>-9e+33</font></div>
<div class="gmail_default" style><br></div><div class="gmail_default" style>Note the list of variables under the "Coordinates:" section. These variables are called "coordinate arrays", and make it easier to plot your data. </div>
<div class="gmail_default" style><br></div><div class="gmail_default" style>However, when you do calculations on these variables, like:</div><div class="gmail_default" style><br></div><div class="gmail_default" style><div class="gmail_default">
<font face="courier new, monospace">Tavg = dim_avg_n(tho,0)</font></div><div class="gmail_default"> </div><div class="gmail_default">
this causes all the metadata to be stripped off, and thus "Tavg" will have no metadata. This makes plotting "Tavg" over a map difficult.</div><div class="gmail_default"><br></div><div class="gmail_default">
Fortunately, we have versions of these functions that will reattach the metadata for you:</div><div class="gmail_default"><br></div><div class="gmail_default"><div class="gmail_default"><font face="courier new, monospace">Tavg = dim_avg_n_Wrap(tho,0) ; Use "_Wrap" to maintain metadata</font></div>
<div class="gmail_default"><br></div><div class="gmail_default">Now, when you do a "printVarSummary(tavg)", you will see that "tavg" also has metadata:</div><div class="gmail_default"><br></div><div class="gmail_default">
<font face="courier new, monospace">Variable: Tavg</font></div><div class="gmail_default"><font face="courier new, monospace">Type: float</font></div><div class="gmail_default"><font face="courier new, monospace">Total Size: 409920 bytes</font></div>
<div class="gmail_default"><font face="courier new, monospace"> 102480 values</font></div><div class="gmail_default"><font face="courier new, monospace">Number of Dimensions: 3</font></div><div class="gmail_default">
<font face="courier new, monospace">Dimensions and sizes:<span class="" style="white-space:pre"> </span>[depth | 20] x [lat | 122] x [lon | 42]</font></div><div class="gmail_default"><font face="courier new, monospace">Coordinates: </font></div>
<div class="gmail_default"><font face="courier new, monospace"> depth: [ 7..2435]</font></div><div class="gmail_default"><font face="courier new, monospace"> lat: [45.375..-45.375]</font></div><div class="gmail_default">
<font face="courier new, monospace"> lon: [-51.375..-20.625]</font></div><div class="gmail_default"><font face="courier new, monospace">Number Of Attributes: 7</font></div><div class="gmail_default"><font face="courier new, monospace"> _FillValue :<span class="" style="white-space:pre"> </span>-9e+33</font></div>
<div class="gmail_default"><font face="courier new, monospace"> standard_name :<span class="" style="white-space:pre"> </span>sea_water_potential_temperature</font></div><div class="gmail_default"><font face="courier new, monospace"> long_name :<span class="" style="white-space:pre"> </span>Sea water potential temperature</font></div>
<div class="gmail_default"><font face="courier new, monospace"> units :<span class="" style="white-space:pre"> </span>C</font></div><div class="gmail_default"><font face="courier new, monospace"> code :<span class="" style="white-space:pre"> </span>2</font></div>
<div class="gmail_default"><font face="courier new, monospace"> missing_value :<span class="" style="white-space:pre"> </span>-9e+33</font></div><div class="gmail_default"><font face="courier new, monospace"> average_op_ncl :<span class="" style="white-space:pre"> </span>dim_avg_n over dimension(s): time </font> </div>
<div>One other issue is that your "pot" variable has coordinate arrays, but the "lon" values go from 308 to 342. I had to subtract 360 from your longitude values so that it would plot correctly.</div><div>
<br></div><div>I cleaned up your script quite a bit, but I'm not sure the plot is correct. I will send you the script and image offline, since I'm not sure you want to share it with ncl-talk.</div><div><br></div>
<div>--Mary</div></div></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Jul 29, 2014 at 9:18 AM, Hemming, Michael <span dir="ltr"><<a href="mailto:michael.hemming@mpimet.mpg.de" target="_blank">michael.hemming@mpimet.mpg.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt"><br>
<div style="font-family:Times New Roman;color:#000000;font-size:16px">
<div></div>
Here is my code, I have attached some data, thanks for your help! This is only 1 year of the data as the file size was too large.<br>
<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt"><br>
<br>
<span style="background-color:rgb(255,255,0)">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
<br>
data12 = addfile ("<a href="http://sector_icon_lr_12_mpiom_tho_ym_0000-0099.nc" target="_blank">sector_icon_lr_12_mpiom_tho_ym_0000-0099.nc</a>", "r")<br>
tho = data12->tho<br>
noshupw = addfile ("<a href="http://con.ym.atlbox.nosh.upw.100y.nc" target="_blank">con.ym.atlbox.nosh.upw.100y.nc</a>", "r")<br>
potT = noshupw->t_acc<br>
<br>
<br>
;----------Surface plot MPIOM<br>
<br>
<br>
wks = gsn_open_wks("pdf", "T_Surface_MPIOM")<br>
gsn_define_colormap(wks,"gui_default") <br>
<br>
res = True<br>
res@tiMainString = "Surface T averaged for years 01-99"<br>
res@gsnSpreadColors = True <br>
res@pmLegendDisplayMode = False<br>
res@cnFillOn = True <br>
<br>
<br>
<br>
Tavg = dim_avg_n(tho,0) <br>
<br>
;plot = gsn_csm_contour(wks,Tavg(1,:,:),res)<br>
<br>
;----------Surface plot ICON<br>
<br>
; Does not work when including all the other plot code....try in new file and works fine (e.g. tester, copy data and plot code parts)<br>
<br>
wks = gsn_open_wks("x11", "T_Surface_ICON")<br>
gsn_define_colormap(wks,"gui_default") <br>
<br>
res = True<br>
res@tiMainString = "Surface T averaged for years 01-99 (ICON)"<br>
res@gsnSpreadColors = True <br>
res@pmLegendDisplayMode = False<br>
res@cnFillOn = True <br>
<br>
<br>
<br>
potTavg = dim_avg_n(potT,0) <br>
<br>
;plot = gsn_csm_contour(wks,potTavg(1,:,:),res)<div class=""><br>
<br>
<br>
;----------MPIOM vs. ICON comparison plots<br>
<br>
<br>
<br></div>
wks = gsn_open_wks("pdf","MPIOM_vs_ICON_panelplot_surface")<div class=""><br>
gsn_define_colormap(wks,"gui_default") <br>
<br>
<br>
plot1 = new(2,graphic) <br>
<br>
res@gsnDraw = False ; don't draw<br>
res@gsnFrame = False ; don't advance frame<br>
resP = True<br>
res@tiYAxisString = "Latitude"<br>
res@tiXAxisString = "Longitude"<br>
resP@txString = "Surface T averaged for years 01-99"<br>
res@lbOrientation = "vertical" ; vertical label bar<br>
<br>
res@trXMinF = 0 ; set minimum X-axis value<br>
res@trXMaxF = 40 ; set maximum X-axis value <br>
res@trYMinF = 0 ; set minimum Y-axis value<br>
res@trYMaxF = 120 ; set maximum Y-axis value<br>
res@tiMainString = "MPIOM exp 12"<br>
<br>
<br></div>
plot1(0) = gsn_csm_contour_map(wks,Tavg(1,:,:),res)<div class=""><br>
<br>
<br>
res@tiMainString = "ICON NOSH.UPW"<br>
<br>
;res@tmXBMode = "Explicit" <br>
;res@tmXBValues = (/-44.5,-27.5,-5.5,14.5,34.5/)<br>
;res@tmXBLabels = "" + res@tmXBValues <br>
;res@tmXBMinorValues = ispan(0,30,1)<br>
<br>
<br>
<br>
<br></div>
plot1(1) = gsn_csm_contour_map(wks,potTavg(1,:,:),res)</span><br>
<br>
<br>
<span style="background-color:rgb(255,255,0)">gsn_panel(wks,plot1,(/1,2/),resP) </span>
<br>
<div style="font-family:Times New Roman;color:#000000;font-size:16px">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>From:</b> Mary Haley [<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>]<br>
<b>Sent:</b> 29 July 2014 05:46<br>
<b>To:</b> Hemming, Michael<br>
<b>Cc:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<b>Subject:</b> Re: [ncl-talk] changing tick marks<br>
</font><br>
</div><div><div class="h5">
<div></div>
<div>
<div dir="ltr">
<div class="gmail_default" style="font-size:small">Hi Michael,</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">It looks like you have a curvilinear grid. The information I gave you was assuming you had a rectilinear grid, which means your data contains 1D lat/lon coordinate arrays attached to it.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">With a curvilinear grid, you have to read the 2D lat/lon coordinate values off the file separately and attach them to your data.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">Can you do a "printVarSummary" on the two variables you are trying to plot, so I can see what kind of metadata they have:</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small"><span style="color:rgb(255,0,0);font-family:Tahoma;font-size:13px;background-color:rgb(255,255,153)">printVarSummary(potTavg)</span><br>
</div>
<div class="gmail_default" style="font-size:small"><span style="color:rgb(255,0,0);font-family:Tahoma;font-size:13px;background-color:rgb(255,255,153)">printVarsummary(Tavg)</span></div>
<div class="gmail_default" style="font-size:small"><span style="font-family:Tahoma;font-size:13px;background-color:rgb(255,255,153)"><font color="#000000"><br>
</font></span></div>
<div class="gmail_default" style="font-size:small"><span style="color:rgb(0,0,0);font-family:Tahoma">This will help me determine what kind of grid you have, and hopefully what the lat/lon arrays associated with it are.</span></div>
<div class="gmail_default" style="font-size:small"><span style="color:rgb(0,0,0);font-family:Tahoma"><br>
</span></div>
<div class="gmail_default" style="font-size:small"><span style="color:rgb(0,0,0);font-family:Tahoma">Better yet, if you can provide me with your full script and data (you can do this offline), then I can probably quickly provide you with an NCL script that
will plot it.</span><br>
</div>
<div class="gmail_default" style="font-size:small"><span style="font-family:Tahoma;font-size:13px;background-color:rgb(255,255,255)"><font color="#000000"><br>
</font></span></div>
<div class="gmail_default" style="font-size:small"><span style="font-family:Tahoma;font-size:13px;background-color:rgb(255,255,255)"><font color="#000000">Thanks,</font></span></div>
<div class="gmail_default" style="font-size:small"><span style="font-family:Tahoma;font-size:13px;background-color:rgb(255,255,255)"><font color="#000000"><br>
</font></span></div>
<div class="gmail_default"><font style="background-color:rgb(255,255,255)" color="#000000" face="Tahoma">--Mary</font></div>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Mon, Jul 28, 2014 at 1:14 AM, Michael Hemming <span dir="ltr">
<<a href="mailto:michael.hemming@mpimet.mpg.de" target="_blank">michael.hemming@mpimet.mpg.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF">Hi Mary,<br>
<br>
Thanks for replying but it doesn't seem to work, it complains:<br>
<br>
'check_for_y_lat_coord: Warning: data does not contain a valid latitude coordinate array or doesn't contain one at all..'<br>
<br>
and the same for the lon_coord<br>
<br>
I'm confused as Ferret plots the data with the correct lat and lon.<br>
<br>
Is there something else I am missing? the only thing I changed within the code was the removal of the two lines you stated below and the addition of '_map' to the gsn command.<br>
<br>
thanks,<br>
Michael
<div>
<div><br>
<br>
<div>On 27/07/2014 19:48, Mary Haley wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div class="gmail_default" style="font-size:small">Hi Michael,</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">Part of the problem is that you've set the minimum of your X axis (the longitude axis to 0):</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small"><span style="color:rgb(255,0,0);font-family:Tahoma;font-size:13px;background-color:rgb(255,255,153)">res@trXMinF = 0 ; set minimum X-axis value</span><br style="color:rgb(255,0,0);font-family:Tahoma;font-size:13px;background-color:rgb(255,255,153)">
<span style="color:rgb(255,0,0);font-family:Tahoma;font-size:13px;background-color:rgb(255,255,153)">res@trXMaxF = 30 ; set maximum X-axis value </span><br>
</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">This is causing your longitudes to only go from 0 to 30, and hence you are not going to see any labels or tickmarks at negative longitude values.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">I suggest not setting the tr*M**F resources at all, unless you need to zoom in on the plot for some reason.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">Also, if you are indeed plotting lat/lon data, you should be calling "gsn_csm_contour_map" and not "gsn_csm_contour". The former will put draw contours over a map, whereas the latter only draws a contour plot.
For examples of using gsn_csm_contour_map, see:</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/cylineq.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/cylineq.shtml</a><br>
</div>
<div class="gmail_default"><br>
</div>
<div class="gmail_default">Note that some of these examples are using "gsn_csm_contour_map_ce", which is the same thing. gsn_csm_contour_map will draw a cylindrical equidistant map plot by default, unless you change the map projection.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">--Mary</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Sat, Jul 26, 2014 at 7:44 AM, Hemming, Michael <span dir="ltr">
<<a href="mailto:michael.hemming@mpimet.mpg.de" target="_blank">michael.hemming@mpimet.mpg.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Hey everyone,<br>
<br>
I am having trouble understanding how to alter tickmarks on a contour plot. I want to change my index numbers on the x and y axis to lat and lon values. I have been trying to use the example given here:
<a href="http://www.ncl.ucar.edu/Applications/Scripts/tm_2.ncl" target="_blank">http://www.ncl.ucar.edu/Applications/Scripts/tm_2.ncl</a>
<br>
<br>
Here is my code ( I have only tried altering the second plot out of the 2 e.g. for plot1(1)..):<br>
<br>
;----------MPIOM vs. ICON comparison plots<br>
<br>
<br>
<span style="background-color:rgb(255,255,153)"><br>
<font color="FF0000">wks = gsn_open_wks("x11","MPIOM_vs_ICON_panelplot_surface")<br>
gsn_define_colormap(wks,"gui_default") <br>
<br>
<br>
plot1 = new(2,graphic) <br>
<br>
res@gsnDraw = False ; don't draw<br>
res@gsnFrame = False ; don't advance frame<br>
resP = True<br>
res@tiYAxisString = "Latitude"<br>
res@tiXAxisString = "Longitude"<br>
resP@txString = "Surface T averaged for years 01-99"<br>
res@lbOrientation = "vertical" ; vertical label bar<br>
<br>
res@trXMinF = 0 ; set minimum X-axis value<br>
res@trXMaxF = 40 ; set maximum X-axis value <br>
res@trYMinF = 0 ; set minimum Y-axis value<br>
res@trYMaxF = 120 ; set maximum Y-axis value<br>
res@tiMainString = "MPIOM exp 12"<br>
<br>
<br>
plot1(0) = gsn_csm_contour(wks,Tavg(1,:,:),res)<br>
<br>
<br>
res@trXMinF = 0 ; set minimum X-axis value<br>
res@trXMaxF = 30 ; set maximum X-axis value <br>
res@trYMinF = 0 ; set minimum Y-axis value<br>
res@trYMaxF = 80 ; set maximum Y-axis value<br>
res@tiMainString = "ICON NOSH.UPW"<br>
<br>
res@tmXBMode = "Explicit" <br>
res@tmXBValues = (/-44.5,-27.5,-5.5,14.5,34.5/)<br>
res@tmXBLabels = "" + res@tmXBValues <br>
res@tmXBMinorValues = ispan(0,30,1)<br>
<br>
<br>
<br>
<br>
plot1(1) = gsn_csm_contour(wks,potTavg(1,:,:),res)<br>
<br>
<br>
gsn_panel(wks,plot1,(/1,2/),resP) </font></span><br>
<br>
<br>
<span style="background-color:rgb(255,255,255)"><font color="000000">I want l<font color="000000">ongitude from -44.5
<font color="000000">E</font> to 34.5 <font color="000000">W, once I know how to do it for longitude, I can then apply the code to the lat.<br>
<br>
<font color="000000">thanks<font color="000000"><font color="000000"> in advance!</font><span><font color="#888888"><br>
<br>
<font color="000000">Micha<font color="000000">el</font></font><br>
</font></span></font></font></font></font></font></span></div>
</div>
<br>
_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</blockquote>
<br>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div></div></div>
</div>
</div>
</div>
</div>
</div>
</blockquote></div><br></div>