Changeset 1067 for LMDZ4/trunk/libf/phylmd/surf_land_mod.F90

Timestamp:

Dec 17, 2008, 2:30:13 PM (15 years ago)

Author:

Laurent Fairhead

Message:

• Modifications lie au premier niveau du modele pour la diffusion turbulent

du vent.

• Preparation pour un couplage des courrant oceaniques.

JG

File:

• LMDZ4/trunk/libf/phylmd/surf_land_mod.F90 (modified) (9 diffs)

LMDZ4/trunk/libf/phylmd/surf_land_mod.F90

@@ -1 @@
-!
-! $Header$
 !
 MODULE surf_land_mod
   
-  USE surface_data, ONLY    : ok_veget
-  USE dimphy
- 
-#ifdef CPP_VEGET
-  USE surf_land_orchidee_mod
-#endif
-  USE surf_land_bucket_mod
-
   IMPLICIT NONE
 
@@ -21 +11 @@
        rlon, rlat, &
        debut, lafin, zlev, ccanopy, swnet, lwnet, albedo, &
-       tsurf, p1lay, tq_cdrag, precip_rain, precip_snow, temp_air, spechum, &
-       petAcoef, peqAcoef, petBcoef, peqBcoef, &
-       pref, u1_lay, v1_lay, rugoro, pctsrf, &
+       tsurf, p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, &
+       AcoefH, AcoefQ, BcoefH, BcoefQ, &
+       AcoefU, AcoefV, BcoefU, BcoefV, & 
+       pref, u1, v1, rugoro, pctsrf, &
        snow, qsol, agesno, tsoil, &
        z0_new, alb1_new, alb2_new, evap, fluxsens, fluxlat, &
        qsurf, tsurf_new, dflux_s, dflux_l, &
+       flux_u1, flux_v1, & 
        lwdown_m)
+
+    USE dimphy
+    USE surface_data, ONLY    : ok_veget
+    USE surf_land_orchidee_mod
+    USE surf_land_bucket_mod
+    USE calcul_fluxs_mod
 
     INCLUDE "indicesol.h"
@@ -46 +44 @@
     REAL, DIMENSION(klon), INTENT(IN)       :: tsurf
     REAL, DIMENSION(klon), INTENT(IN)       :: p1lay
-    REAL, DIMENSION(klon), INTENT(IN)       :: tq_cdrag
+    REAL, DIMENSION(klon), INTENT(IN)       :: cdragh, cdragm
     REAL, DIMENSION(klon), INTENT(IN)       :: precip_rain, precip_snow
     REAL, DIMENSION(klon), INTENT(IN)       :: temp_air, spechum
-    REAL, DIMENSION(klon), INTENT(IN)       :: petAcoef, peqAcoef
-    REAL, DIMENSION(klon), INTENT(IN)       :: petBcoef, peqBcoef
+    REAL, DIMENSION(klon), INTENT(IN)       :: AcoefH, AcoefQ, BcoefH, BcoefQ
+    REAL, DIMENSION(klon), INTENT(IN)       :: AcoefU, AcoefV, BcoefU, BcoefV
     REAL, DIMENSION(klon), INTENT(IN)       :: pref   ! pressure reference
-    REAL, DIMENSION(klon), INTENT(IN)       :: u1_lay, v1_lay
+    REAL, DIMENSION(klon), INTENT(IN)       :: u1, v1
     REAL, DIMENSION(klon), INTENT(IN)       :: rugoro
     REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf
@@ -75 +73 @@
     REAL, DIMENSION(klon), INTENT(OUT)       :: tsurf_new
     REAL, DIMENSION(klon), INTENT(OUT)       :: dflux_s, dflux_l      
+    REAL, DIMENSION(klon), INTENT(OUT)       :: flux_u1, flux_v1  ! flux for U and V at first model level
 
 ! Local variables
@@ -84 +83 @@
     REAL, DIMENSION(klon) :: epot_air           ! potential air temperature
     REAL, DIMENSION(klon) :: tsol_rad, emis_new ! output from interfsol not used
+    REAL, DIMENSION(klon) :: u0, v0     ! surface speed
     INTEGER               :: i
 
@@ -117 +117 @@
        END DO
 
-#ifdef CPP_VEGET
        ! temporary for keeping same results using lwdown_m instead of lwdown
        CALL surf_land_orchidee(itime, dtime, date0, knon, &
             knindex, rlon, rlat, pctsrf, &
             debut, lafin, &
-            zlev,  u1_lay, v1_lay, temp_air, spechum, epot_air, ccanopy, & 
-            tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, &
+            zlev,  u1, v1, temp_air, spechum, epot_air, ccanopy, & 
+            cdragh, AcoefH, AcoefQ, BcoefH, BcoefQ, &
             precip_rain, precip_snow, lwdown_m, swnet, swdown, &
             pref_tmp, &
@@ -129 +128 @@
             tsol_rad, tsurf_new, alb1_new, alb2_new, &
             emis_new, z0_new, qsurf)       
-#endif
 
 !  
@@ -144 +142 @@
 !****************************************************************************************
        CALL surf_land_bucket(itime, jour, knon, knindex, debut, dtime,&
-            tsurf, p1lay, tq_cdrag, precip_rain, precip_snow, temp_air, &
-            spechum, petAcoef, peqAcoef, petBcoef, peqBcoef, pref, &
-            u1_lay, v1_lay, rugoro, swnet, lwnet, &
+            tsurf, p1lay, cdragh, precip_rain, precip_snow, temp_air, &
+            spechum, AcoefH, AcoefQ, BcoefH, BcoefQ, pref, &
+            u1, v1, rugoro, swnet, lwnet, &
             snow, qsol, agesno, tsoil, &
             qsurf, z0_new, alb1_new, alb2_new, evap, &
@@ -153 +151 @@
     ENDIF ! ok_veget
 
+!****************************************************************************************
+! Calculation for all land models
+! - Flux calculation at first modele level for U and V
+!****************************************************************************************
+! Suppose zero surface speed
+    u0(:)=0.0
+    v0(:)=0.0
+    CALL calcul_flux_wind(knon, dtime, &
+         u0, v0, u1, v1, cdragm, &
+         AcoefU, AcoefV, BcoefU, BcoefV, &
+         p1lay, temp_air, &
+         flux_u1, flux_v1)
+    
   END SUBROUTINE surf_land
 !