Last change
on this file since 928 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.