mod_da_admin.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module mod_da_admin
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: da_admin_setup
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  character(len=H_SHORT), public :: da_type = 'OFF'
 
  logical, public :: da_do ! do data assimilation?
 
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine da_admin_setup
    use scale_prc, only: &
       prc_abort
    implicit none
 
    namelist / param_da / &
       da_type
 
    integer :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("DA_ADMIN_setup",*) 'Setup'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_da,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("DA_ADMIN_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("DA_ADMIN_setup",*) 'Not appropriate names in namelist PARAM_DA. Check!'
       call prc_abort
    endif
    log_nml(param_da)
 
    !-----< module component check >-----
 
    log_newline
    log_info("DA_ADMIN_setup",*) 'Data Assimilation components '
 
    if ( da_type /= 'OFF' ) then
       log_info_cont(*) 'Data Assimilation: ON, ', trim(da_type)
       da_do = .true.
    else
       log_info_cont(*) 'Data Assimilation: OFF'
       da_do = .false.
    endif
 
    return
  end subroutine da_admin_setup
 
end module mod_da_admin