New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 13786 – NEMO

Changeset 13786


Ignore:
Timestamp:
2020-11-12T19:13:23+01:00 (4 years ago)
Author:
mathiot
Message:

ticket #1900: cosmetic change

Location:
NEMO/branches/2020/tickets_icb_1900/src/OCE/ICB
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/tickets_icb_1900/src/OCE/ICB/icbdyn.F90

    r13374 r13786  
    216216         ! 
    217217         ! berg reach a new t-cell, but an ocean one 
     218         ! .AND. needed in case berg hit an isf (tmask(ii,ij,1) == 0 and tmask(ii,ij,ikb) /= 0) 
    218219         IF(  tmask(ii,ij,ikb) /= 0._wp .AND. tmask(ii,ij,1) /= 0._wp ) RETURN 
    219220         ! 
  • NEMO/branches/2020/tickets_icb_1900/src/OCE/ICB/icbini.F90

    r13374 r13786  
    7373      ! 
    7474      IF( .NOT. ln_icebergs )   RETURN 
    75       ! 
    76       ! check nn_fsbc 
    77       IF ( ln_icebergs .AND. ln_M2016 .AND. nn_fsbc /= 1) THEN 
    78          IF (lwp) WRITE(numout,*) 'The use of ln_M2016 within the iceberg module is only compatible with nn_fsbc == 1' 
    79          IF (lwp) WRITE(numout,*) 'because it needs the 3d t,s,u,v model data' 
    80          CALL ctl_stop('Compatibility error in the iceberg module') 
    81       END IF 
    8275      ! 
    8376      !                          ! allocate gridded fields 
  • NEMO/branches/2020/tickets_icb_1900/src/OCE/ICB/icbthm.F90

    r13359 r13786  
    102102         zM   = pt%mass 
    103103         zT   = pt%thickness                               ! total thickness 
    104          zD   = rho_berg_1_oce * zT          ! draught (keel depth) 
    105        ! F   = zT - D ! freeboard 
     104         zD   = rho_berg_1_oce * zT                        ! draught (keel depth) 
    106105         zW   = pt%width 
    107106         zL   = pt%length 
     
    125124            CALL icb_utl_getkb( pt%kb, ze3t, zD ) 
    126125 
    127             ikb = MIN(pt%kb,mbkt(ii,ij)) 
    128             ztb = ztoce(ikb)                                                ! basal temperature 
     126            ikb = MIN(pt%kb,mbkt(ii,ij))                             ! limit pt%kb by mbkt  
     127                                                                     ! => bottom temperature used to fill ztoce(mbkt:jpk) 
     128            ztb = ztoce(ikb)                                         ! basal temperature 
    129129            zub = zuoce(ikb) 
    130130            zvb = zvoce(ikb) 
  • NEMO/branches/2020/tickets_icb_1900/src/OCE/ICB/icbutl.F90

    r13374 r13786  
    9595      CALL lbc_lnk_icb( 'icbutl', ssu_e, 'U', -1._wp, 1, 1 ) 
    9696      CALL lbc_lnk_icb( 'icbutl', ssv_e, 'V', -1._wp, 1, 1 ) 
    97       CALL lbc_lnk_icb( 'icbutl', ua_e, 'U', -1._wp, 1, 1 ) 
    98       CALL lbc_lnk_icb( 'icbutl', va_e, 'V', -1._wp, 1, 1 ) 
     97      CALL lbc_lnk_icb( 'icbutl', ua_e , 'U', -1._wp, 1, 1 ) 
     98      CALL lbc_lnk_icb( 'icbutl', va_e , 'V', -1._wp, 1, 1 ) 
    9999#if defined key_si3 
    100100      hi_e(1:jpi, 1:jpj) = hm_i (:,:)   
     
    274274         ! since we're looking for four T points containing quadrant we're in of  
    275275         ! current T cell 
    276          kii = MAX(0, INT( pi     )) 
    277          kij = MAX(0, INT( pj     ))    ! T-point 
     276         kii = MAX(0, INT( pi        )) 
     277         kij = MAX(0, INT( pj        ))    ! T-point 
    278278         zwi = pi - REAL(kii,wp) 
    279279         zwj = pj - REAL(kij,wp) 
    280280      CASE ( 'U' ) 
    281281         kii = MAX(0, INT( pi-0.5_wp )) 
    282          kij = MAX(0, INT( pj     ))    ! U-point 
     282         kij = MAX(0, INT( pj        ))    ! U-point 
    283283         zwi = pi - 0.5_wp - REAL(kii,wp) 
    284284         zwj = pj - REAL(kij,wp) 
    285285      CASE ( 'V' ) 
    286          kii = MAX(0, INT( pi     )) 
     286         kii = MAX(0, INT( pi        )) 
    287287         kij = MAX(0, INT( pj-0.5_wp ))    ! V-point 
    288288         zwi = pi - REAL(kii,wp) 
     
    896896 
    897897   SUBROUTINE test_icb_utl_getkb 
     898      !!---------------------------------------------------------------------- 
     899      !!                 ***  FUNCTION test_icb_utl_getkb  *** 
     900      !! 
     901      !! ** Purpose : Test routine icb_utl_getkb, icb_utl_zavg 
     902      !! ** Methode : Call each subroutine with specific input data 
     903      !!              What should be the output is easy to determined and check  
     904      !!              if NEMO return the correct answer. 
     905      !! ** Comments : not called, if needed a CALL test_icb_utl_getkb need to be added in icb_step 
     906      !!---------------------------------------------------------------------- 
    898907      INTEGER :: ikb 
    899908      REAL(wp) :: zD, zout 
Note: See TracChangeset for help on using the changeset viewer.