budgetmodule Module Reference

This module contains the BudgetModule. More...

Data Types
type	budgettype
	Derived type for the Budget object. More...

Functions/Subroutines
subroutine, public	budget_cr (this, name_model)
	@ brief Create a new budget object More...
subroutine	budget_df (this, maxsize, bdtype, bddim, labeltitle, bdzone)
	@ brief Define information for this object More...
subroutine	value_to_string (val, string, big, small)
	@ brief Convert a number to a string More...
subroutine	budget_ot (this, kstp, kper, iout)
	@ brief Output the budget table More...
subroutine	budget_da (this)
	@ brief Deallocate memory More...
subroutine	reset (this)
	@ brief Reset the budget object More...
subroutine	add_single_entry (this, rin, rout, delt, text, isupress_accumulate, rowlabel)
	@ brief Add a single row of information More...
subroutine	add_multi_entry (this, budterm, delt, budtxt, isupress_accumulate, rowlabel)
	@ brief Add multiple rows of information More...
subroutine	finalize_step (this, delt)
	@ brief Update accumulators More...
subroutine	allocate_scalars (this, name_model)
	@ brief allocate scalar variables More...
subroutine	allocate_arrays (this)
	@ brief allocate array variables More...
subroutine	resize (this, maxsize)
	@ brief Resize the budget object More...
subroutine, public	rate_accumulator (flow, rin, rout)
	@ brief Rate accumulator subroutine More...
subroutine	set_ibudcsv (this, ibudcsv)
	@ brief Set unit number for csv output file More...
subroutine	writecsv (this, totim)
	@ brief Write csv output More...
subroutine	write_csv_header (this)
	@ brief Write csv header More...

Detailed Description

New entries can be added for each time step, however, the same number of entries must be provided, and they must be provided in the same order. If not, the module will terminate with an error.

Maxsize is required as part of the df method and the arrays will be allocated to maxsize. If additional entries beyond maxsize are added, the arrays will dynamically increase in size, however, to avoid allocation and copying, it is best to set maxsize large enough up front.

vbvl(1, :) contains cumulative rate in vbvl(2, :) contains cumulative rate out vbvl(3, :) contains rate in vbvl(4, :) contains rate out vbnm(:) contains a LENBUDTXT character text string for each entry rowlabel(:) contains a LENBUDROWLABEL character text string to write as a label for each entry

Function/Subroutine Documentation

add_multi_entry()

subroutine budgetmodule::add_multi_entry ( class(budgettype)  this, real(dp), dimension(:, :), intent(in)  budterm, real(dp), intent(in)  delt, character(len=lenbudtxt), dimension(:), intent(in)  budtxt, integer(i4b), intent(in), optional  isupress_accumulate, character(len=*), intent(in), optional  rowlabel  )

private

Add information corresponding to one multiple rows in the budget table budterm is an array with inflow in column 1 and outflow in column 2 delt is the time step length budtxt is the name of the entries. It should have one entry for each row in budterm isupress_accumulate is an optional flag. If specified as 1, then the volume is NOT added to the accumulators on vbvl(1, :) and vbvl(2, :). rowlabel is a LENBUDROWLABEL character text entry that is written to the right of the table. It can be used for adding package names to budget entries. For multiple entries, the same rowlabel is used for each entry.

Parameters

	this	BudgetType object
[in]	budterm	array of budget terms
[in]	delt	time step length
[in]	budtxt	name of the entries
[in]	isupress_accumulate	suppress accumulate
[in]	rowlabel	row label

