All_par.f90

Go to the documentation of this file.

!!============================================================================== 
!!
!  Name   :  All_par
!
!  Remarks:  contains the global parameters
!
!==============================================================================
! design          : H. Lischke, N. Zimmermann
! author(s)       : H. Lischke, N. Zimmermann
! implementation  : H. Lischke, N. Zimmermann
! cleaner         : T.J. Loeffler
! copyright       : (c) 1999, 03 by H. Lischke
!==============================================================================
module all_par
 
!****************** CONTROL PARAMETERS ***************************************!
!----- Important values are read in ReadControlPars from file control.pars
 
    !----- Basic setup for simulation
    INTEGER numyrs
    INTEGER simustartyear
    INTEGER, parameter:: climatestartyear = 0
    INTEGER, dimension  (10, 3) :: reportintervals
    INTEGER                     :: nreportintervals
    INTEGER                     :: stateoutputdist
    INTEGER, dimension (10)      :: stateoutputyears
    ! INTEGER, parameter   :: maxArraylat            = 2001! 1800 ! 500! 220 ! 1001! 1132 !2201        !> maximal number of rows, only for declarations
    ! INTEGER, parameter   :: maxArraylon            = 2001! 1800  ! 500! 180! 1834             !> maximal number of colums, only for declarations
 
    INTEGER, parameter :: maxcctimes = 1
    INTEGER, parameter :: maximmigrationsperspec = 1
    ! INTEGER, parameter :: maxKernelDimFine     =  150!1024! 1024 !1601! 300 !150 ! 600 !200!  1601                     !> maximal dimension of coarse kernel
    ! INTEGER, parameter :: maxKernelDim          = 150! 512! 512! 600 !25  ! 600 !200!  1601                     !> maximal dimension of coarse kernel
    ! INTEGER, parameter :: maxSeededCells        = (2*maxKernelDim+1)**2                 !> max. number of cells which receive seeds
    ! INTEGER, parameter :: maxSeededCells         = (2*maxKernelDim*maxArrayLon)                !> max. number of cells which receive seeds
    INTEGER, parameter :: nmeasuredtimes = 200
    INTEGER, parameter :: maxnumbprocs = 64
!     ! INTEGER             :: requests(maxArraylat*maxArraylon)
    INTEGER maxlat
    INTEGER maxlon
  !  INTEGER maxlatEven, maxlonEven                                  !>  chosen number of colums, if necessary made even
!     INTEGER maxlatExt, maxlonExt                                                 !> extended chosen number of colums
!     INTEGER:: maxlatEvenext, maxlonEvenext                           !> next larger power of 2
    REAL ::  latrealstart, latrealend, lonrealstart, lonrealend
    INTEGER ncctimes
    INTEGER startlatstart, startlatend, startlonstart, startlonend
    REAL unitofspatialdata
    REAL unitofspatialdata2
    INTEGER inoculationtime
    LOGICAL readstatefile                                        !! read initial state from file statefile?
    REAL, parameter         :: biothreshold = 0.1
    LOGICAL includeenv
    LOGICAL envfromfile
    
    INTEGER, parameter :: edge = 0
    REAL alpha_all
    REAL cellsidelength
    REAL epskernel
    REAL, parameter::     seednumberthresh = 0.0000
    INTEGER kerneltype
    LOGICAL, parameter :: dodispersal = .true.
    LOGICAL active
    LOGICAL dostochseeddisp
    LOGICAL seedsdifferent
    LOGICAL mathgtdifferent
    LOGICAL dispersaldifferent
    LOGICAL dofft
    LOGICAL mastseeding
    LOGICAL calcseedbank
    LOGICAL readdispersalkernel
    LOGICAL withseedantagonists
    REAL seedantagraz
    REAL seedantaeff
    REAL seedantamort
    REAL seedantarain
    LOGICAL withseedcarrcap
    REAL seedcarrcap
    LOGICAL immigration
    INTEGER immiperiod
    LOGICAL disturbances
    REAL disturbprob
    REAL disturbintensity
    LOGICAL tabsepoutput
    LOGICAL :: estabeqhc0
    LOGICAL :: contlightdepgerm
 
    character(len=:), allocatable :: boundaries
    INTEGER   :: outerseeds
    LOGICAL    :: cellnotyetallocated
    !----- Basic model setup
