[1632] | 1 | ! |
---|
| 2 | ! $Header$ |
---|
| 3 | ! |
---|
[5246] | 4 | SUBROUTINE covnat_loc(klevel,ucov, vcov, unat, vnat ) |
---|
| 5 | USE parallel_lmdz |
---|
[5271] | 6 | USE dimensions_mod, ONLY: iim, jjm, llm, ndm |
---|
[5272] | 7 | USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, & |
---|
| 8 | ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm |
---|
[5271] | 9 | IMPLICIT NONE |
---|
[1632] | 10 | |
---|
[5246] | 11 | !======================================================================= |
---|
| 12 | ! |
---|
| 13 | ! Auteur: F Hourdin Phu LeVan |
---|
| 14 | ! ------- |
---|
| 15 | ! |
---|
| 16 | ! Objet: |
---|
| 17 | ! ------ |
---|
| 18 | ! |
---|
| 19 | ! ********************************************************************* |
---|
| 20 | ! calcul des compos. naturelles a partir des comp.covariantes |
---|
| 21 | ! ******************************************************************** |
---|
| 22 | ! |
---|
| 23 | !======================================================================= |
---|
| 24 | INCLUDE "comgeom.h" |
---|
[1632] | 25 | |
---|
[5246] | 26 | INTEGER :: klevel |
---|
| 27 | REAL :: ucov( ijb_u:ije_u,klevel ), vcov( ijb_v:ije_v,klevel ) |
---|
| 28 | REAL :: unat( ijb_u:ije_u,klevel ), vnat( ijb_v:ije_v,klevel ) |
---|
| 29 | INTEGER :: l,ij |
---|
| 30 | INTEGER :: ijb,ije |
---|
[1632] | 31 | |
---|
[5246] | 32 | |
---|
| 33 | ijb=ij_begin |
---|
| 34 | ije=ij_end |
---|
| 35 | |
---|
| 36 | if (pole_nord) then |
---|
| 37 | |
---|
| 38 | !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) |
---|
| 39 | DO l = 1,klevel |
---|
| 40 | DO ij = 1, iip1 |
---|
| 41 | unat (ij,l) =0. |
---|
| 42 | END DO |
---|
| 43 | ENDDO |
---|
[1632] | 44 | !$OMP ENDDO NOWAIT |
---|
[5246] | 45 | endif |
---|
[1632] | 46 | |
---|
[5246] | 47 | if (pole_sud) then |
---|
| 48 | !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) |
---|
| 49 | DO l = 1,klevel |
---|
| 50 | DO ij = ip1jm+1, ip1jmp1 |
---|
| 51 | unat (ij,l) =0. |
---|
| 52 | END DO |
---|
| 53 | ENDDO |
---|
[1632] | 54 | !$OMP ENDDO NOWAIT |
---|
[5246] | 55 | endif |
---|
[1632] | 56 | |
---|
[5246] | 57 | ijb=ij_begin |
---|
| 58 | ije=ij_end |
---|
| 59 | if (pole_nord) ijb=ij_begin+iip1 |
---|
| 60 | if (pole_sud) ije=ij_end-iip1 |
---|
| 61 | |
---|
| 62 | !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) |
---|
| 63 | DO l = 1,klevel |
---|
| 64 | DO ij = ijb, ije |
---|
| 65 | unat( ij,l ) = ucov( ij,l ) / cu(ij) |
---|
| 66 | ENDDO |
---|
| 67 | END DO |
---|
[1632] | 68 | !$OMP ENDDO NOWAIT |
---|
| 69 | |
---|
[5246] | 70 | ijb=ij_begin-iip1 |
---|
| 71 | ije=ij_end |
---|
| 72 | if (pole_nord) ijb=ij_begin |
---|
| 73 | if (pole_sud) ije=ij_end-iip1 |
---|
| 74 | |
---|
| 75 | !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) |
---|
| 76 | DO l = 1,klevel |
---|
| 77 | DO ij = ijb,ije |
---|
| 78 | vnat( ij,l ) = vcov( ij,l ) / cv(ij) |
---|
| 79 | ENDDO |
---|
| 80 | ENDDO |
---|
[1632] | 81 | !$OMP ENDDO NOWAIT |
---|
[5246] | 82 | |
---|
| 83 | RETURN |
---|
| 84 | END SUBROUTINE covnat_loc |
---|