FileListModule.f90

Go to the documentation of this file.

 ! uses preprocessor
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!==============================================================================
! Module FileListModule  
! ***************************************** 
! contains
!  File declarations
!  subroutine  ReadControlParsAndInitFiles   
!==============================================================================
module filelistmodule
    use filehandlingmodule
!==============================================================================
! Type Declarations 
! *****************  
! for different text and netcdf files. The types are defined in FileHandlingModule
!
! USE  FileHandlingModule   
!==============================================================================
    character(len=:), allocatable :: resultpath, resultpath_h, statepath; ! defined in read control pars
    TYPE(file) :: log_file
    TYPE(file) :: statein_file
    TYPE(file) :: stateout_file
    TYPE(file) :: statebackupout_file
    TYPE(file) :: statebackupout_file_0
    TYPE(file) :: statebackupout_file_1
    TYPE(file) :: dispersalkernel_file
    TYPE(file) :: dispersalkernelnew_file
    TYPE(file) :: controlpars_file
    TYPE(file) :: controlparsout_file
    TYPE(file) :: species_file
    TYPE(file) :: immigration_file
    TYPE(file) :: bioclimate_file
    TYPE(file) :: stockarea_file
    TYPE(file) :: landcoverfilter_file
    TYPE(file) :: immigrationout_file
    TYPE(file) :: speciesout_file
    TYPE(file) :: biomassoutfile
    TYPE(file) :: numberoutfile
    TYPE(file) :: seedoutfile
    TYPE(file) :: antaoutfile
    TYPE(file) :: laioutfile   
    TYPE(file) :: pollenoutfile
    TYPE(file) :: nppoutfile
    TYPE(file) :: basareaoutfile
    TYPE(file) :: ingrowthoutfile
    TYPE(file) :: biodivoutfile
    TYPE(file) :: laipertreehcloutfile
    TYPE(file) :: dbhpertreehcloutfile
    TYPE(file) :: biompertreehcloutfile
    TYPE(file) :: seedprodfractpertreehcloutfile
    TYPE(file) :: heightperhcloutfile
    TYPE(file) :: lightdistrperhcloutfile
    TYPE(file) :: lightandleafareaperlightclass
    
    TYPE(file) :: heightstructoutfiles(16)
    
    TYPE ncvar
        character(len=:), allocatable :: name        
        character(len=:), allocatable :: attr_unit
        INTEGER :: varid = -1
        REAL , allocatable :: values2d(:,:)    
        REAL , allocatable :: values3d(:,:,:)        
        REAL , allocatable :: values4d(:,:,:,:)
    end TYPE
    INTEGER :: nc_file_id
    character(len=:), allocatable :: nc_file_name
    TYPE(ncvar) :: biodiv_nc
    TYPE(ncvar) :: biomass_nc
    TYPE(ncvar) :: number_nc
    TYPE(ncvar) :: seed_nc
    TYPE(ncvar) :: ingrowth_nc
    TYPE(ncvar) :: npp_nc
    TYPE(ncvar) :: basalarea_nc
    TYPE(ncvar) :: antagonist_nc
    TYPE(ncvar) :: pollen_nc
    TYPE(ncvar) :: lai_nc
    TYPE(ncvar) :: light_nc
    TYPE(ncvar) :: heightstruct_nc(16)
contains
 
!==============================================================================
!==============================================================================
    SUBROUTINE readcontrolparsandinitfiles()
        use loggermodule
        use filehandlingmodule, only: relativepathletter, pathletter, openfiler, closefile, file
        use all_par, only: specfilename, immifilename, bioclimfilename, stockareafilename,landcoverfilterfilename, experimentid
        IMPLICIT NONE
        character(len=128) :: controlFileName
        character(len=128) :: workingDir
        character(len=:), allocatable :: FileOutID, FileOutID2
        INTEGER :: i
        LOGICAL :: exists
        character(len=2) :: structID
!---- prepares the logfile and also the logging on the console       
        log_file = file(10,"logfile.txt")  
        call initlogger(log_file);        
        call loginfo("Starting TreeMig v2.0")
        call loginfo("Initializing files...")
