source: codes/icosagcm/devel/src/kernels_unst/div.k90 @ 645

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

devel : added unstructured-mesh kernels
File size: 639 bytes
Line 
1   !--------------------------------------------------------------------------
2   !---------------------------- div ----------------------------------
3   !$OMP DO SCHEDULE(STATIC)
4   DO ij = 1, primal_num
5      DO l = 1, llm
6         div_ij=0.
7         DO iedge = 1, primal_deg(ij)
8            edge = primal_edge(iedge,ij)
9            div_ij = div_ij + primal_ne(iedge,ij)*le_de(edge)*u(l,edge)
10         END DO
11         divu(l,ij) = div_ij / Ai(ij)
12      END DO
13   END DO
14   !$OMP END DO
15   !---------------------------- div ----------------------------------
16   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.