source: codes/icosagcm/devel/src/kernels_unst/pvort_only.k90 @ 614

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

devel : added unstructured-mesh kernels
File size: 1.0 KB
Line 
1   !--------------------------------------------------------------------------
2   !---------------------------- pvort_only ----------------------------------
3   !$OMP DO SCHEDULE(STATIC)
4   DO ij = 1, dual_num
5      DO l = 1, llm
6         etav = 0.d0
7         DO iedge = 1, dual_deg(ij)
8            edge = dual_edge(iedge,ij)
9            etav = etav + dual_ne(iedge,ij)*u(l,edge)
10         END DO
11         hv=0.
12         DO ivertex = 1, dual_deg(ij)
13            vertex = dual_vertex(ivertex,ij)
14            hv = hv + Riv2(ivertex,ij)*rhodz(l,vertex)
15         END DO
16         qv(l,ij) = (etav + fv(ij)*Av(ij) )/(hv*Av(ij))
17      END DO
18   END DO
19   !$OMP END DO
20   !$OMP DO SCHEDULE(STATIC)
21   DO edge = 1, edge_num
22      ij_up = up(edge)
23      ij_down = down(edge)
24      DO l = 1, llm
25         qu(l,edge)=0.5d0*(qv(l,ij_down)+qv(l,ij_up))
26      END DO
27   END DO
28   !$OMP END DO
29   !---------------------------- pvort_only ----------------------------------
30   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.