Changeset 914


Ignore:
Timestamp:
2008-04-30T16:30:23+02:00 (13 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

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

Legend:

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

    r892 r914  
    270270      REAL(wp) ::   zdtetar, ztvmoyr, zlxins, zchcm, zclcm      !    -         - 
    271271      REAL(wp) ::   zmt1, zmt2, zmt3, ztatm3, ztamr, ztaevbk    !    -         - 
    272       REAL(wp) ::   zsst, ztatm, zcco1, zpatm, zinda            !    -         - 
     272      REAL(wp) ::   zsst, ztatm, zcco1, zpatm                   !    -         - 
    273273      REAL(wp) ::   zrhoa, zev, zes, zeso, zqatm, zevsqr        !    -         - 
    274274      !! 
     
    296296       
    297297      CALL blk_clio_qsr_oce( qsr ) 
    298  
    299       ! CAUTION: ocean shortwave radiation sets to zero if more than 50% of sea-ice !!gm to be removed 
    300       DO jj = 1, jpj 
    301          DO ji = 1, jpi 
    302             zinda    = MAX(  0.e0, SIGN(  1.e0, -( -1.5 - freeze(ji,jj) )  )  ) 
    303             qsr(ji,jj) = zinda * qsr(ji,jj) 
    304          END DO 
    305       END DO 
    306  
    307298 
    308299      !------------------------! 
  • trunk/NEMO/OPA_SRC/SBC/sbcice_lim.F90

    r906 r914  
    5252   USE in_out_manager  ! I/O manager 
    5353   USE prtctl          ! Print control 
     54   USE ocfzpt          ! ocean freezing point 
    5455 
    5556   IMPLICIT NONE 
     
    99100      INTEGER  ::   alert_id           ! number of the current alert 
    100101      REAL(wp) ::   ztmelts            ! ice layer melting point 
    101       INTEGER , DIMENSION(20) ::  numal                     ! number of alerts positive 
     102      REAL(wp) ::   zinda     
     103      INTEGER , DIMENSION(20) ::  numal                  ! number of alerts positive 
    102104      CHARACTER (len=30), DIMENSION(20) ::   alname      ! name of alert 
    103105      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   alb_ice_os   ! albedo of the ice under overcast sky 
     
    159161               &                             tprecip  , sprecip   ,                        & 
    160162               &                             fr1_i0   , fr2_i0    , cl_grid  ) 
     163 
     164            ! CAUTION: ocean shortwave radiation sets to zero if more than 50% of sea-ice !!gm to be removed 
     165            DO jj = 1, jpj 
     166               DO ji = 1, jpi 
     167                  zinda    = MAX(  0.e0, SIGN(  1.e0, -( -1.5 - freeze(ji,jj) )  )  ) 
     168                  qsr(ji,jj) = zinda * qsr(ji,jj) 
     169               END DO 
     170            END DO 
     171 
    161172         CASE( 4 )           ! CORE bulk formulation 
    162173            CALL blk_ice_core( t_su , u_ice , v_ice   , alb_ice_cs,                      & 
  • trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

    r900 r914  
    4646   USE in_out_manager  ! I/O manager 
    4747   USE prtctl          ! Print control 
     48   USE ocfzpt          ! ocean freezing point 
    4849 
    4950   IMPLICIT NONE 
     
    8990      !! 
    9091      INTEGER  ::   ji, jj   ! dummy loop indices 
     92      REAL(wp) ::   zinda     
    9193      REAL(wp), DIMENSION(jpi,jpj,1) ::   alb_ice_os   ! albedo of the ice under overcast sky 
    9294      REAL(wp), DIMENSION(jpi,jpj,1) ::   alb_ice_cs   ! albedo of ice under clear sky 
     
    152154               &                               tprecip   , sprecip    ,                          & 
    153155               &                               fr1_i0    , fr2_i0     , cl_grid  ) 
     156 
     157            ! CAUTION: ocean shortwave radiation sets to zero if more than 50% of sea-ice !!gm to be removed 
     158            DO jj = 1, jpj 
     159               DO ji = 1, jpi 
     160                  zinda    = MAX(  0.e0, SIGN(  1.e0, -( -1.5 - freeze(ji,jj) )  )  ) 
     161                  qsr(ji,jj) = zinda * qsr(ji,jj) 
     162               END DO 
     163            END DO 
     164 
    154165         CASE( 4 )           ! CORE bulk formulation 
    155166            CALL blk_ice_core( zsist , ui_ice , vi_ice   , alb_ice_cs ,                         & 
Note: See TracChangeset for help on using the changeset viewer.