mod_atmos_dyn_vars.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "inc_openmp.h"
module mod_atmos_dyn_vars
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_grid_index
  use scale_index
  use scale_tracer
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_dyn_vars_setup
  public :: atmos_dyn_vars_fillhalo
  public :: atmos_dyn_vars_restart_read
  public :: atmos_dyn_vars_restart_write

  public :: atmos_dyn_vars_restart_create
  public :: atmos_dyn_vars_restart_open
  public :: atmos_dyn_vars_restart_def_var
  public :: atmos_dyn_vars_restart_enddef
  public :: atmos_dyn_vars_restart_close

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  logical,               public :: atmos_dyn_restart_output                = .false.             

  character(len=H_LONG),  public :: atmos_dyn_restart_in_basename           = ''
  logical,                public :: atmos_dyn_restart_in_postfix_timelabel  = .false.             
  character(len=H_LONG),  public :: atmos_dyn_restart_out_basename          = ''
  logical,                public :: atmos_dyn_restart_out_postfix_timelabel = .true.              
  character(len=H_MID),   public :: atmos_dyn_restart_out_title             = 'ATMOS_DYN restart'
  character(len=H_SHORT), public :: atmos_dyn_restart_out_dtype             = 'DEFAULT'

  ! prognostic variables
  real(RP), public, allocatable :: prog(:,:,:,:)

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  integer,                private, parameter :: vmax = 100
  character(len=H_SHORT), private            :: var_name(vmax)
  character(len=H_MID),   private            :: var_desc(vmax)
  character(len=H_SHORT), private            :: var_unit(vmax)
  integer,                private            :: var_id(vmax)
  integer,                private            :: restart_fid = -1  ! file ID

  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_vars_setup
    use scale_process, only: &
       prc_mpistop
    use scale_const, only: &
       undef => const_undef
    use mod_atmos_admin, only: &
       atmos_dyn_type
    use scale_atmos_dyn_tstep_short, only: &
       atmos_dyn_tstep_short_regist
    implicit none

    namelist / param_atmos_dyn_vars / &
       atmos_dyn_restart_in_basename,           &
       atmos_dyn_restart_in_postfix_timelabel,  &
       atmos_dyn_restart_output,                &
       atmos_dyn_restart_out_basename,          &
       atmos_dyn_restart_out_postfix_timelabel, &
       atmos_dyn_restart_out_title,             &
       atmos_dyn_restart_out_dtype

    integer :: ierr
    integer :: iv
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[VARS] / Categ[ATMOS DYN] / Origin[SCALE-RM]'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_dyn_vars,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       write(*,*) 'xxx Not appropriate names in namelist PARAM_ATMOS_DYN_VARS. Check!'
       call prc_mpistop
    endif
    if( io_nml ) write(io_fid_nml,nml=param_atmos_dyn_vars)

    call atmos_dyn_tstep_short_regist( atmos_dyn_type, & ! [IN]
                                       va,             & ! [OUT]
                                       var_name,       & ! [OUT]
                                       var_desc,       & ! [OUT]
                                       var_unit        ) ! [OUT]

    if ( va > 0 ) then ! additional prognostic variables

       if ( va > vmax ) then
          write(*,*) 'xxx number of the prognostic variables is exceed the limit', va, ' > ', vmax
          call prc_mpistop
       endif
       allocate( prog(ka,ia,ja,va) )
       prog(:,:,:,:) = undef

       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** [ATMOS_DYN] prognostic/diagnostic variables'
       if( io_l ) write(io_fid_log,'(1x,A,A24,A,A48,A,A12,A)') &
                  '***       |', 'VARNAME                 ','|', &
                  'DESCRIPTION                                     ', '[', 'UNIT        ', ']'
       do iv = 1, va
          if( io_l ) write(io_fid_log,'(1x,A,I3,A,A24,A,A48,A,A12,A)') &
                     '*** NO.',iv,'|',var_name(iv),'|',var_desc(iv),'[',var_unit(iv),']'
       enddo

       if( io_l ) write(io_fid_log,*)
       if ( atmos_dyn_restart_in_basename /= '' ) then
          if( io_l ) write(io_fid_log,*) '*** Restart input?  : YES, file = ', trim(atmos_dyn_restart_in_basename)
          if( io_l ) write(io_fid_log,*) '*** Add timelabel?  : ', atmos_dyn_restart_in_postfix_timelabel
       else
          if( io_l ) write(io_fid_log,*) '*** Restart input?  : NO'
       endif
       if (       atmos_dyn_restart_output             &
            .AND. atmos_dyn_restart_out_basename /= '' ) then
          if( io_l ) write(io_fid_log,*) '*** Restart output? : YES, file = ', trim(atmos_dyn_restart_out_basename)
          if( io_l ) write(io_fid_log,*) '*** Add timelabel?  : ', atmos_dyn_restart_out_postfix_timelabel
       else
          if( io_l ) write(io_fid_log,*) '*** Restart output? : NO'
          atmos_dyn_restart_output = .false.
       endif

    else ! no additional prognostic variables

       allocate( prog(ka,ia,ja,1) ) ! for safety
       prog(:,:,:,:) = undef
       atmos_dyn_restart_output = .false.

    endif

    return
  end subroutine atmos_dyn_vars_setup

  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_vars_fillhalo
    use scale_comm, only: &
       comm_vars8, &
       comm_wait
    implicit none

    integer :: i, j, iv
    !---------------------------------------------------------------------------

    do iv = 1, va
       do j  = js, je
       do i  = is, ie
          prog(   1:ks-1,i,j,iv) = prog(ks,i,j,iv)
          prog(ke+1:ka,  i,j,iv) = prog(ke,i,j,iv)
       enddo
       enddo

       call comm_vars8( prog(:,:,:,iv), iv )
    enddo

    do iv = 1, va
       call comm_wait ( prog(:,:,:,iv), iv )
    enddo

    return
  end subroutine atmos_dyn_vars_fillhalo

  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_vars_restart_open
    use scale_time, only: &
       time_gettimelabel
    use scale_fileio, only: &
       fileio_open
    implicit none

    character(len=19)     :: timelabel
    character(len=H_LONG) :: basename
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '*** Open restart file (ATMOS_DYN) ***'

    if ( atmos_dyn_restart_in_basename /= '' ) then

       if ( atmos_dyn_restart_in_postfix_timelabel ) then
          call time_gettimelabel( timelabel )
          basename = trim(atmos_dyn_restart_in_basename)//'_'//trim(timelabel)
       else
          basename = trim(atmos_dyn_restart_in_basename)
       endif

       if( io_l ) write(io_fid_log,*) '*** basename: ', trim(basename)

       call fileio_open( restart_fid, basename )
    else
       if( io_l ) write(io_fid_log,*) '*** restart file for ATMOS_DYN is not specified.'
    endif

    return
  end subroutine atmos_dyn_vars_restart_open

  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_vars_restart_read
    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_fileio, only: &
       fileio_read, &
       fileio_flush
    implicit none

    real(RP) :: total
    integer  :: i, j, iv
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then
       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** Read from restart file (ATMOS_DYN) ***'

       do iv = 1, va
          call fileio_read( prog(:,:,:,iv),                          & ! [OUT]
                            restart_fid, var_name(iv), 'ZXY', step=1 ) ! [IN]
       enddo

       if ( io_aggregate ) then
          call fileio_flush( restart_fid )
          ! X/Y halos have been read from file

          ! fill K halos
          do iv = 1, va
             do j  = 1, ja
             do i  = 1, ia
                prog(   1:ks-1,i,j,iv) = prog(ks,i,j,iv)
                prog(ke+1:ka,  i,j,iv) = prog(ke,i,j,iv)
             enddo
             enddo
          enddo
       else
          call atmos_dyn_vars_fillhalo
       end if

       if ( statistics_checktotal ) then
          do iv = 1, va
             call stat_total( total, prog(:,:,:,iv), var_name(iv) )
          enddo
       endif
    else
       if( io_l ) write(io_fid_log,*) '*** invalid restart file ID for ATMOS_DYN.'
    endif

    return
  end subroutine atmos_dyn_vars_restart_read

  subroutine atmos_dyn_vars_restart_create
    use scale_time, only: &
       time_gettimelabel
    use scale_fileio, only: &
       fileio_create
    implicit none

    character(len=19)     :: timelabel
    character(len=H_LONG) :: basename
    !---------------------------------------------------------------------------

    if ( atmos_dyn_restart_out_basename /= '' ) then

       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** Create restart file (ATMOS_DYN) ***'

       if ( atmos_dyn_restart_out_postfix_timelabel ) then
          call time_gettimelabel( timelabel )
          basename = trim(atmos_dyn_restart_out_basename)//'_'//trim(timelabel)
       else
          basename = trim(atmos_dyn_restart_out_basename)
       endif

       if( io_l ) write(io_fid_log,*) '*** basename: ', trim(basename)

       call fileio_create( restart_fid,                                                       & ! [OUT]
                           basename, atmos_dyn_restart_out_title, atmos_dyn_restart_out_dtype ) ! [IN]

    endif

    return
  end subroutine atmos_dyn_vars_restart_create

  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_vars_restart_enddef
    use scale_fileio, only: &
       fileio_enddef
    implicit none

    if ( restart_fid /= -1 ) then
       call fileio_enddef( restart_fid ) ! [IN]
    endif

    return
  end subroutine atmos_dyn_vars_restart_enddef

  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_vars_restart_close
    use scale_fileio, only: &
       fileio_close
    implicit none
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then
       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** Close restart file (ATMOS_DYN) ***'

       call fileio_close( restart_fid ) ! [IN]

       restart_fid = -1
    endif

    return
  end subroutine atmos_dyn_vars_restart_close

  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_vars_restart_def_var
    use scale_fileio, only: &
       fileio_def_var
    implicit none

    integer iv
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then

       do iv = 1, va
          call fileio_def_var( restart_fid, var_id(iv), var_name(iv), var_desc(iv), &
                               var_unit(iv), 'ZXY', atmos_dyn_restart_out_dtype     ) ! [IN]
       enddo

    endif

    return
  end subroutine atmos_dyn_vars_restart_def_var

  !-----------------------------------------------------------------------------
  subroutine atmos_dyn_vars_restart_write
    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_fileio, only: &
       fileio_write_var
    implicit none

    real(RP) :: total
    integer  :: iv
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then

       call atmos_dyn_vars_fillhalo

       if ( statistics_checktotal ) then
          do iv = 1, va
             call stat_total( total, prog(:,:,:,iv), var_name(iv) )
          enddo
       endif

       do iv = 1, va
          call fileio_write_var( restart_fid, var_id(iv), prog(:,:,:,iv), var_name(iv), 'ZXY' ) ! [IN]
       enddo

    endif

    return
  end subroutine atmos_dyn_vars_restart_write

end module mod_atmos_dyn_vars