source: codes/icosagcm/devel/src/kernels_hex/compute_pmid.k90 @ 929

Last change on this file since 929 was 929, checked in by dubos, 5 years ago

devel : added missing kernels
File size: 838 bytes
Line 
1   !--------------------------------------------------------------------------
2   !---------------------------- compute_pmid ----------------------------------
3   IF(ap_bp_present) THEN
4      IF(offset>0) THEN
5         DO l = ll_begin, ll_end
6            !DIR$ SIMD
7            DO ij=ij_begin_ext, ij_end_ext
8               pmid(ij,l) = .5*(AP(ij,l)+AP(ij,l+1) + (BP(ij,l)+BP(ij,l+1)) * ps(ij) )
9            END DO
10         END DO
11      ELSE
12         DO l = ll_begin, ll_end
13            !DIR$ SIMD
14            DO ij=ij_begin, ij_end
15               pmid(ij,l) = .5*(AP(ij,l)+AP(ij,l+1) + (BP(ij,l)+BP(ij,l+1)) * ps(ij) )
16            END DO
17         END DO
18      END IF
19   END IF
20   !---------------------------- compute_pmid ----------------------------------
21   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.