allowing to specify multiple tracers in control file and reading tracer names

Author: nmizukami

route/build/src/globalData.f90

@@ -252,10 +252,12 @@ MODULE globalData
   real(dp),           allocatable, public :: basinEvapo_main(:)     ! HRU evaporation array (m/s) for mainstem
   real(dp),           allocatable, public :: basinPrecip_trib(:)    ! HRU precipitation array (m/s) for tributaries
   real(dp),           allocatable, public :: basinPrecip_main(:)    ! HRU precipitation array (m/s) for mainstem
- ! solute (tracer) data
+
+  ! solute (tracer) data
   integer(i4b)                            :: nTracer                ! number of tracers
   real(dp),           allocatable, public :: basinSolute_trib(:,:)  ! HRU constituent array [nhru, nTracer] (mg/s) for tributaries
   real(dp),           allocatable, public :: basinSolute_main(:,:)  ! HRU constituent array [nhru, nTracer] (mg/s) for mainstem
+
   ! seg water management fluxes and target volume
   type(wm),                        public :: wm_data                ! SEG flux and target vol data structure for one time step for river network
   real(dp),           allocatable, public :: flux_wm_trib(:)        ! SEG flux array (m3/s) for tributaries

route/build/src/read_control.f90

@@ -37,6 +37,7 @@ SUBROUTINE read_control(ctl_fname, err, message)
  USE globalData, ONLY: idxSUM,idxIRF,idxKWT, &
                        idxKW,idxMC,idxDW
  USE globalData, ONLY: runMode                 ! mizuRoute run mode: standalone, cesm-coupling
+ USE globalData, ONLY: nTracer                 ! number of active tracers
  ! index of named variables in each structure
  USE var_lookup, ONLY: ixHRU
  USE var_lookup, ONLY: ixHRU2SEG
@@ -48,6 +49,7 @@ SUBROUTINE read_control(ctl_fname, err, message)
  ! external subroutines
  USE ascii_utils, ONLY: file_open        ! open file (performs a few checks as well)
  USE ascii_utils, ONLY: get_vlines       ! get a list of character strings from non-comment lines
+ USE ascii_utils, ONLY: split_line       ! split string into sub-strings with delimiter
  USE ascii_utils, ONLY: lower            ! convert string to lower case
  USE nr_utils,    ONLY: char2int         ! convert integer number to a array containing individual digits
 
@@ -61,7 +63,7 @@ SUBROUTINE read_control(ctl_fname, err, message)
  character(len=strLen)             :: cName,cData             ! name and data from cLines(iLine)
  character(len=strLen)             :: cLength,cTime           ! length and time units
  character(len=strLen)             :: cMass                   ! mass units needed only when tracer is on
- character(len=strLen)             :: vname                   ! multiple variable names with semi-colon delimited
+ character(len=strLen),allocatable :: vname_tmp(:)            ! multiple variable names with semi-colon delimited
  logical(lgt)                      :: isGeneric               ! temporal logical scalar
  logical(lgt)                      :: onlyOneRouting          ! temporal logical scalar
  integer(i4b)                      :: ipos                    ! index of character string
@@ -159,7 +161,7 @@ SUBROUTINE read_control(ctl_fname, err, message)
    case('<vname_qsim>');           vname_qsim   = trim(cData)                          ! name of runoff variable
    case('<vname_evapo>');          vname_evapo  = trim(cData)                          ! name of actual evapoartion variable
    case('<vname_precip>');         vname_precip = trim(cData)                          ! name of precipitation variable
-   case('<vname_solute>');         vname        = trim(cData)                          ! name of solute mass flux variables (up to 10 is ok)
+   case('<vname_solute>');         vname_tmp    = split_line(trim(cData),delim=';,:')  ! name of solute mass flux variables (up to 10 is ok), acceptable delimiter: : , :
    case('<vname_time>');           vname_time   = trim(cData)                          ! name of time variable in the runoff file
    case('<vname_hruid>');          vname_hruid  = trim(cData)                          ! name of the HRU id
    case('<dname_time>');           dname_time   = trim(cData)                          ! name of time variable in the runoff file
@@ -383,6 +385,9 @@ SUBROUTINE read_control(ctl_fname, err, message)
  end do  ! looping through lines in the control file
 
  ! ---------- Perform minor processing and checking control variables ----------------------------------------
+ ! extract tracer names into vname_solute
+ nTracer=size(vname_tmp)
+ vname_solute(1:nTracer)=vname_tmp(:)
 
  ! ---------- directory option  ---------------------------------------------------------------------
  if (trim(restart_dir)==charMissing) then