Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/filtrez/filtreg_mod.F90

Timestamp:

Jul 22, 2024, 9:29:09 PM (4 months ago)

Author:

abarral

Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:

• LMDZ6/branches/Amaury_dev/libf/filtrez/filtreg_mod.F90 (modified) (28 diffs)

LMDZ6/branches/Amaury_dev/libf/filtrez/filtreg_mod.F90

@@ -1 @@
-!
+
 ! $Id$
-!
+
 MODULE filtreg_mod
 
@@ -19 +19 @@
 
     !    ... H. Upadhyaya, O.Sharma   ...
-    !
+
     IMPLICIT NONE
-    !
+
     !     version 3 .....
 
@@ -34 +34 @@
     REAL  dlonu(iim),dlatu(jjm)
     REAL  rlamda( iim ),  eignvl( iim )
-    !
 
     REAL    lamdamax,pi,cof
@@ -47 +46 @@
     INTEGER iymin 
     INTEGER ixmineq
-    !
+
     ! ------------------------------------------------------------
     !   This routine computes the eigenfunctions of the laplacien
     !   on the stretched grid, and the filtering coefficients
-    !      
+
     !  We designate:
     !   eignfn   eigenfunctions of the discrete laplacien
@@ -61 +60 @@
     !   coefil   filtering coefficients ( lamda_max*COS(rlat)/lamda )
     !   sdd      SQRT( dx )
-    !      
+
     !     the modes are filtered from modfrst to modemax
-    !      
+
     !-----------------------------------------------------------
-    !
+
     if ( iim == 1 ) return ! No filtre in 2D y-z
 
@@ -73 +72 @@
        dlonu(i) = xprimu( i )
     ENDDO
-    !
+
     CALL inifgn(eignvl)
-    !
+
     PRINT *,'inifilr: EIGNVL '
     PRINT 250,eignvl
 250 FORMAT( 1x,5e14.6)
-    !
+
     ! compute eigenvalues and eigenfunctions
-    !
-    !
+
+
     !.................................................................
-    !
+
     !  compute the filtering coefficients for scalar lines and 
     !  meridional wind v-lines
-    !
+
     !  we filter all those latitude lines WHERE coefil < 1
     !  NO FILTERING AT POLES
-    !
+
     !  colat0 is to be used  when alpha (stretching coefficient)
     !  is set equal to zero for the regular grid CASE 
-    !
+
     !    .......   Calcul  de  colat0   .........
     !     .....  colat0 = minimum de ( 0.5, min dy/ min dx )   ...
-    !
-    !
+
+
     DO j = 1,jjm
        dlatu( j ) = rlatu( j ) - rlatu( j+1 )
     ENDDO
-    !
 
     dxmin   =  dlonu(1)
@@ -111 +109 @@
        dymin = MIN( dymin,dlatu(j) )
     ENDDO
