Paso 1)  Correr  
              computes_km_mum_from_v33_v55_tordillo_quintuco_VM3_VM4.f  
              (VM4 tiene velocidades mas altas (10% menos que los datos de Ariel) que el background de VM (20% menos que VM4)

Paso 2) Correr 
                 genero_modulos_gassmann_densities_formations.f
        con input file 
                 inp_genero_modulus_Q_T_VM3_VM4
                 
          Genera los archivos     
                p_modulus.freq,  s_modulus.freq,   bulk_density.mat
                
                
Paso 3) Correr                 
                v16_onda_plana_v1.f
                
        con input file
                 input_QT_VM_for_v16
                 