!    LOGICAL, parameter :: slow = .FALSE.                     !> slowing down initial growth (yes/no)
    INTEGER, parameter :: maxspc = 70
    INTEGER, parameter :: maxhc = 15
    INTEGER, parameter :: maxlc = 10
    character(len=:), allocatable:: specfilename
    character(len=:), allocatable:: immifilename
    character(len=:), allocatable:: bioclimfilename
    character(len=:), allocatable:: stockareafilename, landcoverfilterfilename
    character(len=:), allocatable:: experimentid
    INTEGER             ::          ndrivers
    character(len=30)   ::          columnnamesofdrivers (1:50)
    INTEGER, parameter  ::   nalloweddrivernames=10
    character(len=30), dimension(nallowedDrivernames),parameter   ::  alloweddrivernames=&
             ["lon        ","lat        ","year       ",&
              "DD         ","WiT        ","DrStr      ", &
              "Distu      ","browsPress ","germDrought","nutrients  "]
    INTEGER, dimension(nallowedDrivernames)::                         columnlocations
 !   REAL, dimension(nallowedDrivernames), parameter ::   driverDefault=[0,0,0,0,0,0,0,8,0,100] !> default values for the drivers, if no explicit values are read in. The defaults should be read in  
    REAL, dimension(nallowedDrivernames) ::   driverdefault
    !     character(len=:), allocatable,dimension(7),parameter   ::  allowedDriverNames = &
!               ["lon  ","lat  ","year ","DD  ","WiT  ","DrStr","Distu"]  !> Names of allowed drivers, x, y, year, bioclim., + optional  drivers  
  !----- Basic parameters
    REAL, parameter :: pi = 3.141592654
    REAL, parameter :: birht = 1.37
    REAL, parameter :: toplt = 1.0
    REAL facvar
    REAL fcvarinit
    REAL fcvardyn
    REAL            :: pltsiz
    REAL, parameter :: potbirth = 0.5
    REAL            :: seedlcrowndiamincm
 
    !----- Basic variables
    REAL :: a1(5)
    REAL :: a2(5)
    REAL :: c1(2)
    REAL :: c2(2)
    REAL :: kn1(3)
    REAL :: kn2(3)
 
 
    TYPE wroutp
         LOGICAL :: biomass, number, seeds, antagonists, hstruct, pollen, lai, basalarea, &
                npp, ingrowth, biodiv, light, netcdf
    end TYPE wroutp 
 
    TYPE(wroutp) :: writeoutput
 
    !----- Overall result variables
    REAL    :: sumdiversity
    REAL    :: avdiversity80, avdiversity160, avdiversity240, avdiversity500, avdiversityeq
    REAL    :: biocrit(10)
    ! REAL    :: specAppearance(10, maxspc, maxArraylat, maxArraylon)
    ! REAL    :: specMigSpeed(10, maxspc, maxArraylat, maxArraylon)
    INTEGER :: nonzerospecs80, nonzerospecs160, nonzerospecs240, nonzerospecs500, nonzerospecseq
    INTEGER :: domspec80(3), domspec160(3), domspec240(3), domspec500(3), domspeceq(3)
    ! REAL        :: kernel(-maxKernelDim:maxKernelDim,-maxKernelDim:maxKernelDim)
    REAL, allocatable  :: kernel(:, :), kernelfine(:, :), kernel1(:, :), kernel2(:, :)
 
    !natalie:
    REAL    :: stmvol
 