-    !
+
     ! For a regular grid, we want the filter to start at latitudes
     ! corresponding to lengths dx of the same size as dy (in terms
@@ -118 +116 @@
     ! Same idea for the zoomed grid: start filtering polewards as soon
     ! as length dx becomes of the same size as dy 
-    !
+
     ! if maxlatfilter >0, prescribe the colat0 value from the .def files
    
@@ -125 +123 @@
     colat0  =  MIN( 0.5, dymin/dxmin )
     ! colat0  =  1.
-    !
+
     IF( .NOT.fxyhypb.AND.ysinus )  THEN
        colat0 = 0.6
@@ -137 +135 @@
 
     ENDIF
-    
-    
-    
-    !
+
     PRINT 50, colat0,alphax
 50  FORMAT(/15x,' Inifilr colat0 alphax ',2e16.7)
-    !
+
     IF(alphax==1. )  THEN
        PRINT *,' Inifilr  alphax doit etre  <  a 1.  Corriger '
        STOP
     ENDIF
-    !
+
     lamdamax = iim / ( pi * colat0 * ( 1. - alphax ) )
 
     !                        ... Correction  le 28/10/97  ( P.Le Van ) ..
-    !
+
     DO i = 2,iim
        rlamda( i ) = lamdamax/ SQRT( ABS( eignvl(i) ) )
     ENDDO
-    !
 
     DO j = 1,jjm
@@ -167 +161 @@
     ENDDO
 
-    !
     !    ... Determination de jfiltnu,jfiltnv,jfiltsu,jfiltsv ....
     !    .........................................................
-    !
+
     modemax = iim
 
@@ -176 +169 @@
 
     imx  = iim
-    !
+
     PRINT *,'inifilr: TRUNCATION AT ',imx
-    !
+
 ! Ehouarn: set up some defaults
     jfiltnu=2 ! avoid north pole
@@ -200 +193 @@
        ENDIF
     ENDDO
-    !
+
     DO j = 1, jjm/2
        cof = COS( rlatv(j) )/ colat0
@@ -216 +209 @@
        ENDIF
     ENDDO
-    !                                 
 
     IF( jfiltnu> jjm/2 +1 )  THEN
@@ -251 +243 @@
     ENDIF
 
-    !                                 
     !   ... Determination de coefilu,coefilv,n=modfrstu,modfrstv ....
     !................................................................
-    !
-    !
+
+
     DO j = 1,jjm
     !default initialization: all modes are retained (i.e. no filtering)
@@ -261 +252 @@
        modfrstv( j ) = iim
     ENDDO
-    !
+
     DO j = 2,jfiltnu
        DO k = 2,modemax
@@ -269 +260 @@
        GOTO 84
 82     modfrstu( j ) = k
-       !
+
        kf = modfrstu( j )
        DO k = kf , modemax
@@ -278 +269 @@
 84     CONTINUE
     ENDDO
-    !                                 
-    !
+
+
     DO j = 1,jfiltnv
-       !
+
        DO k = 2,modemax
           cof = rlamda(k) * COS( rlatv(j) )
@@ -288 +279 @@
        GOTO 89
 87     modfrstv( j ) = k
-       !
+
        kf = modfrstv( j )
        DO k = kf , modemax
@@ -297 +288 @@
 89     CONTINUE
     ENDDO
-    !
+
     DO j = jfiltsu,jjm
        DO k = 2,modemax
@@ -305 +296 @@
        GOTO 94
 92     modfrstu( j ) = k
-       !
+
        kf = modfrstu( j )
        DO k = kf , modemax
@@ -314 +305 @@
 94     CONTINUE
     ENDDO
-    !                                 
+
     DO j = jfiltsv,jjm
        DO k = 2,modemax
@@ -322 +313 @@
        GOTO 99
 97     modfrstv( j ) = k
-       !
+
        kf = modfrstv( j )
        DO k = kf , modemax
@@ -331 +322 @@
 99     CONTINUE
     ENDDO
-    !
 
     IF(jfiltnv>=jjm/2 .OR. jfiltnu>=jjm/2)THEN
@@ -348 +338 @@
     PRINT 334,modfrstu
 
-    !  
-    !   ...................................................................
-    !
+    !   ...................................................................
+
     !   ... Calcul de la matrice filtre 'matriceu'  pour les champs situes
     !                       sur la grille scalaire                 ........
     !   ...................................................................
-    !
+
     DO j = 2, jfiltnu
 
@@ -410 +399 @@
 
     !   ...................................................................
-    !
+
     !   ... Calcul de la matrice filtre 'matricev'  pour les champs situes
     !                       sur la grille   de V ou de Z           ........
     !   ...................................................................
-    !
+
     DO j = 1, jfiltnv
 
@@ -471 +460 @@
 
     !   ...................................................................
-    !
+
     !   ... Calcul de la matrice filtre 'matrinv'  pour les champs situes
     !              sur la grille scalaire , pour le filtre inverse ........
     !   ...................................................................
-    !
+
     DO j = 2, jfiltnu
 
@@ -539 +528 @@
     !   ...................................................................
 
-    !
 334 FORMAT(1x,24i3)