Re: vertical plot

I was about to send the following when I saw Jamie
Scott's reply

What I sent before was a bit crude.

[1]
Use "pres_hybrid_ccm" to calculate the *actual* pressures at each
hybrid pressure level at each time step.

       PH = pres_hybrid_ccm(...)

[2]
Use the "where" to mask out all grid points and levels
< sfcP

   x = where(PH.gt.PS, x@_FillValue, x)

[3] I assume you have "Z3" ... geopotential height, maybe, also
     PHIS (sfc geopotential)

   ZZ = Z3 ; copy meta data
   ZZ = where(ismissing(x), Z3@_FillValue, Z3)

[4]
Use int2p_n_Wrap .... [similar to linint1_n_Wrap] to interpolate

D

On 5/19/10 8:45 AM, Michael Notaro wrote:
> 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
Received on Wed May 19 09:31:05 2010