Ignore:
Timestamp:
Jul 24, 2024, 6:46:45 PM (2 months ago)
Author:
abarral
Message:

enforce PRIVATE by default in several modules, expose PUBLIC as needed
move eigen.f90 to obsolete/
(lint) aslong the way

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv2.f90

    r5106 r5119  
    1 
    21! $Header$
    32
    4 SUBROUTINE gradiv2(klevel, xcov, ycov, ld, gdx, gdy )
     3SUBROUTINE gradiv2(klevel, xcov, ycov, ld, gdx, gdy)
    54  !
    65  ! P. Le Van
    … …  
    1716  !
    1817  USE lmdz_filtreg, ONLY: filtreg
     18  USE lmdz_ssum_scopy, ONLY: scopy
     19
    1920  IMPLICIT NONE
    2021  !
    … …  
    2728
    2829  INTEGER :: klevel
    29   REAL :: xcov( ip1jmp1,klevel ), ycov( ip1jm,klevel )
    30   REAL :: gdx( ip1jmp1,klevel ),  gdy( ip1jm,klevel )
     30  REAL :: xcov(ip1jmp1, klevel), ycov(ip1jm, klevel)
     31  REAL :: gdx(ip1jmp1, klevel), gdy(ip1jm, klevel)
    3132  !
    3233  ! ........       variables locales       .........
    3334  !
    34   REAL :: div(ip1jmp1,llm)
     35  REAL :: div(ip1jmp1, llm)
    3536  REAL :: signe, nugrads
    36   INTEGER :: l,ij,iter,ld
     37  INTEGER :: l, ij, iter, ld
    3738
    3839  !    ........................................................
    3940  !
    4041  !
    41   CALL SCOPY( ip1jmp1 * klevel, xcov, 1, gdx, 1 )
    42   CALL SCOPY(   ip1jm * klevel, ycov, 1, gdy, 1 )
     42  CALL SCOPY(ip1jmp1 * klevel, xcov, 1, gdx, 1)
     43  CALL SCOPY(ip1jm * klevel, ycov, 1, gdy, 1)
    4344  !
    4445  !
    45   signe   = (-1.)**ld
     46  signe = (-1.)**ld
    4647  nugrads = signe * cdivu
    4748  !
    4849
     50  CALL    divergf(klevel, gdx, gdy, div)
    4951
    50   CALL    divergf( klevel, gdx,   gdy , div )
     52  IF(ld>1)   THEN
    5153
    52   IF( ld>1 )   THEN
     54    CALL laplacien (klevel, div, div)
    5355
    54     CALL laplacien ( klevel, div,  div     )
     56    !    ......  Iteration de l'operateur laplacien_gam   .......
    5557
    56   !    ......  Iteration de l'operateur laplacien_gam   .......
    57 
    58     DO iter = 1, ld -2
    59      CALL laplacien_gam ( klevel,cuvscvgam1,cvuscugam1,unsair_gam1, &
    60            unsapolnga1, unsapolsga1,  div, div       )
     58    DO iter = 1, ld - 2
     59      CALL laplacien_gam (klevel, cuvscvgam1, cvuscugam1, unsair_gam1, &
     60              unsapolnga1, unsapolsga1, div, div)
    6161    ENDDO
    6262
    6363  ENDIF
    6464
    65 
    66    CALL filtreg( div   , jjp1, klevel, 2, 1, .TRUE., 1 )
    67    CALL  grad  ( klevel,  div,   gdx,  gdy             )
     65  CALL filtreg(div, jjp1, klevel, 2, 1, .TRUE., 1)
     66  CALL  grad  (klevel, div, gdx, gdy)
    6867
    6968  !
    70    DO   l = 1, klevel
    71      DO  ij = 1, ip1jmp1
    72       gdx( ij,l ) = gdx( ij,l ) * nugrads
    73      ENDDO
    74      DO  ij = 1, ip1jm
    75       gdy( ij,l ) = gdy( ij,l ) * nugrads
    76      ENDDO
    77    ENDDO
     69  DO   l = 1, klevel
     70    DO  ij = 1, ip1jmp1
     71      gdx(ij, l) = gdx(ij, l) * nugrads
     72    ENDDO
     73    DO  ij = 1, ip1jm
     74      gdy(ij, l) = gdy(ij, l) * nugrads
     75    ENDDO
     76  ENDDO
    7877  !
    79    RETURN
     78  RETURN
    8079END SUBROUTINE gradiv2
Note: See TracChangeset for help on using the changeset viewer.