Re: Trying to calculate vertical wind shear

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Nov 15 2010 - 14:12:45 MST

The function to use is

https://www.ncl.ucar.edu/Document/Functions/Built-in/pres_hybrid_ccm.shtml

This calculates the pressure at every grid point level.
It can be placed in your 'pt' calculation.

*** Please use "printVarSummary" often. ****

On 11/15/2010 11:45 AM, sbasu@iarc.uaf.edu wrote:
> Hello,
>
> I am trying to calculate the vertical gradient of potential temperature
> and zonal wind. I am trying to use 'center_finite_diff_n' method in NCL.So
> far I am able to plot d(pt)/dp, pt-potential temp. Is it possible to
> calculate du/dz and d(pt)/dz by this method ?
> My dataset is standard output from CAM 3.1_p2 output files. I am attaching
> the script I am using.
> I am getting this error:
> fatal:center_finite_diff_n: r must either be a scalar, a 1D array the same
> length as the dim-th dimemsion of q, or the same size as q
> fatal:Execute: Error occurred at or near line 58 in file
> baroclinincity_Xseries5.ncl
>
> Thanks,
> Soumik
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================


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

Received on Mon Nov 15 14:12:54 2010

This archive was generated by hypermail 2.1.8 : Wed Nov 17 2010 - 13:14:26 MST