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

Timestamp:

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

Author:

Laurent Fairhead

Message:

Modifications sur l'albedo JG
LF

File:

• LMDZ4/trunk/libf/phylmd/surf_seaice_mod.F90 (modified) (7 diffs)

LMDZ4/trunk/libf/phylmd/surf_seaice_mod.F90

@@ -16 +16 @@
 !
   SUBROUTINE surf_seaice( & 
-       rlon, rlat, sollw, albedo, &
-       fder, &
+       rlon, rlat, swnet, lwnet, alb1, fder, &
        itime, dtime, jour, knon, knindex, &
-       debut, lafin, swdown, &
+       debut, lafin, &
        tsurf, p1lay, tq_cdrag, precip_rain, precip_snow, temp_air, spechum, &
        petAcoef, peqAcoef, petBcoef, peqBcoef, &
        ps, u1_lay, v1_lay, rugoro, pctsrf, &
-       radsol, snow, qsurf, qsol, agesno, &
-       tsoil, z0_new, alblw, evap, fluxsens, fluxlat, &
-       tsurf_new, alb_new, dflux_s, dflux_l, pctsrf_sic)
+       snow, qsurf, qsol, agesno, tsoil, &
+       z0_new, alb1_new, alb2_new, evap, fluxsens, fluxlat, &
+       tsurf_new, dflux_s, dflux_l, pctsrf_sic)
 !
 ! This subroutine will make a call to ocean_XXX_ice according to the ocean mode (force, 
@@ -41 +40 @@
     REAL, INTENT(IN)                         :: dtime
     REAL, DIMENSION(klon), INTENT(IN)        :: rlon, rlat
-    REAL, DIMENSION(klon), INTENT(IN)        :: sollw
-    REAL, DIMENSION(klon), INTENT(IN)        :: albedo
+    REAL, DIMENSION(klon), INTENT(IN)        :: swnet  ! net shortwave radiation at surface  
+    REAL, DIMENSION(klon), INTENT(IN)        :: lwnet  ! net longwave radiation at surface  
+    REAL, DIMENSION(klon), INTENT(IN)        :: alb1   ! albedo in visible SW interval
     REAL, DIMENSION(klon), INTENT(IN)        :: fder
-    REAL, DIMENSION(klon), INTENT(IN)        :: swdown    
     REAL, DIMENSION(klon), INTENT(IN)        :: tsurf
     REAL, DIMENSION(klon), INTENT(IN)        :: p1lay
@@ -59 +58 @@
 ! In/Output arguments
 !****************************************************************************************
-    REAL, DIMENSION(klon), INTENT(INOUT)          :: radsol
     REAL, DIMENSION(klon), INTENT(INOUT)          :: snow, qsurf, qsol
     REAL, DIMENSION(klon), INTENT(INOUT)          :: agesno
@@ -67 +65 @@
 !****************************************************************************************
     REAL, DIMENSION(klon), INTENT(OUT)       :: z0_new
-    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
 
+! Local arguments
+!****************************************************************************************
+    REAL, DIMENSION(klon)  :: radsol
+!
 ! End definitions
 !****************************************************************************************
+
+
+!****************************************************************************************
+! Calculate total net radiance at surface
+!
+!****************************************************************************************
+    radsol(:) = 0.0
+    radsol(1:knon) = swnet(1:knon) + lwnet(1:knon)
 
 !****************************************************************************************
@@ -83 +94 @@
     CASE('couple')
        CALL ocean_cpl_ice( &
-          rlon, rlat, sollw, albedo, & 
+          rlon, rlat, swnet, lwnet, alb1, & 
           fder, & 
           itime, dtime, knon, knindex, &
           lafin,&
-          swdown, &
           p1lay, tq_cdrag, precip_rain, precip_snow, temp_air, spechum,&
           petAcoef, peqAcoef, petBcoef, peqBcoef, &
           ps, u1_lay, v1_lay, pctsrf, &
           radsol, snow, qsurf, &
-          alblw, evap, fluxsens, fluxlat, &
-          tsurf_new, alb_new, dflux_s, dflux_l, pctsrf_sic)
+          alb1_new, alb2_new, evap, fluxsens, fluxlat, &
+          tsurf_new, dflux_s, dflux_l, pctsrf_sic)
 
     CASE('slab')
@@ -102 +112 @@
           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)
     
     CASE('force')
@@ -113 +122 @@
           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, qsol, agesno, tsoil, &
+          qsurf, alb1_new, alb2_new, evap, fluxsens, fluxlat, &
+          tsurf_new, dflux_s, dflux_l, pctsrf_sic)
     END SELECT