source: codes/icosagcm/devel/src/kernels_hex/compute_rhodz.k90 @ 1007

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

devel : compute_rhodz for unstructured mesh
File size: 938 bytes
Line 
1   !--------------------------------------------------------------------------
2   !---------------------------- compute_rhodz ----------------------------------
3   IF(comp) THEN
4      DO l = ll_begin, ll_end
5         !DIR$ SIMD
6         DO ij=ij_begin, ij_end
7            m = ( AP(ij,l)-AP(ij,l+1) + (BP(ij,l)-BP(ij,l+1))*ps(ij) )/g
8            rhodz(ij,l)=m
9         END DO
10      END DO
11   ELSE
12      err=0.
13      DO l = ll_begin, ll_end
14         !DIR$ SIMD
15         DO ij=ij_begin_ext, ij_end_ext
16            m = ( AP(ij,l)-AP(ij,l+1) + (BP(ij,l)-BP(ij,l+1))*ps(ij) )/g
17            err = MAX(err, ABS(m-rhodz(ij,l)))
18         END DO
19      END DO
20      IF(err>1e-10) THEN
21         PRINT *, 'Discrepancy between ps and rhodz detected', err
22         STOP
23      END IF
24   END IF
25   !---------------------------- compute_rhodz ----------------------------------
26   !--------------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.