Changeset 1109 for trunk/NEMO

Timestamp:

2008-06-13T11:59:54+02:00 (16 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:

• sbcice_lim.F90 (modified) (2 diffs)
• sbcice_lim_2.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/SBC/sbcice_lim.F90

@@ -103 +103 @@
       INTEGER, INTENT(in) ::   kico    ! ice-ocean stress treatment
       !!
+      INTEGER  ::   jl                 ! loop index
       REAL(wp) ::   zcoef              ! temporary scalar
       REAL(wp), DIMENSION(jpi,jpj,jpl) ::   alb_ice_os   ! albedo of the ice under overcast sky
@@ -134 +135 @@
          CALL albedo_ice( t_su, ht_i, ht_s, alb_ice_cs, alb_ice_os )  ! ... ice albedo
 
+         DO jl = 1, jpl
+            t_su(:,:,jl) = t_su(:,:,jl) +  rt0 * ( 1. - tmask(:,:,1) )
+         END DO
                                                      ! Bulk formulea - provides the following fields:
          ! utaui_ice, vtaui_ice : surface ice stress                     (U- & V-points)   [N/m2]

trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

@@ -127 +127 @@
          tfu(:,:) = tfreez( sss_m ) +  rt0 
 
-         zsist (:,:,1) = sist (:,:)
-         zhicif(:,:,1) = hicif(:,:)   ;   zhsnif(:,:,1) = hsnif(:,:)
+         zsist (:,:,1) = sist (:,:) + rt0 * ( 1. - tmask(:,:,1) )
+         zhicif(:,:,1) = hicif(:,:)
+         zhsnif(:,:,1) = hsnif(:,:)
 
          ! ... ice albedo