[ncl-talk] Temp advection, wrf_contour fails
Ryan Connelly
rconne01 at gmail.com
Wed Apr 20 11:49:24 MDT 2016
Hi Dennis,
My file is now in ../incoming/wrf_temp_adv_rconne01.ncl.
It doesn't like dimension reduction because the resulting argument isn't
the right size:
fatal:center_finite_diff_n: r must either be a scalar, a 1D array the same
length as the dim-th dimemsion of q, or the same size as q
It also throws a fatal on the declaration of lon_plane and/or lat_plane:
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading dims
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 224 in
file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 97 in file
wrf_temp_adv.ncl
My best success so far has been with your first suggestion, merely changing
the last arg of the dTDx declaration from 0 to 1, which gave me realistic
values for dTdx instead of inf's:
(0) min=-15.1425 max=18.5633
(0) min=-24.1597 max=15.1911
But that ultimately still hangs up, though interestingly enough, it now
gets one step farther. Before, wrf_contour was where it stopped. Now, it
gets past wrf_contour, but fails on plot=.
Thanks,
Ryan
On Wed, Apr 20, 2016 at 1:19 PM, Ryan Connelly <rconne01 at gmail.com> wrote:
> Hi Mary,
>
> That shows that dTdx is the problem. (dTdx is first, dTdy second:)
>
> (0) min=-inf max=inf
> (0) min=-24.1597 max=15.1911
>
> I'll work on Dennis's suggestion now.
>
> Ryan
>
> On Wed, Apr 20, 2016 at 11:59 AM, Dennis Shea <shea at ucar.edu> wrote:
>
>> NOTE: I have not looked at this in detail.
>>
>> [1]
>> You could not use 'advect_variable' because it requires "[t]he array *must
>> be global and ordered south to north*." Also, this assumes a rectilinear
>> grid (eg: fixed or gaussian)
>>
>> Your data are not global :-(
>> Yourhorizontal grid is not rectilinear.
>>
>> Note: I could change ''advect_variable" which uses spherical harmonics
>> (hence, the global grid requirement) to derive the horizontal and
>> latitudinal gradients tusing the 6.4.0
>>
>> https://www.ncl.ucar.edu/Document/Functions/Contributed/grad_latlon_cfd.shtml
>> BUT that would take a day or two.
>>
>> [2]
>> center_finite_diff "the issue" .. yes and no
>>
>> The error occurs in center_finite_diff but the issue is, I speculate, the
>> arguments. Specifically, the use of lat and lon
>>
>> ====
>> lon = wrf_user_getvar(a,"lon",it)
>> lat = wrf_user_getvar(a,"lat",it)
>>
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>> ; Interpolate to isobaric planes
>>
>> pressure = 850.
>>
>> tc_plane = wrf_user_intrp3d(tc,p,"h",pressure,0.,False)
>> z_plane = wrf_user_intrp3d(z_dec,p,"h",pressure,0.,False)
>> u_plane = wrf_user_intrp3d( u,p,"h",pressure,0.,False)
>> v_plane = wrf_user_intrp3d( v,p,"h",pressure,0.,False)
>>
>>
>> ; Define variables needed to calculate temp advection
>>
>> dTdx = center_finite_diff_n(tc_plane,lon,False,0,0)
>> dTdy = center_finite_diff_n(tc_plane,lat,False,0,0)
>> ===
>>
>> Maybe,
>>
>>
>> dTdx = center_finite_diff_n(tc_plane,lon,False,0,1) ; <==
>> dTdy = center_finite_diff_n(tc_plane,lat,False,0,0)
>>
>> Or [281] x [353]=>nlat,mlon
>>
>> dTdx = center_finite_diff_n(tc_plane,lon(nlat/2,:),False,0,0)
>> dTdy = center_finite_diff_n(tc_plane,lat(:,mlon/2),False,0,0)
>>
>> Note that NCL's 'dimension reduction' (elimination of degenerate
>> dimensions), lon(nlat/2,:) and lat(:,mlon/2) are one dimensional.
>> Use printVarSummary( lon(nlat/2,:) ) to verify
>>
>> [3]
>>
>> [2] is a thumb in the dike approach. I think it would be best (if
>> possible) to calculate
>>
>> lon_plane =
>> lat_plane =
>>
>> dTdx = center_finite_diff_n(tc_plane,lon_plane,False,0,0)
>> dTdy = center_finite_diff_n(tc_plane,lat_plane,False,0,0)
>>
>> [4]
>> Maybe you could send a sample filw (one only) to
>>
>> ftp ftp.cgd.ucar.edu
>> anonymous
>> your_email
>> cd incoming
>> put .....
>> quit
>>
>> then let ncl_talk know the name of the file.
>>
>> On Tue, Apr 19, 2016 at 4:21 PM, Mary Haley <haley at ucar.edu> wrote:
>>
>>> Hi Ryan,
>>>
>>> That makes sense why you can't upgrade.
>>>
>>> As for the "inf", this is definitely an issue. To further determine
>>> where the problem is, I suggest doing a printMinMax on the individual
>>> variables in the calculation:
>>>
>>> printMinMax(dTdx,0)
>>> printMinMax(dTdy,0)
>>>
>>> etc.
>>>
>>> You may want to do this on every variable that you are getting with
>>> "wrf_user_getvar" to make sure those all look okay.
>>>
>>> If it looks like center_finite_diff is the issue, then revisit the
>>> documentation to make sure all the input parameters are being input
>>> correctly:
>>>
>>>
>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/center_finite_diff.shtml
>>>
>>> --Mary
>>>
>>> On Tue, Apr 19, 2016 at 3:53 PM, Ryan Connelly <rconne01 at gmail.com>
>>> wrote:
>>>
>>>> Hi Mary,
>>>>
>>>> I use NCL that's installed on a shared server where WRF is configured.
>>>> I don't suspect I have permissions to remove an old version. I recently
>>>> updated to 6.3.0 on my local machine, but then I have to copy the wrfout
>>>> over, which is a bit time-consuming since I'm running a triple nest down to
>>>> 1.333 km! :O
>>>>
>>>> These extra commands that you suggested certainly give a clue that I'm
>>>> not calculating things right:
>>>>
>>>> Variable: temp_adv
>>>> Type: float
>>>> Total Size: 396772 bytes
>>>> 99193 values
>>>> Number of Dimensions: 2
>>>> Dimensions and sizes: [281] x [353]
>>>> Coordinates:
>>>> Number Of Attributes: 1
>>>> _FillValue : 9.96921e+36
>>>> (0) min=-inf max=inf
>>>>
>>>>
>>>> 281x353 is my grid size, so that's good, but the infinities probably
>>>> are not...
>>>>
>>>> Thanks,
>>>> Ryan
>>>>
>>>> On Tue, Apr 19, 2016 at 10:38 AM, Mary Haley <haley at ucar.edu> wrote:
>>>>
>>>>> Hi Ryan,
>>>>>
>>>>> Is there a reason why you don't want to upgrade from V6.1.2? The
>>>>> version you have is over 3 years old.
>>>>>
>>>>> I can't be sure of why your temp_adv variable is not plotting, without
>>>>> actually being able to run the script.
>>>>>
>>>>> However, rather than printing the whole temp_adv array, what does the
>>>>> following report:
>>>>>
>>>>> printVarSummary(temp_adv)
>>>>> printMinMax(temp_adv,0)
>>>>>
>>>>> 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.
>>>>>
>>>>> You can also try plotting temp_adv with a basic contour plot call:
>>>>>
>>>>> plot = gsn_csm_contour(wks,temp_adv,False)
>>>>>
>>>>> just to make sure that the problem isn't with wrf_contour itself.
>>>>>
>>>>> --Mary
>>>>>
>>>>>
>>>>> On Mon, Apr 18, 2016 at 12:49 PM, Ryan Connelly <rconne01 at gmail.com>
>>>>> wrote:
>>>>>
>>>>>> Hi,
>>>>>>
>>>>>> Running 6.1.2, so can't use advect_variable without upgrading. So
>>>>>> instead I have...
>>>>>>
>>>>>> ;do it = 12,30,3 ; TIME LOOP
>>>>>> do it = 12,13,1
>>>>>>
>>>>>> print("Working on time: " + times(it) )
>>>>>> res at TimeLabel = times(it) ; Set Valid time to use on plots
>>>>>>
>>>>>>
>>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>>> ; First get the variables we will need
>>>>>>
>>>>>> tc = wrf_user_getvar(a,"tc",it) ; 3D tc
>>>>>> td = wrf_user_getvar(a,"td",it) ; 3D td
>>>>>> u = wrf_user_getvar(a,"ua",it) ; 3D U at mass points
>>>>>> v = wrf_user_getvar(a,"va",it) ; 3D V at mass points
>>>>>> p = wrf_user_getvar(a, "pressure",it) ; pressure is our vertical
>>>>>> coordinate
>>>>>> z = wrf_user_getvar(a,"z",it) ; Full model height in meters
>>>>>> z_dec = z/10. ; Height in decameters
>>>>>> lon = wrf_user_getvar(a,"lon",it)
>>>>>> lat = wrf_user_getvar(a,"lat",it)
>>>>>>
>>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>>> ; Interpolate to isobaric planes
>>>>>>
>>>>>> pressure = 850.
>>>>>>
>>>>>> tc_plane = wrf_user_intrp3d(tc,p,"h",pressure,0.,False)
>>>>>> z_plane = wrf_user_intrp3d(z_dec,p,"h",pressure,0.,False)
>>>>>> u_plane = wrf_user_intrp3d( u,p,"h",pressure,0.,False)
>>>>>> v_plane = wrf_user_intrp3d( v,p,"h",pressure,0.,False)
>>>>>>
>>>>>>
>>>>>> ; Define variables needed to calculate temp advection
>>>>>>
>>>>>> dTdx = center_finite_diff_n(tc_plane,lon,False,0,0)
>>>>>> dTdy = center_finite_diff_n(tc_plane,lat,False,0,0)
>>>>>>
>>>>>> temp_adv = u_plane*dTdx + v_plane*dTdy
>>>>>>
>>>>>> print(temp_adv)
>>>>>>
>>>>>>
>>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>>>
>>>>>> ; Plotting options for Tc
>>>>>> opts = res
>>>>>> print("Defined ops")
>>>>>> opts at cnFillOn = True
>>>>>> print("Turned fill on")
>>>>>> opts at cnLabelMasking = True
>>>>>> opts at cnInfoLabelOn = False
>>>>>> opts at cnLineLabelPerimOn = False
>>>>>> opts at ContourParameters = (/ -60., 60., 3./)
>>>>>> print("Accepted contour parameters")
>>>>>> ;opts at gsnSpreadColorEnd = -3 ; End third from the last color
>>>>>> in color map
>>>>>> contour_tc = wrf_contour(a,wks,temp_adv,opts) ; <- Breaks right
>>>>>> here
>>>>>> print("wrf_contour called")
>>>>>> delete(opts)
>>>>>> print("opts deleted")
>>>>>>
>>>>>> print("Got past plotting options for Tc")
>>>>>>
>>>>>>
>>>>>> ; MAKE PLOTS
>>>>>> ;plot =
>>>>>> wrf_map_overlays(a,wks,(/contour_tc,contour_td,contour_z/),pltres,mpres)
>>>>>> plot = wrf_map_overlays(a,wks,(/contour_tc/),pltres,mpres)
>>>>>> ;plot =
>>>>>> wrf_map_overlays(a,wks,(/contour_td,vector/),pltres,mpres)
>>>>>>
>>>>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>>>>>>
>>>>>> end do ; END OF TIME LOOP
>>>>>>
>>>>>>
>>>>>> 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.
>>>>>>
>>>>>> Thanks,
>>>>>> Ryan
>>>>>>
>>>>>> --
>>>>>> Ryan Connelly
>>>>>> M.S. Student in Atmospheric Sciences, Stony Brook University
>>>>>> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
>>>>>> University
>>>>>> rconne01 at gmail.com
>>>>>> ryan.connelly at stonybrook.edu
>>>>>>
>>>>>> _______________________________________________
>>>>>> ncl-talk mailing list
>>>>>> ncl-talk at ucar.edu
>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>
>>>>>>
>>>>>
>>>>
>>>>
>>>> --
>>>> Ryan Connelly
>>>> M.S. Student in Atmospheric Sciences, Stony Brook University
>>>> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
>>>> University
>>>> rconne01 at gmail.com
>>>> ryan.connelly at stonybrook.edu
>>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>
>
> --
> Ryan Connelly
> M.S. Student in Atmospheric Sciences, Stony Brook University
> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> University
> rconne01 at gmail.com
> ryan.connelly at stonybrook.edu
>
--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160420/0b315d77/attachment.html
More information about the ncl-talk
mailing list