Dear all,
I am very new to NCL and I am trying to calculate cape from Japan reanalysis data
(JRA25) using the ncl function wrf_cape_2d but I am getting a constant value. I am
attaching both the data file (test2d.nc) and ncl script. May I know where I went wrong.
The ncl script is :--------
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
a = addfile("test2d.nc","r")
P = a->PL
T = a->TK
Q = a->QG
Z = a->GPH
Z = Z/9.81
Q = Q
ZSFC = a->SPF3(:,0,:,:)
ZSFC = ZSFC/9.81
PSFC = a->PS(:,0,:,:)
PSFC = PSFC/100
cinfo = wrf_cape_2d( P, T, Q, Z, ZSFC, PSFC, False )
cape = cinfo(0,:,:,:)
copy_VarCoords(ZSFC, cape)
filo = "cape2d.nc"
system ("/bin/rm "+filo) ; remove any pre-existing file
fo = addfile(filo , "c") ; open output file
fo_at_title = "CAPE/CIN from JRA-25"
cape_at_long_name = "Convective Available Potential Energy"
cape_at_units = "j/kg"
fo->cape = cape
end
Thanking in advance
Dr. Samir Pokhrel
Climate and Global Modelling Division
Indian Institute of Tropical Meteorology
Dr. Homi Bhabha Road Pashan Pune-411008
Office - 020-25904307
Mobile - 09881745234
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.2.0 : Fri Jun 19 2009 - 13:23:25 MDT