Re: calculating "CAPE"

Dear Wei,   Actually Alan is right. The dimension of "cape_3d" is =
5, as it is shown below. Can you give more suggestion?   Thanks ve=
ry much. Mark   Copyright (C) 1995-2013 - All Rights Reserved =
 University Corporation for Atmospheric Research  NCAR Command =
Language Version 6.1.2  The use of this software is governed by a Li=
cense Agreement.  See http://www.ncl.ucar.edu/ for more details. V=
ariable: cape_3d Type: float Total Size: 2233079200 bytes  Â=
 Â Â Â Â Â Â Â Â Â  558269800 values=

2] x [Time | 121] x [bottom_top | 59] x [south_north | 170] x [west_east | =
230] Coordinates: Number Of Attributes: 6   _FillValue :  9=
.96921e+36   description : cape ; cin   FieldType : Â=
  104   MemoryOrder : XYZ   stagger :   coordinates =
: XLONG XLAT fatal:["Execute.c":7432]:performASSIGN_VAR_VAR_OP: Number of=
 subscripts on rhs do not match       Â=
 Â Â Â Â Â Â Â Â Â Â Â Â =
     number of dimensions of variable, (4) Subscripts u=
sed, (5) (null) fatal:["Execute.c":8128]:Execute: Error occurred at or ne=
ar line 22 in file my3_CAPE.ncl On Monday, June 16, 2014 11:11 =
PM, Wei Huang <huangwei@ucar.edu> wrote: Mark, Add prin=
tVarSummary(cape_3d) after: cape_3d = wrf_user_getvar(a,"cape_3d",time)=

 the next line should be. Wei On Mon, Jun 16, 2014 at 2:46=
 PM, Mark Chan <cym263@yahoo.com> wrote: Dear Alan, > > >Thanks =
a lot for your email! Do you have any suggestion about how to match the dim=
ensions based on my earlier ncl script? > > >Thanks a lot!Mark >=

lbany.edu> wrote: > > > >Mark,  > > >fatal:["Execut=
e.c":7432]:performASSIGN_VAR_VAR_OP: Number of subscripts on rhs do not mat=
ch >>           =
             nu=
mber of dimensions of variable, (4) Subscripts used, (5) â–’ >> >=
>fatal:["Execute.c":8128]:Execute: Error occurred at or near line 22 in fil=
e my3_CAPE.ncl > > >Line 22:      cape = cape_3=
d(0,:,:,:)     > > > > >So in this line you’=
re telling NCL that cape_3d has 4 dimensions "(4) Subscripts used�.=
 Â  NCL is telling you that there are actually 5 dimensions. > >Pu=
t a  printVarSummary(cape_3d) in Line 21 and edit your subscripts acco=
rdingly.   >The 5 dimensions are I imagine, (time, cape/cin, level,=
 lat, lon ) maybe not in that order.  > > >You probably then w=
ant to also check that the output file CAPE3D matches the output dimensions=
 of cape.  > > > > > > > > >Alan.  > =
> > > >############################## >Alan Brammer, >PhD Stude=
nt, > > >Department of Atmospheric and Environmental Sciences,=

ny, NY, 12222 >abrammer@albany.edu >############################## =

 for calculating CAPE and put the result into a file "CAPE_3D.nc". But I go=
t the error message as: >> >> >>-----error message------------ >>=
  >> Copyright (C) 1995-2013 - All Rights Reserved >> Univers=
ity Corporation for Atmospheric Research >> NCAR Command Language Ve=
rsion 6.1.2 >> The use of this software is governed by a License Agr=
eement. >> See http://www.ncl.ucar.edu/ for more details. >>fatal:=
["Execute.c":7432]:performASSIGN_VAR_VAR_OP: Number of subscripts on rhs do=
 not match >>          =
            Â=
  number of dimensions of variable, (4) Subscripts used, (5) â–’ =
>> >>fatal:["Execute.c":8128]:Execute: Error occurred at or near line 22 =
in file my3_CAPE.ncl >>------------------------------------- >> >> =

based on my script here. >> >> >>Thanks so much! >>Mark >> =
>> >> >>===NCL script====== >> >>load "$NCARG_R=
OOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" >>load "$NCARG_ROOT/lib/ncarg/=
nclscripts/wrf/WRFUserARW.ncl" >> >>begin >>; >>; The WRF ARW input=
 file. >> >>   a = addfile("wrfout_d02_1994-01-11_00:00:00.=
nc","r") >> >>   system("/bin/rm -f " + "CAPE_3D.nc") >> =
>>   b = addfile("CAPE_3D.nc","c") >> >>   dim_na=
mes = getvardims(a) >>   dim_sizes = getfiledimsizes(a) >=
>   dimUnlim = (/True, False, False, False, False, False, False=
, False, False/) >>   filedimdef(b, dim_names, dim_sizes, dimUn=
lim) >>   time = -1 >> >>     cape_3d=
 = wrf_user_getvar(a,"cape_3d",time) >> >>     ca=
pe = cape_3d(0,:,:,:) >> >>   lat = a->XLAT >> Â=
  lon = a->XLONG >>   times = a->Times >>   f=
ileattdef(b,a) >> >>   b->Times = times >>   b-=
>XLAT = lat >>   b->XLONG = lon >>   b->CAPE3D =
= cape >> >>end >> >> >> >> ____________________________=
___________________ >>ncl-talk mailing list >>List instructions, subscr=
iber options, unsubscribe: >>http://mailman.ucar.edu/mailman/listinfo/ncl=
-talk >> > > > >__________________________________________=
_____ >ncl-talk mailing list >List instructions, subscriber options, un=
subscribe: >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 Tue Jun 17 05:52:55 2014