source: codes/icosagcm/trunk/src/kernels/div.k90 @ 580

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

trunk : upgrading to devel
File size: 868 bytes
Line 
1   !--------------------------------------------------------------------------
2   !---------------------------- div ----------------------------------
3   DO l = ll_begin, ll_end
4      !DIR$ SIMD
5      DO ij=ij_begin_ext, ij_end_ext
6         div_ij=0.
7         div_ij = div_ij + ne_rup*le_de(ij+u_rup)*u(ij+u_rup,l)
8         div_ij = div_ij + ne_lup*le_de(ij+u_lup)*u(ij+u_lup,l)
9         div_ij = div_ij + ne_left*le_de(ij+u_left)*u(ij+u_left,l)
10         div_ij = div_ij + ne_ldown*le_de(ij+u_ldown)*u(ij+u_ldown,l)
11         div_ij = div_ij + ne_rdown*le_de(ij+u_rdown)*u(ij+u_rdown,l)
12         div_ij = div_ij + ne_right*le_de(ij+u_right)*u(ij+u_right,l)
13         divu(ij,l) = div_ij / Ai(ij)
14      END DO
15   END DO
16   !---------------------------- div ----------------------------------
17   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.