Skip to content
Snippets Groups Projects
Commit 2aa1aa56 authored by Kyle's avatar Kyle
Browse files

copied files from summa

parent eb9a5c90
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
build/source/engine/sundials/type4IDA.f90
+ 61
− 68
View file @ 2aa1aa56
module type4IDA module type4IDA
! data types ! data types
...@@ -16,69 +14,64 @@ USE data_types,only:& ...@@ -16,69 +14,64 @@ USE data_types,only:&
implicit none implicit none
type eqnsData type eqnsData
type(c_ptr) :: ida_mem ! IDA memory type(c_ptr) :: ida_mem ! IDA memory
real(rkind) :: dt ! time step real(rkind) :: dt ! time step
integer(i4b) :: nSnow ! number of snow layers integer(i4b) :: nSnow ! number of snow layers
integer(i4b) :: nSoil ! number of soil layers integer(i4b) :: nSoil ! number of soil layers
integer(i4b) :: nLayers ! total number of layers integer(i4b) :: nLayers ! total number of layers
integer :: nState ! total number of state variables integer :: nState ! total number of state variables
integer(i4b) :: ixMatrix ! form of matrix (dense or banded) integer(i4b) :: ixMatrix ! form of matrix (dense or banded)
logical(lgt) :: firstSubStep ! flag to indicate if we are processing the first sub-step logical(lgt) :: firstSubStep ! flag to indicate if we are processing the first sub-step
logical(lgt) :: firstFluxCall logical(lgt) :: firstFluxCall
logical(lgt) :: firstSplitOper logical(lgt) :: firstSplitOper
logical(lgt) :: computeVegFlux ! flag to indicate if computing fluxes over vegetation logical(lgt) :: computeVegFlux ! flag to indicate if computing fluxes over vegetation
logical(lgt) :: scalarSolution ! flag to denote if implementing the scalar solution logical(lgt) :: scalarSolution ! flag to denote if implementing the scalar solution
type(zLookup) :: lookup_data ! lookup tables type(zLookup) :: lookup_data ! lookup tables
type(var_i) :: type_data ! type of vegetation and soil type(var_i) :: type_data ! type of vegetation and soil
type(var_d) :: attr_data ! spatial attributes type(var_d) :: attr_data ! spatial attributes
type(var_dlength) :: mpar_data ! model parameters type(var_dlength) :: mpar_data ! model parameters
type(var_d) :: forc_data ! model forcing data type(var_d) :: forc_data ! model forcing data
type(var_dlength) :: bvar_data ! model variables for the local basin type(var_dlength) :: bvar_data ! model variables for the local basin
type(var_dlength) :: prog_data ! prognostic variables for a local HRU type(var_dlength) :: prog_data ! prognostic variables for a local HRU
type(var_ilength) :: indx_data ! indices defining model states and layers type(var_ilength) :: indx_data ! indices defining model states and layers
type(var_dlength) :: diag_data ! diagnostic variables for a local HRU type(var_dlength) :: diag_data ! diagnostic variables for a local HRU
type(var_dlength) :: flux_data ! model fluxes for a local HRU type(var_dlength) :: flux_data ! model fluxes for a local HRU
type(var_dlength) :: deriv_data ! derivatives in model fluxes w.r.t. relevant state variables type(var_dlength) :: deriv_data ! derivatives in model fluxes w.r.t. relevant state variables
real(rkind) :: scalarCanopyTempTrial ! trial value of canopy temperature (K) real(rkind) :: scalarCanopyTempTrial ! trial value of canopy temperature (K)
real(rkind) :: scalarCanopyTempPrev ! previous value of canopy temperature (K) real(rkind) :: scalarCanopyTempPrev ! previous value of canopy temperature (K)
real(rkind) :: scalarCanopyIceTrial ! trial value of canopy ice content (kg m-2) real(rkind) :: scalarCanopyIceTrial ! trial value of canopy ice content (kg m-2)
real(rkind) :: scalarCanopyIcePrev ! value of canopy ice content (kg m-2) at previous step real(rkind) :: scalarCanopyIcePrev ! value of canopy ice content (kg m-2) at previous step
real(rkind) :: scalarCanopyLiqTrial ! trial value of canopy ice content (kg m-2) real(rkind) :: scalarCanopyLiqTrial ! trial value of canopy ice content (kg m-2)
real(rkind) :: scalarCanopyLiqPrev ! value of canopy ice content (kg m-2) at previous step real(rkind) :: scalarCanopyLiqPrev ! value of canopy ice content (kg m-2) at previous step
real(rkind) :: scalarCanopyEnthalpyTrial ! trial enthalpy of the vegetation canopy (J m-3) real(rkind) :: scalarCanopyEnthalpyTrial ! trial enthalpy of the vegetation canopy (J m-3)
real(rkind) :: scalarCanopyEnthalpyPrev ! previous enthalpy of the vegetation canopy (J m-3) real(rkind) :: scalarCanopyEnthalpyPrev ! previous enthalpy of the vegetation canopy (J m-3)
real(qp), allocatable :: sMul(:) ! state vector multiplier (used in the residual calculations) real(qp), allocatable :: sMul(:) ! state vector multiplier (used in the residual calculations)
real(rkind), allocatable :: dMat(:) ! diagonal of the Jacobian matrix real(rkind), allocatable :: dMat(:) ! diagonal of the Jacobian matrix
real(rkind), allocatable :: fluxVec(:) ! flux vector real(rkind), allocatable :: fluxVec(:) ! flux vector
real(qp), allocatable :: resSink(:) ! additional (sink) terms on the RHS of the state equation real(qp), allocatable :: resSink(:) ! additional (sink) terms on the RHS of the state equation
real(rkind), allocatable :: atol(:) ! vector of absolute tolerances real(rkind), allocatable :: atol(:) ! vector of absolute tolerances
real(rkind), allocatable :: rtol(:) ! vector of relative tolerances real(rkind), allocatable :: rtol(:) ! vector of relative tolerances
real(rkind), allocatable :: mLayerMatricHeadLiqTrial(:) ! trial value for liquid water matric potential (m) real(rkind), allocatable :: mLayerMatricHeadLiqTrial(:) ! trial value for liquid water matric potential (m)
real(rkind), allocatable :: mLayerMatricHeadTrial(:) ! trial vector of total water matric potential (m) real(rkind), allocatable :: mLayerMatricHeadTrial(:) ! trial vector of total water matric potential (m)
real(rkind), allocatable :: mLayerMatricHeadPrev(:) ! vector of total water matric potential (m) at previous step real(rkind), allocatable :: mLayerMatricHeadPrev(:) ! vector of total water matric potential (m) at previous step
real(rkind), allocatable :: mLayerVolFracWatTrial(:) ! trial value for volumetric fraction of total water (-) real(rkind), allocatable :: mLayerVolFracWatTrial(:) ! trial value for volumetric fraction of total water (-)
real(rkind), allocatable :: mLayerVolFracWatPrev(:) ! value for volumetric fraction of total water (-) at previous step real(rkind), allocatable :: mLayerVolFracWatPrev(:) ! value for volumetric fraction of total water (-) at previous step
real(rkind), allocatable :: mLayerVolFracIceTrial(:) ! trial value for volumetric fraction of ice (-) real(rkind), allocatable :: mLayerVolFracIceTrial(:) ! trial value for volumetric fraction of ice (-)
real(rkind), allocatable :: mLayerVolFracIcePrev(:) ! value for volumetric fraction of ice (-) at previous step real(rkind), allocatable :: mLayerVolFracIcePrev(:) ! value for volumetric fraction of ice (-) at previous step
real(rkind), allocatable :: mLayerVolFracLiqTrial(:) ! trial value for volumetric fraction of liquid water (-) real(rkind), allocatable :: mLayerVolFracLiqTrial(:) ! trial value for volumetric fraction of liquid water (-)
real(rkind), allocatable :: mLayerVolFracLiqPrev(:) ! value for volumetric fraction of liquid water (-) at previous step real(rkind), allocatable :: mLayerVolFracLiqPrev(:) ! value for volumetric fraction of liquid water (-) at previous step
real(rkind) :: scalarAquiferStoragePrev real(rkind) :: scalarAquiferStoragePrev
real(rkind) :: scalarAquiferStorageTrial real(rkind) :: scalarAquiferStorageTrial
real(rkind), allocatable :: mLayerEnthalpyTrial(:) ! trial enthalpy of snow and soil (J m-3) real(rkind), allocatable :: mLayerEnthalpyTrial(:) ! trial enthalpy of snow and soil (J m-3)
real(rkind), allocatable :: mLayerEnthalpyPrev(:) ! enthalpy of snow and soil (J m-3) at previous step real(rkind), allocatable :: mLayerEnthalpyPrev(:) ! enthalpy of snow and soil (J m-3) at previous step
real(rkind), allocatable :: mLayerTempTrial(:) ! trial vector of layer temperature (K) real(rkind), allocatable :: mLayerTempTrial(:) ! trial vector of layer temperature (K)
real(rkind), allocatable :: mLayerTempPrev(:) ! vector of layer temperature (K) at previous step real(rkind), allocatable :: mLayerTempPrev(:) ! vector of layer temperature (K) at previous step
real(rkind), allocatable :: dBaseflow_dMatric(:,:) ! derivative in baseflow w.r.t. matric head (s-1) real(rkind), allocatable :: dBaseflow_dMatric(:,:) ! derivative in baseflow w.r.t. matric head (s-1)
integer :: ixSaturation integer :: ixSaturation
integer(i4b) :: err ! error code integer(i4b) :: err ! error code
character(len = 50) :: message ! error message character(len = 50) :: message ! error message
end type eqnsData end type eqnsData
end module type4IDA end module type4IDA
\ No newline at end of file
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment