-
Notifications
You must be signed in to change notification settings - Fork 0
/
Main.f
68 lines (55 loc) · 2.03 KB
/
Main.f
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
program main
use GlobalSetup
IMPLICIT NONE
character*500 VersionGIT
integer i
real*8 starttime,stoptime
call SetOutputMode(.true.)
call cpu_time(starttime)
call getcommandline()
call GetOutputDir
open(unit=9,file=trim(outputdir) // "log.dat",FORM="FORMATTED",ACCESS="STREAM")
call output("Output dir: " // trim(outputdir))
call output("==================================================================")
call output(" ARtful modelling code for exoplanet Science - ARCiS")
call output("==================================================================")
c terms of use
call output("By using ARCiS you agree to the terms of use.")
call output("It basically means you consult with any of the authors before")
call output("publication if there is any doubt on the outcome of the results.")
call output("==================================================================")
call output("Let's get the show on the road!!")
call output("ARCiS version "//trim(VersionGIT()))
call output("==================================================================")
call Init()
if(nobs.ne.0) call ReadObs()
call cpu_time(stoptime)
call output("Initialisation time: " // trim(dbl2string((stoptime-starttime),'(f10.2)')) // " s")
call output("==================================================================")
if(opacitymode) then
if(dochemistry) then
call OnlyChemCompute
else
do i=1,nmol
mixrat_r(:,i)=mixrat(i)
enddo
endif
call SetupOpacities()
call WriteOpacityFITS()
else if(dopostequalweights) then
call PostEqualWeights()
else if(domakeai) then
call MakeAI()
else if(retrieval) then
call DoRetrieval()
else
call ComputeModel(.true.)
call WriteStructure()
call WriteOutput()
endif
call cpu_time(stoptime)
call output("==================================================================")
call output("Total runtime: " // trim(dbl2string((stoptime-starttime),'(f10.2)')) // " s")
call output("All done!")
call output("==================================================================")
end