/[lmdze]/trunk/Sources/dyn3d/Guide/guide.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/Guide/guide.f

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

revision 108 by guez, Tue Sep 16 14:00:41 2014 UTC revision 109 by guez, Wed Sep 17 10:08:00 2014 UTC
# Line 30  CONTAINS Line 30  CONTAINS
30      USE q_sat_m, ONLY: q_sat      USE q_sat_m, ONLY: q_sat
31      use read_reanalyse_m, only: read_reanalyse      use read_reanalyse_m, only: read_reanalyse
32      USE serre, ONLY: clat, clon      USE serre, ONLY: clat, clon
33      use tau2alpha_m, only: tau2alpha, dxdys      use tau2alpha_m, only: tau2alpha
34      use writefield_m, only: writefield      use writefield_m, only: writefield
35    
36      INTEGER, INTENT(IN):: itau      INTEGER, INTENT(IN):: itau
# Line 71  CONTAINS Line 71  CONTAINS
71    
72      INTEGER ilon, ilat      INTEGER ilon, ilat
73      REAL factt ! pas de temps entre deux appels au guidage, en fraction de jour      REAL factt ! pas de temps entre deux appels au guidage, en fraction de jour
     real ztau(iim + 1, jjm + 1)  
74    
75      INTEGER ij, l      INTEGER ij, l
76      INTEGER ncid, dimid      INTEGER ncid, dimid
# Line 139  CONTAINS Line 138  CONTAINS
138         CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, &         CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, &
139              masserea2, nlev)              masserea2, nlev)
140         qrea2 = max(qrea2, 0.1)         qrea2 = max(qrea2, 0.1)
141    
142           CALL writefield("alpha_u", alpha_u)
143           CALL writefield("alpha_t", alpha_t)
144      END IF first_call      END IF first_call
145    
146      ! IMPORTATION DES VENTS, PRESSION ET TEMPERATURE REELS:      ! IMPORTATION DES VENTS, PRESSION ET TEMPERATURE REELS:
# Line 157  CONTAINS Line 159  CONTAINS
159              masserea2, nlev)              masserea2, nlev)
160         qrea2 = max(qrea2, 0.1)         qrea2 = max(qrea2, 0.1)
161         factt = dtvr * iperiod / daysec         factt = dtvr * iperiod / daysec
        ztau = factt / max(alpha_t, 1E-10)  
        CALL writefield("aire", aire)  
        CALL writefield("dxdys", dxdys)  
        CALL writefield("alpha_u", alpha_u)  
        CALL writefield("alpha_t", alpha_t)  
        CALL writefield("ztau", ztau)  
162         CALL writefield("ucov", ucov)         CALL writefield("ucov", ucov)
163         CALL writefield("ucovrea2", ucovrea2)         CALL writefield("ucovrea2", ucovrea2)
164         CALL writefield("teta", teta)         CALL writefield("teta", teta)
Legend:
Removed from v.108  
changed lines
  Added in v.109
  ViewVC Help
Powered by ViewVC 1.1.21