[ncl-talk] Stability profile

[Dennis Shea]

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/20180623/1bd28474/attachment.html>