Ignore:
Timestamp:
Jul 28, 2024, 4:17:54 PM (8 weeks ago)
Author:
abarral
Message:

Put comgeom.h, comgeom2.h into modules

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r5134 r5136  
    1 
    21! $Header$
    32
    4   SUBROUTINE interpost(q,qppm)
     3SUBROUTINE interpost(q, qppm)
     4  USE lmdz_comgeom2
    55
    6    IMPLICIT NONE
    7 
     6  IMPLICIT NONE
    87
    98  INCLUDE "dimensions.h"
    109  INCLUDE "paramet.h"
    11   INCLUDE "comgeom2.h"
    1210
    1311  ! Arguments
    14   REAL :: q(iip1,jjp1,llm)
    15   REAL :: qppm(iim,jjp1,llm)
     12  REAL :: q(iip1, jjp1, llm)
     13  REAL :: qppm(iim, jjp1, llm)
    1614  ! Local
    17   INTEGER :: l,i,j
     15  INTEGER :: l, i, j
    1816
    1917  ! RE-INVERSION DES NIVEAUX
    … …  
    2220  ! On passe donc des niveaux de Lin à ceux du LMDZ
    2321
    24     do l=1,llm
    25       do j=1,jjp1
    26          do i=1,iim
    27              q(i,j,l)=qppm(i,j,llm-l+1)
    28          enddo
     22  do l = 1, llm
     23    do j = 1, jjp1
     24      do i = 1, iim
     25        q(i, j, l) = qppm(i, j, llm - l + 1)
    2926      enddo
    30      enddo
     27    enddo
     28  enddo
    3129
    3230  ! BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D
    3331
    34      do l=1,llm
    35        do j=1,jjp1
    36         q(iip1,j,l)=q(1,j,l)
    37        enddo
    38      enddo
     32  do l = 1, llm
     33    do j = 1, jjp1
     34      q(iip1, j, l) = q(1, j, l)
     35    enddo
     36  enddo
    3937
    40 
    41    return
     38  return
    4239
    4340END SUBROUTINE interpost
Note: See TracChangeset for help on using the changeset viewer.