
#############################################################
# Running Simulations
#############################################################

#Include KG_Module to run simulation
include("../src/KG_Module.jl")

p0 = Param(

    #Temporal variables
    tmax           = 1000.0     , 
    out_every_t    = 0.1        ,  
    deltat         = 0.05       , 

    #Spacial variables
    xmin           =  -30.0     ,  
    xmax           =  30.0      , 
    xnodes         =  128       ,    

    #Spacial variables
    ymin           =  -30.0     ,  
    ymax           =  30.0      , 
    ynodes         =   128      , 
    
    #Initial Conditions
    A0              =   3.5     , 
    σ               =   1.5     , 
    r0              =   17.0    ,
    ω               =   0.5     , 
    m               =   2       , 
    
    #Klein Gordon Mass
    μ               =   0.0     ,

    #Dissipation terms
    dissipation     =  :true    ,   
    
    R_orbit         = 0.0       ,
    R1              = 1.5       ,
    R2              = 0.0       ,
       
    Ω               = 0.45      , #(3)^(-3/2) ,      ,
    alpha1          = 10.0      ,
    alpha2          = 0.0       ,    

    #Periodic BC?
    Boundaries      = :radial   , # :square :periodic 
    
    #Data variables
    folder          = "/test_runs/"  ,       #must have a bar before the folder
    fname           = "run_004"

)

solve_wave_equation_2D(p0)


#############################################################
# Analysis
#############################################################

#Include plot utils
include("../src/KG_Plot_Module.jl")

#name of file to analyse
#filename = "./test_runs/run_003.h5"

#Slices Observer (also works with filename)
slices_observer(p0)

#Generate animation
#generate_animation_heatmap_makie2( filename , "test_animation.mp4")