Changeset 8544 for branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG/CS15mini/MY_SRC/domvvl.F90

Timestamp:

2017-09-19T17:41:15+02:00 (7 years ago)

Author:

deazer

Message:

Added reference level (above all potential wet points) to avoid negative depth bathymetry as a work around.
Require reference level to be emedded into the domain configuration file
uses ht_0 instead of ht_wd. This si still in the code but should be removed in time.

File:

• branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG/CS15mini/MY_SRC/domvvl.F90 (modified) (1 diff)

branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG/CS15mini/MY_SRC/domvvl.F90

@@ -889 +889 @@
                   vn   (:,:,:)   = vb  (:,:,:)
                ELSEIF( ln_wd .or. ln_rwd   ) THEN ! if not test case
-                  sshb(:,:) = 0.0_wp ! slwa
-                  sshn(:,:) = 0.0_wp ! slwa
+                  sshn(:,:) = -rn_ssh_ref
+                  sshb(:,:) = -rn_ssh_ref
 
                   DO jj = 1, jpj
                      DO ji = 1, jpi
-                        IF( ht_wd(ji,jj) <=  rn_wdmin1 ) THEN ! if total depth is less than min depth
-
-                           sshb(ji,jj) = rn_wdmin1 - ht_wd(ji,jj)       
-                           sshn(ji,jj) = rn_wdmin1 - ht_wd(ji,jj)
-                           ssha(ji,jj) = rn_wdmin1 - ht_wd(ji,jj)
+                        IF( ht_0(ji,jj)-rn_ssh_ref <  rn_wdmin1 ) THEN ! if total depth is less than min depth
+
+                           sshb(ji,jj) = rn_wdmin1 - (ht_0(ji,jj) )
+                           sshn(ji,jj) = rn_wdmin1 - (ht_0(ji,jj) )
+                           ssha(ji,jj) = rn_wdmin1 - (ht_0(ji,jj) )
                         ENDIF
                      ENDDO