Changeset 1144 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Apr 8, 2009, 12:47:23 PM (16 years ago)

Author:

jghattas

Message:

Ajoute des variables q2m et t2m comme argument d'entrees dans ORCHIDEE. Pour cela, ajoute de la calcul de ces variables avec stdlevvar avant calcul sur sous-surface.

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/pbl_surface_mod.F90 (modified) (2 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/pbl_surface_mod.F90

@@ -765 +765 @@
        r_co2_ppm(:) = co2_ppm
 
+
+!****************************************************************************************
+!
+! Calulate t2m and q2m for the case of calculation at land grid points 
+! t2m and q2m are needed as input to ORCHIDEE
+!
+!****************************************************************************************
+       IF (nsrf == is_ter) THEN
+
+          DO i = 1, knon
+             zgeo1(i) = RD * yt(i,1) / (0.5*(ypaprs(i,1)+ypplay(i,1))) &
+                  * (ypaprs(i,1)-ypplay(i,1))
+          END DO
+
+          ! Calculate the temperature et relative humidity at 2m and the wind at 10m 
+          CALL stdlevvar(klon, knon, is_ter, zxli, &
+               yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, &
+               yts, yqsurf, yrugos, ypaprs(:,1), ypplay(:,1), &
+               yt2m, yq2m, yt10m, yq10m, yu10m, yustar)
+          
+       END IF
+
 !****************************************************************************************
 !
@@ -787 +809 @@
                AcoefU, AcoefV, BcoefU, BcoefV, & 
                ypsref, yu1, yv1, yrugoro, pctsrf, &
+               ylwdown, yq2m, yt2m, &
                ysnow, yqsol, yagesno, ytsoil, &
                yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
                yqsurf, ytsurf_new, y_dflux_t, y_dflux_q, &
-               y_flux_u1, y_flux_v1, & 
-               ylwdown)
+               y_flux_u1, y_flux_v1 )
+               
      
        CASE(is_lic)