PIP/main.f90 at master · AstroSnow/PIP · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
program main
  !====================================================================
  ! This is the main program for the PIP2 code. (2014/07/08)
  ! Author A.Hillier
  ! First version - 2013/04/23
  ! Second version - 2013/05/29 AH - Main changes to put the do loop on the top level
  ! Third version -2013/06/04 NN bug fix
  ! Fourth version - 2013/06/28 AH Changed dinmension switch to inside solver routines
  ! Fifth version - 2013/07/10 Moved scheme switch into solver routines
  ! PIP2-----------------------
  ! First version - 2014/07/08 NN
  !====================================================================
  !  use globalvar,only:U_h,U_m,nt,time,dt,flag_stop,tend,nmax
  use globalvar,only:U_h,U_m,nt,time,dt,flag_stop,tend,nmax,nout,my_rank, flag_IR,flag_PIP,&
	esav,emsavtime,U_m_backup,U_h_backup
  use initial_rot,only:prologue
  use io_rot,only:output,stop_sim,epilogue
  use solver_rot,only:run_solver
  use scheme_rot,only:cfl
  use pip_rot,only:set_IR
  implicit none
  double precision::dtnew,dtmod
  dtmod=1.d0
  !----------------------------------------------------------------------|
  ! initial setting
  !----------------------------------------------------------------------|
  call prologue
  !----------------------------------------------------------------------|
  ! Time intergration
  !----------------------------------------------------------------------|
  do while(time .lt. tend)
    if((flag_IR.ge.1) .and. (flag_pip.ge.1)) then
       call set_IR(U_h,U_m)
    endif
     ! calculate cfl condition  -----------------------------------------
     call cfl(U_h,U_m,dt)
     dt=dt*dtmod
     !Back up the array
     U_m_backup=U_m
     U_h_backup=U_h
     ! Time intergration  -----------------------------------------------
     call run_solver(U_h,U_m)
     !Check change in dt
     if((flag_IR.ge.1) .and. (flag_pip.ge.1)) then
       call set_IR(U_h,U_m)
     endif
     call cfl(U_h,U_m,dtnew)
     !print*,maxval(U_m(:,:,:,2)),maxval(U_m_backup(:,:,:,2))
     !print*,dt,dtnew
     !if (dtnew .lt. dt/10.d0) then
     !	dtmod=0.1d0*dtmod
     !   if (my_rank .eq. 0) print *, 'Reducing dt by ',dtmod,time
     !   U_m=U_m_backup
     !   U_h=U_h_backup
     !else
		 ! Output of simulation data  ---------------------------------------
		 if (dt.ne.dt) print *, my_rank,dt
		 time=time+dt !! NB: dt can be modified by STS routines
		 call output(0)
		 dtmod=1.d0
	     ! Error check  -----------------------------------------------------
		 call stop_sim
		 if (flag_stop .eq.1) then
		    call output(1)
		    exit
		 endif
	 !endif
  enddo
  !----------------------------------------------------------------------|
  ! End of simulation
  !----------------------------------------------------------------------|
  print *, my_rank, tend,time
  call epilogue
end program main