Changeset 2626 for trunk/LMDZ.COMMON/libf/filtrez
- Timestamp:
- Feb 23, 2022, 2:22:42 PM (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/filtrez/filtreg.F
r979 r2626 5 5 & griscal ,iter) 6 6 7 USE filtreg_mod 7 USE filtreg_mod, only: matriceun, matricevn, matrinvn, matrinvs 8 USE filtreg_mod, only: matriceus, matricevs 8 9 9 10 IMPLICIT NONE … … 46 47 c 47 48 c 48 #include "dimensions.h" 49 #include "paramet.h" 50 #include "coefils.h" 51 52 INTEGER nlat,nbniv,ifiltre,iter 49 include "dimensions.h" 50 include "paramet.h" 51 include "coefils.h" 52 53 REAL,INTENT(INOUT) :: champ( iip1,nlat,nbniv) 54 INTEGER,INTENT(IN) :: nlat,nbniv,ifiltre,iter 55 INTEGER,INTENT(IN) :: iaire 56 53 57 INTEGER i,j,l,k 54 58 INTEGER iim2,immjm 55 59 INTEGER jdfil1,jdfil2,jffil1,jffil2,jdfil,jffil 56 60 57 REAL champ( iip1,nlat,nbniv)58 61 59 62 REAL eignq(iim,nlat,nbniv), sdd1(iim),sdd2(iim) 60 LOGICAL 61 INTEGER hemisph , iaire63 LOGICAL,INTENT(IN) :: griscal 64 INTEGER hemisph 62 65 63 66 LOGICAL,SAVE :: first=.TRUE. … … 132 135 jffil2 = jjm 133 136 END IF 134 ELSE 137 ELSE ! of IF( griscal ) 135 138 IF( nlat.NE.jjm ) THEN 136 139 PRINT 2222 … … 159 162 jffil2 = jjm 160 163 END IF 161 END IF 164 END IF ! of IF( griscal ) 162 165 163 166 DO hemisph = 1, 2 … … 223 226 END DO 224 227 225 ENDIF 226 227 ELSE 228 ENDIF ! of IF( ifiltre.EQ.-2 ) ... ELSE IF ( griscal ) ... ELSE 229 230 ELSE ! of IF( hemisph. EQ. 1 ) 228 231 229 232 IF( ifiltre. EQ. -2 ) THEN … … 273 276 END DO 274 277 275 ENDIF 276 277 ENDIF 278 ENDIF ! of IF( ifiltre.EQ.-2 ) ... ELSE IF ( griscal ) ... ELSE 279 280 ENDIF ! of IF( hemisph. EQ. 1 ) 278 281 279 282 IF( ifiltre.EQ. 2 ) THEN … … 301 304 END DO 302 305 303 ENDIF 306 ENDIF ! of IF( ifiltre.EQ. 2 ) 304 307 305 308 DO l = 1, nbniv … … 310 313 311 314 312 ENDDO 315 ENDDO ! of DO hemisph = 1, 2 313 316 314 317 1111 FORMAT(//20x,'ERREUR dans le dimensionnement du tableau CHAMP a … … 316 319 2222 FORMAT(//20x,'ERREUR dans le dimensionnement du tableau CHAMP a fi 317 320 & ltrer, sur la grille de V ou de Z'/) 318 RETURN 321 319 322 END
Note: See TracChangeset
for help on using the changeset viewer.