/[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

trunk/libf/dyn3d/Dissipation/inidissip.f90 revision 65 by guez, Thu Sep 20 09:57:03 2012 UTC trunk/dyn3d/Dissipation/inidissip.f revision 107 by guez, Thu Sep 11 15:09:15 2014 UTC
# Line 23  contains Line 23  contains
23      USE comconst, ONLY: dtvr      USE comconst, ONLY: dtvr
24      use comdissnew, only: nitergdiv, nitergrot, niterh, tetagdiv, tetagrot, &      use comdissnew, only: nitergdiv, nitergrot, niterh, tetagdiv, tetagrot, &
25           tetatemp           tetatemp
26      USE comvert, ONLY: preff, presnivs      USE disvert_m, ONLY: preff, presnivs
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
# Line 45  contains Line 45  contains
45    
46      PRINT *, 'Call sequence information: inidissip'      PRINT *, 'Call sequence information: inidissip'
47      call random_seed(size=seed_size)      call random_seed(size=seed_size)
48      call random_seed(put=(/(0, ii = 1, seed_size)/))      call random_seed(put=(/(1, ii = 1, seed_size)/))
49    
50      PRINT *, 'Calcul des valeurs propres de divgrad'      PRINT *, 'Calcul des valeurs propres de divgrad'
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.65  
changed lines
  Added in v.107
  ViewVC Help
Powered by ViewVC 1.1.21