Changeset 8544 for branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG/CS15mini/MY_SRC/usrdef_istate.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/usrdef_istate.F90 (modified) (3 diffs)

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

@@ -14 +14 @@
    !!----------------------------------------------------------------------
    USE par_oce        ! ocean space and time domain
+   USE dom_oce        ! ocean space and time domain
    USE phycst         ! physical constants
    USE wet_dry        ! Wetting and drying
@@ -85 +86 @@
 
 ! subtract the height of z=0 above the geoid (this allows z = 0 to be higher than all points that may become wet)   
-      pssh(:,:) = pssh(:,:) - rn_ssh_ref
+      pssh(:,:) =  -rn_ssh_ref
 
       !
@@ -92 +93 @@
       do jj = 1,jpj
          do ji = 1,jpi
-            IF( ht_wd(ji,jj) + pssh(ji,jj) < rn_wdmin1 ) THEN
-               pssh(ji,jj) = ptmask(ji,jj,1)*( rn_wdmin1 - ht_wd(ji,jj) )
+            IF( ht_0(ji,jj) + pssh(ji,jj)  < rn_wdmin1 ) THEN
+               pssh(ji,jj) = ptmask(ji,jj,1)*( rn_wdmin1 - (ht_0(ji,jj)) )
             ENDIF
          end do