Changeset 914 for trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.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_2.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

@@ -46 +46 @@
    USE in_out_manager  ! I/O manager
    USE prtctl          ! Print control
+   USE ocfzpt          ! ocean freezing point
 
    IMPLICIT NONE
@@ -89 +90 @@
       !!
       INTEGER  ::   ji, jj   ! dummy loop indices
+      REAL(wp) ::   zinda    
       REAL(wp), DIMENSION(jpi,jpj,1) ::   alb_ice_os   ! albedo of the ice under overcast sky
       REAL(wp), DIMENSION(jpi,jpj,1) ::   alb_ice_cs   ! albedo of ice under clear sky
@@ -152 +154 @@
                &                               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( zsist , ui_ice , vi_ice   , alb_ice_cs ,                         &