forked from garrelt/C2-Ray3Dm
-
Notifications
You must be signed in to change notification settings - Fork 0
/
time_ini.F90
100 lines (73 loc) · 2.7 KB
/
time_ini.F90
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
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
!>
!! \brief This module handles the time variables
!!
!! Module for C2-Ray (f90)
!!
!! \b Author: Garrelt Mellema
!!
!! \b Date: 20-Aug-2006
!<
module times
! This module handles the time variables
use precision, only: dp
use my_mpi
use file_admin, only: stdinput, file_input
use astroconstants, only: YEAR
use cosmology, only: t0, zred_t0, zred2time
implicit none
integer :: number_timesteps !< Number of time steps between redshift slices
integer :: number_outputs !< Number of outputs between redshift slices
contains
! =======================================================================
!>
!! Initializes the module variables
!<
subroutine time_ini( )
! Initializes number of time steps per frame (integration and output)
! Author: Garrelt Mellema
! Date: 20-Aug-2006 (f77: 10-Mar-2004)
#ifdef MPI
integer :: ierror
#endif
if (rank == 0) then
! Ask for number of time steps
if (.not.file_input) &
write(*,'(A,$)') 'Enter number of time steps between slices: '
read(stdinput,*) number_timesteps
! Ask for interval between outputs
if (.not.file_input) &
write(*,'(A,$)') 'Enter number of outputs between slices: '
read(stdinput,*) number_outputs
endif
#ifdef MPI
! Distribute the input parameters to the other nodes
call MPI_BCAST(number_timesteps,1,MPI_INTEGER,0,&
MPI_COMM_NEW,ierror)
call MPI_BCAST(number_outputs,1,MPI_INTEGER,0,MPI_COMM_NEW,ierror)
#endif
end subroutine time_ini
!=======================================================================
!>
!! Sets the time steps for the interval zred0-zred_end
!! The time steps are constant, given by dividing the redshift
!! range into equal time intervals (number_timesteps)
!<
subroutine set_timesteps (zred0,zred_end,end_time,dt,output_dt)
! Sets time steps for calculation and outputting
! Author: Garrelt Mellema
! Date: 20-Aug-2006 (f77: 19-May-2005)
real(kind=dp),intent(in) :: zred0 !< starting redshift
real(kind=dp),intent(in) :: zred_end !< ending redshift
real(kind=dp),intent(out) :: end_time !< end time
real(kind=dp),intent(out) :: dt !< time step
real(kind=dp),intent(out) :: output_dt !< output time step
real(kind=dp) :: current_time
! Convert to time (in seconds)
current_time=zred2time(zred0) !t0*( ( (1.0+zred_t0)/(1.0+zred0) )**1.5-1.0)
end_time=zred2time(zred_end) !t0*( ( (1.0+zred_t0)/(1.0+zred_end) )**1.5-1.0)
! Set value of time step
dt=(end_time-current_time)/real(number_timesteps)
! Convert to time
output_dt=(end_time-current_time)/real(number_outputs)
end subroutine set_timesteps
end module times