/[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/dyn3d/fluxstokenc.f revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC trunk/Sources/dyn3d/fluxstokenc.f revision 150 by guez, Thu Jun 18 13:49:26 2015 UTC
# Line 15  contains Line 15  contains
15      USE tracstoke, ONLY: istdyn, istphy      USE tracstoke, ONLY: istdyn, istphy
16    
17      REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)      REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
18      REAL masse(ip1jmp1, llm)      REAL, intent(in):: masse(ip1jmp1, llm)
19      real, intent(in):: phi(ip1jmp1, llm)      real, intent(in):: phi(ip1jmp1, llm)
20      real, intent(in):: teta(ip1jmp1, llm)      real, intent(in):: teta(ip1jmp1, llm)
21      REAL, intent(in):: phis(ip1jmp1)      REAL, intent(in):: phis(ip1jmp1)
# Line 37  contains Line 37  contains
37      !-------------------------------------------------------------      !-------------------------------------------------------------
38    
39      IF (itau == 0) THEN      IF (itau == 0) THEN
40         CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, nqmx, &         CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, &
41              fluxid, fluxvid, fluxdid)              fluxid, fluxvid, fluxdid)
42         CALL histwrite(fluxid, 'phis', 1, phis)         CALL histwrite(fluxid, 'phis', 1, phis)
43         CALL histwrite(fluxid, 'aire', 1, aire)         CALL histwrite(fluxid, 'aire', 1, aire)
# Line 67  contains Line 67  contains
67      END DO      END DO
68    
69      !   selection de la masse instantannee des mailles avant le transport.      !   selection de la masse instantannee des mailles avant le transport.
70      IF (itau == 0) THEN      IF (itau == 0) massem = masse
        CALL scopy(ip1jmp1*llm, masse, 1, massem, 1)  
     END IF  
71    
72      IF (mod(itau + 1, istdyn) == 0) THEN      IF (mod(itau + 1, istdyn) == 0) THEN
73         ! on advecte a ce pas de temps         ! on advecte a ce pas de temps
# Line 89  contains Line 87  contains
87         !     1. calcul de w         !     1. calcul de w
88         !     2. groupement des mailles pres du pole.         !     2. groupement des mailles pres du pole.
89    
90         CALL groupe(massem, pbaruc, pbarvc, pbarug, pbarvg, wg)         CALL groupe(pbaruc, pbarvc, pbarug, pbarvg, wg)
91    
92         CALL histwrite(fluxid, 'masse', itau, massem)         CALL histwrite(fluxid, 'masse', itau, massem)
93         CALL histwrite(fluxid, 'pbaru', itau, pbarug)         CALL histwrite(fluxid, 'pbaru', itau, pbarug)
Legend:
Removed from v.82  
changed lines
  Added in v.150
  ViewVC Help
Powered by ViewVC 1.1.21