- Timestamp:
- Aug 2, 2024, 9:58:25 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_gam_loc.f90
r5140 r5159 1 1 SUBROUTINE diverg_gam_loc(klevel,cuvscvgam,cvuscugam,unsairegam, & 2 2 unsapolnga,unsapolsga, x, y, div ) 3 ! 3 4 4 ! P. Le Van 5 ! 5 6 6 ! ********************************************************************* 7 7 ! ... calcule la divergence a tous les niveaux d'1 vecteur de compos. … … 13 13 USE lmdz_comgeom 14 14 15 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 16 USE lmdz_paramet 15 17 IMPLICIT NONE 16 ! 18 17 19 ! x et y sont des arguments d'entree pour le s-prog 18 20 ! div est un argument de sortie pour le s-prog 19 21 20 INCLUDE "dimensions.h" 21 INCLUDE "paramet.h" 22 ! 22 23 24 23 25 ! .......... variables en arguments ................... 24 ! 26 25 27 INTEGER :: klevel 26 28 REAL :: x( ijb_u:ije_u,klevel ),y( ijb_v:ije_v,klevel ) … … 28 30 REAL :: cuvscvgam(ip1jm),cvuscugam(ip1jmp1),unsairegam(ip1jmp1) 29 31 REAL :: unsapolnga,unsapolsga 30 ! 32 31 33 ! ............... variables locales ......................... 32 34 … … 36 38 ! ................................................................... 37 39 INTEGER :: ijb,ije,jjb,jje 38 ! 39 ! 40 41 40 42 ijb=ij_begin 41 43 ije=ij_end … … 45 47 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 46 48 DO l = 1,klevel 47 ! 49 48 50 DO ij = ijb, ije - 1 49 51 div( ij + 1, l ) = ( & … … 52 54 unsairegam( ij+1 ) 53 55 ENDDO 54 ! 56 55 57 ! .... correction pour div( 1,j,l) ...... 56 58 ! .... div(1,j,l)= div(iip1,j,l) .... 57 ! 59 58 60 !DIR$ IVDEP 59 61 DO ij = ijb,ije,iip1 60 62 div( ij,l ) = div( ij + iim,l ) 61 63 ENDDO 62 ! 64 63 65 ! .... calcul aux poles ..... 64 ! 66 65 67 IF (pole_nord) THEN 66 68 DO ij = 1,iim … … 68 70 ENDDO 69 71 sumypn = SSUM ( iim,aiy1,1 ) * unsapolnga 70 ! 72 71 73 DO ij = 1,iip1 72 74 div( ij , l ) = - sumypn … … 79 81 ENDDO 80 82 sumyps = SSUM ( iim,aiy2,1 ) * unsapolsga 81 ! 83 82 84 DO ij = 1,iip1 83 85 div( ij + ip1jm, l ) = sumyps
Note: See TracChangeset
for help on using the changeset viewer.