problem of PotVortIsobaric

Hi,

I had a problem to use the function of "PotVortIsobaric".
The error message is below:

"Error occurred at or near line 3831 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl
 And, in the line near 162 in file pv.ncl "

The line near 162 in pv.ncl is where I call PotVortIsobaric
Here is the brief description in my scrip:

   p0 = 100000. ; since ps is in Pa
   pm60 = pres_hybrid_ccm(PS60a,p0,hyam60,hybm60)
   pm30 = pres_hybrid_ccm(PS30a,p0,hyam30,hybm30)

   theta60 = ppT60*(p0/pm60)^0.286
   theta30 = ppT30*(p0/pm30)^0.286

  copy_VarCoords(ppT60,theta60)
  copy_VarCoords(ppT30,theta30)

  printVarSummary(theta60)
  printVarSummary(theta30)

  lev30b = lev30*100.0
  lev60b = lev60*100.0

 copy_VarAtts(lev30,lev30b)
 copy_VarAtts(lev60,lev60b)

  printVarSummary(lev30b)
  printVarSummary(lev60b)

  PV30 = PotVortIsobaric(lev30b,ppU30,ppV30,ppT30,lat30, 0, 0)
  PV60 = PotVortIsobaric(lev60b,ppU60,ppV60,ppT60,lat60, 0, 0)

Thank you very much for your help

Sincerely,
Leslie

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Jan 12 19:37:18 2013