TimeLoop.f90

Go to the documentation of this file.

!==============================================================================
!===============================================================
SUBROUTINE timeloop(yearStart, yearEnd, &
                    latStart, latEnd, lonStart, lonEnd, &
                    currLatStart, currLatEnd, currLonStart, currLonEnd, &
                    facvar, &
                    nspc, doDispersal, ncallSpatSimu, doParallel, doFFT)
 
    ! <CALL modules for TYPE definitions and +- fixed variables>---------------------------------------------------------------
    use all_par, only:  simustartyear, stateoutputdist, stateoutputyears, numyrs
!    use TMMPITypes, only: master
    use loggermodule, only: logmessage, loginfo
    IMPLICIT NONE
 
    ! <Passed variables>--------------------------------------------------------------------------------
    INTEGER, INTENT(in)    :: yearStart, yearEnd, &  ! start and end of simulated period
                              latStart, latEnd, lonStart, lonEnd, &  ! start and end of simulated areaINTEGER :: nspc
                              currLatStart, currLatEnd, currLonStart, currLonEnd, &  ! part of total simulation area currently simulated, e.g. for spin-up. Is not subarea!
                              nspc, ncallSpatSimu                                   ! number of species involved [=spec]
    REAL, INTENT(in)       :: facvar
    LOGICAL, INTENT(in)    :: doDispersal, doParallel, doFFT                                  ! flag: simulate dispersal?
 
    ! <Local variables>--------------------------------------------------------------------------------
    INTEGER :: year             ! looping variable for simulated year
 
    ! <Here the SUBROUTINE starts>**********************************************************************
 
    looptime: do year = yearstart, yearend
        write (logmessage, "(A10,I8,A7,L)") "    Year: ", year, " doDispersal: ", dodispersal
        call loginfo(logmessage)
!---- Reads the bioclimate if necessary. In the current version bioclimate is read in each year.        
        call readnewbioclimifnecessary(year)
!---- Run the model for this time step
        if (doparallel) then
! #ifdef WITHMPI
!             call MastersParallelWorkForDistributedAndDispersalSimus( &
!                 year, &
!                 latStart, latEnd, lonStart, lonEnd, &
!                 currLatStart, currLatEnd, currLonStart, currLonEnd, &
!                 facvar, &
!                 nspc, doDispersal)
! #endif
        else
            call spatialdynonetimestep(year, & 
                                       currlatstart, currlatend, currlonstart, currlonend, &
                                       nspc, dodispersal, dofft)
        end if
!---- Write out the current state to the statefile if wished at this year
        if (((mod(year, stateoutputdist) == 0) &
            .or.    any((year + simustartyear) == stateoutputyears)  &
            .or.    (year == numyrs-1)) &
            .and.   (year > 0))&
            call writestate(year, latstart, latend, lonstart, lonend, nspc)  ! Store intermediate results all 400 years
    end do looptime
    
! #ifdef WITHMPI
!     if (doParallel) call BCastTimeLoopIsReady(ncallSpatSimu, thisyear) ! Comment out if run without MPI
! #endif
    !    call measuretime(master,20) ! Comment out if run without MPI
 
end SUBROUTINE timeloop
!done