Changeset 1109


Ignore:
Timestamp:
2008-06-13T11:59:54+02:00 (12 years ago)
Author:
ctlod
Message:

trunk: ensure that surface sea-ice temperature is set to 273 degK over land area, see ticket: #199

Location:
trunk/NEMO/OPA_SRC/SBC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/SBC/sbcice_lim.F90

    r1055 r1109  
    103103      INTEGER, INTENT(in) ::   kico    ! ice-ocean stress treatment 
    104104      !! 
     105      INTEGER  ::   jl                 ! loop index 
    105106      REAL(wp) ::   zcoef              ! temporary scalar 
    106107      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   alb_ice_os   ! albedo of the ice under overcast sky 
     
    134135         CALL albedo_ice( t_su, ht_i, ht_s, alb_ice_cs, alb_ice_os )  ! ... ice albedo 
    135136 
     137         DO jl = 1, jpl 
     138            t_su(:,:,jl) = t_su(:,:,jl) +  rt0 * ( 1. - tmask(:,:,1) ) 
     139         END DO 
    136140                                                     ! Bulk formulea - provides the following fields: 
    137141         ! utaui_ice, vtaui_ice : surface ice stress                     (U- & V-points)   [N/m2] 
  • trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

    r1055 r1109  
    127127         tfu(:,:) = tfreez( sss_m ) +  rt0  
    128128 
    129          zsist (:,:,1) = sist (:,:) 
    130          zhicif(:,:,1) = hicif(:,:)   ;   zhsnif(:,:,1) = hsnif(:,:) 
     129         zsist (:,:,1) = sist (:,:) + rt0 * ( 1. - tmask(:,:,1) ) 
     130         zhicif(:,:,1) = hicif(:,:) 
     131         zhsnif(:,:,1) = hsnif(:,:) 
    131132 
    132133         ! ... ice albedo 
Note: See TracChangeset for help on using the changeset viewer.