!****************** SPECIES TYPES AND PARAMETERS  ***************************************!
!****************** User: do NOT TOUCH!>           ***************************************!
    TYPE idata
        INTEGER :: year
        REAL    :: rellat
        REAL    :: rellon
    end TYPE idata
 
    TYPE immigrationdata
        INTEGER     :: n
        TYPE(idata) :: dat(maximmigrationsperspec)
    end TYPE immigrationdata
 
    TYPE specproperties          ! scalar species properties, and kernels
        INTEGER               :: nb
        INTEGER               :: stn
        INTEGER               :: kdmax
        INTEGER               :: khmax
        INTEGER               :: kamax
        INTEGER               :: kg
        INTEGER               :: kddmin
        INTEGER               :: kddmax
        INTEGER               :: kwit
        INTEGER               :: kntol
        INTEGER               :: kbrow
        INTEGER               :: klighs
        INTEGER               :: kligha
        INTEGER               :: klq
        INTEGER               :: kimm
        INTEGER               :: resld
        INTEGER               :: ploto
        INTEGER               :: par0
        INTEGER               :: par1
        INTEGER               :: par2
        INTEGER               :: par3
        INTEGER               :: par4
        INTEGER               :: par5
        INTEGER               :: red
        INTEGER               :: green
        INTEGER               :: blue
        INTEGER               :: par9
        INTEGER               :: tb
        INTEGER               :: maxhtc
        INTEGER               :: mxseed
        INTEGER               :: rad
        INTEGER               :: maxlatext
        INTEGER               :: maxlonext
        INTEGER               :: maxlatevenext
        INTEGER               :: maxlonevenext
!      INTEGER               :: kdimi_FT, kdimj_FT !>  even dimension of kernel
        REAL                  :: kdrt
        REAL                  :: kgermdrsens
        REAL                  :: par6
        REAL                  :: par7
        REAL                  :: par8
        REAL                  :: maxdhm
        REAL                  :: maxgro
        REAL                  :: cmort
        REAL                  :: seedgerm
        REAL                  :: seedsurv
        REAL                  :: seedloss
        REAL                  :: lthres
        REAL                  :: envdep
        REAL                  :: enveta
        REAL                  :: minmat
        REAL                  :: maxseed
        REAL                  :: dispfac
        REAL                  :: alfa1
        REAL                  :: alfa2
        REAL                  :: birth
        REAL                  :: period
        REAL                  :: mathcp
        REAL                  :: repfac
        REAL                  :: sumtimebio
        ! REAL                  :: kernel(-maxKernelDim:maxKernelDim,-maxKernelDim:maxKernelDim)
        REAL, allocatable     :: kernel(:, :)
        complex, allocatable  :: kernel_ft(:, :)
        character(len=:), allocatable         :: nl
        character(len=:), allocatable         :: ns
        character(1)          :: stb
        TYPE(immigrationdata) :: immi
    end TYPE specproperties
 
    TYPE specpropertiesincl
        REAL :: htmin(0:maxhc + 1)
        REAL :: diain(0:maxhc)
        REAL :: lain(0:maxhc)
        REAL :: bioin(0:maxhc)
        REAL :: hlgro(0:maxhc, maxlc)
        REAL :: envvitality(0:maxhc, maxlc)
        REAL :: growthvitality(0:maxhc, maxlc)
        REAL :: seedprod(0:maxhc)
        REAL :: lightdepa(1:maxlc)
        REAL :: lightdepe(1:maxlc)
    end TYPE specpropertiesincl
 
    TYPE specstate
        REAL    :: sb    = 0
        REAL    :: producedseeds = 0
!        REAL    :: caches         !> cache site area (areaper grid-cell area)         [=..] !>DSP
        REAL    :: htmps
        REAL    :: htdmp
        REAL    :: antagonist
        REAL    :: newseeds
        REAL    :: verynewseeds = 0
        REAL    :: seedbank
        REAL    :: ingrowth   = 0
        REAL    :: spectotalbiomass = 0
        REAL    :: npp
        REAL    :: numin(0:maxhc)
        REAL    :: seeds(0:maxhc) = 0
