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
|
! Copyright (C) 2005 Barbara Ercolano
!
! MoCaSSiNoutput = MOnte CArlo SImulationS of Nebulae
! this is the output-only driver of the simulation
! (requires grid1.out,grid2.out,grid3.out files)
!
! Version 2.02
program MoCaSSiNoutput
use common_mod
use constants_mod
use grid_mod
use iteration_mod
use output_mod
use set_input_mod
use xSec_mod
implicit none
include 'mpif.h'
type(grid_type) :: grid3D(maxGrids) ! the 3D Cartesian grid
integer :: iGrid
call mpi_init(ierr)
call mpi_comm_rank(MPI_COMM_WORLD, taskid, ierr)
call mpi_comm_size(MPI_COMM_WORLD, numtasks, ierr)
! starttime = mpi_wtime()
lgWarm = .true.
lgRecombination = .true.
if (taskid == 0) then
print*, "MOCASSIN 2005 output Version ",VERSION
print*, "Creating output files from current grid*.out files "
print*, " stored in the output/ directory"
print*, " "
end if
if (taskid == 0) then
print*, " "
end if
! reset the 3D cartesian grid
call resetGrid(grid3D)
call setStarPosition(grid3D(1)%xAxis,grid3D(1)%yAxis,grid3D(1)%zAxis,grid3D)
! initialize opacities x sections array
call initXSecArray()
! set the ionzing continuum according to the contShape variable
call setContinuum()
! prepare atomica data stuff
call makeElements()
if (taskid == 0) then
! determine final statistics
if (lgGas) call outputGas(grid3D)
end if
call mpi_barrier(mpi_comm_world, ierr)
! free all space allocated to the 3D grid
do iGrid=1, nGrids
call freeGrid(grid3D(iGrid))
end do
! endtime = mpi_wtime()
if (taskid == 0) then
! print*, "time: ", endtime-starttime
end if
call mpi_finalize(ierr)
stop 'mpi done'
end program MoCaSSiNoutput
|
