66 |
CALL filtreg(zh, jjp1, 1, 2, 1, .TRUE., 1) |
CALL filtreg(zh, jjp1, 1, 2, 1, .TRUE., 1) |
67 |
|
|
68 |
CALL minmax(iip1*jjp1, zh, zhmin, zhmax) |
CALL minmax(iip1*jjp1, zh, zhmin, zhmax) |
69 |
|
IF (zhmin >= zhmax) THEN |
70 |
IF (zhmin>=zhmax) THEN |
PRINT *, 'zhmin zhmax', zhmin, zhmax |
71 |
PRINT *, ' Inidissip zh min max ', zhmin, zhmax |
print *, 'Problème générateur aléatoire dans inidissip' |
72 |
STOP 'probleme generateur alleatoire dans inidissip' |
STOP 1 |
73 |
END IF |
END IF |
74 |
|
|
75 |
zllm = abs(zhmax) |
zllm = abs(zhmax) |
100 |
PRINT *, 'calcul des valeurs propres' |
PRINT *, 'calcul des valeurs propres' |
101 |
|
|
102 |
DO ii = 1, 2 |
DO ii = 1, 2 |
|
|
|
103 |
DO ij = 1, ip1jmp1 |
DO ij = 1, ip1jmp1 |
104 |
zu(ij) = ran1(idum) - .5 |
zu(ij) = ran1(idum) - .5 |
105 |
END DO |
END DO |
159 |
crot = zllm**(-1./nitergrot) |
crot = zllm**(-1./nitergrot) |
160 |
END IF |
END IF |
161 |
END IF |
END IF |
|
|
|
162 |
END DO |
END DO |
163 |
|
|
164 |
PRINT *, 'cdivu = ', cdivu |
PRINT *, 'cdivu = ', cdivu |