heat_capacity_base.F90 Source File


This file depends on

sourcefile~~heat_capacity_base.f90~~EfferentGraph sourcefile~heat_capacity_base.f90 heat_capacity_base.F90 sourcefile~heat_capacity_interface.f90 heat_capacity_interface.F90 sourcefile~heat_capacity_base.f90->sourcefile~heat_capacity_interface.f90 sourcefile~core.f90 core.F90 sourcefile~heat_capacity_interface.f90->sourcefile~core.f90 sourcefile~density_interface.f90 density_interface.F90 sourcefile~heat_capacity_interface.f90->sourcefile~density_interface.f90 sourcefile~input_interface.f90 input_interface.F90 sourcefile~heat_capacity_interface.f90->sourcefile~input_interface.f90 sourcefile~allocate.f90 allocate.F90 sourcefile~core.f90->sourcefile~allocate.f90 sourcefile~check_range.f90 check_range.F90 sourcefile~core.f90->sourcefile~check_range.f90 sourcefile~deallocate.f90 deallocate.F90 sourcefile~core.f90->sourcefile~deallocate.f90 sourcefile~error.f90 error.F90 sourcefile~core.f90->sourcefile~error.f90 sourcefile~fortran_utils.f90 fortran_utils.F90 sourcefile~core.f90->sourcefile~fortran_utils.f90 sourcefile~string_utils.f90 string_utils.F90 sourcefile~core.f90->sourcefile~string_utils.f90 sourcefile~types.f90 types.F90 sourcefile~core.f90->sourcefile~types.f90 sourcefile~unique.f90 unique.F90 sourcefile~core.f90->sourcefile~unique.f90 sourcefile~vtk.f90 vtk.F90 sourcefile~core.f90->sourcefile~vtk.f90 sourcefile~vtk_constants.f90 vtk_constants.F90 sourcefile~core.f90->sourcefile~vtk_constants.f90 sourcefile~density_interface.f90->sourcefile~core.f90 sourcefile~input.f90 input.F90 sourcefile~density_interface.f90->sourcefile~input.f90 sourcefile~input_interface.f90->sourcefile~core.f90 sourcefile~project_settings.f90 project_settings.F90 sourcefile~input_interface.f90->sourcefile~project_settings.f90 sourcefile~allocate.f90->sourcefile~error.f90 sourcefile~deallocate.f90->sourcefile~error.f90 sourcefile~memory_stats_wrapper.f90 memory_stats_wrapper.F90 sourcefile~fortran_utils.f90->sourcefile~memory_stats_wrapper.f90 sourcefile~signal_flag_wrapper.f90 signal_flag_wrapper.F90 sourcefile~fortran_utils.f90->sourcefile~signal_flag_wrapper.f90 sourcefile~system_info_wrapper.f90 system_info_wrapper.F90 sourcefile~fortran_utils.f90->sourcefile~system_info_wrapper.f90 sourcefile~input.f90->sourcefile~input_interface.f90 sourcefile~project_settings.f90->sourcefile~core.f90 sourcefile~string_utils.f90->sourcefile~allocate.f90 sourcefile~array.f90 array.F90 sourcefile~types.f90->sourcefile~array.f90 sourcefile~gauss.f90 gauss.F90 sourcefile~types.f90->sourcefile~gauss.f90 sourcefile~pointer.f90 pointer.F90 sourcefile~types.f90->sourcefile~pointer.f90 sourcefile~variable.f90 variable.F90 sourcefile~types.f90->sourcefile~variable.f90 sourcefile~vector.f90 vector.F90 sourcefile~types.f90->sourcefile~vector.f90 sourcefile~unique.f90->sourcefile~allocate.f90 sourcefile~vtk.f90->sourcefile~allocate.f90 sourcefile~vtk.f90->sourcefile~deallocate.f90 sourcefile~vtk.f90->sourcefile~unique.f90 sourcefile~vtk.f90->sourcefile~vtk_constants.f90 sourcefile~vtk.f90->sourcefile~array.f90 sourcefile~vtk_wrapper.f90 vtk_wrapper.F90 sourcefile~vtk.f90->sourcefile~vtk_wrapper.f90 sourcefile~vtu_wrapper.f90 vtu_wrapper.F90 sourcefile~vtk.f90->sourcefile~vtu_wrapper.f90 sourcefile~array.f90->sourcefile~allocate.f90 sourcefile~array.f90->sourcefile~deallocate.f90 sourcefile~c_utils.f90 c_utils.F90 sourcefile~memory_stats_wrapper.f90->sourcefile~c_utils.f90 sourcefile~signal_flag.f90 signal_flag.F90 sourcefile~signal_flag_wrapper.f90->sourcefile~signal_flag.f90 sourcefile~system_info_wrapper.f90->sourcefile~c_utils.f90 sourcefile~variable.f90->sourcefile~allocate.f90 sourcefile~c_utils.f90->sourcefile~signal_flag.f90 sourcefile~memory_stats.f90 memory_stats.F90 sourcefile~c_utils.f90->sourcefile~memory_stats.f90 sourcefile~system_info.f90 system_info.F90 sourcefile~c_utils.f90->sourcefile~system_info.f90

Source Code

submodule(calculate_volumetric_heat_capacity) Calc_VHC_Base
    implicit none

contains

    module subroutine initialize_holder_vhcs(self, input, i_material)
        implicit none
        class(holder_vhcs), intent(inout) :: self
        type(type_input), intent(in) :: input
        integer(int32), intent(in) :: i_material

        select case (input%basic%materials(i_material)%phase)
        case (3)
            if (input%basic%materials(i_material)%is_frozen) then
                self%p = type_vhc_3phase_apparent(input, i_material)
            else
                self%p = type_vhc_3phase(input, i_material)
            end if
        end select

    end subroutine initialize_holder_vhcs

    module function calc_vhc_3(VHC_soil, phi_soil, &
                               VHC_water, phi_water, &
                               VHC_ice, phi_ice) result(VHC)
        implicit none
        real(real64), intent(in) :: VHC_soil
        real(real64), intent(in) :: phi_soil
        real(real64), intent(in) :: VHC_water
        real(real64), intent(in) :: phi_water
        real(real64), intent(in) :: VHC_ice
        real(real64), intent(in) :: phi_ice
        real(real64) :: VHC

        VHC = VHC_soil * phi_soil + VHC_water * phi_water + VHC_ice * phi_ice

    end function calc_vhc_3

    module function calc_vhc_3a(VHC_soil, phi_soil, VHC_water, phi_water, &
                                VHC_ice, phi_ice, Lf, DEN_ice, dQi_dT) result(VHC)
        implicit none
        real(real64), intent(in) :: VHC_soil
        real(real64), intent(in) :: phi_soil
        real(real64), intent(in) :: VHC_water
        real(real64), intent(in) :: phi_water
        real(real64), intent(in) :: VHC_ice
        real(real64), intent(in) :: phi_ice
        real(real64), intent(in) :: Lf
        real(real64), intent(in) :: DEN_ice
        real(real64), intent(in) :: dQi_dT
        real(real64) :: VHC

        VHC = VHC_soil * phi_soil + VHC_water * phi_water + VHC_ice * phi_ice - Lf * DEN_ice * dQi_dT

    end function calc_vhc_3a

end submodule Calc_VHC_Base