Hi,
I have a question. What is the use of the subroutine
cape = cape_thermo(p,tc,plcl,iprnt)
in skewt_func.ncl ?
It recomputes the plcl, but in a weird way in my sense.
In my case, it recomputes the LCL to be at the first level, which is
not saturated.
The saturation occurs at 951.9 mb, not at 994 mb. Why does this
subroutine recomputes the LCL pressure??
Many thanks
cape.f: LCL pressure: original 951.8922729492188 mb
cape.f: adjusted to 994.0000000000000 mb
cape.f: and lies at level 1
data= 994.0000000000000 287.7500003814697
data= 925.0000000000000 282.7500003814697
data= 907.0000000000000 281.3799995422363
data= 866.0000000000000 278.1500000000000
data= 850.0000000000000 277.5500000953674
data= 811.0000000000000 275.7499999046325
data= 753.0000000000000 273.1500000000000
data= 725.0000000000000 271.0500000953674
data= 700.0000000000000 269.2499999046325
data= 673.0000000000000 267.2499999046325
data= 646.0000000000000 265.0499996185303
data= 596.0000000000000 260.6500000000000
data= 582.0000000000000 260.8499998092651
data= 500.0000000000000 253.6500000000000
data= 492.0000000000000 252.8500007629394
data= 458.0000000000000 248.4499992370605
data= 446.0000000000000 246.7600006103515
data= 400.0000000000000 239.8500007629394
data= 359.0000000000000 233.2499984741211
data= 304.0000000000000 224.2499984741211
data= 300.0000000000000 223.6500000000000
data= 282.0000000000000 220.9999984741211
data= 250.0000000000000 215.8500007629394
data= 245.0000000000000 215.0500015258789
data= 223.0000000000000 216.4400009155273
data= 200.0000000000000 218.0500015258789
data= 167.0000000000000 220.8500007629394
data= 163.0000000000000 220.4000000000000
data= 150.0000000000000 218.8500007629394
data= 134.0000000000000 216.8500007629394
data= 107.0000000000000 218.0500015258789
data= 100.0000000000000 217.8500007629394
cape.f: initial parcel pressure: 994.0000000000000 mb
cape.f: parcel temp is 287.7500003814697 K
cape.f: initial P and T lie at level 1
cape.f: P zlev tpar tenv
1 994.00 0.00 287.75 287.75
2 925.00 0.07 284.82 282.75
3 907.00 0.09 283.95 281.38
4 866.00 0.14 282.05 278.15
5 850.00 0.16 281.25 277.55
6 811.00 0.20 279.25 275.75
7 753.00 0.28 275.95 273.15
8 725.00 0.32 274.25 271.05
9 700.00 0.35 272.65 269.25
10 673.00 0.39 270.75 267.25
11 646.00 0.43 268.85 265.05
12 596.00 0.51 264.85 260.65
13 582.00 0.54 263.65 260.85
14 500.00 0.69 255.65 253.65
15 492.00 0.70 254.75 252.85
16 458.00 0.77 250.65 248.45
17 446.00 0.80 249.15 246.76
18 400.00 0.91 242.75 239.85
19 359.00 1.02 236.15 233.25
20 304.00 1.18 225.85 224.25
21 300.00 1.20 225.05 223.65
22 282.00 1.26 221.25 221.00
23 250.00 1.38 213.95 215.85
24 245.00 1.40 212.75 215.05
25 223.00 1.49 207.15 216.44
26 200.00 1.60 200.85 218.05
27 167.00 1.78 190.75 220.85
28 163.00 1.81 189.45 220.40
29 150.00 1.89 185.05 218.85
30 134.00 2.00 179.15 216.85
31 107.00 2.23 168.05 218.05
32 100.00 2.30 164.85 217.85
LFC pressure= 994.0000000000000 Crossing pressure=
282.0000000000000
jlfc, jcross= 1 22
Env Temp & Parcel Temp
1 994.0 0.00 0.00 -273.15 14.60287.75
2 925.0 21.37 0.07 -273.15 11.67284.82
3 907.0 34.46 0.09 -273.15 10.80283.95
4 866.0 77.40 0.14 -273.15 8.90282.05
5 850.0 97.74 0.16 -273.15 8.10281.25
6 811.0 146.27 0.20 -273.15 6.10279.25
7 753.0 213.35 0.28 -273.15 2.80275.95
8 725.0 245.98 0.32 -273.15 1.10274.25
9 700.0 279.21 0.35 -273.15 -0.50272.65
10 673.0 318.16 0.39 -273.15 -2.40270.75
11 646.0 361.05 0.43 -273.15 -4.30268.85
12 596.0 453.53 0.51 -273.15 -8.30264.85
13 582.0 477.41 0.54 -273.15 -9.50263.65
14 500.0 582.01 0.69 -273.15 -17.50255.65
15 492.0 591.04 0.70 -273.15 -18.40254.75
16 458.0 633.17 0.77 -273.15 -22.50250.65
17 446.0 650.66 0.80 -273.15 -24.00249.15
18 400.0 733.29 0.91 -273.15 -30.40242.75
19 359.0 823.30 1.02 -273.15 -37.00236.15
20 304.0 930.69 1.18 -273.15 -47.30225.85
21 300.0 936.39 1.20 -273.15 -48.10225.05
22 282.0 951.04 1.26 -273.15 -51.90221.25
-- Christian Pagé http://meteocentre.com/ http://meteoalerte.com/ Etudiant au Doctorat en Sciences de l'environnement UQAM_______________________________________________ ncl-talk mailing list ncl-talk AT ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Tue Oct 12 2004 - 17:06:26 MDT