/[lmdze]/trunk/libf/dyn3d/leapfrog.f90
ViewVC logotype

Diff of /trunk/libf/dyn3d/leapfrog.f90

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

revision 10 by guez, Fri Apr 18 14:45:53 2008 UTC revision 22 by guez, Fri Jul 31 15:18:47 2009 UTC
# Line 6  module leapfrog_m Line 6  module leapfrog_m
6    
7  contains  contains
8    
9    SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, nq, q, clesphy0, &    SUBROUTINE leapfrog(ucov, vcov, teta, ps, masse, phis, nq, q, time_0)
        time_0)  
10    
11      ! From dyn3d/leapfrog.F, version 1.6 2005/04/13 08:58:34      ! From dyn3d/leapfrog.F, version 1.6 2005/04/13 08:58:34
12    
# Line 34  contains Line 33  contains
33      use comvert, only: ap, bp      use comvert, only: ap, bp
34      use conf_gcm_m, only: day_step, iconser, idissip, iphysiq, iperiod, nday, &      use conf_gcm_m, only: day_step, iconser, idissip, iphysiq, iperiod, nday, &
35           offline, periodav           offline, periodav
36      use logic, only: ok_guide, apdiss, apphys, conser, forward, iflag_phys, &      use logic, only: ok_guide, iflag_phys
          leapf, statcl  
37      use comgeom      use comgeom
38      use serre      use serre
39      use temps, only: itaufin, day_ini, dt      use temps, only: itaufin, day_ini, dt
# Line 48  contains Line 46  contains
46      use pression_m, only: pression      use pression_m, only: pression
47      use pressure_var, only: p3d      use pressure_var, only: p3d
48    
49      integer nq      integer, intent(in):: nq
     REAL clesphy0(:)  
50    
51      ! Variables dynamiques:      ! Variables dynamiques:
52      REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants      REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants
# Line 95  contains Line 92  contains
92    
93      REAL tppn(iim), tpps(iim), tpn, tps      REAL tppn(iim), tpps(iim), tpn, tps
94    
95      INTEGER itau, itaufinp1      INTEGER itau ! index of the time step of the dynamics, starts at 0
96        integer itaufinp1
97      INTEGER iday ! jour julien      INTEGER iday ! jour julien
98      REAL time ! Heure de la journee en fraction d'1 jour      REAL time ! time of day, as a fraction of day length
99    
100      REAL SSUM      REAL SSUM
101      real finvmaold(ip1jmp1, llm)      real finvmaold(ip1jmp1, llm)
# Line 120  contains Line 118  contains
118    
119      logical:: dissip_conservative = .true.      logical:: dissip_conservative = .true.
120      LOGICAL:: prem = .true.      LOGICAL:: prem = .true.
121        logical forward, leapf, apphys, conser, apdiss
122    
123      !---------------------------------------------------      !---------------------------------------------------
124    
# Line 167  contains Line 166  contains
166            ! gestion des appels de la physique et des dissipations:            ! gestion des appels de la physique et des dissipations:
167    
168            apphys = .FALSE.            apphys = .FALSE.
           statcl = .FALSE.  
169            conser = .FALSE.            conser = .FALSE.
170            apdiss = .FALSE.            apdiss = .FALSE.
171    
# Line 196  contains Line 194  contains
194    
195            ! integrations dynamique et traceurs:            ! integrations dynamique et traceurs:
196            CALL integrd(2, vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, &            CALL integrd(2, vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, &
197                 dteta, dq, dp, vcov, ucov, teta, q, ps, masse, phis, finvmaold)                 dteta, dq, dp, vcov, ucov, teta, q, ps, masse, phis, &
198                   finvmaold, leapf)
199    
200            ! calcul des tendances physiques:            ! calcul des tendances physiques:
201    
# Line 220  contains Line 219  contains
219    
220               CALL calfis(nq, lafin, rdayvrai, time, ucov, vcov, teta, q, &               CALL calfis(nq, lafin, rdayvrai, time, ucov, vcov, teta, q, &
221                    masse, ps, pk, phis, phi, du, dv, dteta, dq, w, &                    masse, ps, pk, phis, phi, du, dv, dteta, dq, w, &
222                    clesphy0, dufi, dvfi, dtetafi, dqfi, dpfi)                    dufi, dvfi, dtetafi, dqfi, dpfi)
223    
224               ! ajout des tendances physiques:               ! ajout des tendances physiques:
225               CALL addfi(nqmx, dtphys, &               CALL addfi(nqmx, dtphys, &
# Line 337  contains Line 336  contains
336                  dt = 2. * dtvr                  dt = 2. * dtvr
337               END IF               END IF
338            ELSE            ELSE
339               ! ...... pas leapfrog .....               ! pas leapfrog
340               leapf = .TRUE.               leapf = .TRUE.
341               dt = 2. * dtvr               dt = 2. * dtvr
342            END IF            END IF
Legend:
Removed from v.10  
changed lines
  Added in v.22
  ViewVC Help
Powered by ViewVC 1.1.21