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 6986 for branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90 – NEMO

Ignore:
Timestamp:
2016-10-04T17:46:55+02:00 (8 years ago)
Author:
acc
Message:

Branch dev_r6393_NOC_WAD. Introduced some WAD test cases, tidied and corrected WAD code

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

    r6604 r6986  
    8888      ENDIF 
    8989      ! 
    90  
     90      z2dt = 2._wp * rdt                          ! set time step size (Euler/Leapfrog) 
     91      IF( neuler == 0 .AND. kt == nit000 )   z2dt = rdt 
     92      zcoef = 0.5_wp * r1_rau0 
     93 
     94      !                                           !------------------------------! 
     95      !                                           !   After Sea Surface Height   ! 
     96      !                                           !------------------------------! 
    9197      IF(ln_wd) THEN 
    92         zcoef = 0.5_wp * r1_rau0 
    93         CALL wad_lmt(sshb, zcoef * (emp_b(:,:) + emp(:,:)), z2dt) 
    94       END IF 
    95  
    96       CALL div_hor( kt )                              ! Horizontal divergence 
    97       ! 
    98       z2dt = 2._wp * rdt                              ! set time step size (Euler/Leapfrog) 
    99       IF( neuler == 0 .AND. kt == nit000 )   z2dt = rdt 
    100  
    101       !                                           !------------------------------! 
    102       !                                           !   After Sea Surface Height   ! 
    103       !                                           !------------------------------! 
     98         CALL wad_lmt(sshb, zcoef * (emp_b(:,:) + emp(:,:)), z2dt) 
     99      ENDIF 
     100 
     101      CALL div_hor( kt )                               ! Horizontal divergence 
     102      ! 
    104103      zhdiv(:,:) = 0._wp 
    105104      DO jk = 1, jpkm1                                 ! Horizontal divergence of barotropic transports 
Note: See TracChangeset for help on using the changeset viewer.