salt_uptake_read.f90 Source File

      subroutine salt_uptake_read
      !read in specified salt uptake mass (kg/ha)
      
      use basin_module
      use input_file_module
      use climate_module
      use time_module
      use maximum_data_module
      use constituent_mass_module
      use hydrograph_module
      use salt_module
      use maximum_data_module
      
      implicit none
      
      character (len=80) :: file = "" !           |filename
      character (len=80) :: header = "" !           |header of file
      character (len=30) :: name = "" !           |header of file
      character (len=30) :: plnt_typ = ""!           |header of file
      logical :: i_exist              !none       |check to determine if file exists
      integer :: eof = 0              !           |end of file
      integer :: i = 0                !           |counters
      integer :: j = 0                !           |counters
      
      
      eof = 0

      !only proceed if there are salt ions in the simulation
      if(cs_db%num_salts > 0) then
      
      salt_uptake_on = 0
      inquire (file='salt_uptake',exist=i_exist)
      if(i_exist) then
        
        !turn on salt uptake feature
        salt_uptake_on = 1
        
        !open and read first lines
        open(5054,file='salt_uptake')
        read(5054,*) header
        read(5054,*) header
        read(5054,*) header
        
        !allocate array
        allocate (salt_uptake_kg(db_mx%plantparm,cs_db%num_salts), source = 0.)
        salt_uptake_kg = 0.
        
        !loop through the plant communities
        do i=1,db_mx%plantparm
          read(5054,*) name,(salt_uptake_kg(i,j),j=1,cs_db%num_salts)
        enddo
        
        close(5054) 
      endif
      
      endif

      return
      end subroutine salt_uptake_read