module INITIAL More...

Functions/Subroutines
subroutine, public	mkinit_setup
	Setup.

subroutine, public	mkinit (output, model_mesh, atm_prgvars_manager, atm_auxvars_manager, atm_trcvars_manager)
	Driver.

Variables
integer, public	mkinit_type = -1

integer, public	i_ignore = 0

Detailed Description

module INITIAL

Description: subroutines for preparing initial data

Author: Yuta Kawai, Team SCALE

NAMELIST

PARAM_MKINIT

name type default value comment

INITNAME character(len=H_SHORT) 'NONE'

History Output: No history output

Function/Subroutine Documentation

◆ mkinit_setup()

subroutine, public mod_mkinit::mkinit_setup

Setup.

Definition at line 72 of file mod_mkinit.F90.

    implicit none
 
    character(len=H_SHORT) :: initname = 'NONE'
 
    namelist / param_mkinit / &
      initname
 
    integer :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("MKINIT_setup",*) 'Setup'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_mkinit,iostat=ierr)
    if( ierr < 0 ) then !--- missing
      log_info("MKINIT_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
      log_error("MKINIT_setup",*) 'Not appropriate names in namelist PARAM_MKINIT. Check!'
      call prc_abort
    endif
    log_nml(param_mkinit)
 
    select case(trim(initname))
    case('NONE')
      mkinit_type = i_ignore
    end select
    
    return

References i_ignore, and mkinit_type.

Referenced by mod_dg_prep::dg_prep().

◆ mkinit()

subroutine, public mod_mkinit::mkinit	(	logical, intent(out)	output,
		class(atmosmesh), intent(in), target	model_mesh,
		class(modelvarmanager), intent(inout)	atm_prgvars_manager,
		class(modelvarmanager), intent(inout)	atm_auxvars_manager,
		class(modelvarmanager), intent(inout)	atm_trcvars_manager )

Driver.

Definition at line 107 of file mod_mkinit.F90.

  
    use scale_model_var_manager, only: modelvarmanager
    use mod_atmos_mesh, only: atmosmesh
    use mod_atmos_vars, only: &
      atmosvars_getlocalmeshprgvars
 
    implicit none
 
    logical, intent(out) :: output
    class(AtmosMesh), target, intent(in) :: model_mesh
    class(ModelVarManager), intent(inout) :: atm_prgvars_manager
    class(ModelVarManager), intent(inout) :: atm_auxvars_manager
    class(ModelVarManager), intent(inout) :: atm_trcvars_manager
 
    class(LocalMeshFieldBase), pointer :: DDENS, MOMX, MOMY, MOMZ, DRHOT
    class(LocalMeshFieldBase), pointer :: DENS_hyd, PRES_hyd
    class(LocalMeshFieldBase), pointer :: Rtot, CPtot, CVtot
 
    integer :: n
    integer :: ke
    class(LocalMesh3D), pointer :: lcmesh3D
    class(MeshBase3D), pointer :: mesh
    !---------------------------------------------------------------------------
 
    mesh => model_mesh%ptr_mesh
 
    if ( mkinit_type == i_ignore ) then
      log_newline
      log_progress(*) 'skip  making initial data'
      output = .false.
    else
      log_newline
      log_progress(*) 'start making initial data'
 
      ! call PROF_rapstart('_MkInit_main',3)   
      
      do n=1, mesh%LOCAL_MESH_NUM
        call atmosvars_getlocalmeshprgvars( n, &
          mesh, atm_prgvars_manager, atm_auxvars_manager, &
          ddens, momx, momy, momz, drhot,                 &
          dens_hyd, pres_hyd, rtot, cvtot, cptot,         &
          lcmesh3d                                        )
 
        !$omp parallel do
        do ke=lcmesh3d%NeS, lcmesh3d%NeE
          momx %val(:,ke) = 0.0_rp
          momy %val(:,ke) = 0.0_rp
          momz %val(:,ke) = 0.0_rp
          ddens%val(:,ke) = 0.0_rp
          drhot%val(:,ke) = 0.0_rp
 
          cptot%val(:,ke) = cpdry
          cvtot%val(:,ke) = cvdry
          rtot%val(:,ke) = rdry
        end do
      end do
      ! call PROF_rapend  ('_MkInit_main',3)
        
      ! LOG_PROGRESS(*) 'end   making initial data'
      output = .true.
    end if
 
    return