source: codes/icosagcm/devel/src/kernels_hex_master/compute_divergence.k90 @ 1052

Last change on this file since 1052 was 1052, checked in by dubos, 4 years ago

devel : diagnose divergence and vorticity
File size: 879 bytes
Line 
1   !--------------------------------------------------------------------------
2   !---------------------------- compute_divergence ----------------------------------
3   DO l = 1, llm
4      !DIR$ SIMD
5      DO ij=ij_begin_ext, ij_end_ext
6         div_ij = 0.d0
7         div_ij = div_ij + ne_rup*ue(ij+u_rup,l)*le(ij+u_rup)
8         div_ij = div_ij + ne_lup*ue(ij+u_lup,l)*le(ij+u_lup)
9         div_ij = div_ij + ne_left*ue(ij+u_left,l)*le(ij+u_left)
10         div_ij = div_ij + ne_ldown*ue(ij+u_ldown,l)*le(ij+u_ldown)
11         div_ij = div_ij + ne_rdown*ue(ij+u_rdown,l)*le(ij+u_rdown)
12         div_ij = div_ij + ne_right*ue(ij+u_right,l)*le(ij+u_right)
13         div(ij,l) = div_ij / Ai(ij)
14      END DO
15   END DO
16   !---------------------------- compute_divergence ----------------------------------
17   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.