salt_urban_read.f90 Source File

      subroutine salt_urban_read !rtb salt

      use maximum_data_module
      use urban_data_module
      use constituent_mass_module
      use salt_module
      
      implicit none
      
      character (len=80) :: titldum = ""!           |title of file
      character (len=80) :: header = "" !           |header of file
      character (len=16) :: urb_type = ""!           |urban land use type
      logical :: i_exist              !none       |check to determine if file exists
      integer :: eof = 0              !           |end of file
      integer :: imax = 0.            !none       |determine max number for array (imax) and total number in file
      integer :: itype = 0            !           |counter
      integer :: iu = 0               !none       |counter 
      integer :: isalt = 0            !none       |salt ion counter
      
      
      !only proceed if there are salt ions in the simulation
      if(cs_db%num_salts > 0) then
      
      inquire (file='salt_urban',exist=i_exist)
      if(i_exist) then
      
        !open file and read first two header lines
        open(5054,file='salt_urban')
        read(5054,*) header
        read(5054,*) header
      
        !read in salt ion concentration for each urban land use type  
        !determine number of land use types listed
        eof = 0
        imax = 0
        do while (eof == 0)
          read (5054,*,iostat=eof) titldum
          if (eof < 0) exit
          imax = imax + 1
        enddo

        !allocate urban salt arrays
        allocate (salt_urban_conc(imax,cs_db%num_salts), source = 0.)
        salt_urban_conc = 0.
        
        !loop through each land use type; verify match with types listed in urban.urb        
        rewind(5054)
        read(5054,*) header
        read(5054,*) header
        do itype=1,imax
          read(5054,*) urb_type
          !see if there is a match (loop through urban land use types)
          do iu=1,db_mx%urban
            if(urb_type == urbdb(iu)%urbnm) then
              backspace(5054)
              read(5054,*) urb_type,(salt_urban_conc(iu,isalt),isalt=1,cs_db%num_salts)
            endif
          enddo
        enddo
        
      endif
      
      endif
      
      close (5054)
  
      return
      end subroutine salt_urban_read