/[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 103 by guez, Wed Mar 26 17:18:58 2014 UTC revision 104 by guez, Thu Sep 4 10:05:52 2014 UTC
# Line 12  contains Line 12  contains
12      ! Objet : calcul des tendances dynamiques      ! Objet : calcul des tendances dynamiques
13    
14      use advect_m, only: advect      use advect_m, only: advect
15        use bernoui_m, only: bernoui
16        USE comconst, ONLY: daysec, dtvr
17      USE comgeom, ONLY: airesurg, constang_2d      USE comgeom, ONLY: airesurg, constang_2d
18        USE conf_gcm_m, ONLY: day_step
19      use convmas_m, only: convmas      use convmas_m, only: convmas
20      USE dimens_m, ONLY: iim, jjm, llm      USE dimens_m, ONLY: iim, jjm, llm
21      USE disvert_m, ONLY: ap, bp      USE disvert_m, ONLY: ap, bp
22      use dteta1_m, only: dteta1      use dteta1_m, only: dteta1
23      use dudv1_m, only: dudv1      use dudv1_m, only: dudv1
24      use dudv2_m, only: dudv2      use dudv2_m, only: dudv2
25        USE dynetat0_m, ONLY: day_ini
26      use flumass_m, only: flumass      use flumass_m, only: flumass
27      use massbar_m, only: massbar      use massbar_m, only: massbar
28      use massbarxy_m, only: massbarxy      use massbarxy_m, only: massbarxy
29      use massdair_m, only: massdair      use massdair_m, only: massdair
30      USE paramet_m, ONLY: iip1, ip1jmp1, jjp1, llmp1      USE paramet_m, ONLY: iip1, ip1jmp1, jjp1, llmp1
31      use sortvarc_m, only: sortvarc      use sortvarc_m, only: sortvarc, ang, etot, ptot, rmsdpdt, rmsv, stot, ztot
32      use tourpot_m, only: tourpot      use tourpot_m, only: tourpot
33      use vitvert_m, only: vitvert      use vitvert_m, only: vitvert
34    
# Line 37  contains Line 41  contains
41      REAL, INTENT(IN):: pk(iip1, jjp1, llm)      REAL, INTENT(IN):: pk(iip1, jjp1, llm)
42      REAL, INTENT(IN):: pkf(ip1jmp1, llm)      REAL, INTENT(IN):: pkf(ip1jmp1, llm)
43      REAL, INTENT(IN):: phis(ip1jmp1)      REAL, INTENT(IN):: phis(ip1jmp1)
44      REAL, INTENT(IN):: phi(ip1jmp1, llm)      REAL, INTENT(IN):: phi(iim + 1, jjm + 1, llm)
45      REAL dudyn(:, :, :) ! (iim + 1, jjm + 1, llm)      REAL dudyn(:, :, :) ! (iim + 1, jjm + 1, llm)
46      real dv((iim + 1) * jjm, llm)      real dv((iim + 1) * jjm, llm)
47      REAL, INTENT(out):: dteta(ip1jmp1, llm)      REAL, INTENT(out):: dteta(ip1jmp1, llm)
# Line 49  contains Line 53  contains
53    
54      ! Local:      ! Local:
55      REAL vcont((iim + 1) * jjm, llm), ucont(ip1jmp1, llm)      REAL vcont((iim + 1) * jjm, llm), ucont(ip1jmp1, llm)
56      REAL ang(iim + 1, jjm + 1, llm), p(ip1jmp1, llmp1)      REAL ang_3d(iim + 1, jjm + 1, llm), p(ip1jmp1, llmp1)
57      REAL massebx(ip1jmp1, llm), masseby((iim + 1) * jjm, llm)      REAL massebx(ip1jmp1, llm), masseby((iim + 1) * jjm, llm)
58      REAL vorpot(iim + 1, jjm, llm)      REAL vorpot(iim + 1, jjm, llm)
59      real ecin(ip1jmp1, llm), convm(ip1jmp1, llm)      real ecin(iim + 1, jjm + 1, llm), convm(ip1jmp1, llm)
60      REAL bern(ip1jmp1, llm)      REAL bern(iim + 1, jjm + 1, llm)
61      REAL massebxy(iim + 1, jjm, llm)      REAL massebxy(iim + 1, jjm, llm)
62      INTEGER ij, l      INTEGER ij, l
63        real heure, time
64    
65      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
66    
# Line 72  contains Line 77  contains
77      CALL tourpot(vcov, ucov, massebxy, vorpot)      CALL tourpot(vcov, ucov, massebxy, vorpot)
78      CALL dudv1(vorpot, pbaru, pbarv, dudyn(:, 2: jjm, :), dv)      CALL dudv1(vorpot, pbaru, pbarv, dudyn(:, 2: jjm, :), dv)
79      CALL enercin(vcov, ucov, vcont, ucont, ecin)      CALL enercin(vcov, ucov, vcont, ucont, ecin)
80      CALL bernoui(ip1jmp1, llm, phi, ecin, bern)      CALL bernoui(phi, ecin, bern)
81      CALL dudv2(teta, pkf, bern, dudyn, dv)      CALL dudv2(teta, pkf, bern, dudyn, dv)
82    
83      forall (l = 1: llm) ang(:, :, l) = ucov(:, :, l) + constang_2d      forall (l = 1: llm) ang_3d(:, :, l) = ucov(:, :, l) + constang_2d
84      CALL advect(ang, vcov, teta, w, massebx, masseby, dudyn, dv, dteta)      CALL advect(ang_3d, vcov, teta, w, massebx, masseby, dudyn, dv, dteta)
85    
86      ! Warning problème de périodicité de dv sur les PC Linux. Problème      ! Warning problème de périodicité de dv sur les PC Linux. Problème
87      ! d'arrondi probablement. Observé sur le code compilé avec pgf90      ! d'arrondi probablement. Observé sur le code compilé avec pgf90
# Line 90  contains Line 95  contains
95      END DO      END DO
96    
97      ! Sorties éventuelles des variables de contrôle :      ! Sorties éventuelles des variables de contrôle :
98      IF (conser) CALL sortvarc(itau, ucov, teta, ps, masse, pk, phis, vorpot, &      IF (conser) then
99           phi, bern, dp, time_0)         CALL sortvarc(ucov, teta, ps, masse, pk, phis, vorpot, phi, bern, dp, &
100                resetvarc = .false.)
101           time = real(itau) / day_step + time_0
102           heure = mod(itau * dtvr / daysec, 1.) * 24.
103           IF (abs(heure-24.) <= 1e-4) heure = 0.
104    
105           PRINT 3500, itau, int(day_ini + time), heure, time
106           PRINT 4000, ptot, rmsdpdt, etot, ztot, stot, rmsv, ang
107        end IF
108    
109    3500 FORMAT (4X, 'pas', I7, 5X, 'jour', i5, 1X, 'heure', F5.1, 4X, 'date', &
110             F10.5)
111    4000 FORMAT (10X, 'masse', 4X, 'rmsdpdt', 7X, 'energie', 2X, 'enstrophie', &
112             2X, 'entropie', 3X, 'rmsv', 4X, 'mt.ang', /, 'GLOB  ', F10.6, &
113             E13.6, 5F10.3/)
114    
115    END SUBROUTINE caldyn    END SUBROUTINE caldyn
116    
Legend:
Removed from v.103  
changed lines
  Added in v.104
  ViewVC Help
Powered by ViewVC 1.1.21