Definition at line 454 of file Budget.f90.

    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    real(DP), dimension(:, :), intent(in) :: budterm !< array of budget terms
    real(DP), intent(in) :: delt !< time step length
    character(len=LENBUDTXT), dimension(:), intent(in) :: budtxt !< name of the entries
    integer(I4B), optional, intent(in) :: isupress_accumulate !< suppress accumulate
    character(len=*), optional, intent(in) :: rowlabel !< row label
    ! -- local
    character(len=LINELENGTH) :: errmsg
    character(len=*), parameter :: fmtbuderr = &
      &"('Error in MODFLOW 6.', 'Entries do not match: ', (a), (a) )"
    integer(i4b) :: iscv, i
    integer(I4B) :: nbudterms, maxsize
    !
    iscv = 0
    if (present(isupress_accumulate)) then
      iscv = isupress_accumulate
    end if
    !
    ! -- ensure budget arrays are large enough
    nbudterms = size(budtxt)
    maxsize = this%msum - 1 + nbudterms
    if (maxsize > this%maxsize) then
      call this%resize(maxsize)
    end if
    !
    ! -- Process each of the multi-entry budget terms
    do i = 1, size(budtxt)
      !
      ! -- If budget has been written at least once, then make sure that the present
      !    text entry matches the last text entry
      if (this%written_once) then
        if (trim(adjustl(this%vbnm(this%msum))) /= &
            trim(adjustl(budtxt(i)))) then
          write (errmsg, fmtbuderr) trim(adjustl(this%vbnm(this%msum))), &
            trim(adjustl(budtxt(i)))
          call store_error(errmsg)
        end if
      end if
      !
      this%vbvl(3, this%msum) = budterm(1, i)
      this%vbvl(4, this%msum) = budterm(2, i)
      this%vbnm(this%msum) = adjustr(budtxt(i))
      if (present(rowlabel)) then
        this%rowlabel(this%msum) = adjustl(rowlabel)
        this%labeled = .true.
      end if
      this%msum = this%msum + 1
      !
    end do
    !
    ! -- Check for errors
    if (count_errors() > 0) then
      call store_error('Could not add multi-entry', terminate=.true.)
    end if
    !
    ! -- Return
    return

Here is the call graph for this function:

Here is the caller graph for this function:

add_single_entry()

subroutine budgetmodule::add_single_entry ( class(budgettype)  this, real(dp), intent(in)  rin, real(dp), intent(in)  rout, real(dp), intent(in)  delt, character(len=lenbudtxt), intent(in)  text, integer(i4b), intent(in), optional  isupress_accumulate, character(len=*), intent(in), optional  rowlabel  )

private

Add information corresponding to one row in the budget table rin the inflow rate rout is the outflow rate delt is the time step length text is the name of the entry isupress_accumulate is an optional flag. If specified as 1, then the volume is NOT added to the accumulators on vbvl(1, :) and vbvl(2, :). rowlabel is a LENBUDROWLABEL character text entry that is written to the right of the table. It can be used for adding package names to budget entries.

Parameters

	this	BudgetType object
[in]	rin	inflow rate
[in]	rout	outflow rate
[in]	delt	time step length
[in]	text	name of the entry
[in]	isupress_accumulate	accumulate flag
[in]	rowlabel	row label

Definition at line 389 of file Budget.f90.

    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    real(DP), intent(in) :: rin !< inflow rate
    real(DP), intent(in) :: rout !< outflow rate
    real(DP), intent(in) :: delt !< time step length
    character(len=LENBUDTXT), intent(in) :: text !< name of the entry
    integer(I4B), optional, intent(in) :: isupress_accumulate !< accumulate flag
    character(len=*), optional, intent(in) :: rowlabel !< row label
    ! -- local
    character(len=LINELENGTH) :: errmsg
    character(len=*), parameter :: fmtbuderr = &
      &"('Error in MODFLOW 6.', 'Entries do not match: ', (a), (a) )"
    integer(i4b) :: iscv
    integer(I4B) :: maxsize
    !
    iscv = 0
    if (present(isupress_accumulate)) then
      iscv = isupress_accumulate
    end if
    !
    ! -- ensure budget arrays are large enough
    maxsize = this%msum
    if (maxsize > this%maxsize) then
      call this%resize(maxsize)
    end if
    !
    ! -- If budget has been written at least once, then make sure that the present
    !    text entry matches the last text entry
    if (this%written_once) then
      if (trim(adjustl(this%vbnm(this%msum))) /= trim(adjustl(text))) then
        write (errmsg, fmtbuderr) trim(adjustl(this%vbnm(this%msum))), &
          trim(adjustl(text))
        call store_error(errmsg, terminate=.true.)
      end if
    end if
    !
    this%vbvl(3, this%msum) = rin
    this%vbvl(4, this%msum) = rout
    this%vbnm(this%msum) = adjustr(text)
    if (present(rowlabel)) then
      this%rowlabel(this%msum) = adjustl(rowlabel)
      this%labeled = .true.
    end if
    this%msum = this%msum + 1
    !
    ! -- Return
    return

Here is the call graph for this function:

