PrepareOutput.f90

Go to the documentation of this file.

!==============================================================================
!
!  Name   :  PrepareOutput
!
!  Remarks:  contains the SUBROUTINEs
!
!            PrepareOutput (<TreeMig)
!            CloseOutputFiles (<TreeMig)
!
!==============================================================================
 
 
!==============================================================================
!==============================================================================
SUBROUTINE prepareoutput(nspc)
    
    use all_par, only: writeoutput
    use filelistmodule, only: laipertreehcloutfile, dbhpertreehcloutfile, biompertreehcloutfile,&
                              seedprodfractpertreehcloutfile, heightperhcloutfile, &
                              lightandleafareaperlightclass,&
                              biomassoutfile, numberoutfile, seedoutfile, antaoutfile, pollenoutfile, &
                              laioutfile, ingrowthoutfile,  nppoutfile, basareaoutfile, lightdistrperhcloutfile,&
                              biodivoutfile, heightstructoutfiles
    IMPLICIT NONE
 
    ! Passed variables>--------------------------------------------------------------------------------
    INTEGER, INTENT(in)       :: nspc              ! site index
 
    ! Local variables>---------------------------------------------------------------------------------
    INTEGER      :: i
    ! prepare constant output files (.txt) for additional information 
    ! species specific, in height classes
    call openandwriteheader(laipertreehcloutfile, nspc)           ! leaf area  
    call openandwriteheader(dbhpertreehcloutfile, nspc)           ! dbh  
    call openandwriteheader(biompertreehcloutfile, nspc)          ! biomass  
    call openandwriteheader(seedprodfractpertreehcloutfile, nspc) ! fraction of total seed production
    ! for all species
    call openandwriteheader(heightperhcloutfile, nspc)            ! height
    ! in light classes
    call openandwriteheaderlight(lightandleafareaperlightclass)     
    !----- 
    if(writeoutput%netcdf) return
    !----- only for .txt output
    !----- open file for biomass output
    if (writeoutput%biomass)     call openandwriteheader(biomassoutfile, nspc)
    !----- open file for number output  
    if (writeoutput%number)      call openandwriteheader(numberoutfile, nspc)
    !----- open file for seedBank output - only works if the seeds are stored in the seedBank variable
    if (writeoutput%seeds)       call openandwriteheader(seedoutfile, nspc)
    !----- open file for antagonist output
    if (writeoutput%antagonists) call openandwriteheader(antaoutfile, nspc)
    !----- open file for pollen output
    if (writeoutput%pollen)      call openandwriteheader(pollenoutfile, nspc)
     !----- open file for LAI output
    if (writeoutput%lai)         call openandwriteheader(laioutfile, nspc)
   !----- open file for ingrowth output
    if (writeoutput%ingrowth)    call openandwriteheader(ingrowthoutfile, nspc)
    !----- open file for NPP output
    if (writeoutput%NPP)         call openandwriteheader(nppoutfile, nspc)
   !----- open file for biodiversity output
    if (writeoutput%biodiv)      call openandwriteheader(biodivoutfile, nspc)
    !---- open file for basal area output
    if (writeoutput%basalArea)   call openandwriteheader(basareaoutfile, nspc)
    !---- open file for height structure output
    if (writeoutput%light)       call openandwriteheaderlight(lightdistrperhcloutfile)
 
    !---- open file for height structure output
    if (writeoutput%hstruct) then
        do i = 1, 16
            call openandwriteheader(heightstructoutfiles(i), nspc)
        end do
    end if
end SUBROUTINE prepareoutput
!==============================================================================
!==============================================================================
SUBROUTINE closeoutputfiles
    use filehandlingmodule,only : closefile
    use filelistmodule, only: laipertreehcloutfile, dbhpertreehcloutfile, biompertreehcloutfile,&
                              seedprodfractpertreehcloutfile, heightperhcloutfile,&
                              biomassoutfile, numberoutfile, seedoutfile, antaoutfile, pollenoutfile, &
                              laioutfile, ingrowthoutfile,  nppoutfile, basareaoutfile, lightdistrperhcloutfile,&
                              biodivoutfile, heightstructoutfiles, nc_file_id
    use netcdf, only: nf90_close
    use all_par, only: writeoutput
    IMPLICIT NONE
    INTEGER :: i 
    
    call closefile(laipertreehcloutfile)
    call closefile(dbhpertreehcloutfile)
    call closefile(biompertreehcloutfile)
    call closefile(seedprodfractpertreehcloutfile)
    call closefile(heightperhcloutfile)
    
    if(writeoutput%netcdf)   then
        call m_check(nf90_close(nc_file_id), __line__, __file__)
        return ! in this case no need to close the other files, they were never opened
    end if
 
    if (writeoutput%biomass)     call closefile(biomassoutfile)
    if (writeoutput%number)      call closefile(numberoutfile)
    if (writeoutput%seeds)       call closefile(seedoutfile)
    if (writeoutput%antagonists) call closefile(antaoutfile)
    if (writeoutput%pollen)      call closefile(pollenoutfile)    
    if (writeoutput%lai)         call closefile(laioutfile)
    if (writeoutput%ingrowth)    call closefile(ingrowthoutfile)    
    if (writeoutput%basalArea)   call closefile(basareaoutfile)
    if (writeoutput%NPP)         call closefile(nppoutfile)
    if (writeoutput%biodiv)      call closefile(biodivoutfile)
    if (writeoutput%light)       call closefile(lightdistrperhcloutfile)
    if (writeoutput%hstruct) then
        do i = 1, 16
            call closefile(heightstructoutfiles(i))
        end do
    end if
end SUBROUTINE
 
!done