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 7022 for branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 – NEMO

Ignore:
Timestamp:
2016-10-11T18:18:38+02:00 (7 years ago)
Author:
acc
Message:

Branch dev_r6393_NOC_WAD. WAD TEST CASES 4 (parabolic bowl) and 5 (double slope) now initialised correctly and tested

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r6986 r7022  
    473473                        ! 
    474474                        DO ji = 1, jpidta 
    475                           zi = MAX(1.0-FLOAT((ji-25)**2)/484.0, 0.0 ) 
     475                          zi = MAX(1.0-FLOAT((ji-25)**2)/484.0, -2.0 ) 
    476476                        DO jj = 1, jpjdta 
    477                           zj = MAX(1.0-FLOAT((jj-17)**2)/196.0, 0.0 ) 
    478                           zdta(ji,jj) = MAX(rn_bathy*zi*zj, 0.0) 
     477                          zj = MAX(1.0-FLOAT((jj-17)**2)/196.0, -2.0 ) 
     478                          zdta(ji,jj) = MAX(rn_bathy*zi*zj, -2.0) 
    479479                        END DO 
    480480                          IF(lwp)write(numout,*) 'ZDTA ',ji,zi,zdta(ji,1) 
     
    521521                          IF(lwp)write(numout,*) 'ZDTA ',ji,zi,zdta(ji,1) 
    522522                        END DO 
     523                        !                                    ! =========================== 
     524                        zdta(1:2,:) = -4._wp 
     525                        zdta(jpidta-1:jpidta,:) = -4._wp 
     526                        zdta(:,1) = -4._wp 
     527                        zdta(:,jpjdta) = -4._wp 
     528                        zdta(:,1:3) = -4._wp 
     529                        zdta(:,jpjdta-2:jpjdta) = -4._wp 
    523530                        !                                    ! =========================== 
    524531                     CASE DEFAULT 
Note: See TracChangeset for help on using the changeset viewer.