Changing the code as you indicate would be difficult.
The code is hard-wired to reproduce the "USAF Skew-t, log p diagram
(form dod-wpc 9-16-1).
As indicated in the code, some intersections are computed
and some are specified.
I think the 'trial and error' approach is the only option.
Unfortunately, I do have the time to halp with the modifed skewT.
Good Luck
===================================================================
[SNIP]
; --- declare isotherm [C] values and pressures [hPa] where isotherms
; --- intersect the edge of the skew-t diagram.
temp = (/-100.,-90.,-80.,-70.,-60.,-50.,-40.,-30. \
, -20.,-10., 0., 10., 20., 30., 40., 50. /)
lendt= (/ 132., 181., 247., 337., 459., 625., 855. \
,1050.,1050.,1050.,1050.,1050.,1050.,1050.,1050.,1050./)
rendt= (/ 100., 100., 100., 100., 100., 100., 100. \
, 135., 185., 251., 342., 430., 500., 580., 730., 993./)
ntemp= dimsizes (temp) ; isotherms
if (dimsizes(lendt).ne.dimsizes(rendt) .or. \
dimsizes(lendt).ne.ntemp) then
print ("Error: dimsizes temp, lendt, rendt do not match")
exit
end if
; --- declare pressure values [hPa] and x coordinates of the endpoints
of each
; --- isobar. these x,y values are computed from the equations in the
; --- transform functions listed at the beginning of this program.
; --- refer to a skew-t diagram for reference if necessary.
pres = (/1050.,1000.,850.,700.,500.,400.,300.,250.,200.,150.,100./)
xpl = (/-19.0,-19.0,-19.0,-19.0,-19.0,-19.0,-19.0,-19.0 \
,-19.0,-19.0,-19.0/)
xpr = (/27.1,27.1,27.1,27.1,22.83,18.6,18.6,18.6,18.6,18.6,18.6/)
npres= dimsizes (pres) ; pressures
if (dimsizes(xpl).ne.dimsizes(xpr) .or. dimsizes(xpl).ne.npres) then
print ("Error: dimsizes pres, xpl, xpr do not match")
exit
end if
; --- declare adiabat values [C] and pressures where adiabats intersect the
; --- edge of the skew-t diagram. refer to a skew-t diagram if necessary.
theta= (/-30.,-20.,-10.,0.,10.,20.,30.,40.,50.,60.,70.,80.,90. \
,100.,110.,120.,130.,140.,150.,160.,170. /)
lendth= (/880.,670.,512.,388.,292.,220.,163.,119.,100.,100.,100. \
,100.,100.,100.,100.,100.,100.,100.,100.,100.,100. /)
rendth= (/1050.,1050.,1050.,1050.,1050.,1050.,1050.,1050. \
,1003.,852.,728.,618.,395.,334.,286., 245., 210. \
[SNIP]
======================================================================
On 12/20/2011 10:21 AM, Manisha Ganeshan wrote:
> Hi,
>
> I am trying to modify the skewt_func.ncl code to be able to plot one
> additional dry adiabat, saturation mixing ratio line and moist adiabat
> from my sounding data.
> I noticed that the skewt_func.ncl file has these values hardwired into
> the code. For example, dry adiabat values are declared as potential
> temperature values at 1000hPa level. The points where these intersect
> the left and right edge of the skewt (y-axis) are also declared. Do you
> know how these values are calculated?
>
> For example, I want to plot one additional adiabat say with potential
> temperature 32 degrees C at 1000hPa, but I don't see how can I determine
> the value of the intersection of this adiabat exactly at the y-axis
> (edge of the skewT).
>
> Kindly help.
>
> Nisha
>
>
> _______________________________________________
> 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 Tue Dec 20 10:47:16 2011
This archive was generated by hypermail 2.1.8 : Wed Dec 21 2011 - 10:44:06 MST