question about simpne

From: James McPhee <jmcphee_at_nyahnyahspammersnyahnyah>
Date: Thu, 27 Jul 2006 15:42:07 -0700

Hi All,

I'm trying to compute CAPE with radiosonde data, so I need to integrate
  a temperature difference term along a vertical profile. I'm trying to
use simpne to achieve this, but my results are confusing. To illustrate
better my pproblem, let me give an example of what I'm seeing. Given the
150-element vectors:

y1 [-] x1 [m]
0 801
0 2901
0 3224
0 3442
0 3640
0 4312
0 5542
0.007306491 6276
0.01352645 7018
0.1340122 12867
0 12934
0 13948
0 15456
0 16732
0 17734
0 19826
0 20522
0 21978
0 23374
0 27944
0 28729
0 30206
0 30605
0 31766
-999 -999
... ...

These are the results I get from simpne:
ncl 24> print(simpne(x1(6:10),y1(6:10)))
(0) 11838.13
ncl 25> print(simpne(x1(7:9),y1(7:9)))
(0) 377.2643
ncl 26> print(simpne(x1(6:9),y1(6:9)))
(0) 380.2147
ncl 27> print(simpne(x1(7:10),y1(7:10)))
(0) 381.7708
ncl 28> print(simpne(x1,y1))
(0) -999
ncl 29> print(simpne(x1(0:23),y1(0:23)))
(0) 11838.13

If I use a simple trapezoidal rule the result I get for the range 0:23
(everything else is _FillValue), is 446.3767. I'll be grateful if
someone could provide some insight into this...

James
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jul 27 2006 - 16:42:07 MDT

This archive was generated by hypermail 2.2.0 : Fri Jul 28 2006 - 16:18:18 MDT