!---- prepares reading of the control file       
!---- gets from the command line (console) the control filename (1st argument) and the working directory (2nd argument)    
        call get_command_argument(1, controlfilename) ! get the name of the control file from command line
        call get_command_argument(2, workingdir)
        if (trim(workingdir) /= '') then ! if a working directory is indicated set it
            call chdir(workingdir)
        end if
 
        if (trim(controlfilename) == '') then  ! if no name for control file is given request one 
            do 
                write(*,*) 'Please enter control-file name:'
                read(*,*) controlfilename
                INQUIRE(file=relativepathletter//"C"//pathletter//trim(controlfilename), exist=exists)
                if(exists) then
                    exit
                else
                    write(*,*) 'Control-file not found:'//relativepathletter//"C"//pathletter//trim(controlfilename)
                end if
            end do
        end if
              
!---- read control pars file, parts of the information is needed for defining the other files!
        controlpars_file        = file(11,relativepathletter//"C"//pathletter//trim(controlfilename))
        controlparsout_file     = file(12,null()) ! path defined in ReadControlPars()
        call readcontrolpars() ! 
 
!---- define (unit, path, name)  other input files       
        statein_file            = file(71, relativepathletter//"InitialValues.txt")
        dispersalkernel_file    = file(100,relativepathletter//"dispKern.txt")
        dispersalkernelnew_file = file(110,relativepathletter//"newdispkernels.txt")
        species_file            = file(20, relativepathletter//'P'//pathletter//specfilename)
        immigration_file        = file(21, relativepathletter//'P'//pathletter//immifilename)
        bioclimate_file         = file(22, relativepathletter//'E'//pathletter//bioclimfilename)
        stockarea_file          = file(23, relativepathletter//'E'//pathletter//stockareafilename)
        landcoverfilter_file    = file(24, relativepathletter//'E'//pathletter//landcoverfilterfilename)
 
!---- define (unit, path, name)  output files       
        fileoutid = experimentid//"_"
        stateout_file           = file(72,statepath//"statefile_0000.txt")
        statebackupout_file     = file(73,statepath//"statefile_backup.txt")
        statebackupout_file_0   = file(74,statepath//"statefile_backup_0.txt")
        statebackupout_file_1   = file(75,statepath//"statefile_backup_1.txt")
 
        immigrationout_file     = file(200, resultpath//fileoutid//"immigration"//".pars")
        speciesout_file         = file(201, resultpath//fileoutid//"species"//".pars")        
 
        biompertreehcloutfile          = file(220, resultpath_h//fileoutid2//"bioPerTreeHCL"//".txt")
        laipertreehcloutfile           = file(221, resultpath_h//fileoutid2//"laiPerTreeHCL"//".txt")
        dbhpertreehcloutfile           = file(222, resultpath_h//fileoutid2//"dbhPerTreeHCL"//".txt")
        seedprodfractpertreehcloutfile = file(223, resultpath_h//fileoutid2//"seedProdPerTreeHCL"//".txt")
        heightperhcloutfile            = file(224, resultpath_h//fileoutid2//"heights"//".txt")
        lightdistrperhcloutfile        = file(225, resultpath_h//fileoutid2//"lightDistr"//".txt")
        lightandleafareaperlightclass  = file(226, resultpath //"lightLAI_LCL"//".txt")
 
        biomassoutfile     = file(300, resultpath//'Biomass'//'.txt')
        numberoutfile      = file(301, resultpath//'Number'//'.txt')
        seedoutfile        = file(302, resultpath//'Seeds'//'.txt')
        antaoutfile        = file(303, resultpath//'Antagonists'//'.txt')
        pollenoutfile      = file(304, resultpath//'Pollen'//'.txt')       !pollen
        nppoutfile         = file(305, resultpath//'NPP'//'.txt')          !NPP
        biodivoutfile      = file(306, resultpath//'Biodiversity'//'.txt') !biodiv
        laioutfile         = file(307, resultpath//'LAI'//'.txt')          !LAI
        basareaoutfile     = file(308, resultpath//'BA'//'.txt')           !basal area
        ingrowthoutfile    = file(309, resultpath//'InGr'//'.txt')         !ingrowth
 
        do i=0,15
            write(structid,"(I2.2)") i
            heightstructoutfiles(i+1)  = file(400+i, resultpath_h//'HeightStruct_'//structid//'.txt')
        end do      
        
!---- Define  (unit, path, name)  netcdf output file and netcdf variables
 
        nc_file_name = resultpath //fileoutid//'OUT'//'.nc'
 
        biodiv_nc       =  ncvar("Biodiversity"    , "1")
        biomass_nc      =  ncvar("Biomass"         , "t/ha")
        number_nc       =  ncvar("Number"          , "1")
        seed_nc         =  ncvar("seedBank"         , "1")
        ingrowth_nc     =  ncvar("Ingrowth"        , "1") 
        npp_nc          =  ncvar("NPP"             , "t/ha") ! net primary productivity
        basalarea_nc    =  ncvar("BasalArea"       , "m^2/h")
        antagonist_nc   =  ncvar("Antagonist"      , "1")
        pollen_nc       =  ncvar("Pollen"          , "percent") ! pollen percent
        lai_nc          =  ncvar("LAI"             , "m^2/m^2") !"leaf_area_index", 
        light_nc        =  ncvar("LightDistr"      , "1") !"light distribution", 
        do i=0,15
            write(structid,"(I2.2)") i
            heightstruct_nc(i+1)  = ncvar("HeightStruct_"//structid, "t/ha")
        end do
    end SUBROUTINE readcontrolparsandinitfiles
    !==============================================================================
 
end module filelistmodule
!==============================================================================