-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathTCs_2a_create_gefs_tracks.ncl
executable file
·116 lines (86 loc) · 2.9 KB
/
TCs_2a_create_gefs_tracks.ncl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
;;;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;a
;;
;; Creates GEFS TC track .csv files from a-deck .dat file
;;
;;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;;
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/ut_string.ncl"
load "./get_atcf_output.ncl"
;************************************************************************
begin
print(" ")
print("**************")
print(" START SCRIPT ")
print("**************")
print(" ")
int_time = "2024070812"
TCname = "beryl"
file_load = "aal022024" ; do not include the .dat at the end
;-------------------------------------------------------------------------
dir_out = "./"
filename = "gefs_"+int_time+"_"+TCname+".csv"
filename_mean = "gefs_det_mean_"+int_time+"_"+TCname+".csv"
print("int_time: "+int_time)
line_num = 0
line_num2 = 0
outlines := new((/3000/),"string")
outlines2 := new((/3000/),"string")
do ens = 0,32,1 ; loop through all ensemble members
if (ens .eq. 0) then
mem = "AC00"
end if
if (ens .ge. 1 .and. ens .le. 10) then
mem = "AP0"+ens
end if
if (ens .ge. 10 .and. ens .le. 30) then
mem = "AP"+ens
end if
if (ens .eq. 31) then
mem = "AEMN"
end if
if (ens .eq. 32) then
mem = "AVNO"
end if
print("mem: "+mem)
atcf_output := get_atcf(file_load,mem,int_time,True)
; printVarSummary(atcf_output)
;print(atcf_output(0,:))
do i = 0,dimsizes(atcf_output(0,:))-1 ; loop through all times in this ensemble member
time = atcf_output(0,i)
time@units = "hours since 1800-01-01 00:00:00"
utc_date = cd_calendar(time, 0)
year = tointeger(utc_date(:,0)) ; Convert to integer for
month = tointeger(utc_date(:,1)) ; use sprinti
day = tointeger(utc_date(:,2))
hour = tointeger(utc_date(:,3))
mslp = atcf_output(1,i)
lat = atcf_output(2,i)
lon = atcf_output(3,i)
wind = atcf_output(4,i)
ensemble = ens
line = year+", "+month+", "+day+", "+hour+", "+mslp+", "+lat+", "+lon+", "+wind+", "+ensemble
; print("line: "+line)
if ens .le. 30 then
print("line_num: "+line_num)
outlines(line_num) = line
line_num = line_num + 1
end if
if ens .ge. 31 then
; print("line_num2: "+line_num2)
outlines2(line_num2) = line
line_num2 = line_num2 + 1
end if
end do
end do
line_num_final = line_num2
print("line_num_final: "+line_num_final)
asciiwrite(dir_out+filename, outlines)
asciiwrite(dir_out+filename_mean, outlines2)
print(" ")
print("**************")
print(" END SCRIPT ")
print("**************")
print(" ")
end