source: codes/icosagcm/devel/src/kernels_unst/caldyn_vert.k90 @ 624

Last change on this file since 624 was 624, checked in by dubos, 7 years ago

devel : introduced kernels for caldyn_vert (optional)
File size: 2.4 KB
Line 
1   !--------------------------------------------------------------------------
2   !---------------------------- caldyn_vert ----------------------------------
3   DO iq=1,nqdyn
4      !$OMP DO SCHEDULE(STATIC)
5      DO ij = 1, primal_num
6         DO l = 2, llm
7            dtheta_rhodz(l,ij,iq) = dtheta_rhodz(l,ij,iq) + 0.5*(theta(l,ij,iq)+theta(l-1,ij,iq))*wflux(l,ij)
8         END DO
9      END DO
10      !$OMP END DO
11      !$OMP DO SCHEDULE(STATIC)
12      DO ij = 1, primal_num
13         DO l = 1, llm-1
14            dtheta_rhodz(l,ij,iq) = dtheta_rhodz(l,ij,iq) - 0.5*(theta(l,ij,iq)+theta(l+1,ij,iq))*wflux(l+1,ij)
15         END DO
16      END DO
17      !$OMP END DO
18   END DO
19   IF(caldyn_vert_variant == caldyn_vert_cons) THEN
20      ! conservative vertical transport of momentum : (F/m)du/deta = 1/m (d/deta(Fu)-u.dF/deta)
21      !$OMP DO SCHEDULE(STATIC)
22      DO edge = 1, edge_num
23         ij_left = left(edge)
24         ij_right = right(edge)
25         DO l = 2, llm
26            wwuu(l,edge) = .25*(wflux(l,ij_left)+wflux(l,ij_right))*(u(l,edge)+u(l-1,edge)) ! Fu
27         END DO
28      END DO
29      !$OMP END DO
30      ! make sure wwuu is up to date
31      !$OMP BARRIER
32      !$OMP DO SCHEDULE(STATIC)
33      DO edge = 1, edge_num
34         ij_left = left(edge)
35         ij_right = right(edge)
36         DO l = 1, llm
37            dFu_deta = wwuu(l+1,edge)-wwuu(l,edge) ! d/deta (F*u)
38            dF_deta = .5*(wflux(l+1,ij_left)+wflux(l+1,ij_right)-(wflux(l,ij_left)+wflux(l,ij_right))) ! d/deta(F)
39            du(l,edge) = du(l,edge) - (dFu_deta-u(l,edge)*dF_deta) / (.5*(rhodz(l,ij_left)+rhodz(l,ij_right))) ! (F/m)du/deta
40         END DO
41      END DO
42      !$OMP END DO
43   ELSE
44      !$OMP DO SCHEDULE(STATIC)
45      DO edge = 1, edge_num
46         ij_left = left(edge)
47         ij_right = right(edge)
48         DO l = 2, llm
49            wwuu(l,edge) = .5*(wflux(l,ij_left)+wflux(l,ij_right))*(u(l,edge)-u(l-1,edge))
50         END DO
51      END DO
52      !$OMP END DO
53      ! make sure wwuu is up to date
54      !$OMP BARRIER
55      !$OMP DO SCHEDULE(STATIC)
56      DO edge = 1, edge_num
57         ij_left = left(edge)
58         ij_right = right(edge)
59         DO l = 1, llm
60            du(l,edge) = du(l,edge) - (wwuu(l,edge)+wwuu(l+1,edge)) / (rhodz(l,ij_left)+rhodz(l,ij_right))
61         END DO
62      END DO
63      !$OMP END DO
64   END IF
65   !---------------------------- caldyn_vert ----------------------------------
66   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.