larmojdsango/main.f90 at main · jafPrudente/larmojdsango · 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
program main

   ! Verde que te quiero verde.

   !===============================================!
   !      Este programa evoluciona estrellas       !
   !===============================================!

   !------------------------------------------------
   ! Usamos el módulo 'arrays' para declarar arreglos
   ! y el módulo 'vars' para las variables de entrada.
   use arrays
   use vars

   implicit none

   ! Llamamos la subrutina que lee la información de entrada.
   call lector

   !------------------------------------------------
   ! Allojamos en memoria los arreglos generales para comenzar la evolución.
   allocate( r(Nr) )
   allocate( m(Nr), s(Nr) )
   allocate( a(Nr), alpha(Nr) )
   allocate( a_p(Nr), alpha_p(Nr) )
   allocate( adot(Nr) )

   !------------------------------------------------
   ! Evolucionamos en sí: tenemos una opción
   ! para cada tipo de materia.
   if(      matterType == 1 ) then
      call evolveEscalar

   else if( matterType == 2 ) then
      call evolveDirac

   else if( matterType == 3 ) then
      call evolveProca
   end if


end program