Changeset 2258 for LMDZ5/branches/testing/libf/phylmd/change_srf_frac_mod.F90

Timestamp:

Apr 13, 2015, 10:21:09 AM (10 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes 2216:2237 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/change_srf_frac_mod.F90 (modified) (5 diffs)

LMDZ5/branches/testing/libf/phylmd/change_srf_frac_mod.F90

@@ -1 +1 @@
 !
-! $Header$
+! $Id$
 !
 MODULE change_srf_frac_mod
@@ -12 +12 @@
 
   SUBROUTINE change_srf_frac(itime, dtime, jour, &
-       pctsrf, alb1, alb2, tsurf, ustar, u10m, v10m, pbl_tke)
+!albedo SB >>>
+!       pctsrf, alb1, alb2, tsurf, ustar, u10m, v10m, pbl_tke)
+        pctsrf, alb_dir, alb_dif, tsurf, ustar, u10m, v10m, pbl_tke)
+!albedo SB <<<
+   
+
+
 !
 ! This subroutine is called from physiq.F at each timestep. 
@@ -32 +38 @@
     INCLUDE "iniprint.h"
     INCLUDE "YOMCST.h"
+!albedo SB >>>
+    include "clesphys.h"
+!albedo SB <<<
+
+
 
 ! Input arguments
@@ -43 +54 @@
    
     REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: pctsrf ! sub-surface fraction
-    REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: alb1   ! albedo first interval in SW spektrum
-    REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: alb2   ! albedo second interval in SW spektrum
+!albedo SB >>>
+!   REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: alb1   ! albedo first interval in SW spektrum
+!   REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: alb2   ! albedo second interval in SW spektrum
+    REAL, DIMENSION(klon,nsw,nbsrf), INTENT(INOUT) :: alb_dir,alb_dif
+!albedo SB <<<
+
     REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: tsurf
     REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: ustar
@@ -160 +175 @@
 ! 
 !****************************************************************************************
-       CALL pbl_surface_newfrac(itime, pctsrf, pctsrf_old, tsurf, alb1, alb2, ustar, u10m, v10m, pbl_tke)
+
+!albedo SB >>>
+! CALL pbl_surface_newfrac(itime, pctsrf, pctsrf_old, tsurf, alb1, alb2, ustar,
+! u10m, v10m, pbl_tke)
+       CALL pbl_surface_newfrac(itime, pctsrf, pctsrf_old, tsurf, alb_dir,alb_dif, ustar, u10m, v10m, pbl_tke)
+!albedo SB <<<
+
+
 
     ELSE