Thank you.
I first followed the prototype, then the list of arguments (leading to the
samples below). When I returned to the proper order of inputs (matching
the prototype) the function returns the value "-0.1" for levels 1-90, and
if I include the level of missing data (level 0), the function returns 0.
I'm not so worried about the level with missing data at the moment.
I played with the units some, and I did not stumble upon any winning match,
rather the given (K, m, kg/kg, hPa) units seemed the best (though still
providing problematic results.
One additional question just to make sure I am thinking correctly: I
currently have pressure decreasing as level index increases. Am I correct
in believing that this is a "top to bottom" ordering of the data?
Any further help on this function (or even code to accomplish the same)
would be greatly appreciated. Thank you for your time and effort already.
-- Jacob Klee
Mary Haley
<haley_at_ucar.edu>
To
04/23/2007 11:07 Jacob.Klee_at_dom.com
PM cc
ncl-talk_at_ucar.edu
Subject
Re: [ncl-talk] rip_cape_3d Question
Hi Jacob,
I think the problem is in your third argument to rip_cape_3d. It
should be units of kg/kg, but you are passing HGT, which has fairly
large values. I think you might have HGT and SPH swapped.
The documentation is a bit confusing, because it lists the arguments
in the correct order in the prototype section, but they are
out-of-order in the "arguments" section. I will reorder the arguments
in the documentation so they are the same as the prototype.
--Mary
On Mon, 23 Apr 2007 Jacob.Klee_at_dom.com wrote:
> I tried 'subscripting out' the missing data (see below) with no change in
> my results.
>
> My data is in hPa, K, m, kg/kg, m, and hPa.... This is what is in the
> documentation, are these truely correct?
>
> -- Jacob
>
> CAPE_CIN
=rip_cape_3d(P(1:90,20,18),TEM(1:90,20,18),HGT(1:90,20,18),SPH(1:90,20,18),HGT_sf
> c(20,18),P_sfc(20,18),False)
>
> -----ncl-talk-bounces_at_ucar.edu wrote: -----
>
> To: Jacob.Klee_at_dom.com
> From: Mary Haley <haley_at_ucar.edu>
> Sent by: ncl-talk-bounces_at_ucar.edu
> Date: 04/23/2007 11:24AM
> cc: ncl-talk_at_ucar.edu
> Subject: Re: rip_cape_3d Question
>
>
> Hi Jacob,
>
> I think the problem is that your data contains missing
> values.
> The rip_cape_3d code is not designed to handle missing
> values, so
> it treats them as real values.
>
> --Mary
>
>
> On Mon, 23 Apr 2007 Jacob.Klee_at_dom.com wrote:
>
> > Hello,
> >
> > I am having some trouble with the rip_cape3d function. I
> am running an
> > early 4.3.0 in CYGWIN.
> >
> > Any thoughts or observations of what is wrong?
> >
> > -- Jacob
> >
> >
> >
> > The error I am encountering is:
> >
> > "capecalc3d: Outside of lookup table bounds. prs,thte=
> 990.01001 INF"
> >
> >
> >
> > My call to the function is:
> >
> > CAPE_CIN
=rip_cape_3d(P(:,20,18),TEM(:,20,18),HGT(:,20,18),SPH(:,20,18)*1000,HGT_sfc(20,
> 18
> > ),P_sfc(20,18),False)
> >
> >
> >
> > The input data is as follows:
> >
> > Variable: P (subsection)
> > Type: float
> > Total Size: 364 bytes
> > 91 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [lv_ISBLx | 91]
> > Coordinates:
> > Number Of Attributes: 1
> > _FillValue : -999
> > (0) 1000.01
> > (1) 990.01
> > (2) 980.01
> > (3) 970.01
> > (4) 960.01
> > (5) 950.01
> > (6) 940.01
> > (7) 930.01
> > (8) 920.01
> > (9) 910.01
> > (10) 900.01
> > (11) 890.01
> > (12) 880.01
> > (13) 870.01
> > (14) 860.01
> > (15) 850.01
> > (16) 840.01
> > (17) 830.01
> > (18) 820.01
> > (19) 810.01
> > (20) 800.01
> > (21) 790.01
> > (22) 780..01
> > (23) 770.01
> > (24) 760.01
> > (25) 750.01
> > (26) 740.01
> > (27) 730.01
> > (28) 720.01
> > (29) 710.01
> > (30) 700.01
> > (31) 690.01
> > (32) 680.01
> > (33) 670.01
> > (34) 660.01
> > (35) 650.01
> > (36) 640.01
> > (37) 630..01
> > (38) 620.01
> > (39) 610.01
> > (40) 600.01
> > (41) 590.01
> > (42) 580.01
> > (43) 570.01
> > (44) 560.01
> > (45) 550.01
> > (46) 540.01
> > (47) 530.01
> > (48) 520.01
> > (49) 510.01
> > (50) 500.01
> > (51) 490.01
> > (52) 480..01
> > (53) 470.01
> > (54) 460.01
> > (55) 450.01
> > (56) 440.01
> > (57) 430.01
> > (58) 420.01
> > (59) 410.01
> > (60) 400.01
> > (61) 390.01
> > (62) 380.01
> > (63) 370.01
> > (64) 360.01
> > (65) 350.01
> > (66) 340.01
> > (67) 330..01
> > (68) 320.01
> > (69) 310.01
> > (70) 300.01
> > (71) 290.01
> > (72) 280.01
> > (73) 270.01
> > (74) 260.01
> > (75) 250.01
> > (76) 240.01
> > (77) 230.01
> > (78) 220.01
> > (79) 210.01
> > (80) 200.01
> > (81) 190.01
> > (82) 180..01
> > (83) 170.01
> > (84) 160.01
> > (85) 150.01
> > (86) 140.01
> > (87) 130.01
> > (88) 120.01
> > (89) 110.01
> > (90) 100.01
> >
> > Variable: TEM (subsection)
> > Type: float
> > Total Size: 364 bytes
> > 91 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [lv_ISBLx | 91]
> > Coordinates:
> > Number Of Attributes: 1
> > _FillValue : -999
> > (0) -999
> > (1) 300.0381
> > (2) 299.6239
> > (3) 299.0023
> > (4) 298.1691
> > (5) 297.3273
> > (6) 296.7134
> > (7) 296.0933
> > (8) 295.412
> > (9) 294.6685
> > (10) 293.9168
> > (11) 293.1544
> > (12) 292.3835
> > (13) 291.6212
> > (14) 290.8677
> > (15) 290.1055
> > (16) 289.405
> > (17) 288.6962
> > (18) 287.986
> > (19) 287.2744
> > (20) 286.554
> > (21) 285.8527
> > (22) 285.1425
> > (23) 284.4469
> > (24) 283.7663
> > (25) 283.0766
> > (26) 282.4055
> > (27) 281.7254
> > (28) 281.0365
> > (29) 280.3386
> > (30) 279.6308
> > (31) 278.8593
> > (32) 278.0764
> > (33) 277.282
> > (34) 276.4756
> > (35) 275.6568
> > (36) 274.7976
> > (37) 273.9248
> > (38) 273.0381
> > (39) 272.137
> > (40) 271.2209
> > (41) 270.4594
> > (42) 269.685
> > (43) 268.8972
> > (44) 268.0954
> > (45) 267.2791
> > (46) 266.6681
> > (47) 266.0459
> > (48) 265.4117
> > (49) 264.7654
> > (50) 264.1061
> > (51) 263.4885
> > (52) 262.8582
> > (53) 262.2145
> > (54) 261.5571
> > (55) 260.8852
> > (56) 259.9771
> > (57) 259.0478
> > (58) 258.0967
> > (59) 257.1227
> > (60) 256.1246
> > (61) 254.7971
> > (62) 253.4348
> > (63) 252.0362
> > (64) 250.5993
> > (65) 249.1219
> > (66) 247.5164
> > (67) 245.863
> > (68) 244.1586
> > (69) 242.4001
> > (70) 240.584
> > (71) 238.5219
> > (72) 236.3871
> > (73) 234.2024
> > (74) 231..9635
> > (75) 229.6367
> > (76) 227.2763
> > (77) 224.8154
> > (78) 222.5061
> > (79) 220.3575
> > (80) 218.1041
> > (81) 215.8187
> > (82) 213.4099
> > (83) 211.5183
> > (84) 210.19
> > (85) 208.7759
> > (86) 208.193
> > (87) 207.5679
> > (88) 207.2869
> > (89) 207.395
> > (90) 207.5133
> >
> > Variable: HGT (subsection)
> > Type: float
> > Total Size: 364 bytes
> > 91 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [lv_ISBLx | 91]
> > Coordinates:
> > Number Of Attributes: 1
> > _FillValue : -999
> > (0) -999
> > (1) 179.4146
> > (2) 269.5749
> > (3) 360.3191
> > (4) 451.6597
> > (5) 543.9567
> > (6) 636.6967
> > (7) 730.428
> > (8) 824.8289
> > (9) 919.9144
> > (10) 1016.051
> > (11) 1112.524
> > (12) 1210.087
> > (13) 1308.552
> > (14) 1407.942
> > (15) 1508.495
> > (16) 1609.129
> > (17) 1710.967
> > (18) 1813.948
> > (19) 1918.101
> > (20) 2023.547
> > (21) 2129.177
> > (22) 2236.15
> > (23) 2344.133
> > (24) 2453.153
> > (25) 2563.616
> > (26) 2675.239
> > (27) 2788.38
> > (28) 2902.593
> > (29) 3017.909
> > (30) 3134.861
> > (31) 3252.454
> > (32) 3371.763
> > (33) 3492.84
> > (34) 3615.737
> > (35) 3740.51
> > (36) 3864.917
> > (37) 3991.279
> > (38) 4119.664
> > (39) 4250.136
> > (40) 4382.765
> > (41) 4515.67
> > (42) 4650.846
> > (43) 4788.373
> > (44) 4928.334
> > (45) 5070.817
> > (46) 5213.727
> > (47) 5359.306
> > (48) 5507.657
> > (49) 5658.89
> > (50) 5813.118
> > (51) 5968.848
> > (52) 6127.788
> > (53) 6290.074
> > (54) 6455.85
> > (55) 6625.27
> > (56) 6795.643
> > (57) 6969.929
> > (58) 7148.317
> > (59) 7331.003
> > (60) 7518.201
> > (61) 7705.781
> > (62) 7898.229
> > (63) 8095.811
> > (64) 8298.806
> > (65) 8507.519
> > (66) 8715.626
> > (67) 8929.939
> > (68) 9150.848
> > (69) 9378.771
> > (70) 9614.167
> > (71) 9850.447
> > (72) 10095.01
> > (73) 10345.53
> > (74) 10602.5
> > (75) 10869.55
> > (76) 11140.67
> > (77) 11423.32
> > (78) 11714.69
> > (79) 12015.64
> > (80) 12331.29
> > (81) 12653.73
> > (82) 12993.6
> > (83) 13349.37
> > (84) 13723.05
> > (85) 14120.85
> > (86) 14540.64
> > (87) 14991.53
> > (88) 15477.61
> > (89) 16005.02
> > (90) 16582.73
> >
> > Variable: SPH (subsection)
> > Type: float
> > Total Size: 364 bytes
> > 91 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [lv_ISBLx | 91]
> > Coordinates:
> > Number Of Attributes: 1
> > _FillValue : -999
> > (0) -999
> > (1) 0.01600164
> > (2) 0.01576904
> > (3) 0.01543304
> > (4) 0.0149916
> > (5) 0.01454555
> > (6) 0.01356957
> > (7) 0.01258263
> > (8) 0.01199793
> > (9) 0.01182368
> > (10) 0.0116475
> > (11) 0.0113606
> > (12) 0.01107035
> > (13) 0.0107751
> > (14) 0.01047475
> > (15) 0.01017089
> > (16) 0.009880443
> > (17) 0.009586534
> > (18) 0.009296402
> > (19) 0.009010126
> > (20) 0.008720292
> > (21) 0.008474605
> > (22) 0.008225836
> > (23) 0.007987712
> > (24) 0.007760477
> > (25) 0.007530232
> > (26) 0.007330647
> > (27) 0.00712838
> > (28) 0.006931528
> > (29) 0.006740219
> > (30) 0.006546196
> > (31) 0.00642047
> > (32) 0.00629298
> > (33) 0.006163602
> > (34) 0.006032278
> > (35) 0.005898949
> > (36) 0.00570709
> > (37) 0.005512153
> > (38) 0.005314096
> > (39) 0.00511282
> > (40) 0.004908216
> > (41) 0.004676931
> > (42) 0.004441668
> > (43) 0.004202314
> > (44) 0.003958724
> > (45) 0.003710744
> > (46) 0.003543196
> > (47) 0.003372603
> > (48) 0.003198761
> > (49) 0.003021543
> > (50) 0.002840815
> > (51) 0.0026752
> > (52) 0.002506188
> > (53) 0.002333619
> > (54) 0.002157338
> > (55) 0.001977182
> > (56) 0.001806408
> > (57) 0.001631721
> > (58) 0.001452924
> > (59) 0.001269818
> > (60) 0.00108219
> > (61) 0.0009804764
> > (62) 0.0008762126
> > (63) 0.0007691681
> > (64) 0.0006591906
> > (65) 0.0005461148
> > (66) 0.000493438
> > (67) 0.0004392533
> > (68) 0.0003834011
> > (69) 0.0003257756
> > (70) 0.0002662606
> > (71) 0.0002394213
> > (72) 0.0002116755
> > (73) 0.0001864293
> > (74) 0.0001638201
> > (75) 0.000140324
> > (76) 0.0001222632
> > (77) 0.0001034403
> > (78) 8.610252e-05
> > (79) 7.034421e-05
> > (80) 5.381693e-05
> > (81) 4.343148e-05
> > (82) 3.24916e-05
> > (83) 2..498029e-05
> > (84) 2.120895e-05
> > (85) 1.719415e-05
> > (86) 1.67442e-05
> > (87) 1.626488e-05
> > (88) 1.481e-05
> > (89) 1.224572e-05
> > (90) 9.436901e-06
> >
> > Variable: HGT_sfc (subsection)
> > Type: float
> > Total Size: 4 bytes
> > 1 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [1]
> > Coordinates:
> > Number Of Attributes: 14
> > initial_time : 07/07/2004 (15:00)
> > forecast_time : 0
> > level : 1
> > model : North American Regional Reanalysis (NARR)
> > parameter_number : 7
> > parameter_table_version : 131
> > gds_grid_type : 3
> > level_indicator : 109
> > coordinates : g3_lat_0 g3_lon_1
> > _FillValue : -999
> > units : gpm
> > long_name : Geopotential height
> > center : US National Weather Service - NCEP (WMC)
> > sub_center : The North American Regional Reanalysis
> (NARR) Project
> > (0) 115.13
> >
> > Variable: P_sfc (subsection)
> > Type: float
> > Total Size: 4 bytes
> > 1 values
> > Number of Dimensions: 1
> > Dimensions and sizes: [1]
> > Coordinates:
> > Number Of Attributes: 13
> > _FillValue : -999
> > initial_time : 07/07/2004 (15:00)
> > forecast_time : 0
> > model : North American Regional Reanalysis (NARR)
> > parameter_number : 1
> > parameter_table_version : 131
> > gds_grid_type : 3
> > level_indicator : 1
> > coordinates : g3_lat_0 g3_lon_1
> > units : Pa
> > long_name : Pressure
> > center : US National Weather Service - NCEP (WMC)
> > sub_center : The North American Regional Reanalysis
> (NARR) Project
> > (0) 993.0554
> >
>
>________________________________________________________________________________
>
> >
> > CONFIDENTIALITY NOTICE: This
> > electronic message contains information which may be
> legally
> > confidential and/or privileged and does not in any case
> represent a
> > firm ENERGY COMMODITY bid or offer relating thereto which
> binds the
> > sender without an additional express written confirmation
> to that
> > effect. The information is intended solely for the
> individual or
> > entity named above and access by anyone else is
> unauthorized. If
> > you are not the intended recipient, any disclosure,
> copying,
> > distribution, or use of the contents of this information is
> > prohibited and may be unlawful. If you have received this
> > electronic transmission in error, please reply immediately
> to the
> > sender that you have received the message in error, and
> delete it.
> > Thank you.
> >
> >
> >
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Apr 24 2007 - 09:29:52 MDT
This archive was generated by hypermail 2.2.0 : Thu Apr 26 2007 - 08:55:37 MDT