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

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

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

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 137 by guez, Wed May 6 15:51:03 2015 UTC
# Line 27  contains Line 27  contains
27      USE conf_gcm_m, ONLY: iperiod      USE conf_gcm_m, ONLY: iperiod
28      USE dimens_m, ONLY: iim, jjm      USE dimens_m, ONLY: iim, jjm
29      use divgrad2_m, only: divgrad2      use divgrad2_m, only: divgrad2
30      use filtreg_m, only: filtreg      use filtreg_scal_m, only: filtreg_scal
31        use filtreg_v_m, only: filtreg_v
32      use gradiv2_m, only: gradiv2      use gradiv2_m, only: gradiv2
33      use jumble, only: new_unit      use jumble, only: new_unit
34      use nxgraro2_m, only: nxgraro2      use nxgraro2_m, only: nxgraro2
# Line 51  contains Line 52  contains
52      deltap = 1.      deltap = 1.
53      call random_number(zh)      call random_number(zh)
54      zh = zh - 0.5      zh = zh - 0.5
55      CALL filtreg(zh, direct = .true., intensive = .true.)      CALL filtreg_scal(zh, direct = .true., intensive = .true.)
56    
57      DO l = 1, 50      DO l = 1, 50
58         CALL divgrad2(1, zh, deltap, niterh, divgra, -1.)         CALL divgrad2(1, zh, deltap, niterh, divgra, -1.)
# Line 65  contains Line 66  contains
66      PRINT *, 'Calcul des valeurs propres de gradiv'      PRINT *, 'Calcul des valeurs propres de gradiv'
67      call random_number(zu)      call random_number(zu)
68      zu = zu - 0.5      zu = zu - 0.5
69      CALL filtreg(zu, direct = .true., intensive = .true.)      CALL filtreg_scal(zu, direct = .true., intensive = .true.)
70      call random_number(zv)      call random_number(zv)
71      zv = zv - 0.5      zv = zv - 0.5
72      CALL filtreg(zv, direct = .true., intensive = .true.)      CALL filtreg_v(zv, intensive = .true.)
73    
74      DO l = 1, 50      DO l = 1, 50
75         CALL gradiv2(zu, zv, nitergdiv, gx, gy, -1.)         CALL gradiv2(zu, zv, nitergdiv, gx, gy, -1.)
# Line 83  contains Line 84  contains
84      PRINT *, 'Calcul des valeurs propres de nxgrarot'      PRINT *, 'Calcul des valeurs propres de nxgrarot'
85      call random_number(zu)      call random_number(zu)
86      zu = zu - 0.5      zu = zu - 0.5
87      CALL filtreg(zu, direct = .true., intensive = .true.)      CALL filtreg_scal(zu, direct = .true., intensive = .true.)
88      call random_number(zv)      call random_number(zv)
89      zv = zv - 0.5      zv = zv - 0.5
90      CALL filtreg(zv, direct = .true., intensive = .true.)      CALL filtreg_v(zv, intensive = .true.)
91    
92      DO l = 1, 50      DO l = 1, 50
93         CALL nxgraro2(zu, zv, nitergrot, gx, gy, -1.)         CALL nxgraro2(zu, zv, nitergrot, gx, gy, -1.)
Legend:
Removed from v.134  
changed lines
  Added in v.137
  ViewVC Help
Powered by ViewVC 1.1.21