Changeset 8404 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

Timestamp:

2017-08-07T15:05:37+02:00 (7 years ago)

Author:

clem

Message:

start changing calls in icestp.F90

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90 (modified) (2 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

@@ -107 +107 @@
       REAL(wp) ::   zqsr             ! New solar flux received by the ocean
       REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zalb_cs, zalb_os     ! 3D workspace
-      REAL(wp), DIMENSION(jpi,jpj)     ::   zalb                 ! 2D workspace
       !!---------------------------------------------------------------------
 
@@ -120 +119 @@
       ENDIF
       
-      ! albedo output
-      zalb(:,:) = 0._wp
-      WHERE     ( at_i_b <= epsi06 )  ;  zalb(:,:) = rn_alb_oce
-      ELSEWHERE                       ;  zalb(:,:) = SUM( alb_ice * a_i_b, dim=3 ) / at_i_b
-      END WHERE
-      IF( iom_use('icealb'  ) )  CALL iom_put( "icealb"   , zalb(:,:) )          ! ice albedo output
-
-      zalb(:,:) = SUM( alb_ice * a_i_b, dim=3 ) + rn_alb_oce * ( 1._wp - at_i_b )      
-      IF( iom_use('albedo'  ) )        CALL iom_put( "albedo"  , zalb(:,:) )           ! surface albedo output
-
-
       DO jj = 1, jpj
          DO ji = 1, jpi