Here is the caller graph for this function:

allocate_arrays()

subroutine budgetmodule::allocate_arrays ( class(budgettype)  this )

private

Allocate array variables of this budget object

Parameters

this	BudgetType object

Definition at line 582 of file Budget.f90.

    ! -- modules
    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    !
    ! -- If redefining, then need to deallocate/reallocate
    if (associated(this%vbvl)) then
      deallocate (this%vbvl)
      nullify (this%vbvl)
    end if
    if (associated(this%vbnm)) then
      deallocate (this%vbnm)
      nullify (this%vbnm)
    end if
    if (associated(this%rowlabel)) then
      deallocate (this%rowlabel)
      nullify (this%rowlabel)
    end if
    !
    ! -- Allocate
    allocate (this%vbvl(4, this%maxsize))
    allocate (this%vbnm(this%maxsize))
    allocate (this%rowlabel(this%maxsize))
    !
    ! -- Initialize values
    this%vbvl(:, :) = dzero
    this%vbnm(:) = ''
    this%rowlabel(:) = ''
    !
    return

allocate_scalars()

subroutine budgetmodule::allocate_scalars ( class(budgettype)  this, character(len=*), intent(in)  name_model  )

private

Allocate scalar variables of this budget object

Parameters

	this	BudgetType object
[in]	name_model	name of the model

Definition at line 544 of file Budget.f90.

    ! -- modules
    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    character(len=*), intent(in) :: name_model !< name of the model
    !
    allocate (this%msum)
    allocate (this%maxsize)
    allocate (this%budperc)
    allocate (this%written_once)
    allocate (this%labeled)
    allocate (this%bdtype)
    allocate (this%bddim)
    allocate (this%labeltitle)
    allocate (this%bdzone)
    allocate (this%ibudcsv)
    allocate (this%icsvheader)
    !
    ! -- Initialize values
    this%msum = 0
    this%maxsize = 0
    this%written_once = .false.
    this%labeled = .false.
    this%bdtype = ''
    this%bddim = ''
    this%labeltitle = ''
    this%bdzone = ''
    this%ibudcsv = 0
    this%icsvheader = 0
    !
    return

budget_cr()

subroutine, public budgetmodule::budget_cr	(	type(budgettype), pointer	this,
		character(len=*), intent(in)	name_model
	)

Create a new budget object.

Parameters

	this	BudgetType object
[in]	name_model	name of the model

Definition at line 83 of file Budget.f90.

    ! -- modules
    ! -- dummy
    type(BudgetType), pointer :: this !< BudgetType object
    character(len=*), intent(in) :: name_model !< name of the model
    !
    ! -- Create the object
    allocate (this)
    !
    ! -- Allocate scalars
    call this%allocate_scalars(name_model)
    !
    ! -- Return
    return

Here is the caller graph for this function:

budget_da()

subroutine budgetmodule::budget_da ( class(budgettype)  this )

private

Deallocate budget memory

Parameters

this	BudgetType object

Definition at line 328 of file Budget.f90.

    class(BudgetType) :: this !< BudgetType object
    !
    ! -- Scalars
    deallocate (this%msum)
    deallocate (this%maxsize)
    deallocate (this%budperc)
    deallocate (this%written_once)
    deallocate (this%labeled)
    deallocate (this%bdtype)
    deallocate (this%bddim)
    deallocate (this%labeltitle)
    deallocate (this%bdzone)
    deallocate (this%ibudcsv)
    deallocate (this%icsvheader)
    !
    ! -- Arrays
    deallocate (this%vbvl)
    deallocate (this%vbnm)
    deallocate (this%rowlabel)
    !
    ! -- Return
    return

budget_df()

subroutine budgetmodule::budget_df ( class(budgettype)  this, integer(i4b), intent(in)  maxsize, character(len=*), optional  bdtype, character(len=*), optional  bddim, character(len=*), optional  labeltitle, character(len=*), optional  bdzone  )

private

Allocate arrays and set member variables

Parameters

	this	BudgetType object
[in]	maxsize	maximum size of budget arrays
	bdtype	type of budget, default is VOLUME
	bddim	dimensions of terms, default is L**3
	labeltitle	budget label, default is PACKAGE NAME
	bdzone	corresponding zone, default is ENTIRE MODEL

