source: LMDZ6/trunk/libf/dyn3dmem/diverg_p.f90 @ 5374

Last change on this file since 5374 was 5298, checked in by abarral, 7 weeks ago

Turn planete.h comsoil.h into module
Remove obsolete message related to /scratch/ common

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 2.4 KB
RevLine 
[5246]1SUBROUTINE diverg_p(klevel,x,y,div)
2  !
3  ! P. Le Van
4  !
5  !  *********************************************************************
6  !  ... calcule la divergence a tous les niveaux d'1 vecteur de compos.
7  ! x et y...
8  !          x et y  etant des composantes covariantes   ...
9  !  *********************************************************************
[5281]10  USE comgeom_mod_h
[5246]11  USE parallel_lmdz
[5271]12  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
[5285]13USE paramet_mod_h
[5271]14IMPLICIT NONE
[5246]15  !
16  !  x  et  y  sont des arguments  d'entree pour le s-prog
17  !    div      est  un argument  de sortie pour le s-prog
18  !
19  !
[5271]20
[5272]21
[5298]22
[5246]23  !
24  !    ..........          variables en arguments    ...................
25  !
26  INTEGER :: klevel
27  REAL :: x( ip1jmp1,klevel ),y( ip1jm,klevel ),div( ip1jmp1,klevel )
28  INTEGER :: l,ij
29  !
30  !    ...............     variables  locales   .........................
[1632]31
[5246]32  REAL :: aiy1( iip1 ) , aiy2( iip1 )
33  REAL :: sumypn,sumyps
34  INTEGER :: ijb,ije
35  !    ...................................................................
36  !
37  EXTERNAL  SSUM
38  REAL :: SSUM
39  !
40  !
41  ijb=ij_begin
42  ije=ij_end
43  if (pole_nord) ijb=ij_begin+iip1
44  if(pole_sud)  ije=ij_end-iip1
[1632]45
[5246]46!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
47  DO l = 1,klevel
48  !
49    DO  ij = ijb, ije - 1
50     div( ij + 1, l )     = &
51           cvusurcu( ij+1 ) * x( ij+1,l ) - cvusurcu( ij ) * x( ij , l) + &
52           cuvsurcv(ij-iim) * y(ij-iim,l) - cuvsurcv(ij+1) * y(ij+1,l)
53    ENDDO
54  !
55  ! ....  correction pour  div( 1,j,l)  ......
56  ! ....   div(1,j,l)= div(iip1,j,l) ....
57  !
58  !DIR$ IVDEP
59    DO  ij = ijb,ije,iip1
60     div( ij,l ) = div( ij + iim,l )
61    ENDDO
62  !
63  ! ....  calcul  aux poles  .....
64  !
65    if (pole_nord) then
66      DO  ij  = 1,iim
67       aiy1(ij) =    cuvsurcv(    ij       ) * y(     ij     , l )
68      ENDDO
69      sumypn = SSUM ( iim,aiy1,1 ) / apoln
70  !
71      DO  ij = 1,iip1
72       div(     ij    , l ) = - sumypn
73      ENDDO
74    endif
[1632]75
[5246]76   if (pole_sud) then
77      DO  ij  = 1,iim
78       aiy2(ij) =    cuvsurcv( ij+ ip1jmi1 ) * y( ij+ ip1jmi1, l )
79      ENDDO
80      sumyps = SSUM ( iim,aiy2,1 ) / apols
81  !
82      DO  ij = 1,iip1
83       div( ij + ip1jm, l ) =   sumyps
84      ENDDO
85    endif
[1632]86
[5246]87
88  END DO
89!$OMP END DO NOWAIT
90  !
91
92  !cc        CALL filtreg( div, jjp1, klevel, 2, 2, .TRUE., 1 )
93
94  !
95!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
96    DO l = 1, klevel
97       DO ij = ijb,ije
98        div(ij,l) = div(ij,l) * unsaire(ij)
99      ENDDO
100    ENDDO
101!$OMP END DO NOWAIT
102  !
103   RETURN
104END SUBROUTINE diverg_p
Note: See TracBrowser for help on using the repository browser.