[5246] | 1 | SUBROUTINE covcont_loc (klevel,ucov, vcov, ucont, vcont ) |
---|
| 2 | USE parallel_lmdz |
---|
[5271] | 3 | USE dimensions_mod, ONLY: iim, jjm, llm, ndm |
---|
[5272] | 4 | USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, & |
---|
| 5 | ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm |
---|
[5271] | 6 | IMPLICIT NONE |
---|
[1632] | 7 | |
---|
[5246] | 8 | !======================================================================= |
---|
| 9 | ! |
---|
| 10 | ! Auteur: P. Le Van |
---|
| 11 | ! ------- |
---|
| 12 | ! |
---|
| 13 | ! Objet: |
---|
| 14 | ! ------ |
---|
| 15 | ! |
---|
| 16 | ! ********************************************************************* |
---|
| 17 | ! calcul des compos. contravariantes a partir des comp.covariantes |
---|
| 18 | ! ******************************************************************** |
---|
| 19 | ! |
---|
| 20 | !======================================================================= |
---|
| 21 | INCLUDE "comgeom.h" |
---|
[1632] | 22 | |
---|
[5246] | 23 | INTEGER :: klevel |
---|
| 24 | REAL :: ucov( ijb_u:ije_u,klevel ), vcov( ijb_v:ije_v,klevel ) |
---|
| 25 | REAL :: ucont( ijb_u:ije_u,klevel ), vcont( ijb_v:ije_v,klevel ) |
---|
| 26 | INTEGER :: l,ij |
---|
| 27 | INTEGER :: ijbu,ijbv,ijeu,ijev |
---|
[1632] | 28 | |
---|
| 29 | |
---|
[5246] | 30 | ijbu=ij_begin-iip1 |
---|
| 31 | ijbv=ij_begin-iip1 |
---|
| 32 | ijeu=ij_end+iip1 |
---|
| 33 | ijev=ij_end+iip1 |
---|
[1632] | 34 | |
---|
[5246] | 35 | if (pole_nord) then |
---|
| 36 | ijbu=ij_begin+iip1 |
---|
| 37 | ijbv=ij_begin |
---|
| 38 | endif |
---|
[1632] | 39 | |
---|
[5246] | 40 | if (pole_sud) then |
---|
| 41 | ijeu=ij_end-iip1 |
---|
| 42 | ijev=ij_end-iip1 |
---|
| 43 | endif |
---|
[1632] | 44 | |
---|
[5246] | 45 | !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) |
---|
| 46 | DO l = 1,klevel |
---|
| 47 | |
---|
| 48 | DO ij = ijb_u,ije_u |
---|
| 49 | ucont( ij,l ) = ucov( ij,l ) * unscu2( ij ) |
---|
| 50 | END DO |
---|
| 51 | |
---|
| 52 | DO ij = ijb_v,ije_v |
---|
| 53 | vcont( ij,l ) = vcov( ij,l ) * unscv2( ij ) |
---|
| 54 | END DO |
---|
| 55 | |
---|
| 56 | END DO |
---|
| 57 | !$OMP END DO NOWAIT |
---|
| 58 | RETURN |
---|
| 59 | END SUBROUTINE covcont_loc |
---|