Re: calculating "CAPE"

From: Mark Chan <cym263_at_nyahnyahspammersnyahnyah>
Date: Tue Jun 17 2014 - 01:10:08 MDT

Dear Wei, Dear Alan,   my issue about "CAPE" has been solve=
d, thanks so much for your kindly help.   Best regards, Mark =
  On Monday, June 16, 2014 11:11 PM, Wei Huang <huangwei@ucar=
.edu> wrote: Mark, Add printVarSummary(cape_3d) after:=

ape_3d look like, and then can figure out what the next line should be. =

oo.com> wrote: Dear Alan, > > >Thanks a lot for your email! Do y=
ou have any suggestion about how to match the dimensions based on my earlie=
r ncl script? > > >Thanks a lot!Mark > > > >On Monday, June =
16, 2014 4:49 PM, "Brammer, Alan P" <abrammer@albany.edu> wrote: > >=

R_VAR_OP: Number of subscripts on rhs do not match >>   =
            =
        number of dimensions of varia=
ble, (4) Subscripts used, (5) ▒ >> >>fatal:["Execute.c":8128]:E=
xecute: Error occurred at or near line 22 in file my3_CAPE.ncl > > =

3d has 4 dimensions "(4) Subscripts used”.   NCL is telling yo=
u that there are actually 5 dimensions. > >Put a  printVarSummary=
(cape_3d) in Line 21 and edit your subscripts accordingly.   >The 5=
 dimensions are I imagine, (time, cape/cin, level, lat, lon ) maybe not in =
that order.  > > >You probably then want to also check that th=
e output file CAPE3D matches the output dimensions of cape.  > > =

################### >Alan Brammer, >PhD Student, > > >Departme=
nt of Atmospheric and Environmental Sciences, >University at Albany,=
 State University of New York, Albany, NY, 12222 >abramme=
r@albany.edu >############################## > >On Jun 15, 2014, =
at 16:16, Mark Chan <cym263@yahoo.com> wrote: > >Dear NCL users, >>=

 put the result into a file "CAPE_3D.nc". But I got the error message as: =

) 1995-2013 - All Rights Reserved >> University Corporation for Atmo=
spheric Research >> NCAR Command Language Version 6.1.2 >> Th=
e use of this software is governed by a License Agreement. >> See ht=
tp://www.ncl.ucar.edu/ for more details. >>fatal:["Execute.c":7432]:perfo=
rmASSIGN_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 =
>>------------------------------------- >> >> >>I would appreciate i=
f you could give more suggestion about the "CAPE" based on my script here. =

 script====== >> >>load "$NCARG_ROOT/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 = add=
file("CAPE_3D.nc","c") >> >>   dim_names = getvardims(a) =
>>   dim_sizes = getfiledimsizes(a) >>   dimUnlim =
= (/True, False, False, False, False, False, False, False, False/) >>=
   filedimdef(b, dim_names, dim_sizes, dimUnlim) >>  =
 time = -1 >> >>     cape_3d = wrf_user_getvar(=
a,"cape_3d",time) >> >>     cape = cape_3d(0,:,:,=
:) >> >>   lat = a->XLAT >>   lon = a->XLONG=

   b->XLONG = lon >>   b->CAPE3D = cape >> >>=
end >> >> >> >> ______________________________________________=
_ >>ncl-talk mailing list >>List instructions, subscriber options, unsu=
bscribe: >>http://mailman.ucar.edu/mailman/listinfo/ncl-talk >> > =

 mailing list >List instructions, subscriber options, unsubscribe: >htt=
p://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 07:10:16 2014

This archive was generated by hypermail 2.1.8 : Wed Jul 23 2014 - 15:33:46 MDT