-
Notifications
You must be signed in to change notification settings - Fork 1
/
drive_ana.pro
47 lines (30 loc) · 855 Bytes
/
drive_ana.pro
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
pro drive_ana,xtag=xtag,tauxy=tauxy,lon=lon,lat=lat
start=[2010L,10L,1,0]
final=[2010L,11L,1,0]
date=start
stopcon=0
dout$ = '/project/amp/juliob/OGW_UnitTest_output/beta07/'
;xtag='AniOutput'
; double check
itime=0
while stopcon eq 0 do begin
y$=padstring( date(0) )
m$=padstring( date(1) , /e1 )
d$=padstring( date(2) , /e1 )
s$=padstring( date(3) , /e4 )
out$ = dout$+xtag+'.'+y$+'-'+m$+'-'+d$+'-'+s$
rdgw3d,fn=out$,tauf=tauf,utgw=utgw,vtgw=vtgw,lon=lon,lat=lat
print,' Got ',out$
if itime eq 0 then begin
s=size(tauf)
tauxy = dblarr( s(1),s(2), 31L*4 )
endif
tauxy(*,*,itime) = tauf(*,*,46,0)
date = bump_time_vec(date,dtvc=[0,0,0,6*3600L],/sssss)
xx=total( abs(date-final ) )
if xx eq 0 then stopcon=1
itime=itime+1
endwhile
STOP
return
end