Changeset 914 for trunk/NEMO/OPA_SRC/SBC/sbcice_lim.F90

Timestamp:

2008-04-30T16:30:23+02:00 (16 years ago)

Author:

ctlod

Message:

move the set to zero of qsr(:,:) (when more than 50% of sea-ice) from sbcblk_clio.F90 into sbcice_lim*, see ticket: #131

File:

• trunk/NEMO/OPA_SRC/SBC/sbcice_lim.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcice_lim.F90

@@ -52 +52 @@
    USE in_out_manager  ! I/O manager
    USE prtctl          ! Print control
+   USE ocfzpt          ! ocean freezing point
 
    IMPLICIT NONE
@@ -99 +100 @@
       INTEGER  ::   alert_id           ! number of the current alert
       REAL(wp) ::   ztmelts            ! ice layer melting point
-      INTEGER , DIMENSION(20) ::  numal                     ! number of alerts positive
+      REAL(wp) ::   zinda    
+      INTEGER , DIMENSION(20) ::  numal                  ! number of alerts positive
       CHARACTER (len=30), DIMENSION(20) ::   alname      ! name of alert
       REAL(wp), DIMENSION(jpi,jpj,jpl) ::   alb_ice_os   ! albedo of the ice under overcast sky
@@ -159 +161 @@
                &                             tprecip  , sprecip   ,                        &
                &                             fr1_i0   , fr2_i0    , cl_grid  )
+
+            ! CAUTION: ocean shortwave radiation sets to zero if more than 50% of sea-ice !!gm to be removed
+            DO jj = 1, jpj
+               DO ji = 1, jpi
+                  zinda    = MAX(  0.e0, SIGN(  1.e0, -( -1.5 - freeze(ji,jj) )  )  )
+                  qsr(ji,jj) = zinda * qsr(ji,jj)
+               END DO
+            END DO
+
          CASE( 4 )           ! CORE bulk formulation
             CALL blk_ice_core( t_su , u_ice , v_ice   , alb_ice_cs,                      &