Re: rip_cape_3d Question

From: <Jacob.Klee_at_nyahnyahspammersnyahnyah>
Date: Mon, 23 Apr 2007 11:39:49 -0400
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_sfc(20,18),P_sfc(20,18),False)

-----ncl-talk-bounces@ucar.edu wrote: -----

To: Jacob.Klee@dom.com
From: Mary Haley <haley@ucar.edu>
Sent by: ncl-talk-bounces@ucar.edu
Date: 04/23/2007 11:24AM
cc: ncl-talk@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@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@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Mon Apr 23 2007 - 09:39:49 MDT

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