Definition at line 104 of file Budget.f90.

    class(BudgetType) :: this !< BudgetType object
    integer(I4B), intent(in) :: maxsize !< maximum size of budget arrays
    character(len=*), optional :: bdtype !< type of budget, default is VOLUME
    character(len=*), optional :: bddim !< dimensions of terms, default is L**3
    character(len=*), optional :: labeltitle !< budget label, default is PACKAGE NAME
    character(len=*), optional :: bdzone !< corresponding zone, default is ENTIRE MODEL
    !
    ! -- Set values
    this%maxsize = maxsize
    !
    ! -- Allocate arrays
    call this%allocate_arrays()
    !
    ! -- Set the budget type
    if (present(bdtype)) then
      this%bdtype = bdtype
    else
      this%bdtype = 'VOLUME'
    end if
    !
    ! -- Set the budget dimension
    if (present(bddim)) then
      this%bddim = bddim
    else
      this%bddim = 'L**3'
    end if
    !
    ! -- Set the budget zone
    if (present(bdzone)) then
      this%bdzone = bdzone
    else
      this%bdzone = 'ENTIRE MODEL'
    end if
    !
    ! -- Set the label title
    if (present(labeltitle)) then
      this%labeltitle = labeltitle
    else
      this%labeltitle = 'PACKAGE NAME'
    end if
    !
    ! -- Return
    return

budget_ot()

subroutine budgetmodule::budget_ot ( class(budgettype)  this, integer(i4b), intent(in)  kstp, integer(i4b), intent(in)  kper, integer(i4b), intent(in)  iout  )

private

Write the budget table for the current set of budget information.

Parameters

	this	BudgetType object
[in]	kstp	time step
[in]	kper	stress period
[in]	iout	output unit number

Definition at line 185 of file Budget.f90.

    class(BudgetType) :: this !< BudgetType object
    integer(I4B), intent(in) :: kstp !< time step
    integer(I4B), intent(in) :: kper !< stress period
    integer(I4B), intent(in) :: iout !< output unit number
    character(len=17) :: val1, val2
    integer(I4B) :: msum1, l
    real(DP) :: two, hund, bigvl1, bigvl2, small, &
                totrin, totrot, totvin, totvot, diffr, adiffr, &
                pdiffr, pdiffv, avgrat, diffv, adiffv, avgvol
    !
    ! -- Set constants
    two = 2.d0
    hund = 100.d0
    bigvl1 = 9.99999d11
    bigvl2 = 9.99999d10
    small = 0.1d0
    !
    ! -- Determine number of individual budget entries.
    this%budperc = dzero
    msum1 = this%msum - 1
    if (msum1 <= 0) return
    !
    ! -- Clear rate and volume accumulators.
    totrin = dzero
    totrot = dzero
    totvin = dzero
    totvot = dzero
    !
    ! -- Add rates and volumes (in and out) to accumulators.
    do l = 1, msum1
      totrin = totrin + this%vbvl(3, l)
      totrot = totrot + this%vbvl(4, l)
      totvin = totvin + this%vbvl(1, l)
      totvot = totvot + this%vbvl(2, l)
    end do
    !
    ! -- Print time step number and stress period number.
    if (this%labeled) then
      write (iout, 261) trim(adjustl(this%bdtype)), trim(adjustl(this%bdzone)), &
        kstp, kper
      write (iout, 266) trim(adjustl(this%bdtype)), trim(adjustl(this%bddim)), &
        trim(adjustl(this%bddim)), this%labeltitle
    else
      write (iout, 260) trim(adjustl(this%bdtype)), trim(adjustl(this%bdzone)), &
        kstp, kper
      write (iout, 265) trim(adjustl(this%bdtype)), trim(adjustl(this%bddim)), &
        trim(adjustl(this%bddim))
    end if
    !
    ! -- Print individual inflow rates and volumes and their totals.
    do l = 1, msum1
      call value_to_string(this%vbvl(1, l), val1, bigvl1, small)
      call value_to_string(this%vbvl(3, l), val2, bigvl1, small)
      if (this%labeled) then
        write (iout, 276) this%vbnm(l), val1, this%vbnm(l), val2, this%rowlabel(l)
      else
        write (iout, 275) this%vbnm(l), val1, this%vbnm(l), val2
      end if
    end do
    call value_to_string(totvin, val1, bigvl1, small)
    call value_to_string(totrin, val2, bigvl1, small)
    write (iout, 286) val1, val2
    !
    ! -- Print individual outflow rates and volumes and their totals.
    write (iout, 287)
    do l = 1, msum1
      call value_to_string(this%vbvl(2, l), val1, bigvl1, small)
      call value_to_string(this%vbvl(4, l), val2, bigvl1, small)
      if (this%labeled) then
        write (iout, 276) this%vbnm(l), val1, this%vbnm(l), val2, this%rowlabel(l)
      else
        write (iout, 275) this%vbnm(l), val1, this%vbnm(l), val2
      end if
    end do
    call value_to_string(totvot, val1, bigvl1, small)
    call value_to_string(totrot, val2, bigvl1, small)
    write (iout, 298) val1, val2
    !
    ! -- Calculate the difference between inflow and outflow.
    !
    ! -- Calculate difference between rate in and rate out.
    diffr = totrin - totrot
    adiffr = abs(diffr)
    !
    ! -- Calculate percent difference between rate in and rate out.
    pdiffr = dzero
    avgrat = (totrin + totrot) / two
    if (avgrat /= dzero) pdiffr = hund * diffr / avgrat
    this%budperc = pdiffr
    !
    ! -- Calculate difference between volume in and volume out.
    diffv = totvin - totvot
    adiffv = abs(diffv)
    !
    ! -- Get percent difference between volume in and volume out.
    pdiffv = dzero
    avgvol = (totvin + totvot) / two
    if (avgvol /= dzero) pdiffv = hund * diffv / avgvol
    !
    ! -- Print differences and percent differences between input
    ! -- and output rates and volumes.
    call value_to_string(diffv, val1, bigvl2, small)
    call value_to_string(diffr, val2, bigvl2, small)
    write (iout, 299) val1, val2
    write (iout, 300) pdiffv, pdiffr
    !
    ! -- flush the file
    flush (iout)
    !
    ! -- set written_once to .true.
    this%written_once = .true.
    !
    ! -- formats
