Ignore:
Timestamp:
Aug 2, 2024, 9:58:25 PM (3 months ago)
Author:
abarral
Message:

Put dimensions.h and paramet.h into modules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_gam_loc.f90

    r5140 r5159  
    11SUBROUTINE diverg_gam_loc(klevel,cuvscvgam,cvuscugam,unsairegam, &
    22        unsapolnga,unsapolsga,  x, y,  div )
    3   !
     3
    44  ! P. Le Van
    5   !
     5
    66  !  *********************************************************************
    77  !  ... calcule la divergence a tous les niveaux d'1 vecteur de compos.
     
    1313  USE lmdz_comgeom
    1414
     15  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
     16  USE lmdz_paramet
    1517  IMPLICIT NONE
    16   !
     18
    1719  !  x  et  y  sont des arguments  d'entree pour le s-prog
    1820  !    div      est  un argument  de sortie pour le s-prog
    1921
    20   INCLUDE "dimensions.h"
    21   INCLUDE "paramet.h"
    22   !
     22
     23
     24
    2325  !    ..........          variables en arguments    ...................
    24   !
     26
    2527  INTEGER :: klevel
    2628  REAL :: x( ijb_u:ije_u,klevel ),y( ijb_v:ije_v,klevel )
     
    2830  REAL :: cuvscvgam(ip1jm),cvuscugam(ip1jmp1),unsairegam(ip1jmp1)
    2931  REAL :: unsapolnga,unsapolsga
    30   !
     32
    3133  !    ...............     variables  locales   .........................
    3234
     
    3638  !    ...................................................................
    3739  INTEGER :: ijb,ije,jjb,jje
    38   !
    39   !
     40
     41
    4042  ijb=ij_begin
    4143  ije=ij_end
     
    4547!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    4648  DO l = 1,klevel
    47   !
     49
    4850    DO  ij = ijb, ije - 1
    4951     div( ij + 1, l )     = ( &
     
    5254           unsairegam( ij+1 )
    5355    ENDDO
    54   !
     56
    5557  ! ....  correction pour  div( 1,j,l)  ......
    5658  ! ....   div(1,j,l)= div(iip1,j,l) ....
    57   !
     59
    5860  !DIR$ IVDEP
    5961    DO  ij = ijb,ije,iip1
    6062     div( ij,l ) = div( ij + iim,l )
    6163    ENDDO
    62   !
     64
    6365  ! ....  calcul  aux poles  .....
    64   !
     66
    6567   IF (pole_nord) THEN
    6668      DO  ij  = 1,iim
     
    6870      ENDDO
    6971      sumypn = SSUM ( iim,aiy1,1 ) * unsapolnga
    70   !
     72
    7173      DO  ij = 1,iip1
    7274       div(     ij    , l ) = - sumypn
     
    7981      ENDDO
    8082      sumyps = SSUM ( iim,aiy2,1 ) * unsapolsga
    81   !
     83
    8284      DO  ij = 1,iip1
    8385       div( ij + ip1jm, l ) =   sumyps
Note: See TracChangeset for help on using the changeset viewer.