<div dir="ltr">Hi,<div><br></div><div>Running 6.1.2, so can'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("Working on time: " + 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,"tc",it) ; 3D tc</font></div><div><font face="monospace, monospace"> td = wrf_user_getvar(a,"td",it) ; 3D td</font></div><div><font face="monospace, monospace"> u = wrf_user_getvar(a,"ua",it) ; 3D U at mass points</font></div><div><font face="monospace, monospace"> v = wrf_user_getvar(a,"va",it) ; 3D V at mass points</font></div><div><font face="monospace, monospace"> p = wrf_user_getvar(a, "pressure",it) ; pressure is our vertical coordinate</font></div><div><font face="monospace, monospace"> z = wrf_user_getvar(a,"z",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,"lon",it)</span></div><div><font face="monospace, monospace"> lat = wrf_user_getvar(a,"lat",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,"h",pressure,0.,False)</div><div> z_plane = wrf_user_intrp3d(z_dec,p,"h",pressure,0.,False)</div><div> u_plane = wrf_user_intrp3d( u,p,"h",pressure,0.,False)</div><div> v_plane = wrf_user_intrp3d( v,p,"h",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("Defined ops")</div><div> opts@cnFillOn = True</div><div> print("Turned fill on")</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("Accepted contour parameters")</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) ; <- Breaks right here</div><div> print("wrf_contour called")</div><div> delete(opts)</div><div> print("opts deleted")</div><div><br></div><div> print("Got past plotting options for Tc")</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<br clear="all"><div><br></div>-- <br><div class="gmail_signature"><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>
</div></div>