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 5066 for branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90 – NEMO

Ignore:
Timestamp:
2015-02-06T17:02:20+01:00 (9 years ago)
Author:
rfurner
Message:

added current state of wetting and drying code to test...note it does not work

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

    r4486 r5066  
    99   !!             -   !  2010-09  (D.Storkey and E.O'Dea) bug fixes for BDY module 
    1010   !!            3.3  !  2011-10  (M. Leclair) split former ssh_wzv routine and remove all vvl related work 
     11   !!            3.?  !  2014-09  (H. Liu) add wetting and drying   
    1112   !!---------------------------------------------------------------------- 
    1213 
     
    4243   USE wrk_nemo        ! Memory Allocation 
    4344   USE timing          ! Timing 
     45   USE wadlmt          ! Wetting/Drying flux limting 
    4446 
    4547   IMPLICIT NONE 
     
    9496      ENDIF 
    9597      ! 
    96       CALL div_cur( kt )                              ! Horizontal divergence & Relative vorticity 
    97       ! 
    9898      z2dt = 2._wp * rdt                              ! set time step size (Euler/Leapfrog) 
    9999      IF( neuler == 0 .AND. kt == nit000 )   z2dt = rdt 
    100  
     100      ! 
     101      z1_rau0 = 0.5_wp * r1_rau0 
     102      ! 
     103      IF(ln_wd) CALL wad_lmt(sshb, z1_rau0 * (emp_b(:,:) + emp(:,:)), z2dt) 
     104 
     105      CALL div_cur( kt )                              ! Horizontal divergence & Relative vorticity 
     106      ! 
    101107      !                                           !------------------------------! 
    102108      !                                           !   After Sea Surface Height   ! 
     
    110116      ! compute the vertical velocity which can be used to compute the non-linear terms of the momentum equations. 
    111117      !  
    112       z1_rau0 = 0.5_wp * r1_rau0 
    113118      ssha(:,:) = (  sshb(:,:) - z2dt * ( z1_rau0 * ( emp_b(:,:) + emp(:,:) ) + zhdiv(:,:) )  ) * tmask(:,:,1) 
    114119 
Note: See TracChangeset for help on using the changeset viewer.