Changeset 7030 for branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

Timestamp:

2016-10-14T14:24:48+02:00 (8 years ago)

Author:

acc

Message:

Branch dev_r6393_NOC_WAD. Added over-topping test case to WAD_TEST_CASES (jp_cfg=6). Also removed last references to redundant wduflt and wdvflt arrays

File:

• branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (1 diff)

branches/2016/dev_r6393_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -529 +529 @@
                         zdta(:,jpjdta-2:jpjdta) = -4._wp
                         !                                    ! ===========================
+                     CASE ( 6 )                              ! WAD 6 configuration
+                        !                                    ! ====================
+                        !
+                        IF(lwp) WRITE(numout,*)
+                        IF(lwp) WRITE(numout,*) 'zgr_bat : Parabolic channel with gaussian ridge'
+                        IF(lwp) WRITE(numout,*) '~~~~~~~~~~'
+                        !
+                        DO ji = 1, jpidta
+                          zi = MAX(1.0-FLOAT((ji-25)**2)/484.0, -2.0 )
+                          zj = 0.95*MAX(EXP(-1.0*FLOAT((ji-25)**2)/32.0) , 0.0 )
+                          zdta(ji,:) = MAX(rn_bathy*(zi-zj), -2.0)
+                          IF(lwp)write(numout,*) 'ZDTA ',ji,zi,zdta(ji,1)
+                        END DO
+                        zdta(1:2,:) = -4._wp
+                        zdta(jpidta-1:jpidta,:) = -4._wp
+                        zdta(:,1) = -4._wp
+                        zdta(:,jpjdta) = -4._wp
+                        zdta(:,1:3) = -4._wp
+                        zdta(:,jpjdta-2:jpjdta) = -4._wp
+                        !                                    ! ===========================
                      CASE DEFAULT
                         !                                    ! ===========================