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

Diff of /trunk/dyn3d/fluxstokenc.f

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

revision 31 by guez, Thu Apr 1 14:59:19 2010 UTC revision 32 by guez, Tue Apr 6 17:52:58 2010 UTC
# Line 23  SUBROUTINE fluxstokenc(pbaru, pbarv, mas Line 23  SUBROUTINE fluxstokenc(pbaru, pbarv, mas
23    real, SAVE:: tetac(ip1jmp1, llm), phic(ip1jmp1, llm)    real, SAVE:: tetac(ip1jmp1, llm), phic(ip1jmp1, llm)
24    
25    REAL pbarug(ip1jmp1, llm), pbarvg(iip1, jjm, llm), wg(ip1jmp1, llm)    REAL pbarug(ip1jmp1, llm), pbarvg(iip1, jjm, llm), wg(ip1jmp1, llm)
   INTEGER:: iadvtr = 0  
26    REAL tst(1), ist(1), istp(1)    REAL tst(1), ist(1), istp(1)
27    INTEGER ij, l    INTEGER ij, l
28    INTEGER, save:: fluxid, fluxvid    INTEGER, save:: fluxid, fluxvid
29    integer fluxdid    integer fluxdid
   LOGICAL:: first = .TRUE.  
30    
31    !-------------------------------------------------------------    !-------------------------------------------------------------
32    
33    IF (first) THEN    IF (itau == 0) THEN
34       CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, nqmx, &       CALL initfluxsto(time_step, istdyn*time_step, istdyn*time_step, nqmx, &
35            fluxid, fluxvid, fluxdid)            fluxid, fluxvid, fluxdid)
36       CALL histwrite(fluxid, 'phis', 1, phis)       CALL histwrite(fluxid, 'phis', 1, phis)
# Line 44  SUBROUTINE fluxstokenc(pbaru, pbarv, mas Line 42  SUBROUTINE fluxstokenc(pbaru, pbarv, mas
42       istp(1) = istphy       istp(1) = istphy
43       CALL histwrite(fluxdid, 'istphy', 1, istp)       CALL histwrite(fluxdid, 'istphy', 1, istp)
44    
      first = .FALSE.  
   END IF  
   
   IF (itau == 0) THEN  
45       CALL initial0(ijp1llm, phic)       CALL initial0(ijp1llm, phic)
46       CALL initial0(ijp1llm, tetac)       CALL initial0(ijp1llm, tetac)
47       CALL initial0(ijp1llm, pbaruc)       CALL initial0(ijp1llm, pbaruc)
# Line 71  SUBROUTINE fluxstokenc(pbaru, pbarv, mas Line 65  SUBROUTINE fluxstokenc(pbaru, pbarv, mas
65       CALL scopy(ip1jmp1*llm, masse, 1, massem, 1)       CALL scopy(ip1jmp1*llm, masse, 1, massem, 1)
66    END IF    END IF
67    
68    iadvtr = iadvtr + 1    IF (mod(itau + 1, istdyn) == 0) THEN
69         ! on advecte a ce pas de temps
   !   Test pour savoir si on advecte a ce pas de temps  
   IF (iadvtr == istdyn) THEN  
70       !    normalisation       !    normalisation
71       DO l = 1, llm       DO l = 1, llm
72          DO ij = 1, ip1jmp1          DO ij = 1, ip1jmp1
# Line 92  SUBROUTINE fluxstokenc(pbaru, pbarv, mas Line 84  SUBROUTINE fluxstokenc(pbaru, pbarv, mas
84       !     2. groupement des mailles pres du pole.       !     2. groupement des mailles pres du pole.
85    
86       CALL groupe(massem, pbaruc, pbarvc, pbarug, pbarvg, wg)       CALL groupe(massem, pbaruc, pbarvc, pbarug, pbarvg, wg)
      iadvtr = 0  
      PRINT *, 'ITAU auqel on stoke les fluxmasses', itau  
87    
88       CALL histwrite(fluxid, 'masse', itau, massem)       CALL histwrite(fluxid, 'masse', itau, massem)
89       CALL histwrite(fluxid, 'pbaru', itau, pbarug)       CALL histwrite(fluxid, 'pbaru', itau, pbarug)
Legend:
Removed from v.31  
changed lines
  Added in v.32
  ViewVC Help
Powered by ViewVC 1.1.21