260 FORMAT(//2x, a, ' BUDGET FOR ', a, ' AT END OF' &
            , ' TIME STEP', i5, ', STRESS PERIOD', i4 / 2x, 78('-'))
261 FORMAT(//2x, a, ' BUDGET FOR ', a, ' AT END OF' &
            , ' TIME STEP', i5, ', STRESS PERIOD', i4 / 2x, 99('-'))
265 FORMAT(1x, /5x, 'CUMULATIVE ', a, 6x, a, 7x &
           , 'RATES FOR THIS TIME STEP', 6x, a, '/T'/5x, 18('-'), 17x, 24('-') &
           //11x, 'IN:', 38x, 'IN:'/11x, '---', 38x, '---')
266 FORMAT(1x, /5x, 'CUMULATIVE ', a, 6x, a, 7x &
           , 'RATES FOR THIS TIME STEP', 6x, a, '/T', 10x, a16, &
           /5x, 18('-'), 17x, 24('-'), 21x, 16('-') &
           //11x, 'IN:', 38x, 'IN:'/11x, '---', 38x, '---')
275 FORMAT(1x, 3x, a16, ' =', a17, 6x, a16, ' =', a17)
276 FORMAT(1x, 3x, a16, ' =', a17, 6x, a16, ' =', a17, 5x, a)
286 FORMAT(1x, /12x, 'TOTAL IN =', a, 14x, 'TOTAL IN =', a)
287 FORMAT(1x, /10x, 'OUT:', 37x, 'OUT:'/10x, 4('-'), 37x, 4('-'))
298 FORMAT(1x, /11x, 'TOTAL OUT =', a, 13x, 'TOTAL OUT =', a)
299 FORMAT(1x, /12x, 'IN - OUT =', a, 14x, 'IN - OUT =', a)
300 FORMAT(1x, /1x, 'PERCENT DISCREPANCY =', f15.2 &
           , 5x, 'PERCENT DISCREPANCY =', f15.2/)
    !
    ! -- Return
    return

Here is the call graph for this function:

finalize_step()

subroutine budgetmodule::finalize_step ( class(budgettype)  this, real(dp), intent(in)  delt  )

private

This must be called before any output is written in order to update the accumulators in vbvl(1,:) and vbl(2,:).

Parameters

this	BudgetType object

Definition at line 522 of file Budget.f90.

    ! -- modules
    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    real(DP), intent(in) :: delt
    ! -- local
    integer(I4B) :: i
    !
    do i = 1, this%msum - 1
      this%vbvl(1, i) = this%vbvl(1, i) + this%vbvl(3, i) * delt
      this%vbvl(2, i) = this%vbvl(2, i) + this%vbvl(4, i) * delt
    end do
    !
    ! -- Return
    return

rate_accumulator()

subroutine, public budgetmodule::rate_accumulator	(	real(dp), dimension(:), intent(in), contiguous	flow,
		real(dp), intent(out)	rin,
		real(dp), intent(out)	rout
	)

Routing for tallying inflows and outflows of an array

Parameters

[in]	flow	array of flows
[out]	rin	calculated sum of inflows
[out]	rout	calculated sum of outflows

Definition at line 663 of file Budget.f90.

    ! -- modules
    ! -- dummy
    real(DP), dimension(:), contiguous, intent(in) :: flow !< array of flows
    real(DP), intent(out) :: rin !< calculated sum of inflows
    real(DP), intent(out) :: rout !< calculated sum of outflows
    integer(I4B) :: n
    !
    rin = dzero
    rout = dzero
    do n = 1, size(flow)
      if (flow(n) < dzero) then
        rout = rout - flow(n)
      else
        rin = rin + flow(n)
      end if
    end do
    return

Here is the caller graph for this function:

reset()

subroutine budgetmodule::reset ( class(budgettype)  this )

private

Reset the budget object in preparation for next set of entries

Parameters

this	BudgetType object

Definition at line 358 of file Budget.f90.

    ! -- modules
    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    ! -- local
    integer(I4B) :: i
    !
    this%msum = 1
    do i = 1, this%maxsize
      this%vbvl(3, i) = dzero
      this%vbvl(4, i) = dzero
    end do
 
    ! -- Return
    return

resize()

subroutine budgetmodule::resize ( class(budgettype)  this, integer(i4b), intent(in)  maxsize  )

private

If the size wasn't allocated to be large enough, then the budget object we reallocate itself to a larger size.

Parameters

	this	BudgetType object
[in]	maxsize	maximum size

Definition at line 620 of file Budget.f90.

    ! -- modules
    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    integer(I4B), intent(in) :: maxsize !< maximum size
    ! -- local
    real(DP), dimension(:, :), allocatable :: vbvl
    character(len=LENBUDTXT), dimension(:), allocatable :: vbnm
    character(len=LENBUDROWLABEL), dimension(:), allocatable :: rowlabel
    integer(I4B) :: maxsizeold
    !
    ! -- allocate and copy into local storage
    maxsizeold = this%maxsize
    allocate (vbvl(4, maxsizeold))
    allocate (vbnm(maxsizeold))
    allocate (rowlabel(maxsizeold))
    vbvl(:, :) = this%vbvl(:, :)
    vbnm(:) = this%vbnm(:)
    rowlabel(:) = this%rowlabel(:)
    !
    ! -- Set new size and reallocate
    this%maxsize = maxsize
    call this%allocate_arrays()
    !
    ! -- Copy from local back into member variables
    this%vbvl(:, 1:maxsizeold) = vbvl(:, 1:maxsizeold)
    this%vbnm(1:maxsizeold) = vbnm(1:maxsizeold)
    this%rowlabel(1:maxsizeold) = rowlabel(1:maxsizeold)
    !
    ! - deallocate local copies
    deallocate (vbvl)
    deallocate (vbnm)
    deallocate (rowlabel)
    !
    ! -- return
    return

set_ibudcsv()

subroutine budgetmodule::set_ibudcsv ( class(budgettype)  this, integer(i4b), intent(in)  ibudcsv  )

private

This routine can be used to activate csv output by passing in a valid unit number opened for output

Parameters

	this	BudgetType object
[in]	ibudcsv	unit number for csv budget output

Definition at line 689 of file Budget.f90.

    ! -- modules
    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    integer(I4B), intent(in) :: ibudcsv !< unit number for csv budget output
    this%ibudcsv = ibudcsv
    return

value_to_string()

subroutine budgetmodule::value_to_string ( real(dp), intent(in)  val, character(len=*), intent(out)  string, real(dp), intent(in)  big, real(dp), intent(in)  small  )

private

This is sometimes needed to avoid numbers that do not fit correctly into a text string

Parameters

[in]	val	value to convert
[out]	string	string to fill
[in]	big	big value
[in]	small	small value

Definition at line 156 of file Budget.f90.

    real(DP), intent(in) :: val !< value to convert
    character(len=*), intent(out) :: string !< string to fill
    real(DP), intent(in) :: big !< big value
    real(DP), intent(in) :: small !< small value
    real(DP) :: absval
    !
    absval = abs(val)
    if (val /= dzero .and. (absval >= big .or. absval < small)) then
      if (absval >= 1.d100 .or. absval <= 1.d-100) then
        ! -- if exponent has 3 digits, then need to explicitly use the ES
        !    format to force writing the E character
        write (string, '(es17.4E3)') val
      else
        write (string, '(1pe17.4)') val
      end if
    else
      ! -- value is within range where number looks good with F format
      write (string, '(f17.4)') val
    end if
    return

Here is the caller graph for this function:

write_csv_header()

subroutine budgetmodule::write_csv_header ( class(budgettype)  this )

private

This routine will write the csv header based on the names in vbnm

Parameters

this	BudgetType object

Definition at line 763 of file Budget.f90.

    ! -- modules
    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    ! -- local
    integer(I4B) :: l
    character(len=LINELENGTH) :: txt, txtl
    write (this%ibudcsv, '(a)', advance='NO') 'time,'
    !
    ! -- first write IN
    do l = 1, this%msum - 1
      txt = this%vbnm(l)
      txtl = ''
      if (this%labeled) then
        txtl = '('//trim(adjustl(this%rowlabel(l)))//')'
      end if
      txt = trim(adjustl(txt))//trim(adjustl(txtl))//'_IN,'
      write (this%ibudcsv, '(a)', advance='NO') trim(adjustl(txt))
    end do
    !
    ! -- then write OUT
    do l = 1, this%msum - 1
      txt = this%vbnm(l)
      txtl = ''
      if (this%labeled) then
        txtl = '('//trim(adjustl(this%rowlabel(l)))//')'
      end if
      txt = trim(adjustl(txt))//trim(adjustl(txtl))//'_OUT,'
      write (this%ibudcsv, '(a)', advance='NO') trim(adjustl(txt))
    end do
    write (this%ibudcsv, '(a)') 'TOTAL_IN,TOTAL_OUT,PERCENT_DIFFERENCE'
    !
    ! -- return
    return

writecsv()

subroutine budgetmodule::writecsv ( class(budgettype)  this, real(dp), intent(in)  totim  )

private

This routine will write a row of output to the csv file, if it is available for output. Upon first call, it will write the csv header.

Parameters

	this	BudgetType object
[in]	totim	time corresponding to this data

Definition at line 705 of file Budget.f90.

    ! -- modules
    ! -- dummy
    class(BudgetType) :: this !< BudgetType object
    real(DP), intent(in) :: totim !< time corresponding to this data
    ! -- local
    integer(I4B) :: i
    real(DP) :: totrin
    real(DP) :: totrout
    real(DP) :: diffr
    real(DP) :: pdiffr
    real(DP) :: avgrat
    !
    if (this%ibudcsv > 0) then
      !
      ! -- write header
      if (this%icsvheader == 0) then
        call this%write_csv_header()
        this%icsvheader = 1
      end if
      !
      ! -- Calculate in and out
      totrin = dzero
      totrout = dzero
      do i = 1, this%msum - 1
        totrin = totrin + this%vbvl(3, i)
        totrout = totrout + this%vbvl(4, i)
      end do
      !
      ! -- calculate percent difference
      diffr = totrin - totrout
      pdiffr = dzero
      avgrat = (totrin + totrout) / dtwo
      if (avgrat /= dzero) then
        pdiffr = dhundred * diffr / avgrat
      end if
      !
      ! -- write data
      write (this%ibudcsv, '(*(G0,:,","))') &
        totim, &
        (this%vbvl(3, i), i=1, this%msum - 1), &
        (this%vbvl(4, i), i=1, this%msum - 1), &
        totrin, totrout, pdiffr
      !
      ! -- flush the file
      flush (this%ibudcsv)
    end if
    !
    ! -- return
    return