Last change
on this file since 1037 was
906,
checked in by dubos, 5 years ago
|
devel : compute_rhodz for unstructured mesh
|
File size:
1010 bytes
|
Rev | Line | |
---|
[906] | 1 | !-------------------------------------------------------------------------- |
---|
| 2 | !---------------------------- compute_rhodz ---------------------------------- |
---|
| 3 | IF(comp) THEN |
---|
| 4 | !$OMP DO SCHEDULE(STATIC) |
---|
| 5 | DO ij = 1, primal_num |
---|
| 6 | !DIR$ SIMD |
---|
| 7 | DO l = 1, llm |
---|
| 8 | m = ( AP(l,ij)-AP(l+1,ij) + (BP(l,ij)-BP(l+1,ij))*ps(ij) )/g |
---|
| 9 | rhodz(l,ij)=m |
---|
| 10 | END DO |
---|
| 11 | END DO |
---|
| 12 | !$OMP END DO |
---|
| 13 | ELSE |
---|
| 14 | err=0. |
---|
| 15 | !$OMP DO SCHEDULE(STATIC) |
---|
| 16 | DO ij = 1, primal_num |
---|
| 17 | !DIR$ SIMD |
---|
| 18 | DO l = 1, llm |
---|
| 19 | m = ( AP(l,ij)-AP(l+1,ij) + (BP(l,ij)-BP(l+1,ij))*ps(ij) )/g |
---|
| 20 | err = MAX(err, ABS(m-rhodz(l,ij))) |
---|
| 21 | END DO |
---|
| 22 | END DO |
---|
| 23 | !$OMP END DO |
---|
| 24 | IF(err>1e-10) THEN |
---|
| 25 | PRINT *, 'Discrepancy between ps and rhodz detected', err |
---|
| 26 | STOP |
---|
| 27 | END IF |
---|
| 28 | END IF |
---|
| 29 | !---------------------------- compute_rhodz ---------------------------------- |
---|
| 30 | !-------------------------------------------------------------------------- |
---|
Note: See
TracBrowser
for help on using the repository browser.