/[lmdze]/trunk/Sources/dyn3d/fluxstokenc.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/fluxstokenc.f

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

trunk/libf/dyn3d/fluxstokenc.f90 revision 69 by guez, Mon Feb 18 16:33:12 2013 UTC trunk/Sources/dyn3d/fluxstokenc.f revision 178 by guez, Fri Mar 11 18:47:26 2016 UTC
# Line 9  contains Line 9  contains
9      ! Author: F. Hourdin      ! Author: F. Hourdin
10    
11      USE histwrite_m, ONLY: histwrite      USE histwrite_m, ONLY: histwrite
12      USE dimens_m, ONLY: jjm, llm, nqmx      use initfluxsto_m, only: initfluxsto
13      USE paramet_m, ONLY: iip1, ijmllm, ijp1llm, ip1jm, ip1jmp1, jjp1      USE dimens_m, ONLY: jjm, llm
14        USE paramet_m, ONLY: iip1, ijmllm, ijp1llm, ip1jm, ip1jmp1
15      USE comgeom, ONLY: aire      USE comgeom, ONLY: aire
16      USE tracstoke, ONLY: istdyn, istphy      USE tracstoke, ONLY: istdyn, istphy
17    
18      REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)      REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
19      REAL masse(ip1jmp1, llm)      REAL, intent(in):: masse(ip1jmp1, llm)
20      real, intent(in):: phi(ip1jmp1, llm)      real, intent(in):: phi(ip1jmp1, llm)
21      real, intent(in):: teta(ip1jmp1, llm)      real, intent(in):: teta(ip1jmp1, llm)
22      REAL, intent(in):: phis(ip1jmp1)      REAL, intent(in):: phis(ip1jmp1)
# Line 37  contains Line 38  contains
38      !-------------------------------------------------------------      !-------------------------------------------------------------
39    
40      IF (itau == 0) THEN      IF (itau == 0) THEN
41         CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, nqmx, &         CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, &
42              fluxid, fluxvid, fluxdid)              fluxid, fluxvid, fluxdid)
43         CALL histwrite(fluxid, 'phis', 1, phis)         CALL histwrite(fluxid, 'phis', 1, phis)
44         CALL histwrite(fluxid, 'aire', 1, aire)         CALL histwrite(fluxid, 'aire', 1, aire)
# Line 67  contains Line 68  contains
68      END DO      END DO
69    
70      !   selection de la masse instantannee des mailles avant le transport.      !   selection de la masse instantannee des mailles avant le transport.
71      IF (itau == 0) THEN      IF (itau == 0) massem = masse
        CALL scopy(ip1jmp1*llm, masse, 1, massem, 1)  
     END IF  
72    
73      IF (mod(itau + 1, istdyn) == 0) THEN      IF (mod(itau + 1, istdyn) == 0) THEN
74         ! on advecte a ce pas de temps         ! on advecte a ce pas de temps
# Line 89  contains Line 88  contains
88         !     1. calcul de w         !     1. calcul de w
89         !     2. groupement des mailles pres du pole.         !     2. groupement des mailles pres du pole.
90    
91         CALL groupe(massem, pbaruc, pbarvc, pbarug, pbarvg, wg)         CALL groupe(pbaruc, pbarvc, pbarug, pbarvg, wg)
92    
93         CALL histwrite(fluxid, 'masse', itau, massem)         CALL histwrite(fluxid, 'masse', itau, massem)
94         CALL histwrite(fluxid, 'pbaru', itau, pbarug)         CALL histwrite(fluxid, 'pbaru', itau, pbarug)
Legend:
Removed from v.69  
changed lines
  Added in v.178
  ViewVC Help
Powered by ViewVC 1.1.21