Re: vertical plot

I basically want to make the same plot as attached but with height above the ground as the
y axis, not pressure. Right now I am showing hybrid sigma pressure versus month.
I computed the average temperature in a region for each month and each hybrid sigma pressure
level and plotted in the attached figure.

I altered the land surface and so the strongest response is near the ground. This is true
both for low and high elevation regions. That is why I want the y axis to show height above the ground.

Mike

On May 19, 2010, at 9:35 AM, Dennis Shea wrote:

> 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