[gradsusr] Computation of thermodynamic Indices

[Geeta Geeta]

Dear All, 
This problem is not particularly pertaining to Grads but  I am still putting it across in anticipation of a possible help.


It goes like this, I want to compute all the Standard Indices  
( like Showalter, K, SWEAT, CAPE, CINE  etc)
from the Rs data of a station. 


I am able to compute only LCL and then the moist adiabat (using the conservation of MOIST static energy) 
and I have computed the LCL from the RS data and then constructed the moist adiabat of the parcel till 200hPa. Part of the output of moist adiabat looks like this, 

(I)
INPUT LCL

 P T Z
 871.0 294.6 312.0
 MOIST ADIABAT
 P Z T Q NO OF ITRN
 871.0 312.0 294.6 .1952E-01 .3476E+06 0
 870.0 322.0 294.5 .1950E-01 .3476E+06 39
 869.0 332.1 294.5 .1947E-01 .3476E+06 39
 868.0 342.1 294.4 .1945E-01 .3476E+06 40
 867.0 352.2 294.4 .1942E-01 .3476E+06 39
 866.0 362.2 294.4 .1940E-01 .3476E+06 39
 865.0 372.3 294.3 .1938E-01 .3476E+06 40
 864.0 382.4 294.3 .1935E-01 .3476E+06 39

"MOIST.DAT" 676 lines, 44480 characters

(II)
My RS data of a station is at an interval of 50hPa and looks like this (Ht, TT, Td, WS, WD at different levels)


1981 03 10 00 2 1026 20.6 16.8 155 02 1520 18.6 11.2 201 04 ^M
1981 03 10 00 3 2038 15.6 8.2 195 03 2584 12.4 3.8 079 02 
 3159 8.4 -5.3 105 04 3768 4.8 044 03 




Kindly help me reg the computation of CAPE and CINE after I construct the MOIST adiabat for a station. 
geeta 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110711/1c9a68ba/attachment-0003.html