# Re: Small/Medium Circles..

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Thu Sep 29 2011 - 01:25:52 MDT

[1]
PLEASE look at the printed output. It should be for *your*
information as part of *your* debugging efforts. It should
not be done and just sent off to ncl-talk. Examine this
to find problems. Only send if *you* have looked 1st
and can not determine the error.

For example:
print(time(ifs)+" : "+min(slp2d)+" ("+imin(ifs)+","+jmin(ifs)+")")
yields
(0) missing : 1003.26 (29,149)

Ummm, why is time(ifs) missing? Well if you look at the code
you can see that you commented out the line that creates time(ifs)
;time(ifs) = wrf_user_list_times(a)
so, uncomment the line.

;time(ifs)= chartostring( time(ifs)) ; why is this here? it is wrong

[2]
The error message says the error is at line 120.
Well, as part of *your* debugging, add print statements
*before* line 120. Look at
do i=0,ndate-1
You have references to i+1, well when
i=ndate-1
the (i+1) refers to (ndate-1)+1=ndate.
This subscript is out of range. Is it not?.

print("debug: i ="+i)
print("debug: imin(i)="+imin(i)+" jmin(i)="+jmin(i))
print("debug: i+1="+(i+1))
print("debug: imin(i+1)="+imin(i+1)+" jmin(i+1)="+jmin(i+1))

xx=(/lon2d(imin(i),jmin(i)),lon2d(imin(i+1),jmin(i+1))/)

so is (i+1) the problem? If so, you may have to change the
loop to
do i = 0,ndate-2

Dees the subscript exceed the bounds of lon2d???

Variable: lon2d
Dimensions and sizes: [south_north | 117] x [west_east | 150]

if imin(i)>116 or jmin(i)>149

or

if imin(i+1)>116 or jmin(i+1)>149

you have a subscript error.

