Changeset 2258 for LMDZ5/branches/testing/libf/phylmd/surf_landice_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/surf_landice_mod.F90 (modified) (4 diffs)

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

@@ -17 +17 @@
        ps, u1, v1, rugoro, pctsrf, &
        snow, qsurf, qsol, agesno, &
-       tsoil, z0_new, alb1, alb2, evap, fluxsens, fluxlat, &
+!albedo SB >>>
+!      tsoil, z0_new, alb1, alb2, evap, fluxsens, fluxlat, &
+       tsoil, z0_new, SFRWL, alb_dir, alb_dif, evap, fluxsens, fluxlat, &
+!albedo SB <<<
        tsurf_new, dflux_s, dflux_l, &
        slope, cloudf, &
@@ -80 +83 @@
     REAL, DIMENSION(klon), INTENT(OUT)            :: qsurf
     REAL, DIMENSION(klon), INTENT(OUT)            :: z0_new
-    REAL, DIMENSION(klon), INTENT(OUT)            :: alb1  ! new albedo in visible SW interval
-    REAL, DIMENSION(klon), INTENT(OUT)            :: alb2  ! new albedo in near IR interval
+!albedo SB >>>
+!    REAL, DIMENSION(klon), INTENT(OUT)            :: alb1  ! new albedo in visible SW interval
+!    REAL, DIMENSION(klon), INTENT(OUT)            :: alb2  ! new albedo in near IR interval
+    REAL, DIMENSION(6), INTENT(IN)              ::SFRWL
+    REAL, DIMENSION(klon,nsw), INTENT(OUT)        ::alb_dir,alb_dif
+!albedo SB <<<
     REAL, DIMENSION(klon), INTENT(OUT)            :: evap, fluxsens, fluxlat
     REAL, DIMENSION(klon), INTENT(OUT)            :: tsurf_new
@@ -116 +123 @@
     CHARACTER (len = 20)                      :: modname = 'surf_landice'
     CHARACTER (len = 80)                      :: abort_message
+
+!albedo SB >>>
+    real,dimension(klon) :: alb1,alb2
+!albedo SB <<<
 
 ! End definition
@@ -315 +326 @@
 
 
+!albedo SB >>>
+     select case(NSW)
+     case(2)
+       alb_dir(1:knon,1)=alb1(1:knon)
+       alb_dir(1:knon,2)=alb2(1:knon)
+     case(4)
+       alb_dir(1:knon,1)=alb1(1:knon)
+       alb_dir(1:knon,2)=alb2(1:knon)
+       alb_dir(1:knon,3)=alb2(1:knon)
+       alb_dir(1:knon,4)=alb2(1:knon)
+     case(6)
+       alb_dir(1:knon,1)=alb1(1:knon)
+       alb_dir(1:knon,2)=alb1(1:knon)
+       alb_dir(1:knon,3)=alb1(1:knon)
+       alb_dir(1:knon,4)=alb2(1:knon)
+       alb_dir(1:knon,5)=alb2(1:knon)
+       alb_dir(1:knon,6)=alb2(1:knon)
+     end select
+alb_dif=alb_dir
+!albedo SB <<<
+
+
+
+
   END SUBROUTINE surf_landice
 !