<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hello,<br>

<br>

I am creating panel plots using gsn_csm_contour. The data I use has been chosen from a 4d matrix prior to plotting. The metadata is correct for longitude but not for the time axis. I would like to show the time in years, 1-100.

<br>

<br>

time metadata:<br>

<br>

tho = hours, 365_day<br>

potT = days, preleptic_gregorian<br>

<br>

<br>

Here is my code:<br>

<br>

<br>

<span style="background-color: rgb(255, 255, 153);">load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>

load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br>

load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<br>

<br>

begin<br>

<br>

;##########################################################################<br>

;##########################################################################<br>

;-------------------------Temperature Profile------------------------------<br>

;--------------------------------------------------------------------------<br>

;##########################################################################<br>

;##########################################################################<br>

<br>

;----------Load in data<br>

<br>

;MPIOM<br>

<br>

data12 &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; = addfile (&quot;sector.nc&quot;, &quot;r&quot;)<br>

tho &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; = data12-&gt;tho<br>

<br>

<br>

variables &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; = getfilevarnames(data12)<br>

<br>

print(&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;&nbsp;&nbsp;&nbsp; &quot;)<br>

print(&quot; MPIOM temperature data Information:&quot;)<br>

print(&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;&nbsp;&nbsp;&nbsp; &quot;)<br>

print(variables)<br>

printVarSummary(tho)<br>

<br>

;ICON<br>

<br>

noshupw&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; = addfile (&quot;con.nc&quot;, &quot;r&quot;)<br>

potT&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; = noshupw-&gt;t_acc<br>

potT&amp;lon&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = potT&amp;lon-360&nbsp;&nbsp;&nbsp;&nbsp; ; Fix the longitudes to be &lt; 360<br>

<br>

<br>

;variables &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; = getfilevarnames(noshupw)<br>

print(&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;&nbsp;&nbsp;&nbsp; &quot;)<br>

print(&quot;ICON temperature data information:&quot;)<br>

print(&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;&nbsp;&nbsp;&nbsp; &quot;)<br>

print(variables)<br>

printVarSummary(potT)<br>

<br>

<br>

;--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%<br>

;--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%<br>

<br>

;--------------Hovmuller T comparison plots<br>

<br>

;--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%<br>

;--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%--%<br>

<br>

;http://www.ncl.ucar.edu/Applications/Scripts/lat_time_3.ncl<br>

<br>

;*********************************************<br>

; manipulate data for plotting<br>

;*********************************************<br>

<br>

dm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = rm_single_dims(tho({depth|7},{lon|-30:-30},lat|:,time|:))&nbsp;&nbsp;&nbsp; <br>

di&nbsp;&nbsp;&nbsp; = rm_single_dims(potT({depth|7},{lon|-30:-30},lat|:,time|:))&nbsp;&nbsp; <br>

<br>

dm30S&nbsp;&nbsp; = rm_single_dims(tho({depth|7},lon|:,{lat|-30:-30},time|:))&nbsp;&nbsp;&nbsp; <br>

di30S&nbsp;&nbsp;&nbsp; = rm_single_dims(potT({depth|7},lon|:,{lat|-30:-30},time|:))&nbsp;&nbsp; <br>

<br>

dmeq&nbsp;&nbsp; = rm_single_dims(tho({depth|7},lon|:,{lat|0:0},time|:))&nbsp;&nbsp;&nbsp; <br>

dieq&nbsp;&nbsp;&nbsp; = rm_single_dims(potT({depth|7},lon|:,{lat|0:0},time|:))&nbsp;&nbsp; <br>

<br>

dm30N&nbsp;&nbsp; = rm_single_dims(tho({depth|7},lon|:,{lat|30:30},time|:))&nbsp;&nbsp;&nbsp; <br>

di30N&nbsp;&nbsp;&nbsp; = rm_single_dims(potT({depth|7},lon|:,{lat|30:30},time|:))&nbsp;&nbsp; <br>

<br>

<br>

printVarSummary(dm)<br>

printVarSummary(di)<br>

printVarSummary(dm30S)<br>

printVarSummary(di30S)<br>

<br>

;------Plotting<br>

<br>

wks = gsn_open_wks(&quot;x11&quot;, &quot;T)<br>

<br>

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<br>

res@cnFillOn&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; True<br>

<br>

plot&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; new(2,graphic)<br>

<br>

res@gsnDraw&nbsp;&nbsp;&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;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; don't draw<br>

res@gsnFrame&nbsp;&nbsp; &nbsp;&nbsp;&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;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; don't advance frame<br>

<br>

resP&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; = True&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; ; modify the panel plot<br>

resP@gsnMaximize&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp; <br>

res@tiMainFontHeightF &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.02<br>

resP@txFontHeightF&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.015<br>

resP@txString&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;Temperature at the surface with time at a longitude of 30W&quot;<br>

<br>

res@tiMainString&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; &quot;M&quot;&nbsp; <br>

plot(0) = gsn_csm_contour(wks,dm,res)<br>

<br>

res@tiMainString&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; &quot;I&quot;&nbsp;&nbsp;&nbsp;&nbsp; <br>

plot(1) = gsn_csm_contour(wks,di,res)<br>

<br>

;gsn_panel(wks,plot,(/1,2/),resP) <br>

<br>

;----------------------------surface plot<br>

<br>

<br>

wks = gsn_open_wks(&quot;x11&quot;, &quot;T_HOVMULLER_MPIOM_ICON_surface_30N&quot;)<br>

<br>

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<br>

res@cnFillOn&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; True<br>

<br>

plot1&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; new(2,graphic)<br>

<br>

res@gsnDraw&nbsp;&nbsp;&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;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; don't draw<br>

res@gsnFrame&nbsp;&nbsp; &nbsp;&nbsp;&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;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; don't advance frame<br>

<br>

<br>

resP@gsnMaximize&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp; <br>

res@tiMainFontHeightF &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.02<br>

resP@txFontHeightF&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.015<br>

resP@txString&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;Temperature at the surface with time at a latitude of 30N&quot;<br>

res@lbLabelBarOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False <br>

res@cnLevelSelectionMode =&nbsp; &quot;ManualLevels&quot;&nbsp;&nbsp; <br>

res@cnMinLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 3.<br>

res@cnMaxLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 7.<br>

res@cnLevelSpacingF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; 0.2 <br>

res@cnSpanFillPalette&nbsp;&nbsp;&nbsp; &nbsp;=&nbsp;&nbsp; True<br>

resP@gsnPanelLabelBar&nbsp;&nbsp;&nbsp; = True&nbsp; <br>

<br>

res@tiMainString&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; &quot;M&nbsp; <br>

plot1(0) = gsn_csm_contour(wks,dm30N,res)<br>

<br>

res@tiMainString&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; &quot;I&nbsp;&nbsp;&nbsp;&nbsp; <br>

plot1(1) = gsn_csm_contour(wks,di30N,res)<br>

<br>

gsn_panel(wks,plot1,(/1,2/),resP) <br>

<br>

<br>

<br>

<br>

<span style="background-color: rgb(255, 255, 255);">I have been told&nbsp; I can use cd_calendar but this involves creating a new variable 'time' which I would like to avoid. Another option was to change the tick labels but I am not sure how to do this using gsn_csm_contour.

<br>

<br>

thanks,<br>

Michael<br>

</span></span><br>

</div>

</body>

</html>