calculating "CAPE"

Dear NCL users, In the below is an ncl script, which is for calculatin=
g CAPE and put the result into a file "CAPE_3D.nc". But I got the error mes=
sage as: -----error message------------  Copyright (C) 1995-=
2013 - All Rights Reserved  University Corporation for Atmospheric R=
esearch  NCAR Command Language Version 6.1.2  The use of this=
 software is governed by a License Agreement.  See http://www.ncl.uc=
ar.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 -----------------------=
-------------- I would appreciate if you could give more suggestion=
 about the "CAPE" based on my script here. Thanks so much! Mark =

g/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("C=
APE_3D.nc","c")    dim_names = getvardims(a)  Â=
  dim_sizes = getfiledimsizes(a)    dimUnlim = (/True, Fa=
lse, False, False, False, False, False, False, False/)    filed=
imdef(b, dim_names, dim_sizes, dimUnlim)    time = -1 =
     cape_3d = wrf_user_getvar(a,"cape_3d",time) =

 = a->XLAT    lon = a->XLONG    times = a->Ti=
mes    fileattdef(b,a)    b->Times = times =
   b->XLAT = lat    b->XLONG = lon   =
 b->CAPE3D = cape end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk


Received on Sun Jun 15 08:16:14 2014