-
Notifications
You must be signed in to change notification settings - Fork 2
/
stat_2Dinit.f90
60 lines (47 loc) · 1.15 KB
/
stat_2Dinit.f90
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
! Initialize 2D output
subroutine stat_2Dinit(ResetStorage)
use vars
implicit none
integer :: ResetStorage
prec_xy(:,:) = 0.
shf_xy(:,:) = 0.
lhf_xy(:,:) = 0.
lwnt_xy(:,:) = 0.
swnt_xy(:,:) = 0.
lwntc_xy(:,:) = 0.
swntc_xy(:,:) = 0.
pw_xy(:,:) = 0.
cw_xy(:,:) = 0.
iw_xy(:,:) = 0.
cld_xy(:,:) = 0.
u200_xy(:,:) = 0.
v200_xy(:,:) = 0.
usfc_xy(:,:) = 0.
vsfc_xy(:,:) = 0.
w500_xy = 0.
lwns_xy(:,:) = 0.
swns_xy(:,:) = 0.
solin_xy(:,:) = 0.
lwnsc_xy(:,:) = 0.
swnsc_xy(:,:) = 0.
qocean_xy(:,:) = 0.
!===================================
! UW ADDITIONS: MOSTLY 2D STATISTICS
if(ResetStorage.eq.1) then
!bloss: store initial profiles for computation of storage terms in budgets
ustor(:) = u0(1:nzm)
vstor(:) = v0(1:nzm)
tstor(:) = t0(1:nzm)
qstor(:) = q0(1:nzm)
utendcor(:) = 0.
vtendcor(:) = 0.
end if
psfc_xy(:,:) = 0.
u850_xy(:,:) = 0.
v850_xy(:,:) = 0.
swvp_xy(:,:) = 0.
cloudtopheight(:,:) = 0.
echotopheight(:,:) = 0.
cloudtoptemp(:,:) = 0.
! END UW ADDITIONS
end