[ncl-talk] Stability profile
Soma Roy
somaroy892 at gmail.com
Mon Jun 25 08:37:26 MDT 2018
Thanks a lot for your valuable opinions.
It worked!
On Mon, Jun 25, 2018, 19:12 Dennis Shea <shea at ucar.edu> wrote:
> [1] s(3,:,:,{70})
>
> Presumably, s(time,level,lat,lon)
>
> *3* => fourth time step [NCL indices are 0-based]
>
> *:,:* => all levels, all latitudes
>
> *{*70*} *=> 70 longitude
> * {...} *please read about coordinate subscripting
>
>
> ===
> [2]
> re: Plot=*gsn_csm_pres_hgt_vector*(wks,s(3,:,:,{70}),u(3,:,:,{70}),v
> (3,:,:,{70}),res)
> Is this averaging vertical column at 70E?
>
> No. This is a plot function, it is plotting level x latitude. You
> can do:
>
> *printVarSummary*( u(3,:,:,{70}) )
>
> Also, see:
>
> http://www.ncl.ucar.edu/Applications/height_lat.shtml#ex7
>
> ===
> [3]
> re: Secondly, if I want to plot 4 time steps in a panel (1x4) then how to
> change the above lines?
>
>
> u=*short2flt*(in->t(15:18,:,::-1,20:80))
> *printVarSummary*(t) ; what do you see?
>
> dimu = *dimsizes*(u)
> ntim = dimu(0)
>
> Generate one plot for each time step
>
> There are many panel examples at:
> http://www.ncl.ucar.edu/Applications/panel.shtml
> Please examine this carefully.
>
> plot = *new*(ntim,"graphic")
> do *nt*=0,ntim-1
> plot(nt) = *gsn_csm_pres_hgt_vector*(wks,s(*nt*,:,:,{70}),u(
> *nt*,:,:,{70}),v(*nt*,:,:,{70}),res)
> end do
> ===============================
>
> Please spend some time learning NCL: A good source is:
> http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/
>
>
>
> On Sun, Jun 24, 2018 at 1:33 AM, Soma Roy <somaroy892 at gmail.com> wrote:
>
>> Thank you!
>>
>> It worked. I want to know two things in this case.
>>
>> What this line suggests;
>> Plot=gsn_csm_pres_hgt_vector(wks,s(3,:,:,{70}),v(3,:,:,{70}),u
>> (3,:,:,{70}),res)
>>
>> And if it is like below then what it says?
>>
>> Plot=gsn_csm_pres_hgt_vector(wks,s(3,:,:,{70}),u(3,:,:,{70}),v
>> (3,:,:,{70}),res)
>>
>> Is this averaging vertical column at 70E?
>>
>> Secondly, if I want to plot 4 time steps in a panel (1x4) then how to
>> change the above lines?
>>
>> For that case my input is like below,
>>
>> t=short2flt(in->t(15:18,:,::-1,20:80))
>> u=short2flt(in->u(15:18,:,::-1,20:80))
>> v=short2flt(in->v(15:18,:,::-1,20:80))
>>
>> Thanks,
>> Soma
>>
>> On Sun, Jun 24, 2018, 02:04 Dennis Shea <shea at ucar.edu> wrote:
>>
>>> The following indicates that you do not have enough memory on yor
>>> machine.
>>>
>>> fatal:*NclMalloc *Failed:[errno=12]
>>> Segmentation fault (core dumped)
>>> ============================
>>> Total Size: *209943360 bytes*
>>> 52485840 values
>>> Number of Dimensions: 4
>>> Dimensions and sizes: [time | 120] x [level | 22] x [latitude | 141] x
>>> [longitude | 141]
>>>
>>> Each variable is ~210MB
>>> ===========================
>>> W
>>> For testing the code, try with a smaller time sample
>>>
>>> tStrt = 0
>>> tLast = 9
>>>
>>> ;---Read needed variables from file
>>> t = short2flt(in->t(tStrt:tLast,:,:,:)) ; select variable
>>> v = short2flt(in->v(tStrt:tLast,:,:,:))
>>> u = short2flt(in->u(tStrt:tLast,:,:,:)))
>>>
>>> ;---Better and more efficient to reorder *after *input
>>>
>>> t = t(:,:,::-1,:) ; reorder South-to-North
>>> u = u(:,:,::-1,:)
>>> v = v(:,:,::-1,:)
>>>
>>> ===========
>>> Within 'pot_vort_isobaric', other variables of size 210MB are calculated
>>>
>>> pv = -G*(vr + con*(dvdp*dtdx-dudp*dtdy))*dthdp
>>>
>>> EG:
>>>
>>> S = static_stability(p,t,npr,1) ; variable of type list
>>>
>>> s = S[0] ; [...] is list syntax
>>> theta = S[1]
>>> dthdp = S[2]
>>>
>>> also, gradients
>>>
>>> dudp
>>> dvdp
>>> dtdx
>>> dtdy
>>>
>>> vorticity
>>>
>>> vr [vorticity]
>>>
>>> On Sat, Jun 23, 2018 at 8:54 AM, Soma Roy <somaroy892 at gmail.com> wrote:
>>>
>>>> Hello,
>>>>
>>>> I am trying to calculate static stability and potential temperature in
>>>> vertical profile using ERA-Interim pressure level data.
>>>>
>>>> Summary of the input data used in this process is like below;
>>>>
>>>> Variable: pnew
>>>> Type: double
>>>> Total Size: 176 bytes
>>>> 22 values
>>>> Number of Dimensions: 1
>>>> Dimensions and sizes: [level | 22]
>>>> Coordinates:
>>>> level: [ 100..1000]
>>>> Number Of Attributes: 5
>>>> standard_name : air_pressure
>>>> long_name : pressure_level
>>>> units : Pa
>>>> positive : down
>>>> axis : Z
>>>>
>>>> Variable: t
>>>> Type: float
>>>> Total Size: 209943360 bytes
>>>> 52485840 values
>>>> Number of Dimensions: 4
>>>> Dimensions and sizes: [time | 120] x [level | 22] x [latitude | 141]
>>>> x [longitude | 141]
>>>> Coordinates:
>>>> time: [994176..994890]
>>>> level: [ 100..1000]
>>>> latitude: [ 5..40]
>>>> longitude: [65..100]
>>>> Number Of Attributes: 7
>>>> standard_name : air_temperature
>>>> long_name : Temperature
>>>> units : K
>>>> _FillValue_original : -32767
>>>> _FillValue : -32767
>>>> missing_value_original : -32767
>>>> missing_value : -32767
>>>>
>>>> Variable: u
>>>> Type: float
>>>> Total Size: 209943360 bytes
>>>> 52485840 values
>>>> Number of Dimensions: 4
>>>> Dimensions and sizes: [time | 120] x [level | 22] x [latitude | 141]
>>>> x [longitude | 141]
>>>> Coordinates:
>>>> time: [994176..994890]
>>>> level: [ 100..1000]
>>>> latitude: [ 5..40]
>>>> longitude: [65..100]
>>>> Number Of Attributes: 7
>>>> standard_name : eastward_wind
>>>> long_name : U component of wind
>>>> units : m s**-1
>>>> _FillValue_original : -32767
>>>> _FillValue : -32767
>>>> missing_value_original : -32767
>>>> missing_value : -32767
>>>>
>>>> Variable: v
>>>> Type: float
>>>> Total Size: 209943360 bytes
>>>> 52485840 values
>>>> Number of Dimensions: 4
>>>> Dimensions and sizes: [time | 120] x [level | 22] x [latitude | 141]
>>>> x [longitude | 141]
>>>> Coordinates:
>>>> time: [994176..994890]
>>>> level: [ 100..1000]
>>>> latitude: [ 5..40]
>>>> longitude: [65..100]
>>>> Number Of Attributes: 7
>>>> standard_name : northward_wind
>>>> long_name : V component of wind
>>>> units : m s**-1
>>>> _FillValue_original : -32767
>>>> _FillValue : -32767
>>>> missing_value_original : -32767
>>>> missing_value : -32767
>>>>
>>>> My code for this is attached below.
>>>>
>>>> At run time it shows the following;
>>>>
>>>> fatal:NclMalloc Failed:[errno=12]
>>>> Segmentation fault (core dumped)
>>>>
>>>>
>>>> Please inform me how to resolve the issue.
>>>>
>>>> Thanks,
>>>> Soma
>>>>
>>>> _______________________________________________
>>>> 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/20180625/62d29704/attachment.html>
More information about the ncl-talk
mailing list