[gradsusr] vint ?

[Li Dong]

Dear list,

I need to calculate mass-weighted mean values of a variable, say air
temperature, from 1000 to 100 mb. The air temperature is at levels of
1000, 925, 800, 700, 600, 500, 400, 300, 200 and 100mb. I am thinking of
using function "vint" to perform this calculation. I'd appreciate if you
could verify if my script below is right.

'open ps.ctl'        (surface pressure control file, ps in millibar)
'open T.ctl'          (air temperature file)
'd vint(ps.1,T.2,150)/((ps.1-100)*(100/9.8))'

Note that the denominator is the scale factor (f/g) and the entire
column between surface level and top layer in millibar.

Thanks so much in advance!

Li