Changeset 878 for LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Jan 14, 2008, 1:03:39 PM (16 years ago)

Author:

Laurent Fairhead

Message:

Bascule de la physique du LMD vers la physique avec thermiques
LF

File:

• LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (10 diffs)

LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -197 +197 @@
        wfbils,    wfbilo,    flux_t,   flux_u, flux_v,&
        dflux_t,   dflux_q,   zxsnow,                  &
-       zxfluxt,   zxfluxq,   q2m,      flux_q )
+       zxfluxt,   zxfluxq,   q2m,      flux_q, tke    )
 !****************************************************************************************
 ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
@@ -240 +240 @@
 ! flux_t---output-R- flux de chaleur sensible (CpT) J/m**2/s (W/m**2)
 !                    (orientation positive vers le bas)
+! tke---input/output-R- tke (kg/m**2/s)
 ! flux_q---output-R- flux de vapeur d'eau (kg/m**2/s)
 ! flux_u---output-R- tension du vent X: (kg m/s)/(m**2 s) ou Pascal
@@ -348 +349 @@
     REAL, DIMENSION(klon, nbsrf),INTENT(OUT)        :: q2m
     REAL, DIMENSION(klon, klev, nbsrf), INTENT(OUT) :: flux_q
+
+! Input/output
+    REAL, DIMENSION(klon, klev+1, nbsrf), INTENT(INOUT) :: tke
 
 
@@ -421 +425 @@
     REAL, DIMENSION(klon,klev)         :: delp
     REAL, DIMENSION(klon,klev+1)       :: ypaprs
+    REAL, DIMENSION(klon,klev+1)       :: ytke
     REAL, DIMENSION(klon,nsoilmx)      :: ytsoil
     REAL, DIMENSION(klon,nbsrf)        :: pctsrf_pot
@@ -451 +456 @@
 
 
+!jg+ temporary test
+    REAL, DIMENSION(klon,klev)         :: y_flux_u_old, y_flux_v_old
+    REAL, DIMENSION(klon,klev)         :: y_d_u_old, y_d_v_old
+!jg-
+    
 !****************************************************************************************
 ! End of declarations
@@ -514 +524 @@
     flux_u = 0.0  ; flux_v = 0.0     ; d_t = 0.0       ; d_q = 0.0      
     d_u = 0.0     ; d_v = 0.0        ; zcoefh = 0.0    ; yqsol = 0.0    
-    ytherm = 0.0
+    ytherm = 0.0  ; ytke=0.
      
     ytsoil = 999999. 
@@ -657 +667 @@
              ypplay(j,k) = pplay(i,k)
              ydelp(j,k) = delp(i,k)
+             ytke(j,k)=tke(i,k,nsrf)
              yu(j,k) = u(i,k)
              yv(j,k) = v(i,k)
@@ -687 +698 @@
        CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
             ypaprs, ypplay, yu, yv, yq, yt, yts, yrugos, yqsurf,  &
-            ycoefm, ycoefh)
+            ycoefm, ycoefh,ytke)
        
 !****************************************************************************************
@@ -839 +850 @@
 !****************************************************************************************
 
+       tke(:,:,nsrf)=0.
        DO k = 1, klev
           DO j = 1, knon
@@ -853 +865 @@
              flux_u(i,k,nsrf) = y_flux_u(j,k)
              flux_v(i,k,nsrf) = y_flux_v(j,k)
+
+             tke(i,k,nsrf)=ytke(j,k)
+
           ENDDO
        ENDDO