TSSolids.f

cc**********************Start Subroutine for TSS Solids*****************************************
    
      Subroutine TSSOLIDS(mm,iab,jnb,j,k,dz,dzh,dref,sedclass)		!face densities from node densities

      use params
      
 
!>> Check if link is marsh overland flow link, if so use CSS in marsh and add sediment to marsh cumulative sediment flux term
      if (linkt(iab) == 8) then
!>> Only calculate marsh sediment exchange if depth in marsh is greater than 0.3 m (assumed surge conditions)
          if (Eh(j,1)-Bedm(j) > 0.3) then
!>> calculate average concentration at interface between marshes (directionality of flow is needed to apply upwind factor correctly)
              if (Q(iab,1) >= 0.0) then
                  Cssface = CSSh(jus(abs(icc(j,k))),1,sedclass)
              else    
                  Cssface = CSSh(jds(abs(icc(j,k))),1,sedclass)
              endif
!>> add sediment flux from this marsh link to the cumulative sediment flux for marsh in this comparment (QSsumh)
              QSSumh(sedclass) = QSsumh(sedclass)
     &            + sicc(j,k)*(Q(abs(icc(j,k)),1))*Cssface
          endif
          
!>> If link is not marsh overland flow link, use CSS in open water and add sediment flux to open water cumulative sediment flux term
      else
!>> calculate concentration at interface between compartments (directionality of flow is needed)
          if (Q(iab,1) >= 0.0) then
              Cssface = CSS(jus(abs(icc(j,k))),1,sedclass)
          else
              Cssface = CSS(jds(abs(icc(j,k))),1,sedclass)
          endif
          
          !! update cssface for Atchafalya diversion links to have a SWR of 0.5 of sand - MP project # 03b.DI.04
          !if (iab == 3859) then
          !    Cssface_new = 0.5*Cssface
          !    Cssface = Cssface_new
          !endif
          !
          
!>> Calculate incoming flux for sand
!>> van Rijn suspended sand flux (in g/sec) entering compartment from current link inflows and their respective van Rijn equilibrium flow CSS from previous timestep
!>> sand flux out of compartment is set by equilibrium rate in vanRijn subroutine and this will only collect incoming sand fluxes
          if (sedclass == 1) then
              QSsum(sedclass) =  QSsum(sedclass) 
     &            + sicc(j,k)*Q(abs(icc(j,k)),1)*Cssface
              
!>> Overland marsh links do not transport sand - cancels out previous calculation for marsh links         
              QSsumh(sedclass) = 0.0 ! do not transport sand across marsh links
          else
!>> add sediment flux from this link to the cumulative sediment flux for the compartment (QSsum)			    
              QSSum(sedclass) = QSSum(sedclass) 
     &            + sicc(j,k)*(Q(abs(icc(j,k)),1))*Cssface				!face exchange
     &            + fe*EAOL(iab)*(CSS(j,1,sedclass)-CSS(jnb,1,sedclass))			!Diffusion
              
          endif
      endif
      
              
      
	
      return
	end

c***********************End Subroutine for TSS SOLIDS*******************************************