Changeset 888 for LMDZ4/trunk/libf/phylmd/ocean_cpl_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/ocean_cpl_mod.F90 (modified) (10 diffs)

LMDZ4/trunk/libf/phylmd/ocean_cpl_mod.F90

@@ -65 +65 @@
 !
   SUBROUTINE ocean_cpl_noice( &
-       sollw, albedo, &
+       swnet, lwnet, alb1, &
        windsp, &
        fder_old, &
        itime, dtime, knon, knindex, &
-       swdown, &
        p1lay, tq_cdrag, precip_rain, precip_snow, temp_air, spechum, &
        petAcoef, peqAcoef, petBcoef, peqBcoef, &
        ps, u1_lay, v1_lay, pctsrf_in, &
-       radsol, snow, qsurf, agesno, &
-       evap, fluxsens, fluxlat, &
+       radsol, snow, agesno, &
+       qsurf, evap, fluxsens, fluxlat, &
        tsurf_new, dflux_s, dflux_l, pctsrf_oce)
 !
@@ -89 +88 @@
     INTEGER, DIMENSION(klon), INTENT(IN)     :: knindex
     REAL, INTENT(IN)                         :: dtime
-    REAL, DIMENSION(klon), INTENT(IN)        :: sollw
-    REAL, DIMENSION(klon), INTENT(IN)        :: albedo
+    REAL, DIMENSION(klon), INTENT(IN)        :: swnet
+    REAL, DIMENSION(klon), INTENT(IN)        :: lwnet
+    REAL, DIMENSION(klon), INTENT(IN)        :: alb1   ! albedo in visible SW interval
     REAL, DIMENSION(klon), INTENT(IN)        :: windsp
     REAL, DIMENSION(klon), INTENT(IN)        :: fder_old
-    REAL, DIMENSION(klon), INTENT(IN)        :: swdown    
     REAL, DIMENSION(klon), INTENT(IN)        :: p1lay
     REAL, DIMENSION(klon), INTENT(IN)        :: tq_cdrag
@@ -199 +198 @@
 
     CALL cpl_send_ocean_fields(itime, knon, knindex, &
-         swdown, sollw, fluxlat, fluxsens, &
-         precip_rain, precip_snow, evap, tsurf_new, fder_new, albedo, taux, tauy, windsp)
+         swnet, lwnet, fluxlat, fluxsens, &
+         precip_rain, precip_snow, evap, tsurf_new, fder_new, alb1, taux, tauy, windsp)
     
 
@@ -208 +207 @@
 !
   SUBROUTINE ocean_cpl_ice( &
-       rlon, rlat, sollw, albedo, &
+       rlon, rlat, swnet, lwnet, alb1, &
        fder_old, &
        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_in, &
        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)
 !
 ! This subroutine treats the ocean where there is ice. The subroutine first receives 
@@ -234 +232 @@
     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
+    REAL, DIMENSION(klon), INTENT(IN)        :: lwnet
+    REAL, DIMENSION(klon), INTENT(IN)        :: alb1   ! albedo in visible SW interval
     REAL, DIMENSION(klon), INTENT(IN)        :: fder_old
-    REAL, DIMENSION(klon), INTENT(IN)        :: swdown
     REAL, DIMENSION(klon), INTENT(IN)        :: p1lay
     REAL, DIMENSION(klon), INTENT(IN)        :: tq_cdrag
@@ -256 +254 @@
 !****************************************************************************************
     REAL, DIMENSION(klon), INTENT(OUT)       :: qsurf
-    REAL, DIMENSION(klon), INTENT(OUT)       :: alblw
+    REAL, DIMENSION(klon), INTENT(OUT)       :: alb1_new, alb2_new
     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
@@ -271 +269 @@
     REAL, DIMENSION(klon)   :: tsurf_cpl, fder_new
     REAL, DIMENSION(klon)   :: taux, tauy
+    REAL, DIMENSION(klon)   :: alb_cpl
 
 ! End definitions
@@ -278 +277 @@
 
 !****************************************************************************************
-! Receive ocean temperature(tsurf_cpl), albedo(alb_new) and new fraction of 
+! Receive ocean temperature(tsurf_cpl), albedo(alb_cpl) and new fraction of 
 ! seaice(pctsrf_sic) from coupler
 !
@@ -284 +283 @@
 
     CALL cpl_receive_seaice_fields(knon, knindex, &
-         tsurf_cpl, alb_new, pctsrf_sic)
+         tsurf_cpl, alb_cpl, pctsrf_sic)
+
+    alb1_new(1:knon) = alb_cpl(1:knon)
+    alb2_new(1:knon) = alb_cpl(1:knon)    
+
     
 !****************************************************************************************
@@ -342 +345 @@
     CALL cpl_send_seaice_fields(itime, dtime, knon, knindex, &
        pctsrf_in, lafin, rlon, rlat, &
-       swdown, sollw, fluxlat, fluxsens, &
-       precip_rain, precip_snow, evap, tsurf_new, fder_new, albedo, taux, tauy)
+       swnet, lwnet, fluxlat, fluxsens, &
+       precip_rain, precip_snow, evap, tsurf_new, fder_new, alb1, taux, tauy)
  
-
-    alblw(1:knon) = alb_new(1:knon)    
 
   END SUBROUTINE ocean_cpl_ice