Skip to content
Snippets Groups Projects
Commit 5767ca55 authored by ashleymedin's avatar ashleymedin
Browse files

removing more files shouldn't need, moved them into summa code with compiler... 

removing more files shouldn't need, moved them into summa code with compiler def. need to update cmake files to reflect this
parent e944a604
No related branches found
No related tags found
2 merge requests!5add in compiler def ACTORS_ACTIVE to turn on some stuff so can move...,!4Removing repeated code between this and the inner Summa the user clones
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 0 additions and 6429 deletions
build/source/dshare/globalData.f90 deleted 100755 → 0
+ 0
− 337
View file @ e944a604
! SUMMA - Structure for Unifying Multiple Modeling Alternatives
! Copyright (C) 2014-2020 NCAR/RAL; University of Saskatchewan; University of Washington
!
! This file is part of SUMMA
!
! For more information see: http://www.ral.ucar.edu/projects/summa
!
! This program is free software: you can redistribute it and/or modify
! it under the terms of the GNU General Public License as published by
! the Free Software Foundation, either version 3 of the License, or
! (at your option) any later version.
!
! This program is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with this program. If not, see <http://www.gnu.org/licenses/>.
! ----------------------------------------------------------------------------------------------------------------
! * part 1: parameters that are fixed across multiple instantiations
! ----------------------------------------------------------------------------------------------------------------
MODULE globalData
! data types
USE nrtype
USE netcdf
USE,intrinsic :: ieee_arithmetic ! IEEE arithmetic
USE data_types,only:gru2hru_map ! mapping between the GRUs and HRUs
USE data_types,only:hru2gru_map ! mapping between the GRUs and HRUs
USE data_types,only:model_options ! the model decision structure
USE data_types,only:file_info ! metadata for model forcing datafile
USE data_types,only:par_info ! default parameter values and parameter bounds
USE data_types,only:var_info ! metadata for variables in each model structure
USE data_types,only:flux2state ! extended metadata to define flux-to-state mapping
USE data_types,only:extended_info ! extended metadata for variables in each model structure
USE data_types,only:struct_info ! summary information on all data structures
USE data_types,only:var_i ! vector of integers
#ifdef ACTORS_ACTIVE
USE data_types,only:var_forc ! for Actors
USE data_types,only:dlength ! for Actors
USE data_types,only:ilength ! for Actors
USE data_types,only:init_cond ! for Actors
#endif
! number of variables in each data structure
USE var_lookup,only:maxvarTime ! time: maximum number variables
USE var_lookup,only:maxvarForc ! forcing data: maximum number variables
USE var_lookup,only:maxvarAttr ! attributes: maximum number variables
USE var_lookup,only:maxvarType ! type index: maximum number variables
USE var_lookup,only:maxvarId ! IDs index: maximum number variables
USE var_lookup,only:maxvarProg ! prognostic variables: maximum number variables
USE var_lookup,only:maxvarDiag ! diagnostic variables: maximum number variables
USE var_lookup,only:maxvarFlux ! model fluxes: maximum number variables
USE var_lookup,only:maxvarDeriv ! model derivatives: maximum number variables
USE var_lookup,only:maxvarIndx ! model indices: maximum number variables
USE var_lookup,only:maxvarMpar ! model parameters: maximum number variables
USE var_lookup,only:maxvarBvar ! basin-average variables: maximum number variables
USE var_lookup,only:maxvarBpar ! basin-average parameters: maximum number variables
USE var_lookup,only:maxvarDecisions ! maximum number of decisions
USE var_lookup,only:maxvarFreq ! maximum number of output files
USE var_lookup,only:maxvarLookup ! maximum number of variables in the lookup
implicit none
private
! ----------------------------------------------------------------------------------------------------------------
! * part 1: parameters that are fixed across multiple instantiations
! ----------------------------------------------------------------------------------------------------------------
! define missing values
real(rkind),parameter,public :: quadMissing = nr_quadMissing ! (from nrtype) missing quadruple precision number
real(rkind),parameter,public :: realMissing = nr_realMissing ! (from nrtype) missing double precision number
integer(i4b),parameter,public :: integerMissing = nr_integerMissing ! (from nrtype) missing integer
! define run modes
integer(i4b),parameter,public :: iRunModeFull=1 ! named variable defining running mode as full run (all GRUs)
integer(i4b),parameter,public :: iRunModeGRU=2 ! named variable defining running mode as GRU-parallelization run (GRU subset)
integer(i4b),parameter,public :: iRunModeHRU=3 ! named variable defining running mode as single-HRU run (ONE HRU)
! define progress modes
integer(i4b),parameter,public :: ixProgress_im=1000 ! named variable to print progress once per month
integer(i4b),parameter,public :: ixProgress_id=1001 ! named variable to print progress once per day
integer(i4b),parameter,public :: ixProgress_ih=1002 ! named variable to print progress once per hour
integer(i4b),parameter,public :: ixProgress_never=1003 ! named variable to print progress never
integer(i4b),parameter,public :: ixProgress_it=1004 ! named variable to print progress every timestep
! define restart frequency
integer(i4b),parameter,public :: ixRestart_iy=1000 ! named variable to print a re-start file once per year
integer(i4b),parameter,public :: ixRestart_im=1001 ! named variable to print a re-start file once per month
integer(i4b),parameter,public :: ixRestart_id=1002 ! named variable to print a re-start file once per day
integer(i4b),parameter,public :: ixRestart_end=1003 ! named variable to print a re-start file at the end of a run
integer(i4b),parameter,public :: ixRestart_never=1004 ! named variable to print a re-start file never
! define output file frequency
integer(i4b),parameter,public :: noNewFiles=1001 ! no new output files
integer(i4b),parameter,public :: newFileEveryOct1=1002 ! create a new file on Oct 1 every year (start of the USA water year)
! define named variables for "yes" and "no"
integer(i4b),parameter,public :: no=0 ! .false.
integer(i4b),parameter,public :: yes=1 ! .true.
! define named variables to describe the domain type
integer(i4b),parameter,public :: iname_cas =1000 ! named variable to denote a canopy air space state variable
integer(i4b),parameter,public :: iname_veg =1001 ! named variable to denote a vegetation state variable
integer(i4b),parameter,public :: iname_soil=1002 ! named variable to denote a soil layer
integer(i4b),parameter,public :: iname_snow=1003 ! named variable to denote a snow layer
integer(i4b),parameter,public :: iname_aquifer=1004 ! named variable to denote a snow layer
! define named variables to describe the state variable type
integer(i4b),parameter,public :: iname_nrgCanair=2001 ! named variable defining the energy of the canopy air space
integer(i4b),parameter,public :: iname_nrgCanopy=2002 ! named variable defining the energy of the vegetation canopy
integer(i4b),parameter,public :: iname_watCanopy=2003 ! named variable defining the mass of total water on the vegetation canopy
integer(i4b),parameter,public :: iname_liqCanopy=2004 ! named variable defining the mass of liquid water on the vegetation canopy
integer(i4b),parameter,public :: iname_nrgLayer=3001 ! named variable defining the energy state variable for snow+soil layers
integer(i4b),parameter,public :: iname_watLayer=3002 ! named variable defining the total water state variable for snow+soil layers
integer(i4b),parameter,public :: iname_liqLayer=3003 ! named variable defining the liquid water state variable for snow+soil layers
integer(i4b),parameter,public :: iname_matLayer=3004 ! named variable defining the matric head state variable for soil layers
integer(i4b),parameter,public :: iname_lmpLayer=3005 ! named variable defining the liquid matric potential state variable for soil layers
integer(i4b),parameter,public :: iname_watAquifer=3006 ! named variable defining the water storage in the aquifer
! define named variables to describe the form and structure of the band-diagonal matrices used in the numerical solver
! NOTE: This indexing scheme provides the matrix structure expected by lapack and sundials. Specifically, they require kl extra rows for additional storage.
! Consequently, all indices are offset by kl and the total number of bands for storage is 2*kl+ku+1 instead of kl+ku+1.
integer(i4b),parameter,public :: nRHS=1 ! number of unknown variables on the RHS of the linear system A.X=B
integer(i4b),parameter,public :: ku=3 ! number of super-diagonal bands, ku>=3 to accommodate coupled layer above
integer(i4b),parameter,public :: kl=4 ! number of sub-diagonal bands, kl>=4 to accommodate vegetation
integer(i4b),parameter,public :: ixDiag=kl+ku+1 ! index for the diagonal band
integer(i4b),parameter,public :: nBands=2*kl+ku+1 ! length of the leading dimension of the band diagonal matrix
! define named variables for the type of matrix used in the numerical solution.
integer(i4b),parameter,public :: ixFullMatrix=1001 ! named variable for the full Jacobian matrix
integer(i4b),parameter,public :: ixBandMatrix=1002 ! named variable for the band diagonal matrix
! define indices describing the first and last layers of the Jacobian to print (for debugging)
integer(i4b),parameter,public :: iJac1=16 ! first layer of the Jacobian to print
integer(i4b),parameter,public :: iJac2=20 ! last layer of the Jacobian to print
! define limit checks
real(rkind),parameter,public :: verySmall=tiny(1.0_rkind) ! a very small number
real(rkind),parameter,public :: veryBig=1.e+20_rkind ! a very big number
! define algorithmic control parameters
real(rkind),parameter,public :: dx = 1.e-8_rkind ! finite difference increment
! define summary information on all data structures
integer(i4b),parameter :: nStruct=14 ! number of data structures
type(struct_info),parameter,public,dimension(nStruct) :: structInfo=(/&
struct_info('time', 'TIME' , maxvarTime ), & ! the time data structure
struct_info('forc', 'FORCE', maxvarForc ), & ! the forcing data structure
struct_info('attr', 'ATTR' , maxvarAttr ), & ! the attribute data structure
struct_info('type', 'TYPE' , maxvarType ), & ! the type data structure
struct_info('id' , 'ID' , maxvarId ), & ! the type data structure
struct_info('mpar', 'PARAM', maxvarMpar ), & ! the model parameter data structure
struct_info('bpar', 'BPAR' , maxvarBpar ), & ! the basin parameter data structure
struct_info('bvar', 'BVAR' , maxvarBvar ), & ! the basin variable data structure
struct_info('indx', 'INDEX', maxvarIndx ), & ! the model index data structure
struct_info('prog', 'PROG', maxvarProg ), & ! the prognostic (state) variable data structure
struct_info('diag', 'DIAG' , maxvarDiag ), & ! the diagnostic variable data structure
struct_info('flux', 'FLUX' , maxvarFlux ), & ! the flux data structure
struct_info('deriv', 'DERIV', maxvarDeriv), & ! the model derivative data structure
struct_info('lookup','LOOKUP',maxvarLookup) /) ! the lookup table data structure
! fixed model decisions
logical(lgt) , parameter, public :: overwriteRSMIN=.false. ! flag to overwrite RSMIN
integer(i4b) , parameter, public :: maxSoilLayers=10000 ! Maximum Number of Soil Layers
! ----------------------------------------------------------------------------------------------------------------
! * part 2: globally constant variables/structures that require initialization
! ----------------------------------------------------------------------------------------------------------------
! define Not-a-Number (NaN)
real(rkind),save,public :: dNaN
! define default parameter values and parameter bounds
type(par_info),save,public :: localParFallback(maxvarMpar) ! local column default parameters
type(par_info),save,public :: basinParFallback(maxvarBpar) ! basin-average default parameters
! define vectors of metadata
type(var_info),save,public :: time_meta(maxvarTime) ! model time information
type(var_info),save,public :: forc_meta(maxvarForc) ! model forcing data
type(var_info),save,public :: attr_meta(maxvarAttr) ! local attributes
type(var_info),save,public :: type_meta(maxvarType) ! local classification of veg, soil, etc.
type(var_info),save,public :: id_meta(maxvarId) ! local classification of veg, soil, etc.
type(var_info),save,public :: mpar_meta(maxvarMpar) ! local model parameters for each HRU
type(var_info),save,public :: indx_meta(maxvarIndx) ! local model indices for each HRU
type(var_info),save,public :: prog_meta(maxvarProg) ! local state variables for each HRU
type(var_info),save,public :: diag_meta(maxvarDiag) ! local diagnostic variables for each HRU
type(var_info),save,public :: flux_meta(maxvarFlux) ! local model fluxes for each HRU
type(var_info),save,public :: deriv_meta(maxvarDeriv) ! local model derivatives for each HRU
type(var_info),save,public :: lookup_meta(maxvarLookup) ! local lookup tables for each HRU
type(var_info),save,public :: bpar_meta(maxvarBpar) ! basin parameters for aggregated processes
type(var_info),save,public :: bvar_meta(maxvarBvar) ! basin variables for aggregated processes
! ancillary metadata structures
type(flux2state), save,public :: flux2state_orig(maxvarFlux) ! named variables for the states affected by each flux (original)
type(flux2state), save,public :: flux2state_liq(maxvarFlux) ! named variables for the states affected by each flux (liquid water)
type(extended_info),save,public,allocatable :: averageFlux_meta(:) ! timestep-average model fluxes
! mapping from original to child structures
integer(i4b),save,public,allocatable :: forcChild_map(:) ! index of the child data structure: stats forc
integer(i4b),save,public,allocatable :: progChild_map(:) ! index of the child data structure: stats prog
integer(i4b),save,public,allocatable :: diagChild_map(:) ! index of the child data structure: stats diag
integer(i4b),save,public,allocatable :: fluxChild_map(:) ! index of the child data structure: stats flux
integer(i4b),save,public,allocatable :: indxChild_map(:) ! index of the child data structure: stats indx
integer(i4b),save,public,allocatable :: bvarChild_map(:) ! index of the child data structure: stats bvar
! child metadata structures
type(extended_info),save,public,allocatable :: statForc_meta(:) ! child metadata for stats
type(extended_info),save,public,allocatable :: statProg_meta(:) ! child metadata for stats
type(extended_info),save,public,allocatable :: statDiag_meta(:) ! child metadata for stats
type(extended_info),save,public,allocatable :: statFlux_meta(:) ! child metadata for stats
type(extended_info),save,public,allocatable :: statIndx_meta(:) ! child metadata for stats
type(extended_info),save,public,allocatable :: statBvar_meta(:) ! child metadata for stats
! ----------------------------------------------------------------------------------------------------------------
! * part 3: run time variables
! ----------------------------------------------------------------------------------------------------------------
! define the model decisions
type(model_options),save,public :: model_decisions(maxvarDecisions) ! the model decision structure
! define index variables describing the indices of the first and last HRUs in the forcing file
integer(i4b),save,public :: ixHRUfile_min ! minimum index
integer(i4b),save,public :: ixHRUfile_max ! maximum index
! define mapping structures
type(gru2hru_map),allocatable,save,public :: gru_struc(:) ! gru2hru map
type(hru2gru_map),allocatable,save,public :: index_map(:) ! hru2gru map
! define variables used for the vegetation phenology
real(rkind),dimension(12),save,public :: greenVegFrac_monthly ! fraction of green vegetation in each month (0-1)
! define the model output file
character(len=256),save,public :: fileout='' ! output filename
character(len=256),save,public :: output_fileSuffix='' ! suffix for the output file
! define controls on model output
logical(lgt),dimension(maxvarFreq),save,public :: finalizeStats=.false. ! flags to reset statistics
integer(i4b),save,public :: maxLayers ! maximum number of layers
integer(i4b),save,public :: maxSnowLayers ! maximum number of snow layers
! define control variables
integer(i4b),save,public :: startGRU ! index of the starting GRU for parallelization run
integer(i4b),save,public :: checkHRU ! index of the HRU for a single HRU run
integer(i4b),save,public :: iRunMode ! define the current running mode
integer(i4b),save,public :: nThreads=1 ! number of threads
integer(i4b),save,public :: ixProgress=ixProgress_id ! define frequency to write progress
integer(i4b),save,public :: ixRestart=ixRestart_never ! define frequency to write restart files
integer(i4b),save,public :: newOutputFile=noNewFiles ! define option for new output files
! define common variables
integer(i4b),save,public :: numtim ! number of time steps
integer(i4b),save,public :: nHRUrun ! number of HRUs in the run domain
integer(i4b),save,public :: nGRUrun ! number of GRUs in the run domain
real(rkind),save,public :: data_step ! length of the time_step
real(rkind),save,public :: refJulday ! reference time in fractional julian days
real(rkind),save,public :: refJulday_data ! reference time in fractional julian days (data files)
real(rkind),save,public :: dJulianStart ! julian day of start time of simulation
real(rkind),save,public :: dJulianFinsh ! julian day of end time of simulation
integer(i4b),save,public :: nHRUfile ! number of HRUs in the file
integer(i4b),save,public :: urbanVegCategory ! vegetation category for urban areas
logical(lgt),save,public :: doJacobian=.false. ! flag to compute the Jacobian
logical(lgt),save,public :: globalPrintFlag=.false. ! flag to compute the Jacobian
integer(i4b),save,public :: chunksize=1024 ! chunk size for the netcdf read/write
integer(i4b),save,public :: outputPrecision=nf90_double ! variable type
integer(i4b),save,public :: outputCompressionLevel=4 ! output netcdf file deflate level: 0-9. 0 is no compression.
! define result from the time calls
integer(i4b),dimension(8),save,public :: startInit,endInit ! date/time for the start and end of the initialization
integer(i4b),dimension(8),save,public :: startSetup,endSetup ! date/time for the start and end of the parameter setup
integer(i4b),dimension(8),save,public :: startRestart,endRestart ! date/time for the start and end to read restart data
integer(i4b),dimension(8),save,public :: startRead,endRead ! date/time for the start and end of the data read
integer(i4b),dimension(8),save,public :: startWrite,endWrite ! date/time for the start and end of the stats/write
integer(i4b),dimension(8),save,public :: startPhysics,endPhysics ! date/time for the start and end of the physics
! define elapsed time
real(rkind),save,public :: elapsedInit ! elapsed time for the initialization
real(rkind),save,public :: elapsedSetup ! elapsed time for the parameter setup
real(rkind),save,public :: elapsedRestart ! elapsed time to read restart data
real(rkind),save,public :: elapsedRead ! elapsed time for the data read
real(rkind),save,public :: elapsedWrite ! elapsed time for the stats/write
real(rkind),save,public :: elapsedPhysics ! elapsed time for the physics
! define ancillary data structures
type(var_i),save,public :: startTime ! start time for the model simulation
type(var_i),save,public :: finshTime ! end time for the model simulation
type(var_i),save,public :: refTime ! reference time for the model simulation
type(var_i),save,public :: oldTime ! time for the previous model time step
! output file information
logical(lgt),dimension(maxvarFreq),save,public :: outFreq ! true if the output frequency is desired
integer(i4b),dimension(maxvarFreq),save,public :: ncid ! netcdf output file id
! look-up values for the choice of the time zone information (formerly in modelDecisions module)
integer(i4b),parameter,public :: ncTime=1 ! time zone information from NetCDF file (timeOffset = longitude/15. - ncTimeOffset)
integer(i4b),parameter,public :: utcTime=2 ! all times in UTC (timeOffset = longitude/15. hours)
integer(i4b),parameter,public :: localTime=3 ! all times local (timeOffset = 0)
#ifdef ACTORS_ACTIVE
! global data structures are managed by FileAccessActor
type(var_forc),allocatable,save,public :: forcingDataStruct(:) ! forcingDataStruct(:)%var(:)%dataFromFile(:,:)
type(dlength),allocatable,save,public :: vecTime(:)
logical(lgt),allocatable,save,public :: failedHRUs(:) ! list of true and false values to indicate if an HRU has failed
type(ilength),allocatable,save,public :: outputTimeStep(:) ! timestep in output files
! inital conditions for Actors
type(init_cond),allocatable,save,public :: init_cond_prog(:) ! variable data for initial conditions
type(init_cond),allocatable,save,public :: init_cond_bvar(:) ! variable data for initial conditions
#else
! define metadata for model forcing datafile non-Actors
type(file_info),save,public,allocatable :: forcFileInfo(:) ! file info for model forcing data
! define indices in the forcing data files non-Actors
integer(i4b),save,public :: iFile=1 ! index of current forcing file from forcing file list
integer(i4b),save,public :: forcingStep=integerMissing ! index of current time step in current forcing file
integer(i4b),save,public :: forcNcid=integerMissing ! netcdf id for current netcdf forcing file
! define controls on model output non-Actors
integer(i4b),dimension(maxvarFreq),save,public :: statCounter=0 ! time counter for stats
integer(i4b),dimension(maxvarFreq),save,public :: outputTimeStep=0 ! timestep in output files
logical(lgt),dimension(maxvarFreq),save,public :: resetStats=.true. ! flags to reset statistics
! define common variables non-Actors
real(rkind),save,public :: fracJulday ! fractional julian days since the start of year
real(rkind),save,public :: tmZoneOffsetFracDay ! time zone offset in fractional days
integer(i4b),save,public :: yearLength ! number of days in the current year
#endif
! define fixed dimensions
integer(i4b),parameter,public :: nBand=2 ! number of spectral bands
integer(i4b),parameter,public :: nTimeDelay=2000 ! number of time steps in the time delay histogram (default: ~1 season = 24*365/4)
! printing step frequency
integer(i4b),parameter,public :: print_step_freq = 1000
character(len=1024),public,save :: fname ! temporary filename
END MODULE globalData
build/source/dshare/var_lookup.f90 deleted 100644 → 0
+ 0
− 953
View file @ e944a604
This diff is collapsed.
build/source/engine/coupled_em.f90 deleted 100644 → 0
+ 0
− 1421
View file @ e944a604
This diff is collapsed.
build/source/engine/derivforce.f90 deleted 100755 → 0
+ 0
− 353
View file @ e944a604
This diff is collapsed.
build/source/engine/mDecisions.f90 deleted 100755 → 0
+ 0
− 753
View file @ e944a604
This diff is collapsed.
build/source/engine/run_oneHRU.f90 deleted 100644 → 0
+ 0
− 245
View file @ e944a604
! SUMMA - Structure for Unifying Multiple Modeling Alternatives
! Copyright (C) 2014-2020 NCAR/RAL; University of Saskatchewan; University of Washington
!
! This file is part of SUMMA
!
! For more information see: http://www.ral.ucar.edu/projects/summa
!
! This program is free software: you can redistribute it and/or modify
! it under the terms of the GNU General Public License as published by
! the Free Software Foundation, either version 3 of the License, or
! (at your option) any later version.
!
! This program is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with this program. If not, see <http://www.gnu.org/licenses/>.
module run_oneHRU_module
! numerical recipes data types
USE nrtype
! data types
USE data_types,only:&
var_i, & ! x%var(:) (i4b)
var_d, & ! x%var(:) (rkind)
var_ilength, & ! x%var(:)%dat (i4b)
var_dlength, & ! x%var(:)%dat (rkind)
zLookup ! x%z(:)%var(:)%lookup(:) (rkind)
! access vegetation data
USE globalData,only:greenVegFrac_monthly ! fraction of green vegetation in each month (0-1)
USE globalData,only:overwriteRSMIN ! flag to overwrite RSMIN
USE globalData,only:maxSoilLayers ! Maximum Number of Soil Layers
! provide access to Noah-MP constants
USE module_sf_noahmplsm,only:isWater ! parameter for water land cover type
! provide access to the named variables that describe elements of parameter structures
USE var_lookup,only:iLookTYPE ! look-up values for classification of veg, soils etc.
USE var_lookup,only:iLookATTR ! look-up values for local attributes
USE var_lookup,only:iLookPARAM ! look-up values for local column model parameters
! provide access to the named variables that describe elements of variable structures
USE var_lookup,only:iLookPROG ! look-up values for local column model prognostic (state) variables
USE var_lookup,only:iLookDIAG ! look-up values for local column model diagnostic variables
USE var_lookup,only:iLookINDEX ! look-up values for local column index variables
! provide access to model decisions
USE globalData,only:model_decisions ! model decision structure
USE var_lookup,only:iLookDECISIONS ! look-up values for model decisions
! these are needed because we cannot access them in modules locally if we might use those modules with Actors
USE globalData,only:fracJulday ! fractional julian days since the start of year
USE globalData,only:yearLength ! number of days in the current year
USE globalData,only:tmZoneOffsetFracDay ! time zone offset in fractional days
! provide access to the named variables that describe model decisions
USE mDecisions_module,only: & ! look-up values for LAI decisions
monthlyTable,& ! LAI/SAI taken directly from a monthly table for different vegetation classes
specified ! LAI/SAI computed from green vegetation fraction and winterSAI and summerLAI parameters
! ----- global variables that are modified ------------------------------------------------------------------------------------------
! Noah-MP parameters
USE NOAHMP_VEG_PARAMETERS,only:SAIM,LAIM ! 2-d tables for stem area index and leaf area index (vegType,month)
USE NOAHMP_VEG_PARAMETERS,only:HVT,HVB ! height at the top and bottom of vegetation (vegType)
USE noahmp_globals,only:RSMIN ! minimum stomatal resistance (vegType)
! urban vegetation category (could be local)
USE globalData,only:urbanVegCategory ! vegetation category for urban areas
implicit none
private
public::run_oneHRU
contains
! ************************************************************************************************
! public subroutine run_oneGRU: simulation for a single GRU
! ************************************************************************************************
! simulation for a single HRU
subroutine run_oneHRU(&
! model control
hru_nc, & ! intent(in): hru index in netcdf
hruId, & ! intent(in): hruId
dt_init, & ! intent(inout): used to initialize the length of the sub-step for each HRU
computeVegFlux, & ! intent(inout): flag to indicate if we are computing fluxes over vegetation (false=no, true=yes)
nSnow,nSoil,nLayers, & ! intent(inout): number of snow and soil layers
! data structures (input)
timeVec, & ! intent(in): model time data
typeData, & ! intent(in): local classification of soil veg etc. for each HRU
attrData, & ! intent(in): local attributes for each HRU
lookupData, & ! intent(in): local lookup tables for each HRU
bvarData, & ! intent(in): basin-average variables
! data structures (input-output)
mparData, & ! intent(inout): local model parameters
indxData, & ! intent(inout): model indices
forcData, & ! intent(inout): model forcing data
progData, & ! intent(inout): prognostic variables for a local HRU
diagData, & ! intent(inout): diagnostic variables for a local HRU
fluxData, & ! intent(inout): model fluxes for a local HRU
! error control
err,message) ! intent(out): error control
! ----- define downstream subroutines -----------------------------------------------------------------------------------
USE module_sf_noahmplsm,only:redprm ! module to assign more Noah-MP parameters
USE derivforce_module,only:derivforce ! module to compute derived forcing data
USE coupled_em_module,only:coupled_em ! module to run the coupled energy and mass model
implicit none
! ----- define dummy variables ------------------------------------------------------------------------------------------
! model control
integer(i4b) , intent(in) :: hru_nc ! hru index in netcdf
integer(8) , intent(in) :: hruId ! hruId
real(rkind) , intent(inout) :: dt_init ! used to initialize the length of the sub-step for each HRU
logical(lgt) , intent(inout) :: computeVegFlux ! flag to indicate if we are computing fluxes over vegetation (false=no, true=yes)
integer(i4b) , intent(inout) :: nSnow,nSoil,nLayers ! number of snow and soil layers
! data structures (input)
integer(i4b) , intent(in) :: timeVec(:) ! int vector -- model time data
type(var_i) , intent(in) :: typeData ! x%var(:) -- local classification of soil veg etc. for each HRU
type(var_d) , intent(in) :: attrData ! x%var(:) -- local attributes for each HRU
type(zLookup) , intent(in) :: lookupData ! x%z(:)%var(:)%lookup(:) -- local lookup tables for each HRU
type(var_dlength) , intent(in) :: bvarData ! x%var(:)%dat -- basin-average variables
! data structures (input-output)
type(var_dlength) , intent(inout) :: mparData ! x%var(:)%dat -- local (HRU) model parameters
type(var_ilength) , intent(inout) :: indxData ! x%var(:)%dat -- model indices
type(var_d) , intent(inout) :: forcData ! x%var(:) -- model forcing data
type(var_dlength) , intent(inout) :: progData ! x%var(:)%dat -- model prognostic (state) variables
type(var_dlength) , intent(inout) :: diagData ! x%var(:)%dat -- model diagnostic variables
type(var_dlength) , intent(inout) :: fluxData ! x%var(:)%dat -- model fluxes
! error control
integer(i4b) , intent(out) :: err ! error code
character(*) , intent(out) :: message ! error message
! ----- define local variables ------------------------------------------------------------------------------------------
! local variables
character(len=256) :: cmessage ! error message
real(rkind) , allocatable :: zSoilReverseSign(:) ! height at bottom of each soil layer, negative downwards (m)
! initialize error control
err=0; write(message, '(A21,I0,A10,I0,A2)' ) 'run_oneHRU (hru nc = ',hru_nc -1 ,', hruId = ',hruId,')/' !netcdf index starts with 0 if want to subset
! ----- hru initialization ---------------------------------------------------------------------------------------------
! water pixel: do nothing
if (typeData%var(iLookTYPE%vegTypeIndex) == isWater) return
! get height at bottom of each soil layer, negative downwards (used in Noah MP)
allocate(zSoilReverseSign(nSoil),stat=err)
if(err/=0)then
message=trim(message)//'problem allocating space for zSoilReverseSign'
err=20; return
endif
zSoilReverseSign(:) = -progData%var(iLookPROG%iLayerHeight)%dat(nSnow+1:nLayers)
! populate parameters in Noah-MP modules
! Passing a maxSoilLayer in order to pass the check for NROOT, that is done to avoid making any changes to Noah-MP code.
! --> NROOT from Noah-MP veg tables (as read here) is not used in SUMMA
call REDPRM(typeData%var(iLookTYPE%vegTypeIndex), & ! vegetation type index
typeData%var(iLookTYPE%soilTypeIndex), & ! soil type
typeData%var(iLookTYPE%slopeTypeIndex), & ! slope type index
zSoilReverseSign, & ! * not used: height at bottom of each layer [NOTE: negative] (m)
maxSoilLayers, & ! number of soil layers
urbanVegCategory) ! vegetation category for urban areas
! deallocate height at bottom of each soil layer(used in Noah MP)
deallocate(zSoilReverseSign,stat=err)
if(err/=0)then
message=trim(message)//'problem deallocating space for zSoilReverseSign'
err=20; return
endif
! overwrite the minimum resistance
if(overwriteRSMIN) RSMIN = mparData%var(iLookPARAM%minStomatalResistance)%dat(1)
! overwrite the vegetation height
HVT(typeData%var(iLookTYPE%vegTypeIndex)) = mparData%var(iLookPARAM%heightCanopyTop)%dat(1)
HVB(typeData%var(iLookTYPE%vegTypeIndex)) = mparData%var(iLookPARAM%heightCanopyBottom)%dat(1)
! overwrite the tables for LAI and SAI
if(model_decisions(iLookDECISIONS%LAI_method)%iDecision == specified)then
SAIM(typeData%var(iLookTYPE%vegTypeIndex),:) = mparData%var(iLookPARAM%winterSAI)%dat(1)
LAIM(typeData%var(iLookTYPE%vegTypeIndex),:) = mparData%var(iLookPARAM%summerLAI)%dat(1)*greenVegFrac_monthly
end if
! ----- hru forcing ----------------------------------------------------------------------------------------------------
! compute derived forcing variables
call derivforce(timeVec, & ! vector of time information
forcData%var, & ! vector of model forcing data
attrData%var, & ! vector of model attributes
mparData, & ! data structure of model parameters
progData, & ! data structure of model prognostic variables
diagData, & ! data structure of model diagnostic variables
fluxData, & ! data structure of model fluxes
tmZoneOffsetFracDay & ! time zone offset in fractional days
err,cmessage) ! error control
if(err/=0)then; err=20; message=trim(message)//trim(cmessage); return; endif
! ----- run the model --------------------------------------------------------------------------------------------------
! initialize the number of flux calls
diagData%var(iLookDIAG%numFluxCalls)%dat(1) = 0._rkind
! run the model for a single HRU
call coupled_em(&
! model control
hruId, & ! intent(in): hruId
dt_init, & ! intent(inout): initial time step
1._rkind, & ! intent(in): used to adjust the length of the timestep with failure in Actors (non-Actors here, always 1)
computeVegFlux, & ! intent(inout): flag to indicate if we are computing fluxes over vegetation
! data structures (input)
typeData, & ! intent(in): local classification of soil veg etc. for each HRU
attrData, & ! intent(in): local attributes for each HRU
forcData, & ! intent(in): model forcing data
mparData, & ! intent(in): model parameters
bvarData, & ! intent(in): basin-average model variables
lookupData, & ! intent(in): lookup tables
! data structures (input-output)
indxData, & ! intent(inout): model indices
progData, & ! intent(inout): model prognostic variables for a local HRU
diagData, & ! intent(inout): model diagnostic variables for a local HRU
fluxData, & ! intent(inout): model fluxes for a local HRU
fracJulDay, & ! intent(in): fractional julian days since the start of year
yearLength, & ! intent(in): number of days in the current year
! error control
err,cmessage) ! intent(out): error control
if(err/=0)then; err=20; message=trim(message)//trim(cmessage); return; endif
! update the number of layers
nSnow = indxData%var(iLookINDEX%nSnow)%dat(1) ! number of snow layers
nSoil = indxData%var(iLookINDEX%nSoil)%dat(1) ! number of soil layers
nLayers = indxData%var(iLookINDEX%nLayers)%dat(1) ! total number of layers
end subroutine run_oneHRU
end module run_oneHRU_module
build/source/engine/sundials/coupled_em.f90 deleted 100644 → 0
+ 0
− 1424
View file @ e944a604
This diff is collapsed.
build/source/engine/sundials/mDecisions.f90 deleted 100755 → 0
+ 0
− 745
View file @ e944a604
This diff is collapsed.
build/source/engine/vegPhenlgy.f90 deleted 100755 → 0
+ 0
− 198
View file @ e944a604
This diff is collapsed.
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