Changeset 8870 for branches/UKMO/ROMS_WAD_7832/NEMOGCM/TOOLS/DOMAINcfg/src/domzgr.f90

Timestamp:

2017-12-01T10:48:12+01:00 (6 years ago)

Author:

deazer

Message:

Changed WAD option names to Iterative and Directional
Removed old Diagnostics
Updated Domain CFG to allow domain generation with ref height for wad cases
Cleaned up TEST_CASES/cfg.txt file (need to not include WAD2 etc)
TEST caaes run ok
SETTE runs OK
AMM15 5 level runs OK

File:

• branches/UKMO/ROMS_WAD_7832/NEMOGCM/TOOLS/DOMAINcfg/src/domzgr.f90 (modified) (1 diff)

branches/UKMO/ROMS_WAD_7832/NEMOGCM/TOOLS/DOMAINcfg/src/domzgr.f90

@@ -1957 +1957 @@
          WRITE(numout,*) '        Bottom cell (Zb) (m) = H*rn_zb_a + rn_zb_b'
       ENDIF
-
+      ! CEOD
+      ! Add on some reference level to make nagative bathy positive to create
+      ! sensible e3t_0's
+! Can used commented code here to force bathy near open bdy not to have any WAD
+! points
+
+!CEOD MAKE points on the perimeter assuming a bondary is here on shallower than
+!10m so that tide can be specified here.
+!      DO jj = 1, jpj
+!         DO ji = 1, jpi
+!            !CEODIF( mig(ji) <= 2 .or. mjg(jj) <= 2  .or. mig(ji) >= jpiglo-2 .or.  mjg(jj) >=jpjglo-2) then
+!            IF( mig(ji) <= 10   .or. mjg(jj) <=10  .or.  mjg(ji) >=jpiglo-10+1 .or. mjg(jj) >=jpjglo-10+1   ) then
+!               IF( bathy(ji,jj) + rn_wd_ref_depth  > 0.0)   then
+!                  bathy(ji,jj) = max(5.0,bathy(ji,jj)) !
+!               ENDIF
+!               IF( bathy(ji,jj)  + rn_wd_ref_depth   <= 0.0)   then
+!                  bathy(ji,jj) =  -rn_wd_ref_depth*2._wp ! Make sure it is land
+!               ENDIF
+!            ENDIF
+!         END DO
+!      END DO
+
+      IF( ln_wd) THEN
+         bathy(:,:) = bathy(:,:) + rn_wd_ref_depth !  Arbitrary reference level for WAD, This needs
+      ELSE
+         bathy(:,:) = bathy(:,:)  !  Arbitrary reference level for WAD, This needs
+      ENDIF
+                                                !  to be accounted for in istate in the actual model and bdys etc.
       hift(:,:) = rn_sbot_min                     ! set the minimum depth for the s-coordinate
       hifu(:,:) = rn_sbot_min