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

Timestamp:

2017-07-11T20:33:18+02:00 (7 years ago)

Author:

clem

Message:

STEP3 (1): clean separation between sea-ice and ocean

File:

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

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

@@ -31 +31 @@
    USE sbc_oce , ONLY : nn_fsbc, ln_ice_embd, sfx, fr_i, qsr_tot, qns, qsr, fmmflx, emp, taum, utau, vtau
    USE sbccpl         ! Surface boundary condition: coupled interface
-   USE albedo         ! albedo parameters
+   USE albedoice      ! albedo parameters
    USE traqsr         ! add penetration of solar flux in the calculation of heat budget
    USE domvvl         ! Variable volume
@@ -48 +48 @@
    PRIVATE
 
-   PUBLIC   lim_sbc_init   ! called by sbc_lim_init
-   PUBLIC   lim_sbc_flx    ! called by sbc_ice_lim
-   PUBLIC   lim_sbc_tau    ! called by sbc_ice_lim
+   PUBLIC   lim_sbc_init   ! called by ice_init
+   PUBLIC   lim_sbc_flx    ! called by ice_stp
+   PUBLIC   lim_sbc_tau    ! called by ice_stp
 
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   utau_oce, vtau_oce   ! air-ocean surface i- & j-stress     [N/m2]
@@ -124 +124 @@
 
       zalb(:,:) = 0._wp
-      WHERE     ( at_i_b <= epsi06 )  ;  zalb(:,:) = 0.066_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 ) + 0.066_wp * ( 1._wp - at_i_b )      
+      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