/[lmdze]/trunk/Sources/dyn3d/Dissipation/inidissip.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/Dissipation/inidissip.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC revision 107 by guez, Thu Sep 11 15:09:15 2014 UTC
# Line 51  contains Line 51  contains
51      deltap = 1.      deltap = 1.
52      call random_number(zh)      call random_number(zh)
53      zh = zh - 0.5      zh = zh - 0.5
54      CALL filtreg(zh, jjp1, 1, 2, 1, .TRUE.)      CALL filtreg(zh, direct = .true., intensive = .true.)
55    
56      DO l = 1, 50      DO l = 1, 50
57         CALL divgrad2(1, zh, deltap, niterh, divgra, -1.)         CALL divgrad2(1, zh, deltap, niterh, divgra, -1.)
# Line 65  contains Line 65  contains
65      PRINT *, 'Calcul des valeurs propres de gradiv'      PRINT *, 'Calcul des valeurs propres de gradiv'
66      call random_number(zu)      call random_number(zu)
67      zu = zu - 0.5      zu = zu - 0.5
68      CALL filtreg(zu, jjp1, 1, 2, 1, .TRUE.)      CALL filtreg(zu, direct = .true., intensive = .true.)
69      call random_number(zv)      call random_number(zv)
70      zv = zv - 0.5      zv = zv - 0.5
71      CALL filtreg(zv, jjm, 1, 2, 1, .FALSE.)      CALL filtreg(zv, direct = .true., intensive = .true.)
72    
73      DO l = 1, 50      DO l = 1, 50
74         CALL gradiv2(zu, zv, nitergdiv, gx, gy, -1.)         CALL gradiv2(zu, zv, nitergdiv, gx, gy, -1.)
# Line 83  contains Line 83  contains
83      PRINT *, 'Calcul des valeurs propres de nxgrarot'      PRINT *, 'Calcul des valeurs propres de nxgrarot'
84      call random_number(zu)      call random_number(zu)
85      zu = zu - 0.5      zu = zu - 0.5
86      CALL filtreg(zu, jjp1, 1, 2, 1, .TRUE.)      CALL filtreg(zu, direct = .true., intensive = .true.)
87      call random_number(zv)      call random_number(zv)
88      zv = zv - 0.5      zv = zv - 0.5
89      CALL filtreg(zv, jjm, 1, 2, 1, .FALSE.)      CALL filtreg(zv, direct = .true., intensive = .true.)
90    
91      DO l = 1, 50      DO l = 1, 50
92         CALL nxgraro2(zu, zv, nitergrot, gx, gy, -1.)         CALL nxgraro2(zu, zv, nitergrot, gx, gy, -1.)
Legend:
Removed from v.82  
changed lines
  Added in v.107
  ViewVC Help
Powered by ViewVC 1.1.21