Re: rip_cape_3d Question

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Mon, 23 Apr 2007 14:05:25 -0600 (MDT)

Jacob,

Are the subsetted values the ones that you printed out below? If so,
you still have missing values in them. TEM, HGT, and SPH all have -999
as their first value. It looks like you need to change your subscripts to
be one index value higher.

However, even when I fixed this, there was still an "out-of-bounds"
error. I'm going to get into the Fortran code and see what's going
on. I think there are some limits to the ranges of data that you can
pass in, but these limits have never been very clear to me.

--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 Mon Apr 23 2007 - 14:05:25 MDT

This archive was generated by hypermail 2.2.0 : Thu Apr 26 2007 - 08:55:37 MDT