constant output with function wrf_cape_2d

From: samirpokhrel <samir_at_nyahnyahspammersnyahnyah>
Date: Tue, 16 Jun 2009 18:44:34 +0630

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 ( 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"


  a = addfile("","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 = ""
  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


 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:

Received on Tue Jun 16 2009 - 06:14:34 MDT

This archive was generated by hypermail 2.2.0 : Fri Jun 19 2009 - 13:23:25 MDT