source: codes/icosagcm/devel/src/kernels_hex_master/compute_vorticity.k90

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

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