/[lmdze]/trunk/dyn3d/caldyn.f
ViewVC logotype

Diff of /trunk/dyn3d/caldyn.f

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

revision 254 by guez, Mon Feb 5 10:39:38 2018 UTC revision 259 by guez, Tue Mar 6 16:19:52 2018 UTC
# Line 23  contains Line 23  contains
23      use dteta1_m, only: dteta1      use dteta1_m, only: dteta1
24      use dudv1_m, only: dudv1      use dudv1_m, only: dudv1
25      use dudv2_m, only: dudv2      use dudv2_m, only: dudv2
26      USE dynetat0_m, ONLY: day_ini      USE dynetat0_m, ONLY: day_ini, ang0, etot0, ptot0, stot0, ztot0
27      use enercin_m, only: enercin      use enercin_m, only: enercin
28      use flumass_m, only: flumass      use flumass_m, only: flumass
29      use massbar_m, only: massbar      use massbar_m, only: massbar
30      use massbarxy_m, only: massbarxy      use massbarxy_m, only: massbarxy
31      use massdair_m, only: massdair      use massdair_m, only: massdair
32      USE paramet_m, ONLY: iip1, ip1jmp1, jjp1, llmp1      USE paramet_m, ONLY: iip1, ip1jmp1, jjp1, llmp1
33      use sortvarc_m, only: sortvarc, ang, etot, ptot, rmsdpdt, rmsv, stot, ztot      use sortvarc_m, only: sortvarc
34      use tourpot_m, only: tourpot      use tourpot_m, only: tourpot
35      use vitvert_m, only: vitvert      use vitvert_m, only: vitvert
36    
# Line 63  contains Line 63  contains
63      REAL massebxy(iim + 1, jjm, llm)      REAL massebxy(iim + 1, jjm, llm)
64      INTEGER ij, l      INTEGER ij, l
65      real heure, time      real heure, time
66        real ang, etot, ptot, ztot, stot, rmsdpdt, rmsv
67    
68      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
69    
# Line 99  contains Line 100  contains
100      ! Sorties éventuelles des variables de contrôle :      ! Sorties éventuelles des variables de contrôle :
101      IF (conser) then      IF (conser) then
102         CALL sortvarc(ucov, teta, ps, masse, pk, phis, vorpot, phi, bern, dp, &         CALL sortvarc(ucov, teta, ps, masse, pk, phis, vorpot, phi, bern, dp, &
103              resetvarc = .false.)              ang, etot, ptot, ztot, stot, rmsdpdt, rmsv)
104    
105         time = real(itau) / day_step         time = real(itau) / day_step
106         heure = mod(itau * dtvr / daysec, 1.) * 24.         heure = mod(itau * dtvr / daysec, 1.) * 24.
107         IF (abs(heure-24.) <= 1e-4) heure = 0.         IF (abs(heure-24.) <= 1e-4) heure = 0.
108    
109         PRINT 3500, itau, int(day_ini + time), heure, time         PRINT 3500, itau, int(day_ini + time), heure, time
110         PRINT 4000, ptot, rmsdpdt, etot, ztot, stot, rmsv, ang         PRINT 4000, ptot / ptot0, rmsdpdt, etot / etot0, ztot / ztot0, &
111                stot / stot0, sqrt(rmsv / ptot), ang / ang0
112      end IF      end IF
113    
114  3500 FORMAT (4X, 'pas', I7, 5X, 'jour', i5, 1X, 'heure', F5.1, 4X, 'date', &  3500 FORMAT (4X, 'pas', I7, 5X, 'jour', i5, 1X, 'heure', F5.1, 4X, 'date', &
Legend:
Removed from v.254  
changed lines
  Added in v.259
  ViewVC Help
Powered by ViewVC 1.1.21