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
This archive was generated by hypermail 2.1.8 : Wed Jul 23 2014 - 15:33:46 MDT