Re: wrf_user_getvar segfault - for cape_2d

The original code was from wrfhelp [cc'd here].

I ran you script with a standard WRF output file.
I did add two printVarSummary statements.
No problems. The output follows

%> ncl cape_tst.ncl

  NCAR Command Language Version 6.0.0

(0) Get dbz

Variable: dbz
Type: float
Total Size: 1018464 bytes
             254616 values
Number of Dimensions: 3
Dimensions and sizes: [Time | 24] x [south_north | 103] x [west_east | 103]
Coordinates:
Number Of Attributes: 6
   description : Max Reflectivity
   units : dBZ
   FieldType : 104
   MemoryOrder : XYZ
   stagger :
   coordinates : XLONG XLAT
(0) Get CAPE

Variable: cape
Type: float
Total Size: 4073856 bytes
             1018464 values
Number of Dimensions: 4
Dimensions and sizes: [mcape_mcin_lcl_lfc | 4] x [Time | 24] x
[south_north | 103] x [west_east | 103]
Coordinates:
Number Of Attributes: 5
   MemoryOrder : XY
   description : mcape ; mcin ; lcl ; lfc
   FieldType : 104
   stagger :
   coordinates : XLONG XLAT

Try your script with (say) 2 time steps.

Also, for debugging you may wish to do something like the following

$> ncks -O -d Time,0,1 -d west_east,,,2 -d south_north,,,2 -d
south_north_stag,,,2 -d west_east_stag,,,2
wrfout_d03_2009-05-18_12:00:00.nc WRF_TST.nc

On 03/02/2012 02:07 PM, Will Thurston wrote:
> Hi there,
>
> I've had problems with ncl seg faulting when using wrf_user_getvar to
> calculate/get cape from some wrfout files, attached is a pared down
> script.
>
> I've put the get dbz and print statements in to show (myself) that ncl
> can calculate other diagnostic from 3D fields, and it is the cape part
> it falls over on.
>
> The input file is about 22 Gb and has dimensions of 24 x 27 x 723 x
> 732. This guess this is quite large so I also tried getting cape for
> just 1 timestep, rather than the whole file at once, but it still seg
> faults.
>
> ncl -V
> 6.0.0
>
> Cheers,
> Will.
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Mar 2 15:29:33 2012