--- trunk/libf/filtrez/inifilr.f90 2011/12/06 15:07:04 54 +++ trunk/Sources/filtrez/inifilr.f 2015/04/30 18:35:49 136 @@ -1,17 +1,20 @@ module inifilr_m - use dimens_m, only: iim - IMPLICIT NONE INTEGER jfiltnu, jfiltsu, jfiltnv, jfiltsv - INTEGER, PARAMETER:: nfilun=3, nfilus=2, nfilvn=2, nfilvs=2 - real matriceun(iim,iim,nfilun), matriceus(iim,iim,nfilus) - real matricevn(iim,iim,nfilvn), matricevs(iim,iim,nfilvs) - real matrinvn(iim,iim,nfilun), matrinvs(iim,iim,nfilus) + ! North: + real, allocatable:: matriceun(:, :, :), matrinvn(:, :, :) + ! (iim, iim, 2:jfiltnu) + + real, allocatable:: matricevn(:, :, :) ! (iim, iim, jfiltnv) - private iim, nfilun, nfilus, nfilvn, nfilvs + ! South: + real, allocatable:: matriceus(:, :, :), matrinvs(:, :, :) + ! (iim, iim, jfiltsu:jjm) + + real, allocatable:: matricevs(:, :, :) ! (iim, iim, jfiltsv:jjm) contains @@ -34,32 +37,27 @@ ! The modes are filtered from modfrst to modemax. - USE dimens_m, ONLY : iim, jjm - USE logic, ONLY : fxyhypb, ysinus - USE comgeom, ONLY : rlatu, rlatv, xprimu - use nr_util, only: pi - USE serre, ONLY : alphax USE coefils, ONLY : coefilu, coefilu2, coefilv, coefilv2, eignfnu, & eignfnv, modfrstu, modfrstv + USE comgeom, ONLY : rlatu, rlatv, xprimu + USE dimens_m, ONLY : iim, jjm + use inifgn_m, only: inifgn + use nr_util, only: pi + USE serre, ONLY : grossismx ! Local: - REAL dlonu(iim), dlatu(jjm) + REAL dlatu(jjm) REAL rlamda(2: iim), eignvl(iim) REAL lamdamax, cof INTEGER i, j, modemax, imx, k, kf - REAL dymin, dxmin, colat0 + REAL dymin, colat0 REAL eignft(iim, iim), coff - EXTERNAL inifgn !----------------------------------------------------------- print *, "Call sequence information: inifilr" - DO i = 1, iim - dlonu(i) = xprimu(i) - END DO - CALL inifgn(eignvl) PRINT *, 'EIGNVL ' @@ -79,32 +77,16 @@ dlatu(j) = rlatu(j) - rlatu(j+1) END DO - dxmin = dlonu(1) - DO i = 2, iim - dxmin = min(dxmin, dlonu(i)) - END DO dymin = dlatu(1) DO j = 2, jjm dymin = min(dymin, dlatu(j)) END DO - colat0 = min(0.5, dymin/dxmin) - - IF (.NOT. fxyhypb .AND. ysinus) THEN - colat0 = 0.6 - ! À revoir pour ysinus - alphax = 0. - END IF + colat0 = min(0.5, dymin / minval(xprimu(:iim))) PRINT *, 'colat0 = ', colat0 - PRINT *, 'alphax = ', alphax - IF (alphax == 1.) THEN - PRINT *, 'alphax doit etre < a 1. Corriger ' - STOP 1 - END IF - - lamdamax = iim / (pi * colat0 * (1. - alphax)) + lamdamax = iim / (pi * colat0 / grossismx) rlamda = lamdamax / sqrt(abs(eignvl(2: iim))) DO j = 1, jjm @@ -255,78 +237,10 @@ PRINT *, 'Modes premiers u ' PRINT 334, modfrstu - IF (nfilun < jfiltnu) THEN - PRINT *, 'le parametre nfilun utilise pour la matrice ', & - 'matriceun est trop petit ! ' - PRINT *, 'Le changer dans parafilt.h et le mettre a ', jfiltnu - PRINT *, 'Pour information, nfilun, nfilus, nfilvn, nfilvs ', & - 'doivent etre egaux successivement a ', jfiltnu, & - jjm - jfiltsu + 1, jfiltnv, jjm - jfiltsv + 1 - STOP 1 - END IF - IF (nfilun > jfiltnu+2) THEN - PRINT *, 'le parametre nfilun utilise pour la matrice ', & - 'matriceun est trop grand ! Gachis de memoire ! ' - PRINT *, 'Le changer dans parafilt.h et le mettre a ', jfiltnu - PRINT *, 'Pour information, nfilun, nfilus, nfilvn, nfilvs ', & - 'doivent etre egaux successivement a ', & - jfiltnu, jjm - jfiltsu + 1, jfiltnv, jjm - jfiltsv + 1 - END IF - IF (nfilus < jjm-jfiltsu+1) THEN - PRINT *, 'le parametre nfilus utilise pour la matrice ', & - 'matriceus est trop petit ! ' - PRINT *, 'Le changer dans parafilt.h et le mettre a ', & - jjm - jfiltsu + 1 - PRINT *, 'Pour information , nfilun, nfilus, nfilvn, nfilvs ', & - 'doivent etre egaux successivement a ', & - jfiltnu, jjm - jfiltsu + 1, jfiltnv, jjm - jfiltsv + 1 - STOP 1 - END IF - IF (nfilus > jjm-jfiltsu+3) THEN - PRINT *, 'le parametre nfilus utilise pour la matrice ', & - 'matriceus est trop grand ! ' - PRINT *, 'Le changer dans parafilt.h et le mettre a ', & - jjm - jfiltsu + 1 - PRINT *, 'Pour information , nfilun, nfilus, nfilvn, nfilvs ', & - 'doivent etre egaux successivement a ', & - jfiltnu, jjm - jfiltsu + 1, jfiltnv, jjm - jfiltsv + 1 - END IF - IF (nfilvn < jfiltnv) THEN - PRINT *, 'le parametre nfilvn utilise pour la matrice ', & - 'matricevn est trop petit ! ' - PRINT *, 'Le changer dans parafilt.h et le mettre a ', jfiltnv - PRINT *, 'Pour information , nfilun, nfilus, nfilvn, nfilvs ', & - 'doivent etre egaux successivement a ', & - jfiltnu, jjm - jfiltsu + 1, jfiltnv, jjm - jfiltsv + 1 - STOP 1 - END IF - IF (nfilvn > jfiltnv+2) THEN - PRINT *, 'le parametre nfilvn utilise pour la matrice ', & - 'matricevn est trop grand ! Gachis de memoire ! ' - PRINT *, 'Le changer dans parafilt.h et le mettre a ', jfiltnv - PRINT *, 'Pour information , nfilun, nfilus, nfilvn, nfilvs ', & - 'doivent etre egaux successivement a ', & - jfiltnu, jjm - jfiltsu + 1, jfiltnv, jjm - jfiltsv + 1 - END IF - IF (nfilvs < jjm-jfiltsv+1) THEN - PRINT *, 'le parametre nfilvs utilise pour la matrice ', & - 'matricevs est trop petit ! Le changer dans parafilt.h ' - PRINT *, 'Le changer dans parafilt.h et le mettre a ', & - jjm - jfiltsv + 1 - PRINT *, 'Pour information , nfilun, nfilus, nfilvn, nfilvs ', & - 'doivent etre egaux successivement a ', & - jfiltnu, jjm - jfiltsu + 1, jfiltnv, jjm - jfiltsv + 1 - STOP 1 - END IF - IF (nfilvs > jjm-jfiltsv+3) THEN - PRINT *, 'le parametre nfilvs utilise pour la matrice ', & - 'matricevs est trop grand ! Gachis de memoire ! ' - PRINT *, 'Le changer dans parafilt.h et le mettre a ', & - jjm - jfiltsv + 1 - PRINT *, 'Pour information , nfilun, nfilus, nfilvn, nfilvs ', & - 'doivent etre egaux successivement a ', & - jfiltnu, jjm - jfiltsu + 1, jfiltnv, jjm - jfiltsv + 1 - END IF + allocate(matriceun(iim, iim, 2:jfiltnu), matrinvn(iim, iim, 2:jfiltnu)) + allocate(matricevn(iim, iim, jfiltnv)) + allocate(matricevs(iim, iim, jfiltsv:jjm)) + allocate(matriceus(iim, iim, jfiltsu:jjm), matrinvs(iim, iim, jfiltsu:jjm)) ! Calcul de la matrice filtre 'matriceu' pour les champs situes ! sur la grille scalaire @@ -352,7 +266,7 @@ end IF eignft(i, :) = eignfnv(:, i) * coff END DO - matriceus(:, :, j - jfiltsu + 1) = matmul(eignfnv, eignft) + matriceus(:, :, j) = matmul(eignfnv, eignft) END DO ! Calcul de la matrice filtre 'matricev' pour les champs situes @@ -379,7 +293,7 @@ end IF eignft(i, :) = eignfnu(:, i)*coff END DO - matricevs(:, :, j-jfiltsv+1) = matmul(eignfnu, eignft) + matricevs(:, :, j) = matmul(eignfnu, eignft) END DO ! Calcul de la matrice filtre 'matrinv' pour les champs situes @@ -406,7 +320,7 @@ end IF eignft(i, :) = eignfnv(:, i)*coff END DO - matrinvs(:, :, j-jfiltsu+1) = matmul(eignfnv, eignft) + matrinvs(:, :, j) = matmul(eignfnv, eignft) END DO 334 FORMAT (1X, 24I3)