solucion=solucion[:n+1,:,:]
answer=input('Quieres guardar la solucion para generar una animación? [y/n]\n')
if answer =='y':
Datos['solucion_animada']=solucion
Datos['Zcambio']=Zcambio
# Animacion_Contorno(xg,yg,Zcambio,n)
# Animacion_Superficie(xg,yg,Zcambio,n)
elif answer =='n':
pass
else:
pass
Datos['error']=errores
Datos['xg']=xg
Datos['yg']=yg
Datos['solucion']=u
Datos['Nt']=n
Datos['ht']=ht
hdf5.saveParametros(out_file_name,Datos)
# =============================================================================
# Generación de la gráfica
# =============================================================================
#figura 1: mapa de contornos
GrafContornos(xg, yg, u, 8, 0.75, 'inferno')
#figura 2: superficie
Graf3D(xg, yg, u, 'inferno')
#figura 3: error
GrafError(errores,'C2-')
print(mensaje)
sys.exit(1)
#Generamos parametros
Datos = {
'ax': 0, #Punto inicial en dirección x
'bx': 1, #Punto final en dirección x
'ay': 0, #Punto inicial en dirección y
'by': 1, #Punto final en dirección y
'Nx': 60, # Número de nodos en x
'Ny': 60, # Número de nodos en y
'Tx1': 40, # Condición de frontera 1 en x
'Tx2': -20, #Condición de frontera 2 en x
'Ty1': 3000, # Condición de frontera 1 en y
'Ty2': -100, # Condición de frontera 2 en x
'Tini': 120, # Condición inicial
'ht': 0.01, # Espaciamiento en tiempo
'Tmax': 3, # Tiempo máximo
'Tolerancia': 1E-4, # Tolerancia
'kappa_x': 1,
'kappa_y': 2,
'fuente': 0,
'c_p': 1,
'rho': 2,
'vel_x': 1,
'vel_y': 1,
}
hdf5.saveParametros(in_file_name, Datos)