circulargeometrymodule Module Reference

Data Types
type	circulargeometrytype

Functions/Subroutines
real(dp) function	area_sat (this)
	Return area as if geometry is fully saturated. More...
real(dp) function	perimeter_sat (this)
	Return perimeter as if geometry is fully saturated. More...
real(dp) function	area_wet (this, depth)
	Return wetted area. More...
real(dp) function	perimeter_wet (this, depth)
	Return wetted perimeter. More...
subroutine	set_attribute (this, line)
	Set a parameter for this circular object. More...
subroutine	print_attributes (this, iout)
	Print the attributes for this object. More...

Function/Subroutine Documentation

area_sat()

real(dp) function circulargeometrymodule::area_sat ( class(circulargeometrytype)  this )

private

Definition at line 28 of file CircularGeometry.f90.

    ! -- modules
    use constantsmodule, only: dtwo, dpi
    ! -- return
    real(DP) :: area_sat
    ! -- dummy
    class(CircularGeometryType) :: this
    !
    ! -- Calculate area
    area_sat = dpi * this%radius**dtwo

area_wet()

real(dp) function circulargeometrymodule::area_wet	(	class(circulargeometrytype)	this,
		real(dp), intent(in)	depth
	)

Definition at line 56 of file CircularGeometry.f90.

    ! -- modules
    use constantsmodule, only: dtwo, dpi, dzero
    ! -- return
    real(DP) :: area_wet
    ! -- dummy
    class(CircularGeometryType) :: this
    real(DP), intent(in) :: depth
    !
    ! -- Calculate area
    if (depth <= dzero) then
      area_wet = dzero
    elseif (depth <= this%radius) then
      area_wet = this%radius * this%radius * &
                 acos((this%radius - depth) / this%radius) - &
                 (this%radius - depth) * &
                 sqrt(this%radius * this%radius - (this%radius - depth)**dtwo)
    elseif (depth <= dtwo * this%radius) then
      area_wet = this%radius * this%radius * &
                 (dpi - acos((depth - this%radius) / this%radius)) - &
                 (this%radius - depth) * &
                 sqrt(this%radius * this%radius - (this%radius - depth)**dtwo)
    else
      area_wet = dpi * this%radius * this%radius
    end if

perimeter_sat()

real(dp) function circulargeometrymodule::perimeter_sat

(

class(circulargeometrytype)

this

)

Definition at line 42 of file CircularGeometry.f90.

    ! -- modules
    use constantsmodule, only: dtwo, dpi
    ! -- return
    real(DP) :: perimeter_sat
    ! -- dummy
    class(CircularGeometryType) :: this
    !
    ! -- Calculate area
    perimeter_sat = dtwo * dpi * this%radius

perimeter_wet()

real(dp) function circulargeometrymodule::perimeter_wet	(	class(circulargeometrytype)	this,
		real(dp), intent(in)	depth
	)

Definition at line 85 of file CircularGeometry.f90.

    ! -- modules
    use constantsmodule, only: dtwo, dpi
    ! -- return
    real(DP) :: perimeter_wet
    ! -- dummy
    class(CircularGeometryType) :: this
    real(DP), intent(in) :: depth
    !
    ! -- Calculate area
    if (depth <= dzero) then
      perimeter_wet = dzero
    elseif (depth <= this%radius) then
      perimeter_wet = dtwo * this%radius * acos((this%radius - depth) / &
                                                this%radius)
    elseif (depth <= dtwo * this%radius) then
      perimeter_wet = dtwo * this%radius * (dpi - acos((depth - this%radius) / &
                                                       this%radius))
    else
      perimeter_wet = dtwo * dpi * this%radius
    end if

print_attributes()

subroutine circulargeometrymodule::print_attributes	(	class(circulargeometrytype)	this,
		integer(i4b), intent(in)	iout
	)

Definition at line 146 of file CircularGeometry.f90.

    ! -- dummy
    class(CircularGeometryType) :: this
    ! -- local
    integer(I4B), intent(in) :: iout
    ! -- formats
    character(len=*), parameter :: fmtnm = "(4x,a,a)"
    character(len=*), parameter :: fmttd = "(4x,a,1(1PG15.6))"
    !
    ! -- call parent to print parent attributes
    call this%BaseGeometryType%print_attributes(iout)
    !
    ! -- Print specifics of this geometry type
    write (iout, fmttd) 'RADIUS = ', this%radius
    write (iout, fmttd) 'SATURATED AREA = ', this%area_sat()
    write (iout, fmttd) 'SATURATED WETTED PERIMETER = ', this%perimeter_sat()

set_attribute()

subroutine circulargeometrymodule::set_attribute	(	class(circulargeometrytype)	this,
		character(len=*), intent(inout)	line
	)

Definition at line 110 of file CircularGeometry.f90.

    ! -- module
    use inputoutputmodule, only: urword
    use constantsmodule, only: linelength
    use simmodule, only: store_error, count_errors
    ! -- dummy
    class(CircularGeometryType) :: this
    character(len=LINELENGTH) :: errmsg
    character(len=*), intent(inout) :: line
    ! -- local
    integer(I4B) :: lloc, istart, istop, ival
    real(DP) :: rval
    !
    ! -- should change this and set id if uninitialized or store it
    lloc = 1
    call urword(line, lloc, istart, istop, 2, ival, rval, 0, 0)
    this%id = ival
 
    ! -- Parse the attribute
    call urword(line, lloc, istart, istop, 1, ival, rval, 0, 0)
    select case (line(istart:istop))
    case ('NAME')
      call urword(line, lloc, istart, istop, 1, ival, rval, 0, 0)
      this%name = line(istart:istop)
    case ('RADIUS')
      call urword(line, lloc, istart, istop, 3, ival, rval, 0, 0)
      this%radius = rval
    case default
      write (errmsg, '(a,a)') &
        'Unknown circular geometry attribute: ', line(istart:istop)
      call store_error(errmsg, terminate=.true.)
    end select

Here is the call graph for this function: