Add a new option for net sink treatment: if the elem is dry with a net

sink, the vsource is reset to 0.

Author: josephzhang8

cmake/SCHISM.local.build

@@ -20,7 +20,7 @@ set( USE_MARSH OFF CACHE BOOLEAN "Use marsh module")
 
 #   Enable/Disable Modules
 set( USE_SED2D OFF CACHE BOOLEAN "Use 2D sediment module")
-set( USE_WWM ON CACHE BOOLEAN "Use wind-wave module")
+set( USE_WWM OFF CACHE BOOLEAN "Use wind-wave module")
 set( USE_ICE OFF CACHE BOOLEAN "Use ICE module")
 
 #Tracer models

src/Core/schism_glbl.F90

@@ -92,7 +92,7 @@ module schism_glbl
                   &moitn0,mxitn0,nchi,ibtrack_test,nramp_elev,islip,ibtp,inunfl,shorewafo, &
                   &inv_atm_bnd,ieos_type,ieos_pres,iupwind_mom,inter_mom,ishapiro,isav, &
                   &nstep_ice,niter_shap,iunder_deep,ibtrack_openbnd,flag_fib,ielm_transport,max_subcyc, &
-                  &itransport_only
+                  &itransport_only,meth_sink
   integer,save :: ntrs(natrm),nnu_pts(natrm),mnu_pts
   integer,save,dimension(:),allocatable :: iof_hydro,iof_wwm,iof_gen,iof_age,iof_sed,iof_eco, &
      &iof_icm,iof_cos,iof_fib,iof_sed2d,iof_ice,iof_ana,iof_marsh,iof_dvd,iadjust_mass_consv

src/Hydro/schism_init.F90

@@ -184,7 +184,7 @@ subroutine schism_init(iorder,indir,iths,ntime)
      &iwbl,cur_wwm,if_source,nramp_ss,dramp_ss,ieos_type,ieos_pres,eos_a,eos_b,slr_rate, &
      &rho0,shw,isav,nstep_ice,iunder_deep,h1_bcc,h2_bcc,hw_depth,hw_ratio, &
      &ibtrack_openbnd,level_age,vclose_surf_frac,iadjust_mass_consv0,ipre2, &
-     &ielm_transport,max_subcyc,i_hmin_airsea_ex,hmin_airsea_ex,itransport_only
+     &ielm_transport,max_subcyc,i_hmin_airsea_ex,hmin_airsea_ex,itransport_only,meth_sink
 
      namelist /SCHOUT/iof_hydro,iof_wwm,iof_gen,iof_age,iof_sed,iof_eco,iof_icm,iof_cos,iof_fib, &
      &iof_sed2d,iof_ice,iof_ana,iof_marsh,iof_dvd, &
@@ -463,7 +463,7 @@ subroutine schism_init(iorder,indir,iths,ntime)
       ipre2=0
       ielm_transport=0; max_subcyc=10
       hmin_airsea_ex=0.2_rkind
-      itransport_only=0
+      itransport_only=0; meth_sink=0
 
       !Output elev, hvel by detault
       iof_hydro=0; iof_wwm=0; iof_gen=0; iof_age=0; iof_sed=0; iof_eco=0; iof_dvd=0

src/Hydro/schism_step.F90

@@ -1457,6 +1457,11 @@ subroutine schism_step(it)
         endif !impose_net_flux
       endif !isconsv/=0
 
+!...  Option to zero out net sink @dry elem
+      if(meth_sink/=0) then
+        where(idry_e==1.and.vsource<0.d0) vsource=0.d0       
+      endif !meth_sink
+
 !     Calculation of cross-section areas and length for flow b.c.
       if(lflbc) then
         allocate(buf1(nope_global,2),buf2(nope_global,2)); buf1=0.d0;