Changeset 13606


Ignore:
Timestamp:
2020-10-14T18:23:00+02:00 (3 months ago)
Author:
techene
Message:

#2385 adapted for SWE

Location:
NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/src/OCE/DOM
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/src/OCE/DOM/domain.F90

    r13427 r13606  
    1515   !!            3.7  !  2015-11  (G. Madec, A. Coward)  time varying zgr by default 
    1616   !!            4.0  !  2016-10  (G. Madec, S. Flavoni)  domain configuration / user defined interface 
    17    !!            4.x  ! 2020-02  (G. Madec, S. Techene) introduce ssh to h0 ratio 
     17   !!            4.0  ! 2020-02  (G. Madec, S. Techene) introduce ssh to h0 ratio 
    1818   !!---------------------------------------------------------------------- 
    1919    
     
    153153      hv_0(:,:) = 0._wp 
    154154      hf_0(:,:) = 0._wp 
    155       DO jk = 1, jpk 
     155      DO jk = 1, jpkm1 
    156156         ht_0(:,:) = ht_0(:,:) + e3t_0(:,:,jk) * tmask(:,:,jk) 
    157157         hu_0(:,:) = hu_0(:,:) + e3u_0(:,:,jk) * umask(:,:,jk) 
     
    160160      END DO 
    161161      ! 
     162      IF( lk_SWE ) THEN      ! SWE case redefine hf_0 
     163         hf_0(:,:) = hf_0(:,:) + e3f_0(:,:,1) * ssfmask(:,:) 
     164      ENDIF 
     165      ! 
    162166      r1_ht_0(:,:) = ssmask (:,:) / ( ht_0(:,:) + 1._wp -  ssmask (:,:) ) 
    163167      r1_hu_0(:,:) = ssumask(:,:) / ( hu_0(:,:) + 1._wp -  ssumask(:,:) ) 
    164168      r1_hv_0(:,:) = ssvmask(:,:) / ( hv_0(:,:) + 1._wp -  ssvmask(:,:) ) 
    165169      r1_hf_0(:,:) = ssfmask(:,:) / ( hf_0(:,:) + 1._wp -  ssfmask(:,:) ) 
    166  
    167170      ! 
    168171#if defined key_qco 
     
    204207      ELSE                       != time varying : initialize before/now/after variables 
    205208         ! 
    206          IF( .NOT.l_offline )  CALL dom_vvl_init( Kbb, Kmm, Kaa ) 
     209         IF( .NOT.l_offline )   CALL dom_vvl_init( Kbb, Kmm, Kaa ) 
    207210         ! 
    208211      ENDIF 
  • NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/src/OCE/DOM/dommsk.F90

    r13305 r13606  
    182182      ssvmask(:,:) = MAXVAL( vmask(:,:,:), DIM=3 ) 
    183183      ssfmask(:,:) = MAXVAL( fmask(:,:,:), DIM=3 ) 
     184      IF( lk_SWE ) THEN      ! Shallow Water Eq. case : redefine ssfmask 
     185         DO_2D( 0,0 , 0,0 ) 
     186            ssfmask(ji,jj) = MAX(  ssmask(ji,jj+1), ssmask(ji+1,jj+1),  &  
     187               &                   ssmask(ji,jj  ), ssmask(ji+1,jj  )   ) 
     188         END_2D 
     189         CALL lbc_lnk( 'dommsk', ssfmask, 'F', 1.0_wp ) 
     190      ENDIF 
    184191 
    185192      ! Interior domain mask  (used for global sum) 
Note: See TracChangeset for help on using the changeset viewer.