Changeset 1555

Timestamp:

Jul 13, 2011, 4:05:38 PM (13 years ago)

Author:

jghattas

Message:

Ajout du calcul de la temperature sol pour le 1D si flux forces(Nicolas Rochetin).

File:

• LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -385 +385 @@
     REAL, DIMENSION(klon)              :: yts, yrugos, ypct, yz0_new
     REAL, DIMENSION(klon)              :: yalb, yalb1, yalb2
-    REAL, DIMENSION(klon)              :: yu1, yv1
+    REAL, DIMENSION(klon)              :: yu1, yv1,ytoto
     REAL, DIMENSION(klon)              :: ysnow, yqsurf, yagesno, yqsol
     REAL, DIMENSION(klon)              :: yrain_f, ysnow_f
@@ -439 +439 @@
     LOGICAL, PARAMETER                 :: zxli=.FALSE. ! utiliser un jeu de fonctions simples
     LOGICAL, PARAMETER                 :: check=.FALSE.
+    REAL, DIMENSION(klon)              :: Kech_h       ! Coefficient d'echange pour l'energie
 
 ! For debugging with IOIPSL
@@ -900 +901 @@
           y_flux_q1(:) =  flat/RLVTT
           yfluxlat(:) =  flat
+
+          Kech_h(:) = ycdragh(:) * (1.0+SQRT(yu(:,1)**2+yv(:,1)**2)) * &
+               ypplay(:,1)/(RD*yt(:,1))
+          ytoto(:)=(1./RCPD)*(AcoefH(:)+BcoefH(:)*y_flux_t1(:)*dtime)
+          ytsurf_new(:)=ytoto(:)-y_flux_t1(:)/(Kech_h(:)*RCPD)
+          y_d_ts(:) = ytsurf_new(:) - yts(:) 
+
        ELSE
           y_flux_t1(:) =  yfluxsens(:)