<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Ryan,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Is there a reason why you don&#39;t want to upgrade from V6.1.2? The version you have is over 3 years old.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I can&#39;t be sure of why your temp_adv variable is not plotting, without actually being able to run the script.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">However, rather than printing the whole temp_adv array, what does the following report:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">printVarSummary(temp_adv)</div><div class="gmail_default" style="font-size:small">printMinMax(temp_adv,0)</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Sometimes this will give you a clue where the problem is, for example, your min/max values are off scale, or the size of the array is not what you were expecting.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You can also try plotting temp_adv with a basic contour plot call:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">    plot = gsn_csm_contour(wks,temp_adv,False)</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">just to make sure that the problem isn&#39;t with wrf_contour itself.</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><div class="gmail_quote">On Mon, Apr 18, 2016 at 12:49 PM, Ryan Connelly <span dir="ltr">&lt;<a href="mailto:rconne01@gmail.com" target="_blank">rconne01@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>Running 6.1.2, so can&#39;t use advect_variable without upgrading.  So instead I have...</div><div><br></div><div><div><font face="monospace, monospace">;do it = 12,30,3             ; TIME LOOP</font></div><div><font face="monospace, monospace">  do it = 12,13,1</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">    print(&quot;Working on time: &quot; + times(it) )</font></div><div><font face="monospace, monospace">    res@TimeLabel = times(it)   ; Set Valid time to use on plots</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</font></div><div><font face="monospace, monospace">; First get the variables we will need</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">    tc = wrf_user_getvar(a,&quot;tc&quot;,it)      ; 3D tc</font></div><div><font face="monospace, monospace">    td = wrf_user_getvar(a,&quot;td&quot;,it)      ; 3D td</font></div><div><font face="monospace, monospace">    u  = wrf_user_getvar(a,&quot;ua&quot;,it)      ; 3D U at mass points</font></div><div><font face="monospace, monospace">    v  = wrf_user_getvar(a,&quot;va&quot;,it)      ; 3D V at mass points</font></div><div><font face="monospace, monospace">    p  = wrf_user_getvar(a, &quot;pressure&quot;,it) ; pressure is our vertical coordinate</font></div><div><font face="monospace, monospace">    z  = wrf_user_getvar(a,&quot;z&quot;,it)       ; Full model height in meters</font></div><div><font face="monospace, monospace">       z_dec = z/10.                     ; Height in decameters</font></div><div><font face="monospace, monospace">    </font><span style="font-family:monospace,monospace">lon = wrf_user_getvar(a,&quot;lon&quot;,it)</span></div><div><font face="monospace, monospace">    lat = wrf_user_getvar(a,&quot;lat&quot;,it)</font></div></div><div><font face="monospace, monospace"><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div>; Interpolate to isobaric planes</div><div><br></div><div>      pressure = 850.</div><div><br></div><div>      tc_plane = wrf_user_intrp3d(tc,p,&quot;h&quot;,pressure,0.,False)</div><div>      z_plane  = wrf_user_intrp3d(z_dec,p,&quot;h&quot;,pressure,0.,False)</div><div>      u_plane  = wrf_user_intrp3d( u,p,&quot;h&quot;,pressure,0.,False)</div><div>      v_plane  = wrf_user_intrp3d( v,p,&quot;h&quot;,pressure,0.,False)</div><div><br></div><div><br></div><div>; Define variables needed to calculate temp advection<br></div><div><br></div><div>      dTdx = center_finite_diff_n(tc_plane,lon,False,0,0)</div><div>      dTdy = center_finite_diff_n(tc_plane,lat,False,0,0)</div><div><br></div><div>      temp_adv = u_plane*dTdx + v_plane*dTdy<br></div><div><br></div><div>      print(temp_adv)</div><div><br></div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>    ; Plotting options for Tc</div><div>      opts = res</div><div>      print(&quot;Defined ops&quot;)</div><div>      opts@cnFillOn = True</div><div>      print(&quot;Turned fill on&quot;)</div><div>      opts@cnLabelMasking = True</div><div>      opts@cnInfoLabelOn = False</div><div>      opts@cnLineLabelPerimOn = False</div><div>      opts@ContourParameters = (/ -60., 60., 3./)</div><div>      print(&quot;Accepted contour parameters&quot;)</div><div>      ;opts@gsnSpreadColorEnd = -3  ; End third from the last color in color map</div><div>      contour_tc = wrf_contour(a,wks,temp_adv,opts) ; &lt;- Breaks right here</div><div>      print(&quot;wrf_contour called&quot;)</div><div>      delete(opts)</div><div>      print(&quot;opts deleted&quot;)</div><div><br></div><div>      print(&quot;Got past plotting options for Tc&quot;)</div><div><br></div><div><br></div><div><div>; MAKE PLOTS</div><div>      ;plot = wrf_map_overlays(a,wks,(/contour_tc,contour_td,contour_z/),pltres,mpres)</div><div>      plot = wrf_map_overlays(a,wks,(/contour_tc/),pltres,mpres)</div><div>      ;plot = wrf_map_overlays(a,wks,(/contour_td,vector/),pltres,mpres)</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div>  end do        ; END OF TIME LOOP</div></div><div><br></div><div><br></div></font></div><div>When it prints out the full grid of temp_adv, or random indexes of the array, I get reasonable values and no reason to suspect the grid is undefined anywhere.  Yet as you can see from my print statements, the call to wrf_contour fails.  It does not throw an error; it just hangs up.  I checked with just tc_plane instead, and that completed fine, so it has to be an issue with the temp_adv grid I calculated.</div><div><br></div><div>Thanks,</div><div>Ryan<span class="HOEnZb"><font color="#888888"><br clear="all"><div><br></div>-- <br><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font size="1">Ryan Connelly</font></div><div dir="ltr"><div dir="ltr" style="font-size:12.8px"><font color="#000000"><font size="1">M.S. Student in Atmospheric Sciences, </font><span style="font-size:x-small">Stony Brook University</span><br></font></div><div dir="ltr" style="font-size:12.8px"><div><font color="#000000"><font size="1">B.S. in Meteorology with Minors in Mathematics and GIS, </font><span style="font-size:x-small">Valparaiso University</span></font></div></div><div><a href="mailto:rconne01@gmail.com" style="font-size:x-small" target="_blank">rconne01@gmail.com</a><br></div><div><font size="1"><a href="mailto:ryan.connelly@stonybrook.edu" target="_blank">ryan.connelly@stonybrook.edu</a></font></div></div></div></div></div></div></div></div></div></div>
</font></span></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>