!        LOGICAL :: there          !> LOGICAL array -> seeds arehere?            [~immiYear]
    end TYPE specstate
 
    TYPE seedcspec
        INTEGER :: lat
        INTEGER :: lon
        REAL    :: seeds = 0
    end TYPE seedcspec
 
    TYPE seedc
        TYPE(seedcspec) ::sp(maxspc)
    end TYPE seedc
 
    TYPE seedsspec
        REAL    :: newseeds = 0
        REAL    :: verynewseeds = 0
    end TYPE seedsspec
 
    TYPE seedspecarr
        TYPE(seedsspec) ::sp(maxspc)
    end TYPE seedspecarr
 
    TYPE specsstateincell
        TYPE(specstate) :: sp(maxspc)
    end TYPE specsstateincell
 
    TYPE biocliminonecell
        REAL :: mudd, stddd, muwit, stdwit, prop0drstr, mudrstr, stddrstr, disturb, avnit, brwpr, germdrought
    end TYPE biocliminonecell
 
    TYPE bioclimdata
        INTEGER                :: cctime
        TYPE(biocliminonecell), dimension(:, :), allocatable:: bc
    end TYPE bioclimdata
 
    ! TYPE BIOCLIMDATA
    !   INTEGER                :: cctime
    !   TYPE(BioClimInOneCell) :: bc(1:maxArraylat,1:maxArraylon)
    ! end TYPE BIOCLIMDATA
 
    TYPE currstateincell
        TYPE(specstate)       :: sp(maxspc)
        INTEGER               :: lat, lon
        INTEGER               :: maxcarpetlength
        REAL                  :: ddegs, wtemp, drstr, avnit, brwpr, disturb, germdrought
    end TYPE currstateincell
 
    TYPE currenvincell
        REAL                  :: ddegs, wtemp, drstr, avnit, brwpr
    end TYPE currenvincell
 
    TYPE newseedsincell
        TYPE(seedsspec)                 :: sp(maxspc)
    end TYPE newseedsincell
 
    TYPE yearcharacteristics
        INTEGER   :: year
        REAL      :: facvar
        Logical   :: dodispersal
    end TYPE yearcharacteristics
 
    TYPE outputline
        INTEGER   :: iam, lat, lon, year
        REAL      :: num(maxspc)
        REAL      :: sum
    end TYPE outputline
 
    TYPE(currstateincell):: thiscell
 
    ! TYPE(CurrStateInCell):: stateGrid(1:maxArraylat,1:maxArraylon) 
    TYPE(currstateincell), dimension(:, :), allocatable :: stategrid
 
    TYPE(newseedsincell), dimension(:, :), allocatable:: seedrain
 
 
    TYPE(specproperties)      :: spec(maxspc)
    TYPE(specpropertiesincl)  :: speccl(maxspc)
    TYPE(yearcharacteristics) :: thisyear
 
    TYPE(outputline):: outputl
 
    TYPE(bioclimdata) oldbioclim, newbioclim
 
    ! REAL    :: stockability(1:maxArraylat,1:maxArraylon), landcoverfilter(1:maxArraylat,1:maxArraylon)
    ! REAL    :: stocksum(1:maxArraylat ) , stocksumcum(1:maxArraylat ) !> summed over lons, used for load balancing
    REAL, dimension(:, :), allocatable     :: stockability, landcoverfilter
    REAL, dimension(:), allocatable     :: stocksum, stocksumcum
    REAL    :: hcwidth
    INTEGER :: htmax
    ! REAL    ::   avbio(1:maxArraylat,1:maxArraylon,1:maxspc,0:3)
    REAL, allocatable, dimension(:, :, :, :)  ::   avbio
 
    double precision, dimension(nmeasuredtimes)::  oldtime, sumtime
    double precision::  measuredinoculationtime
!****************** STAND TYPES AND PARAMETERS ***************************************!
!****************** User: do NOT TOUCH!>           ***************************************!
 
    REAL    :: ltfrq(maxlc)
    REAL    :: ltcum(maxlc)
    REAL    :: ltprf(0:maxhc, maxlc)
    REAL    :: ltval(maxlc)
    REAL    :: la(maxlc)
 
    INTEGER :: htcl
    INTEGER :: ltcl
 
end module all_par