Ignore:
Timestamp:
2020-05-25T08:46:52+02:00 (14 months ago)
Author:
jchanut
Message:

Change place for declaration of ice-ocean drag, #2468

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/MOI/NEMO_4.03_IODRAG/src/OCE/ZDF/zdfdrg.F90

    r12953 r12964  
    3232   USE lib_mpp        ! distributed memory computing 
    3333   USE prtctl         ! Print control 
     34   USE sbc_oce, ONLY: nn_ice  
    3435 
    3536   IMPLICIT NONE 
     
    4647   LOGICAL          ::   ln_loglayer  ! logarithmic drag: Cd = vkarmn/log(z/z0) 
    4748   LOGICAL , PUBLIC ::   ln_drgimp    ! implicit top/bottom friction flag 
    48    LOGICAL , PUBLIC ::   ln_drgice_imp ! implicit ice-ocean drag 
    49  
     49   LOGICAL , PUBLIC ::   ln_drgice_imp ! implicit ice-ocean drag  
    5050   !                                 !!* Namelist namdrg_top & _bot: TOP or BOTTOM coefficient namelist * 
    5151   REAL(wp)         ::   rn_Cd0       !: drag coefficient                                           [ - ] 
     
    7272   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC ::   rCd0_top, rCd0_bot   !: precomputed top/bottom drag coeff. at t-point (>0) 
    7373   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC ::   rCdU_top, rCdU_bot   !: top/bottom drag coeff. at t-point (<0)  [m/s] 
    74    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:), PUBLIC ::   rCdU_ice             !: ice-ocean drag (used if ln_drgice_imp=T) 
    7574 
    7675   !! * Substitutions 
     
    270269         &                CALL ctl_stop( 'zdf_drg_init: ln_drgice_imp=T requires ln_drgimp=T' ) 
    271270      ! 
     271      IF ( ln_drgice_imp.AND.( nn_ice /=2 ) ) & 
     272         &  CALL ctl_stop( 'zdf_drg_init: ln_drgice_imp=T requires si3' ) 
     273      ! 
    272274      !                     !==  BOTTOM drag setting  ==!   (applied at seafloor) 
    273275      ! 
     
    281283      IF( ln_isfcav.OR.ln_drgice_imp ) THEN              ! Ocean cavities: top friction setting 
    282284         ALLOCATE( rCdU_top(jpi,jpj) ) 
    283          rCdU_top(:,:) = 0._wp 
    284       ENDIF 
    285       ! 
    286       IF (ln_drgice_imp) ALLOCATE( rCdU_ice(jpi,jpj) ) 
     285      ENDIF 
    287286      ! 
    288287      IF( ln_isfcav ) THEN 
Note: See TracChangeset for help on using the changeset viewer.