InitAnalyze Subroutine

public subroutine InitAnalyze()

Uses

  • proc~~initanalyze~~UsesGraph proc~initanalyze InitAnalyze module~modgvec_analyze_vars MODgvec_Analyze_Vars proc~initanalyze->module~modgvec_analyze_vars module~modgvec_globals MODgvec_Globals proc~initanalyze->module~modgvec_globals module~modgvec_mpi MODgvec_MPI proc~initanalyze->module~modgvec_mpi module~modgvec_readintools MODgvec_ReadInTools proc~initanalyze->module~modgvec_readintools module~modgvec_analyze_vars->module~modgvec_globals iso_fortran_env iso_fortran_env module~modgvec_globals->iso_fortran_env module~modgvec_readintools->module~modgvec_globals

Initialize Module

Arguments

None

Calls

proc~~initanalyze~~CallsGraph proc~initanalyze InitAnalyze getint getint proc~initanalyze->getint getintarray getintarray proc~initanalyze->getintarray getlogical getlogical proc~initanalyze->getlogical getrealarray getrealarray proc~initanalyze->getrealarray proc~getrealallocarray GETREALALLOCARRAY proc~initanalyze->proc~getrealallocarray proc~par_barrier par_Barrier proc~initanalyze->proc~par_barrier swrite swrite proc~initanalyze->swrite proc~getrealallocarray->swrite interface~findstr FindStr proc~getrealallocarray->interface~findstr proc~converttoproposalstr ConvertToProposalStr proc~getrealallocarray->proc~converttoproposalstr proc~count_sep count_sep proc~getrealallocarray->proc~count_sep mpi_barrier mpi_barrier proc~par_barrier->mpi_barrier interface~findstr->interface~findstr proc~remove_blanks remove_blanks proc~converttoproposalstr->proc~remove_blanks

Source Code

SUBROUTINE InitAnalyze
! MODULES
USE MODgvec_Globals,ONLY:UNIT_stdOut,fmt_sep
USE MODgvec_MPI,ONLY:par_Barrier
USE MODgvec_Analyze_Vars
USE MODgvec_ReadInTools,ONLY:GETINT,GETINTARRAY,GETREALARRAY,GETLOGICAL,GETREALALLOCARRAY
IMPLICIT NONE
!-----------------------------------------------------------------------------------------------------------------------------------
! INPUT VARIABLES
!-----------------------------------------------------------------------------------------------------------------------------------
! INPUT/OUTPUT VARIABLES
!-----------------------------------------------------------------------------------------------------------------------------------
! OUTPUT VARIABLES
!-----------------------------------------------------------------------------------------------------------------------------------
! LOCAL VARIABLES
REAL(wp):: visu_minmax(3,0:1)
!===================================================================================================================================
  CALL par_Barrier(beforeScreenOut='INIT ANALYZE ...')
  visu1D    = GETINT('visu1D',Proposal=-1)
  visu2D    = GETINT('visu2D',Proposal=-1)
  visu3D    = GETINT('visu3D',Proposal=-1)

  outFileType  = GETINT('outfileType',Proposal=1)

  visu_minmax(1:3,0)=GETREALARRAY("visu_min",3,Proposal=(/0.0_wp,0.0_wp,0.0_wp/))
  visu_minmax(1:3,1)=GETREALARRAY("visu_max",3,Proposal=(/1.0_wp,1.0_wp,1.0_wp/))

  iAnalyze=-1

  IF(visu1D.GT.0)THEN
    np_1d          = GETINT(     "np_1d",Proposal=5)
    IF(np_1d.LE.1) CALL abort(__STAMP__,&
       "np_1d must be >1")
  END IF
  IF(visu2D.GT.0)THEN
    np_visu_BC     = GETINTARRAY("np_visu_BC",2,Proposal=(/20,30/))
    IF(any(np_visu_BC.LE.1)) CALL abort(__STAMP__,&
       "all point numbers in np_visu_BC must be >1")
    visu_BC_minmax(2:3,0)=GETREALARRAY("visu_BC_min",2,Proposal=visu_minmax(2:3,0),quiet_def_in=.TRUE.)
    visu_BC_minmax(2:3,1)=GETREALARRAY("visu_BC_max",2,Proposal=visu_minmax(2:3,1),quiet_def_in=.TRUE.)
    np_visu_planes = GETINTARRAY("np_visu_planes",3, (/5,12,10/))
    IF(any(np_visu_planes.LE.1)) CALL abort(__STAMP__,&
       "all point numbers in np_visu_planes must be >1")
    visu_planes_minmax(1:3,0)=GETREALARRAY("visu_planes_min",3,Proposal=visu_minmax(1:3,0),quiet_def_in=.TRUE.)
    visu_planes_minmax(1:3,1)=GETREALARRAY("visu_planes_max",3,Proposal=visu_minmax(1:3,1),quiet_def_in=.TRUE.)
  END IF
  IF(visu3D.GT.0)THEN
    np_visu_3D     = GETINTARRAY("np_visu_3D",3,Proposal=(/5,12,10/))
    IF(any(np_visu_3D.LE.1)) CALL abort(__STAMP__,&
       "all point numbers in np_visu_3D must be >1")
    visu_3D_minmax(1:3,0)=GETREALARRAY("visu_3D_min",3,Proposal=visu_minmax(1:3,0),quiet_def_in=.TRUE.)
    visu_3D_minmax(1:3,1)=GETREALARRAY("visu_3D_max",3,Proposal=visu_minmax(1:3,1),quiet_def_in=.TRUE.)
  END IF
  SFLout    = GETINT('SFLout',Proposal=-1)
  IF(SFLout.GT.-1)THEN
    SFLout_mn_max = GETINTARRAY("SFLout_mn_max",2,Proposal=(/-1,-1/))
    SFLout_mn_pts = GETINTARRAY("SFLout_mn_pts",2,Proposal=(/40,40/)) !off by default
    SFLout_endpoint=GETLOGICAL("SFLout_endpoint",Proposal=.FALSE.)
    SFLout_relambda=GETLOGICAL("SFLout_relambda",Proposal=.TRUE.)
    CALL GETREALALLOCARRAY("SFLout_radialPos",SFLout_radialpos,SFLout_nrp,Proposal=(/1.0_wp/))
  END IF !SFLout
  CALL par_Barrier(afterScreenOut='... DONE')
  SWRITE(UNIT_stdOut,fmt_sep)
END SUBROUTINE InitAnalyze