[ncl-talk] Eady _growth_rate

George Vandenberghe - NOAA Affiliate george.vandenberghe at noaa.gov
Fri Aug 17 09:33:11 MDT 2018


There was a tangential thread about the relevance of Eady Growth Rates in a
tropical cyclone environment.   The function can give you
a number but is it realistic with the approximations made for
quasi-geostrophic baroclinic theory?

Apart from that though the problem is very relevant when considering a
tropical cyclone as a source of a perturbation
in a baroclinic zone and that situation is common to general with recurving
tropical cyclones.


On Fri, Aug 17, 2018 at 10:28 AM, Dennis Shea <shea at ucar.edu> wrote:

>
>    printVarSummary(T41)          ; ; [Time|12]x[ilev |* 26*] x [lat | 96]
> x [lon | 144]
>    printVarSummary(T411)         ; ; [Time|12]x[ilev | *26*] x [lat | 96]
> x [lon | 144]
>    printVarSummary(T412)         ; ; [Time|12]x[ilev | *27*] x [lat | 96]
> x [lon | 144]
>
> Variable T412 is the issue. The 'ilev' dimension is different.
>
> You must interpolate it to the same dimensions as T41, T411
>
> pressure: int2p_n
> <http://www.ncl.ucar.edu/Document/Functions/Built-in/int2p_n.shtml> or,
> better, *int2p_n_Wrap
> <http://www.ncl.ucar.edu/Document/Functions/Contributed/int2p_n_Wrap.shtml>*
> ====
>
> hybrid vertical coordinates: vinth2p
> <http://www.ncl.ucar.edu/Document/Functions/Built-in/vinth2p.shtml>
>
> On Fri, Aug 17, 2018 at 8:07 AM, Sri Nandini <snandini at marum.de> wrote:
>
>> ​Hello
>>
>> I am trying to plot the eady growth rate as a measure of baroclinic
>> instability. I cannot due to:
>> fatal:conform: The array to be conformed must have the same number of
>> dimensions as indicated by the length of the last argument
>> I cannot conform the latitude to same dimensions as the th variable.
>>
>> Below is my script, would be grateful for any suggestions.
>>
>>
>> load "lat_lon_masked_lambert.ncl"
>>
>> ; ==============================================================
>> ; Open the file: Read only the user specified period first observations
>> then model
>> ; ==============================================================
>>    f= addfile("SLP_PHIS_Z3_PI.nc", "r") ;
>>    T41    = f->Z3(:,:,:,:)
>>
>>    printVarSummary(T41)                                 ; ;
>> [Time|12]x[ilev | 26] x [lat | 96] x [lon | 144]
>>    T41 at _FillValue = -9.96921e+36
>>
>> ;==============================================================
>>
>>    f1= addfile("totalwinds_PI.nc", "r") ;
>>    T411    = f1->U(:,:,:,:)
>>
>>    printVarSummary(T411)                                 ; ;
>> [Time|12]x[ilev | 26] x [lat | 96] x [lon | 144]
>>    T411 at _FillValue = -9.96921e+36
>>
>>
>> ;==============================================================
>>    f2= addfile("th_PI.nc", "r") ;
>>    T412    = f2->TH(:,:,:,:)
>>
>>    printVarSummary(T412)                                 ; ;
>> [Time|12]x[ilev | 27] x [lat | 96] x [lon | 144]
>>    T412 at _FillValue = -9.96921e+36
>>
>> ;==============================================================
>>
>> ;    Read latitudes
>> ;    The 'eady_growth_rate' function requires that 'lat' and 'th' agree
>> ;    Use 'conform' the propogate the lat values
>> ;==============================================================
>>
>>    xlat = f->lat                               ;  [lat | 96]
>>    printVarSummary(xlat)
>>
>>    XLAT = conform(T412, xlat, (/0,2,3/))     ; problem here
>>    printVarSummary(XLAT)
>>
>>    egr = eady_growth_rate(aveX2, aveX1, aveX, XLAT, 0,  1)
>>    printVarSummary(egr)
>>    printMinMax(egr, 0)
>>
>> ;==============================================================
>>    wks = gsn_open_wks("pdf","Eady")              ; send graphics to PNG
>> file
>>
>> ;---Set some basic plot options
>>
>>    res               = True
>>    res at gsnMaximize   = True       ; maximize plot in frame
>>    res at gsnAddCyclic  = False
>>
>>    res at cnFillOn      = True
>>    res at cnLinesOn     = False
>>   ;res at cnFillMode    = "RasterFill"                 ; slow here
>>    res at cnFillMode    = "CellFill"                   ; faster
>>
>>
>>   minlat = 25.                          ; min lat to mask
>>   maxlat = 80.                          ; max lat to mask
>>   minlon = -10.                          ; min lon to mask
>>   maxlon =  110.                          ; max lon to mask
>>   res at mpProjection = "LambertConformal"            ; choose projection
>>
>> ;---masked plot
>>   res at gsnAddCyclic = True                ; regional plot
>>
>>   res at mpMinLatF = minlat              ; min lat to mask
>>   res at mpMaxLatF = maxlat              ; max lat to mask
>>   res at mpMinLonF = minlon              ; min lon to mask
>>   res at mpMaxLonF = maxlon              ; max lon to mask
>>
>>   res at gsnMaskLambertConformal = True                ; turn on lc masking
>>
>> ; specify a level or levels ... within boundary layer
>>
>>    klStrt = 5
>>    klLast = 5
>>    nt     = 0
>>
>> ;--- Eady growth rate (1/day)
>>
>>    egr        = egr*86400
>>    egr at units  = "1/day"
>>
>>    res at cnFillPalette        = "precip2_17lev"
>>    res at cnLevelSelectionMode = "ManualLevels"     ; set manual contour
>> levels
>>    res at cnMinLevelValF       =  0.5               ; set min contour level
>>    res at cnMaxLevelValF       =  4.0               ; set max contour level
>>    res at cnLevelSpacingF      =  0.25              ; set contour spacing
>>
>>    do kl=klStrt,klLast
>>       res at gsnCenterString = "znu="+znu(nt,kl)
>>       res at gsnRightString  =  egr at units
>>       contour = gsn_csm_contour_map(wks,  egr(nt,kl,:,:),res)
>>    end do
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180817/0df237b6/attachment.html>


More information about the ncl-talk mailing list