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 5930 for trunk/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90 – NEMO

Ignore:
Timestamp:
2015-11-26T17:07:10+01:00 (8 years ago)
Author:
jchanut
Message:

#1620 Merge free surface simplification into trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

    r5836 r5930  
    106106      ssha(:,:) = (  sshb(:,:) - z2dt * ( zcoef * ( emp_b(:,:) + emp(:,:) ) + zhdiv(:,:) )  ) * ssmask(:,:) 
    107107 
    108 #if ! defined key_dynspg_ts 
    109       ! These lines are not necessary with time splitting since 
    110       ! boundary condition on sea level is set during ts loop 
     108      IF ( .NOT.ln_dynspg_ts ) THEN 
     109         ! These lines are not necessary with time splitting since 
     110         ! boundary condition on sea level is set during ts loop 
    111111# if defined key_agrif 
    112       CALL agrif_ssh( kt ) 
     112         CALL agrif_ssh( kt ) 
    113113# endif 
    114114# if defined key_bdy 
    115       IF( lk_bdy ) THEN 
    116          CALL lbc_lnk( ssha, 'T', 1. )    ! Not sure that's necessary 
    117          CALL bdy_ssh( ssha )             ! Duplicate sea level across open boundaries 
    118       ENDIF 
     115         IF( lk_bdy ) THEN 
     116            CALL lbc_lnk( ssha, 'T', 1. )    ! Not sure that's necessary 
     117            CALL bdy_ssh( ssha )             ! Duplicate sea level across open boundaries 
     118         ENDIF 
    119119# endif 
    120 #endif 
     120      ENDIF 
    121121 
    122122#if defined key_asminc 
     
    250250      ENDIF 
    251251 
    252 # if defined key_dynspg_ts 
    253       IF( ( neuler == 0 .AND. kt == nit000 ) .OR. ln_bt_fw ) THEN    !** Euler time-stepping: no filter 
    254 # else 
    255       IF ( neuler == 0 .AND. kt == nit000 ) THEN   !** Euler time-stepping at first time-step : no filter 
    256 #endif 
     252      IF( ( neuler == 0 .AND. kt == nit000 ) .OR. ( ln_bt_fw .AND. ln_dynspg_ts ) ) THEN  
     253                                                   !** Euler time-stepping: no filter 
    257254         sshb(:,:) = sshn(:,:)                           ! before <-- now 
    258255         sshn(:,:) = ssha(:,:)                           ! now    <-- after  (before already = now) 
Note: See TracChangeset for help on using the changeset viewer.