# SI HAY ALGUN PROBLEMA CON LOS VALORES, ES DE ESCALA NO DE UNIDADES.
# usamos metros para calcular beta --> eps por lo tanto usamos TODO en metros.
Lx = 4000000 # dimenciones de la cuenca en metros
Ly = 2000000
nx = 200 # puntos de grilla
ny = 100
beta = 10**-11 # en metros
D = 2000 # metros
#arrays (y,x) !
from Cargar import cargar
psi_temp1, vort_temp1, psiF1, vortF1, QG_diag1, QG_curlw1, X1, Y1, dx1, dy1 = cargar(
'C:/Users/Jimena/Desktop/FCEN/circulación/prácticas/TP1/out_tmp1/', Lx, Ly,
nx, ny)
psi_temp2, vort_temp2, psiF2, vortF2, QG_diag2, QG_curlw2, X2, Y2, dx2, dy2 = cargar(
'C:/Users/Jimena/Desktop/FCEN/circulación/prácticas/TP1/out_tmp2/', Lx, Ly,
nx, ny)
psi_temp3, vort_temp3, psiF3, vortF3, QG_diag3, QG_curlw3, X3, Y3, dx3, dy3 = cargar(
'C:/Users/Jimena/Desktop/FCEN/circulación/prácticas/TP1/out_tmp3/', Lx, Ly,
nx, ny)
###----------------1----------------###
# dimensionalizando
Tau = 0.25
U = (2 * np.pi * Tau) / (1025 * D * (beta) * (Lx))
corriente1 = psiF1 * U * Lx
l = 4000000 #tamaño de la cuenca 4000km
D = 2000 #profundidad en metros
b = 0.00000000001 #funcion beta
ro = 1025 # densidad
tau = 0.3 #esfuerzo del viento
U = ((2 * (np.pi) * tau) / (ro * D * b * l)) #velocidad U
# Defino las variables
# psi_temp= Funcion corriente para cada t
# vort_temp= vorticidad para cada t
# psiF=funcion corriente final (estacionario)
# vortF=vorticidad final (estacionaria)
# Para EPS=0.29
psi_temp1, vort_temp1, psiF1, vortF1, QG_diag1, QG_curlw1, X1, Y1, dx1, dy1 = cargar(
dir_salida1, 4000, 2000, 200, 100)
# Para EPS=0.75
psi_temp2, vort_temp2, psiF2, vortF2, QG_diag2, QG_curlw2, X2, Y2, dx2, dy2 = cargar(
dir_salida2, 4000, 2000, 200, 100)
# Para EPS=1.16
psi_temp3, vort_temp3, psiF3, vortF3, QG_diag3, QG_curlw3, X3, Y3, dx3, dy3 = cargar(
dir_salida3, 4000, 2000, 200, 100)
# Paso a dimensionar la funcion corriente y la vorticidad para cada corrida.
# Para EPS=0.29
psiF1_dim = psiF1 * U * l
vortF1_dim = vortF1 * (U / l)
# Para EPS=0.79
psiF2_dim = psiF2 * U * l
vortF2_dim = vortF2 * (U / l)
# Para EPS=1.29
"""
Editor de Spyder
Este es un archivo temporal
"""
import numpy as np
from matplotlib import pyplot as plt
dir_salida = 'C:/Users/Jimena/Desktop/FCEN/circulación/prácticas/TP1/out_tmp3/'
Lx = 4000000
Ly = 2000000
nx = 200
ny = 100
beta = 10e-11
from Cargar import cargar
psi_temp, vort_temp, psiF, vortF, QG_diag, QG_curlw, X, Y, dx, dy = cargar(
dir_salida, Lx, Ly, nx, ny)
L = Lx
Tau = 0.25
U = (2 * np.pi * Tau) / (1025 * 2000 * (beta) * (L))
corriente = psiF * U * L
vort = vortF * (U / L) #chequear lo de la escala
plt.plot(QG_diag[:, 3], 'c', label='Energia cinetica')
plt.xlabel('Energia cinetica')
plt.ylabel('Tiempo')
plt.legend()
plt.savefig('Energiacin3.png')
plt.show()
plt.contourf(corriente)
from matplotlib import pyplot as plt
# SI HAY ALGUN PROBLEMA CON LOS VALORES, ES DE ESCALA NO DE UNIDADES.
# usamos metros para calcular beta --> eps por lo tanto usamos TODO en metros.
Lx = 4000000 # dimenciones de la cuenca en metros
Ly = 2000000
nx = 200 # puntos de grilla
ny = 100
beta = 10**-11 # en metros
D = 2000 # metros
#arrays (y,x) !
from Cargar import cargar
psi_temp1,vort_temp1,psiF1,vortF1,QG_diag1,QG_curlw1,X1,Y1,dx1,dy1=cargar("D:/Facultad/Materias/Circulacion/TP1/out_tmp1/",Lx,Ly,nx,ny)
psi_temp2,vort_temp2,psiF2,vortF2,QG_diag2,QG_curlw2,X2,Y2,dx2,dy2=cargar('D:/Facultad/Materias/Circulacion/TP1/out_tmp2/',Lx,Ly,nx,ny)
psi_temp3,vort_temp3,psiF3,vortF3,QG_diag3,QG_curlw3,X3,Y3,dx3,dy3=cargar('D:/Facultad/Materias/Circulacion/TP1/out_tmp3/',Lx,Ly,nx,ny)
###----------------1----------------###
# dimensionalizando
Tau = 0.25
U = (2*np.pi*Tau)/(1025*D*(beta)*(Lx))
corriente1 = psiF1*U*Lx
corriente2 = psiF2*U*Lx
corriente3 = psiF3*U*Lx
vort1 = vortF1*(U/Lx)
vort2 = vortF2*(U/Lx)
"""
### TP3 ###
import numpy as np
from matplotlib import pyplot as plt
Lx = 1500000 # dimensiones de la cuenca en metros
Ly = 1500000
nx = 150 # puntos de grilla
ny = 150
beta = 2*10**-11 # en metros
D = 2000 # metros
from Cargar import cargar
psi_temp1,vort_temp1,psiF1,vortF1,QG_diag1,QG_curlw1,X1,Y1,dx1,dy1=cargar("D:/Facultad/Materias/Circulacion/TP3/out_tmp1/",Lx,Ly,nx,ny)
psi_temp2,vort_temp2,psiF2,vortF2,QG_diag2,QG_curlw2,X2,Y2,dx2,dy2=cargar('D:/Facultad/Materias/Circulacion/TP3/out_tmp2/',Lx,Ly,nx,ny)
psi_temp3,vort_temp3,psiF3,vortF3,QG_diag3,QG_curlw3,X3,Y3,dx3,dy3=cargar('D:/Facultad/Materias/Circulacion/TP3/out_tmp3/',Lx,Ly,nx,ny)
psi_temp4,vort_temp4,psiF4,vortF4,QG_diag4,QG_curlw4,X4,Y4,dx4,dy4=cargar('D:/Facultad/Materias/Circulacion/TP3/out_tmp4/',Lx,Ly,nx,ny)
psi_temp5,vort_temp5,psiF5,vortF5,QG_diag5,QG_curlw5,X5,Y5,dx5,dy5=cargar('D:/Facultad/Materias/Circulacion/TP3/out_tmp5/',Lx,Ly,nx,ny)
### Escape inercial ###
# 1 (despues pruebo otro)
# R = 0.01, eps = 0.025 da exactamente 4 el cociente.
# simulacion x..(la x es de extremo....)
psi_tempx,vort_tempx,psiFx,vortFx,QG_diagx,QG_curlwx,X_x,Yx,dxx,dyx=cargar("D:/Facultad/Materias/Circulacion/TP3/out_tmp_x/",Lx,Ly,nx,ny)
# 2
# R = 0.005, eps = 0.0175
psi_tempxx,vort_tempxx,psiFxx,vortFxx,QG_diagxx,QG_curlwxx,Xxx,Yxx,dxxx,dyxx=cargar("D:/Facultad/Materias/Circulacion/TP3/out_tmp_xx/",Lx,Ly,nx,ny)
import numpy as np
import matplotlib.pyplot as plt
#%%
dir_salida1 = "C:/Users\Alicia\Desktop\circulacion\k1/" #de esta manera buscamos por donde salen los datos y si o si va la barra al final para que me identificar que debe entrar en la ultima carpeta
dir_salida2 = "C:/Users\Alicia\Desktop\circulacion\k2/"
dir_salida3 = "C:/Users\Alicia\Desktop\circulacion\k3/"
Lx = 4e6
Ly = 2e6
nx = 200 #grillas en x
ny = 100 #grillas en y
from Cargar import cargar #llamamos a la funcion Cargar para cargar los datos
psi_temp1, vort_temp1, psiF1, vortF1, QG_diag1, QG_curlw1, X, Y, dx, dy = cargar(
dir_salida1, Lx, Ly, nx, ny)
from Cargar import cargar
psi_temp2, vort_temp2, psiF2, vortF2, QG_diag2, QG_curlw2, X, Y, dx, dy = cargar(
dir_salida2, Lx, Ly, nx, ny)
from Cargar import cargar
psi_temp3, vort_temp3, psiF3, vortF3, QG_diag3, QG_curlw3, X, Y, dx, dy = cargar(
dir_salida3, Lx, Ly, nx, ny)
#%% Transformo X e Y de m a Km para poder graficar mas prolijo
X1 = X / 1000
Y1 = Y / 1000
#%%dimensionalizamos primero para luego graficar
"""
dimensionalizamos primero para luego graficar
"""
#funcion corriente
beta = 1.9e-11 #definimo beta
tau = 0.3 #tension del viento
import numpy as np
from matplotlib import pyplot as plt
# SI HAY ALGUN PROBLEMA CON LOS VALORES, ES DE ESCALA NO DE UNIDADES.
# usamos metros para calcular beta, usamos TODO en metros.
#ejercicio 1
Lx = 1000000 # dimensiones de la cuenca en metros
Ly = 500000
nx = 100 # puntos de grilla
ny = 50
beta = 2 * 10**-11 # en metros
D = 2000 # metros
from Cargar import cargar
psi_temp1, vort_temp1, psiF1, vortF1, QG_diag1, QG_curlw1, X1, Y1, dx1, dy1 = cargar(
'/home/daniu/Documentos/materias/circulacion_2019/repos_alumnes/CG_mica_luciano/out_tmp1/',
Lx, Ly, nx, ny)
psi_temp2, vort_temp2, psiF2, vortF2, QG_diag2, QG_curlw2, X2, Y2, dx2, dy2 = cargar(
'/home/daniu/Documentos/materias/circulacion_2019/repos_alumnes/CG_mica_luciano/out_tmp2/',
Lx, Ly, nx, ny)
psi_temp3, vort_temp3, psiF3, vortF3, QG_diag3, QG_curlw3, X3, Y3, dx3, dy3 = cargar(
'/home/daniu/Documentos/materias/circulacion_2019/repos_alumnes/CG_mica_luciano/out_tmp3/',
Lx, Ly, nx, ny)
psi_temp4, vort_temp4, psiF4, vortF4, QG_diag4, QG_curlw4, X4, Y4, dx4, dy4 = cargar(
'/home/daniu/Documentos/materias/circulacion_2019/repos_alumnes/CG_mica_luciano/out_tmp4/',
Lx, Ly, nx, ny)
psi_temp5, vort_temp5, psiF5, vortF5, QG_diag5, QG_curlw5, X5, Y5, dx5, dy5 = cargar(
'/home/daniu/Documentos/materias/circulacion_2019/repos_alumnes/CG_mica_luciano/out_tmp5/',
Lx, Ly, nx, ny)
Tau = 0.25
import matplotlib.pyplot as plt
# Directorio madre
os.chdir('/home/kevin/Documentos/MATERIAS_UBA/CG/P/P_02/oceano/salidas_5m')
# Parametros de funcion Cargar
dir_salida1 = '/home/kevin/Documentos/MATERIAS_UBA/CG/P/P_02/oceano/salidas_5m/out_tmp1/'
dir_salida2 = '/home/kevin/Documentos/MATERIAS_UBA/CG/P/P_02/oceano/salidas_5m/out_tmp2/'
dir_salida3 = '/home/kevin/Documentos/MATERIAS_UBA/CG/P/P_02/oceano/salidas_5m/out_tmp3/'
dir_salida4 = '/home/kevin/Documentos/MATERIAS_UBA/CG/P/P_02/oceano/salidas_5m/out_tmp4/'
dir_salida5 = '/home/kevin/Documentos/MATERIAS_UBA/CG/P/P_02/oceano/salidas_5m/out_tmp5/'
Lx = 1000000
Ly = 500000
nx = 100
ny = 50
# Funcion Cargar
from Cargar import cargar
psi_temp1, vort_temp1, psiF1, vortF1, QG_diag1, QG_curlw1, X1, Y1, dx1, dy1 = cargar(
dir_salida1, Lx, Ly, nx, ny)
psi_temp2, vort_temp2, psiF2, vortF2, QG_diag2, QG_curlw2, X2, Y2, dx2, dy2 = cargar(
dir_salida2, Lx, Ly, nx, ny)
psi_temp3, vort_temp3, psiF3, vortF3, QG_diag3, QG_curlw3, X3, Y3, dx3, dy3 = cargar(
dir_salida3, Lx, Ly, nx, ny)
psi_temp4, vort_temp4, psiF4, vortF4, QG_diag4, QG_curlw4, X4, Y4, dx4, dy4 = cargar(
dir_salida4, Lx, Ly, nx, ny)
psi_temp5, vort_temp5, psiF5, vortF5, QG_diag5, QG_curlw5, X5, Y5, dx5, dy5 = cargar(
dir_salida5, Lx, Ly, nx, ny)
#%%
#1
plt.plot(np.array(QG_diag1[:, 3]), label="Ev1=0.025")
plt.plot(np.array(QG_diag2[:, 3]), label="Ev2=0.015")
plt.plot(np.array(QG_diag3[:, 3]), label="Ev3=0.010")
plt.plot(np.array(QG_diag4[:, 3]), label="Ev4=0.005")
plt.plot(np.array(QG_diag5[:, 3]), label="Ev5=0.001")
import numpy as np
from matplotlib import pyplot as plt
# SI HAY ALGUN PROBLEMA CON LOS VALORES, ES DE ESCALA NO DE UNIDADES.
# usamos metros para calcular beta, usamos TODO en metros.
#ejercicio 1
Lx = 1000000 # dimensiones de la cuenca en metros
Ly = 500000
nx = 100 # puntos de grilla
ny = 50
beta = 10**-11 # en metros
D = 2000 # metros
from Cargar import cargar
psi_temp1, vort_temp1, psiF1, vortF1, QG_diag1, QG_curlw1, X1, Y1, dx1, dy1 = cargar(
'C:/Users/Jimena/Desktop/FCEN/circulación/prácticas/TP2/out_tmp1/', Lx, Ly,
nx, ny)
psi_temp2, vort_temp2, psiF2, vortF2, QG_diag2, QG_curlw2, X2, Y2, dx2, dy2 = cargar(
'C:/Users/Jimena/Desktop/FCEN/circulación/prácticas/TP2/out_tmp2/', Lx, Ly,
nx, ny)
psi_temp3, vort_temp3, psiF3, vortF3, QG_diag3, QG_curlw3, X3, Y3, dx3, dy3 = cargar(
'C:/Users/Jimena/Desktop/FCEN/circulación/prácticas/TP2/out_tmp3/', Lx, Ly,
nx, ny)
psi_temp4, vort_temp4, psiF4, vortF4, QG_diag4, QG_curlw4, X4, Y4, dx4, dy4 = cargar(
'C:/Users/Jimena/Desktop/FCEN/circulación/prácticas/TP2/out_tmp4/', Lx, Ly,
nx, ny)
psi_temp5, vort_temp5, psiF5, vortF5, QG_diag5, QG_curlw5, X5, Y5, dx5, dy5 = cargar(
'C:/Users/Jimena/Desktop/FCEN/circulación/prácticas/TP2/out_tmp5/', Lx, Ly,
nx, ny)
Tau = 0.25