Changeset 888 for LMDZ4/trunk/libf/phylmd/ocean_slab_mod.F90

Timestamp:

Feb 4, 2008, 5:24:28 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Modifications sur l'albedo JG
LF

File:

• LMDZ4/trunk/libf/phylmd/ocean_slab_mod.F90 (modified) (8 diffs)

LMDZ4/trunk/libf/phylmd/ocean_slab_mod.F90

@@ -139 +139 @@
        petAcoef, peqAcoef, petBcoef, peqBcoef, &
        ps, u1_lay, v1_lay, &
-       radsol, snow, qsurf, agesno, &
-       evap, fluxsens, fluxlat, &
-       tsurf_new, &
-       dflux_s, dflux_l, pctsrf_oce)
+       radsol, snow, agesno, &
+       qsurf, evap, fluxsens, fluxlat, &
+       tsurf_new, dflux_s, dflux_l, pctsrf_oce)
 
     INCLUDE "indicesol.h"
@@ -179 +178 @@
     INTEGER                :: i
     REAL, DIMENSION(klon)  :: cal, beta, dif_grnd
-    REAL, DIMENSION(klon)  :: alb_neig, tsurf_temp
+    REAL, DIMENSION(klon)  :: tsurf_temp
 
 !****************************************************************************************
@@ -191 +190 @@
     beta = 1.
     dif_grnd = 0.
-    alb_neig(:) = 0.
     agesno(:) = 0.
     
@@ -221 +219 @@
        petAcoef, peqAcoef, petBcoef, peqBcoef, &
        ps, u1_lay, v1_lay, &
-       radsol, snow, qsurf, qsol, agesno, &
-       tsoil, &
-       alblw, evap, fluxsens, fluxlat, &
-       tsurf_new, alb_new, dflux_s, dflux_l, pctsrf_sic)
+       radsol, snow, qsurf, qsol, agesno, tsoil, &
+       alb1_new, alb2_new, evap, fluxsens, fluxlat, &
+       tsurf_new, dflux_s, dflux_l, pctsrf_sic)
 
     INCLUDE "indicesol.h"
@@ -258 +255 @@
 !****************************************************************************************
     REAL, DIMENSION(klon), INTENT(OUT)            :: qsurf
-    REAL, DIMENSION(klon), INTENT(OUT)            :: alblw
+    REAL, DIMENSION(klon), INTENT(OUT)            :: alb1_new  ! new albedo in visible SW interval
+    REAL, DIMENSION(klon), INTENT(OUT)            :: alb2_new  ! new albedo in near IR interval
     REAL, DIMENSION(klon), INTENT(OUT)            :: evap, fluxsens, fluxlat
-    REAL, DIMENSION(klon), INTENT(OUT)            :: tsurf_new, alb_new
+    REAL, DIMENSION(klon), INTENT(OUT)            :: tsurf_new
     REAL, DIMENSION(klon), INTENT(OUT)            :: dflux_s, dflux_l      
     REAL, DIMENSION(klon), INTENT(OUT)            :: pctsrf_sic
@@ -282 +280 @@
 
 ! Initialization of output variables
-    alb_new(:) = 0.0
+    alb1_new(:) = 0.0
 
 !****************************************************************************************
@@ -347 +345 @@
     WHERE (snow(1 : knon) .LT. 0.0001) agesno(1 : knon) = 0.
     zfra(1:knon) = MAX(0.0,MIN(1.0,snow(1:knon)/(snow(1:knon)+10.0)))
-    alb_new(1 : knon) = alb_neig(1 : knon) *zfra(1:knon) + & 
+    alb1_new(1 : knon) = alb_neig(1 : knon) *zfra(1:knon) + & 
          0.6 * (1.0-zfra(1:knon))
-    !!      alb_new(1 : knon) = 0.6
+    
+    alb2_new(:) = alb1_new(:)
 
 !
@@ -397 +396 @@
     ENDDO
      
-!!$    z0_new = 0.002
-!!$    z0_new = SQRT(z0_new**2+rugoro**2)
-    alblw(1:knon) = alb_new(1:knon)
         
 !****************************************************************************************