Re: vertical plot

Hi Mike

It is not clear to me what you are looking for ...
"vertical profile" or "vertical cross-section"

    x = f->X(:,:,{lats:latN},{lonL:lonR)) ; (time,lev,lat,lon)
    ps = f->PS(:,{lats:latN},{lonL:lonR)) ; (time,lat,lon)

    lev = x&lev ; nominal pressure levels (:)
    LEV = conform(x,lev,1) ; (ntim,klev,nlat,mlon)
    PS = conform(x, ps, (/2,3/))

    x = where(PS.lt.LEV, x@_FillValue, x)

The above will set all levels less than the surface pressure to _fillValue

Cross sections can be plotted for each time.
You could also average acroos times for climatologies.

Vertical profiles: Perform areal averages at each time step.

D

On 5/19/10 8:05 AM, Michael Notaro wrote:
> Does anyone have a script or recommended approach for plotting a
> vertical profile of temperature for a region in CCSM3, where the y axis
> is height above the ground? The model uses hybrid sigma pressure coordinates.
>
> Thanks, Michael
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed May 19 08:35:10 2010