Good luck
On 9/28/11 11:36 PM, ramesh karanam wrote:
> Sir,
>
> I am runing WRF model for singal domain case with 15msX15kms
> resolution. I am having single wrfout file with 127 time steps.I am
> getting the following error when I am using the attached script. iam
> getting the following error.
>
> ***********************************************************
> linux:~/rmcncl/wind # ncl WRF_track_1.ncl
> University Corporation for Atmospheric Research
> NCAR Command Language Version 5.2.1
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
>
>
> Variable: nfs
> Type: integer
> Total Size: 4 bytes
> 1 values
> Number of Dimensions: 1
> Dimensions and sizes: [1]
> Coordinates:
> (0) 1
> (0) Check input data:1 .ne. 8
>
>
> Variable: Time
> Type: double
> Total Size: 1016 bytes
> 127 values
> Number of Dimensions: 1
> Dimensions and sizes: [Time | 127]
> Coordinates:
> Time: [962856..962982]
> Number Of Attributes: 4
> units : hours since 1901-01-01 00:00:00
> calendar : standard
> long_name : Time
> description : Time
> (0) 962856
> (1) 962857
> (2) 962858
> (3) 962859
> (4) 962860
> (5) 962861
> (6) 962862
> (7) 962863
> (8) 962864
> (9) 962865
> (10) 962866
> (11) 962867
> (12) 962868
> (13) 962869
> (14) 962870
> (15) 962871
> (16) 962872
> (17) 962873
> (18) 962874
> (19) 962875
> (20) 962876
> (21) 962877
> (22) 962878
> (23) 962879
> (24) 962880
> (25) 962881
> (26) 962882
> (27) 962883
> (28) 962884
> (29) 962885
> (30) 962886
> (31) 962887
> (32) 962888
> (33) 962889
> (34) 962890
> (35) 962891
> (36) 962892
> (37) 962893
> (38) 962894
> (39) 962895
> (40) 962896
> (41) 962897
> (42) 962898
> (43) 962899
> (44) 962900
> (45) 962901
> (46) 962902
> (47) 962903
> (48) 962904
> (49) 962905
> (50) 962906
> (51) 962907
> (52) 962908
> (53) 962909
> (54) 962910
> (55) 962911
> (56) 962912
> (57) 962913
> (58) 962914
> (59) 962915
> (60) 962916
> (61) 962917
> (62) 962918
> (63) 962919
> (64) 962920
> (65) 962921
> (66) 962922
> (67) 962923
> (68) 962924
> (69) 962925
> (70) 962926
> (71) 962927
> (72) 962928
> (73) 962929
> (74) 962930
> (75) 962931
> (76) 962932
> (77) 962933
> (78) 962934
> (79) 962935
> (80) 962936
> (81) 962937
> (82) 962938
> (83) 962939
> (84) 962940
> (85) 962941
> (86) 962942
> (87) 962943
> (88) 962944
> (89) 962945
> (90) 962946
> (91) 962947
> (92) 962948
> (93) 962949
> (94) 962950
> (95) 962951
> (96) 962952
> (97) 962953
> (98) 962954
> (99) 962955
> (100) 962956
> (101) 962957
> (102) 962958
> (103) 962959
> (104) 962960
> (105) 962961
> (106) 962962
> (107) 962963
> (108) 962964
> (109) 962965
> (110) 962966
> (111) 962967
> (112) 962968
> (113) 962969
> (114) 962970
> (115) 962971
> (116) 962972
> (117) 962973
> (118) 962974
> (119) 962975
> (120) 962976
> (121) 962977
> (122) 962978
> (123) 962979
> (124) 962980
> (125) 962981
> (126) 962982
>
>
> Variable: imin (subsection)
> Type: integer
> Total Size: 4 bytes
> 1 values
> Number of Dimensions: 1
> Dimensions and sizes: [1]
> Coordinates:
> Number Of Attributes: 1
> _FillValue : -999
> (0) 29
>
>
> Variable: jmin (subsection)
> Type: integer
> Total Size: 4 bytes
> 1 values
> Number of Dimensions: 1
> Dimensions and sizes: [1]
> Coordinates:
> Number Of Attributes: 1
> _FillValue : -999
> (0) 149
> (0) missing : 1003.26 (29,149)
> warning:gsnAddCyclic is not a valid resource in map at this time
>
>
> Variable: lon2d
> Type: float
> Total Size: 70200 bytes
> 17550 values
> Number of Dimensions: 2
> Dimensions and sizes: [south_north | 117] x [west_east | 150]
> Coordinates:
> Number Of Attributes: 5
> FieldType : 104
> MemoryOrder : XY
> description : LONGITUDE, WEST IS NEGATIVE
> units : degree_east
> stagger :
>
>
> Variable: ndate
> Type: integer
> Total Size: 4 bytes
> 1 values
> Number of Dimensions: 1
> Dimensions and sizes: [1]
> Coordinates:
> fatal:Subscript out of range, error in subscript #0
> fatal:An error occurred reading lon2d
> fatal:Execute: Error occurred at or near line 120 in file WRF_track_1.ncl
> *********************************************************************************
>
> I am new to NCL, kindly help me to debug it or if don't mind kindly
> provide new script whic is suitable for my case.
>
>
> With regards
> Ramesh k
> On 9/29/11, David Brown<dbrown@ucar.edu> wrote:
>>
>> I see that Adam has sent an answer concerning how to get a realistic size
>> for each filled circle. I have just created an example that has a very crude
>> representation of a possible hurricane track -- 3 locations only. It
>> modifies the example polyg_8.ncl. Your situation is actually simpler than
>> polyg_8, because you basically just need to draw one circle at a time. So
>> you can see that mostly I just removed code from the original example.
>> -dave
>>
>>
>>
>> On Sep 28, 2011, at 7:48 AM, Erick Rivera Lemus wrote:
>>
>>> Hi David,
>>>
>>> Yes I see example 8. But the data that I have is
>>> lat lon and radius of isobar. The position changes as well as the radius.
>>> I am not sure how to use the example since we are trying to show the
>>> actual radius of the storm moving along the northeast. I want to show the
>>> size of the storm. Do you have any suggestions on how to do it using the
>>> example you sent me?
>>>
>>> Thanks.
>>>
>>>
>>> Cheers,
>>> Erick Rafael Rivera Lemus
>>> Software Technician
>>> Institute of Marine and Coastal Sciences
>>> Coastal Ocean Observation Laboratory
>>>
>>> Rutgers, The State University of New Jersey
>>> 71 Dudley Road New Brunswick, NJ 08901
>>>
>>> Office: 848-932-3261
>>> Fax: 732-932-8578
>>> Cell: 617-447-5106
>>>
>>>
>>>
>>>
>>> <OI_MARA_2011_08_28_0300.jpg>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> On Sep 27, 2011, at 4:53 PM, David Brown wrote:
>>>
>>>>
>>>> For filled circles you can use poly marker index 16, or for an empty
>>>> circle you can use index 4. You can control the size using
>>>> gsMarkerSizeF. See http://www.ncl.ucar.edu/Applications/polyg.shtml,
>>>> example 8.
>>>> -dave
>>>>
>>>>
>>>> On Sep 27, 2011, at 2:27 PM, Erick Rivera Lemus wrote:
>>>>
>>>>> Hello,
>>>>>
>>>>> I Have a question about circles. Is there a function on NCL that is able
>>>>> to create circles based upon the lat, lon and radius that you specify?
>>>>>
>>>>> In my example I have a map and I have a marker with constant width for
>>>>> each hour. I have data that represent the range in width due to pressure
>>>>> and winds. So for example when is closer to the water I expect the width
>>>>> of the marker to be bigger and smaller as it approaches closer to land.
>>>>> Is there any function that shapes markers according to lat lon and
>>>>>
>>>>> I would like to plot small or medium circles on a map that would give us
>>>>> a better representation of Hurricane Irene as it made his travel along
>>>>> the NE.
>>>>>
>>>>>
>>>>> PS. I sent an image to previous email but It was hold by moderator due
>>>>> to image size.
>>>>>
>>>>>
>>>>> Cheers,
>>>>> Erick Rafael Rivera Lemus
>>>>> Software Technician
>>>>> Institute of Marine and Coastal Sciences
>>>>> Coastal Ocean Observation Laboratory
>>>>>
>>>>> Rutgers, The State University of New Jersey
>>>>> 71 Dudley Road New Brunswick, NJ 08901
>>>>>
>>>>> Office: 848-932-3261
>>>>> Fax: 732-932-8578
>>>>> Cell: 617-447-5106
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>
>>
>>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 29 01:25:58 2011

This archive was generated by hypermail 2.1.8 : Sun Oct 09 2011 - 13:05:26 MDT