output_empire.f90

Go to the documentation of this file.

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Time-stamp: <2016-10-17 12:19:06 pbrowne>
!!!
!!!    Module that stores information about outputting from empire
!!!    Copyright (C) 2015  Philip A. Browne
!!!
!!!    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/>.
!!!
!!!    Email: p.browne @ reading.ac.uk
!!!    Mail:  School of Mathematical and Physical Sciences,
!!!           University of Reading,
!!!       Reading, UK
!!!       RG6 6BB
!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Module output_empire
  use, intrinsic :: iso_fortran_env, only : stderr=>error_unit,stdout=>output_unit
  implicit none
  integer,parameter :: emp_e=stderr
  integer,parameter :: emp_o=stdout
  integer,parameter :: unit_nml=10
  integer,parameter :: unit_obs=11
  integer,parameter :: unit_truth=12
  integer,parameter :: unit_weight=13
  integer,parameter :: unit_mean=14
  integer,parameter :: unit_state=15
  integer,parameter :: unit_ens_rmse=16
  integer,parameter :: unit_mat_tri=17
  integer,parameter :: unit_spatial_rmse=18
  integer,parameter :: unit_variance=19
  integer,parameter :: unit_hist_read=20
  integer,parameter :: unit_hist_write=21
  integer,parameter :: unit_hist_readt=22
  integer,parameter :: unit_hist_readp=23
  integer,parameter :: unit_traj_read=24
  integer,parameter :: unit_traj_write=24
  integer,parameter :: unit_vardata=25
  
contains

  subroutine open_emp_o(id_num)
    implicit none
    integer, intent(in) :: id_num
    character(14) :: filename
    character(10) :: basename
    logical :: opend
    character(3), parameter :: msnul='nul'
    character(9), parameter :: unixnul='/dev/null'
    character(10), parameter :: emp_out_name='empire.nml'
    logical :: file_exists
    integer :: ios
    
    namelist/empire_output/basename

    inquire(file=emp_out_name,exist=file_exists)
    if(file_exists) then
       open(unit_nml,file=emp_out_name,iostat=ios,action='read'&
            &,status='old',form='formatted')
       if(ios .ne. 0) then
          print*,'Cannot open ',emp_out_name
          stop 'open_emp_o ERROR' 
       end if
       read(unit_nml,nml=empire_output,iostat=ios)
       close(unit_nml)
       if(ios .ne. 0) basename ='emp.out'
    else
       basename = 'emp.out'
    end if
    
    select case(basename)
    case(msnul)
       filename = msnul
    case(unixnul)
       filename = unixnul
    case default
       !define filename appropriately
       write(filename,'(A,A,i0)') trim(basename),'.',id_num
    end select

    print*,"this process' output will henceforth be directed to ",filename

    inquire(unit=emp_o,opened=opend)
    if(opend) then
       close(emp_o)
    end if

    open(emp_o,file=filename,action='write',status='replace'&
         &,form='formatted')

    inquire(unit=emp_e,opened=opend)
    if(opend) then
       close(emp_e)
    end if
    write(filename,'(A,i0)') 'emp.err.',id_num
    open(emp_e,file=filename,action='write',status='replace'&
         &,form='formatted')
    
  end subroutine open_emp_o


  subroutine close_emp_o()
    close(emp_o)
    close(emp_e)
  end subroutine close_emp_o
End Module output_empire