CURSO : MECÁNICA DE FLUIDOS

TEMA: MODELO COMPUTACIONAL ELCOM.

(SEGUNDA PARTE)
DOCENTE: AGUIRRE CESPEDES CESAR AUGUSTO
ESTUDIANTES : Bastidas villarreal Andrea katherinne 19170307
Arellano Huertas Jordan Moises 19170103

Salas Galarza André Javier 19170118

Tineo Huamanguilla Edith Milagros 19170309

Velasquez Ayala, Eduardo Alberto 18170339

CICLO :​2020-2
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 2​ de ​28
PARTE 2 Fluidos

1. OBJETIVOS
● Organizar por tercera vez los datos meteorológicos tales como la magnitud
del viento, dirección del viento, temperatura del aire, humedad relativa,
presión atmosférica y radiación solar, para poder sacar nuevas gráficas que
nos ayudarán a diferenciarlas.
● Usar el Matlab para poder pasar los datos que sacamos en el excel y
llevarlos al software para poder convertirlo al wordpad, luego usar el subplot
y el plot para gráficas los datos metereológicos.

2. MARCO TEÓRICO

El preprocesador de ELCOM en el capítulo anterior observamos que se ubica en la

carpeta PRE. Esta contiene una carpeta bathymetry, la aplicación pre_elcom.exe
para correr, y un archivo de configuración run_pre.dat

La configuración para que funcione el preprocesador viene de un archivo

run_pre.dat. Es un archivo de texto que contiene nombres y localización de los
ficheros de entrada y de los archivos de salida del preprocesador.

La salida del preprocesador consiste en dos archivos (binarios) sin formato del
FORTRAN: sparsedata.unf y usedata.unf. El primero contiene la información sobre
espacio tridimensional, y el segundo contiene las relaciones espaciales
unidimensionales requeridas para la simulación. Estos archivos pueden ser
retitulados, pero el usuario debe no perder de vista qué archivos son archivos del
“sparsedata” y qué archivos son archivos del “usedata”.

La carpeta bathymetry presentaba tres archivos: dos archivos de texto

indispensables: bathymetry.dat, bc.dat, y dos opcionales levee.dat, y update.dat.

ELCOM también realiza simulaciones de la hidrodinámica de una laguna cuando

se configura el archivo run_pre.dat que se encuentra en la carpeta Run.
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 3​ de ​28
PARTE 2 Fluidos

En esta carpeta se encuentran además del ejecutable, la carpeta infiles, donde se

ubica los datos meteorológicos y los perfiles de temperatura que sirven como
forzantes externos e internos de la laguna.

También se encuentra el archivo Datablock, que configura la salida del modelo.

Está preparado en bloques. Nos dará resultados de la velocidad de la corriente del
agua, como la estructura térmica de la laguna en tres capas, en la superficie, fondo
y promedio.

Los resultados son almacenados en archivos con extensión ****.unf que no se

pueden abrir, hay que convertir a archivo en formato Netcdf que puede leer Matlab.
De esta forma podremos visualizar los resultados del modelo.

3. PROCEDIMIENTO

Procedimiento:

1-En primera instancia, deberemos utilizar el excel que contenía los datos meteorológicos
(buscar en su ordenador)

2-Estando en el excel, ya que tendremos que pasar su contenido a otro excel que contenga
cierto orden y condiciones horarias, haremos lo siguiente:

-Seleccionar alguna celda dentro de los datos meteorológicos

-Presionar la combinación de teclas CTRL+* y posteriormente CTRL+T, esto le dará formato

de tabla a nuestros datos
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 4​ de ​28
PARTE 2 Fluidos

3-Ahora que nuestra tabla posee formato, podremos filtrar los datos en base a los
parámetros presentados, dejamos abierto este excel para usarlo más tarde.

4-Abrimos un nuevo excel, lo llamaremos ​DatMetElcom​, estando dentro del archivo,

digitamos los datos de tiempo en formato Juliano (Año/Día/Hora), ya que nuestros datos
son a partir del 11 de marzo del 2019 a las 0 horas de dicho día, entonces empezamos con
el siguiente dato: 2019070.00, ahora, en base a la frecuencia (en este caso de 1 hora), le
sumaremos 0.0.42 (tasa a la que aumenta las horas en 1 día de 24 horas). Y sumamos
tantas veces sea necesario hasta llegar la fecha del 19 de marzo de 2019 a las 3:00 pm
(según nuestros Datos Metereológicos en general), es decir, hasta el número 2019078.69
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 5​ de ​28
PARTE 2 Fluidos

Para ello, digitamos en la celda A1 lo siguiente: 2019070.000

En la celda A2 digitamos la siguiente fórmula: =A1+0.042

Ahora, seleccionamos la celda A2, y manteniendo el click izquierdo en la parte inferior derecha de
dicha celda (cuadrito verde), arrastramos hacia abajo hasta que el valor llegue hasta el 2019078.69
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 6​ de ​28
PARTE 2 Fluidos

5-Ahora, regresamos a nuestro excel que contiene los Datos Metereológicos, y filtramos
los datos de tiempo, haciendo click en la casilla que contiene una flecha en dirección hacia
abajo, esta nos desplegará un cuadro con los distintos horarios que se muestran en la
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 7​ de ​28
PARTE 2 Fluidos

tabla, dándonos opción a esconder algunos y mostrar otros, entonces, ya que

trabajaremos (en este caso) a una frecuencia de 1 hora, deseleccionamos todo y
marcamos, a partir de las 12:00 am, los intervalos de 1 hora, es decir 1am, 2am,
3am,...12pm, 1pm, 2pm,..etc. Terminado de seleccionar los horarios damos click en
aceptar
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 8​ de ​28
PARTE 2 Fluidos
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 9​ de ​28
PARTE 2 Fluidos

Nos quedaría algo así:

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 10​ de ​28
PARTE 2 Fluidos

Listo para copiar los datos a nuestro Excel ​DatMetElcom.

6-Ya que cada dato de fecha está asociado a su respectivo valor meteorológico, ahora nos
será fácil trasladar los datos del Excel que contiene los DM a nuestro excel ​DatMetElcom​,
solo hace falta copiar cada columna hacia el otro excel (Nota: Tener en cuenta de qué
parámetro meteorológico se está trabajando, en este caso, el excel con los DM original
seguía el siguiente orden:
RapidezViento/DireccionViento/TemperaturaAire/HumedadRelativa/PresionAtm/Radiacio
nSolar/Precipitacion.

Tenemos que adaptarlo al siguiente formato u orden:

RapidezViento/DireccionViento/TemperaturaAire/HumedadRelativa/RadiacionSolar/Pr
esionAtm/Precipitación.

Nota: Al nombre de la hoja en cuestión, lo llamaremos ​DataElcom​, puede ser cualquier

otro nombre, pero asegurarse de digitarlo bien en el MatLab

7-Terminado el proceso de copiado y ordenamiento de los datos, guardamos nuestro excel

en una carpeta aparte, para su posterior uso en MatLab.
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 11​ de ​28
PARTE 2 Fluidos
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 12​ de ​28
PARTE 2 Fluidos

8-Abrimos MatLab, seleccionamos la carpeta de nuestro excel, y escribimos el siguiente

script:

(Script Personalizado para la frecuencia de 1 hora, pero de aplicación general)

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 13​ de ​28
PARTE 2 Fluidos

9-Hecho esto, nos dirigimos a la carpeta donde se creó nuestro archivo de extensión .dat,
y lo abrimos usando WordPad. Nos quedaría de la siguiente manera, lo cual quiere decir
que el proceso se realizó con éxito.
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 14​ de ​28
PARTE 2 Fluidos

10-Para hacer las gráficas, escribimos en otro script el siguiente código:

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 15​ de ​28
PARTE 2 Fluidos

GRÁFICAS DE LOS DM (Frecuencia 1 hora/Demostración):

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 16​ de ​28
PARTE 2 Fluidos
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 17​ de ​28
PARTE 2 Fluidos

4. RESULTADO
● Subir junto con el informe, el archivo, *****.dat de cada alumno
considerando diferentes frecuencias. El responsable del grupo asigna
una frecuencia a cada integrante, esto permite que cada uno realice
su tarea y aprende a generar el archivo DM. Hacerlo personalizado.

FRECUENCIAS ESTUDIANTES

Cada 60 min SALAS

Cada 45 min BASTIDAS

cada 30 min VELASQUEZ

cada 20 minutos ARELLANO

cada 10 minutos TINEO

● Presentar el gráfico de los DM de cada alumno. Editado.

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 18​ de ​28
PARTE 2 Fluidos

➢ SALAS GALARZA ANDRE

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 19​ de ​28
PARTE 2 Fluidos

➢ BASTIDAS VILLARREAL ANDREA

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 20​ de ​28
PARTE 2 Fluidos

➢ VELASQUEZ AYALA EDUARDO

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 21​ de ​28
PARTE 2 Fluidos

➢ ARELLANO HUERTAS JORDAN

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 22​ de ​28
PARTE 2 Fluidos
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 23​ de ​28
PARTE 2 Fluidos

➢ TINEO HUAMANGUILLA EDITH

 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 24​ de ​28
PARTE 2 Fluidos

5. DISCUSIÓN DE RESULTADO
● Enumerar las dificultades encontradas y la posible explicación.
➢ La principal dificultad que encontré al realizar el trabajo fue poder
ordenar los datos de mi excel; porque cuando los lleve al matlab me
salia como NAN y era porque algunas celdas de mi excel estaban
como textos y lo que hice fue darle click a una llavecita que salia ahi
en la celda y lo cambio a dato número y de ahí normal me corrió
todo el programa.
➢ La principal dificultad que tuve al realizar todo este proceso fue el
ordenamiento y traslado de los datos de los DM originales a los
usados para el ELCOM, para su solución empleé herramienta
FILTROS que ofrece excel, para una mejor gestión de los datos, sin
necesidad de pasar los datos uno por uno a través del MatLab,
mejorando así la eficiencia de mi trabajo.
➢ Mi dificultad más notable fue al momento de ordenar los datos de
excel, ya que me pedían que fuese con una frecuencia de 30
minutos, por lo que organizarlos de forma rápida y eficiente fue una
tarea algo complicada, aunque tras conseguir ordenar los datos,
finalmente se pudo continuar sin problemas el trabajo solicitado.
➢ Mi primera dificultad es hacer el fitro de los minutos que me tocó en
el excel ,y de ahí hacer el ordenamiento según la guía , después en
mi matlab en el momento de correr el archivo me aparece dos
documento uno en formato DAT donde sale solo el título del
documento nada más y el otro documento con formato desconocido
me aparecia todos mis datos del excel , tuve que realizar dos veces
el programa para ver mi error y recien me aparecio el documento
DAT pero ya completo con todos mis datos.
➢ La dificultad que tuve al realizar este informe y/o trabajo fue poder
ordenar los datos en el excel de la manera que el profesor lo hizo,
para luego poder llevarlo al matlab y poder graficarlos en el
wordpad y en el mismo matlab los datos meteorológicos
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 25​ de ​28
PARTE 2 Fluidos

presentados, de ahí al ponerlo de manera ordenada me pudo correr

bien el programa y tratando de que no me salga error.
● Comentar los gráficos de los DM encontrando diferencias al
graficarlas en diferentes frecuencias.
➢ Atendiendo a las gráficas de rapidez de viento en las diferentes
frecuencias (1h, 45 min, 30 min, 20 min, 10 min) podemos observar que la
curva que la compone obtiene mayor agudeza, es decir, existe una mayor
cantidad de datos o que la densidad de datos por periodo ha aumentado,
entonces podemos decir que, a mayor cantidad de datos muestreados, la
gráfica arrojará valores más precisos, ya que podemos identificar de qué
manera se comporta la rapidez del viento en intervalos más cortos.
➢ La gráfica de la dirección del viento en las diferentes frecuencias, a medida
que estas van haciéndose menores, ofrecen una curva que posee altibajos
más abruptos, esto es lógico ya que los datos con los que se trabaja son en
grados (rosa de viento) y estos pueden variar rápidamente en extremos,
es decir, de norte a sur, de este a oeste, de noreste a suroeste, entre otras
combinaciones. Por tanto, a menor frecuencia, la gráfica mostrará crestas
y valles más agudos, con menor espaciado.
➢ En la tercera gráfica, la de temperatura del aire, se aprecia una
periodicidad ligeramente marcada, con algunos cambios en ciertos
intervalos de tiempo, pero en líneas generales, se conserva. Vemos que a
medida que la frecuencia disminuye, las crestas presentan una ligera
variación, es decir, que los cambios más notorios de temperatura a
medida que cambia la frecuencia se dan entre mediodía y las tardes de
cada día.
➢ La variación de frecuencia de los tiempos en la magnitud de la Humedad
Relativa genera que la curva presente una línea gráfica más accidentada,
siguiendo el patrón general visto a frecuencia de 1 hora, pero con mayor
agudeza, osea, mayores picos altos y bajos, aunque esta no difiere mucho,
por lo que podemos decir que la humedad relativa, a menores
frecuencias, se conserva medianamente a lo largo del tiempo, es decir, la
variación es tenue o mínima.
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 26​ de ​28
PARTE 2 Fluidos

➢ En las gráficas de la magnitud de la Radiación Solar, a frecuencia de una

hora, vemos una curva de “suave relieve”, pero, a medida que
disminuimos la frecuencia, notamos un cambio considerable entre los días
14-15, 16-17 y 18-19 de marzo, aquí podemos apreciar que a un menor
intervalo de tiempo, los datos recolectados a esa frecuencia ofrecen
variaciones abruptas pero bastante notorias a mediodía, por lo que
podemos decir que la temperatura, llegando a sus máximos, comprende
cambios importantes de temperatura.
➢ La gráfica de la presión atmosférica no cambia mucho a medida que
disminuimos la frecuencia, siendo similar las gráficas de 1hora con las de
20min, por ejemplo, aunque vale decir que se aprecian ligeras variaciones
de presión a lo largo del tiempo, mas solo eso, pequeñas fluctuaciones de
presión.
➢ La gráfica de la lluvia o precipitación se mantiene constante en todo
momento a un valor de 0, es decir, en estos 8 días, no ha llovido y este un
gráfico dato adicional con respecto a los datos meteorológicos del primer
laboratorio.
6. CUESTIONARIO

1. ¿Cómo cambiar el espacio entre columnas del archivo DM generado?

Podemos corregir el espacio entre columnas mientras digitamos nuestro programa

en MatLab, de modo que todos los datos y los nombres de cada columna queden
ordenados correctamente.

2. ¿Cómo insertar una línea en el archivo DM generado?

Podemos añadir una línea de dos maneras: Añadiendo en el Excel trabajado, o

directamente en Wordpad, yendo a la parte en la que queramos que esté dicha
línea, presionamos la tecla Enter y el texto se moverá una línea, dejando un
espacio para poder añadir la línea deseada.
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 27​ de ​28
PARTE 2 Fluidos

3. ¿Cómo alinear el cabecero de cada columna con los datos debajo?

Tenemos que cuadrar en el matlab los espacio donde queremos que entre nuestra
columna como ejemplo .

4. ¿Cómo escribirías solo una columna de datos?

En los resultados de los datos del matlab podemos obtener los resultados por cada
columna y también en general , abrimos esa nueva pestaña del dato que
deseamos ver y solo sale ese dato , copiamos y lo pasamos al documento DAT , al
igual que también copiarlo en el excel.

7. CONCLUSIONES
● Los datos obtenidos después de la organización nos sirvió para graficar
según el tiempo juliano, como también hay diferencias muy notables con
respecto a los anteriores gráficos que tenía una frecuencia de cinco
minutos, con la ultima grafica que sus frecuencias eran muy variadas
motivo por el cual se diferenciaba en la amplitud de cada onda de la
gráfica, así mismo también sus valores máximos y mínimos varían
considerablemente.
● Usando una frecuencia de media hora, podemos notar, además de que la
cantidad de medidas es el doble que las obtenidas cada hora y por ende,
más exacta, que los datos que manejamos mantienen una correlación unos
con otros, además de mostrar cierto patrón y variación casi predecible, así
como provocar la creación de gráficos más exactos y con más puntos e
información.
● En los datos obtenidos del tiempo cada 10 minutos se puede observar que
en el caso de la magnitud la frecuencia es relativamente variada , en la
dirección del viento al inicio una ligera variación y a la mitad de la semana
se mantiene constante, su temperatura está relativamente aproximada a 25
a 30°C , en la humedad relativa ahí si vemos durante toda la semana un
aumento constante , al observar las demás gráficas de los diferentes
tiempos vemos un ligero cambio en lo que es su magnitud.
● Ya obtenidos los datos del tiempo de cada 20 minutos contando a partir de
las 12:00 am del día de 11 del mes de marzo se puede observar que la
 INFORME N°7
CURSO
Página
MODELO COMPUTACIONAL ELCOM. PRIMERA Mecánica de 28​ de ​28
PARTE 2 Fluidos

magnitud de la frecuencia es variada, varía en la dirección del viento tan

igual como en la temperatura que se mantiene constante en 23 a 25 ºC, en
la humedad relativa aumentó mínimamente con el pasar de los días. en las
demás gráficas solo se tiene ligeramente una diferencia.

8. REFERENCIAS
● Reyes,Carlos, Microcontroladores PIC Programación en basic,Quito-Ecuador,
RISPERGRAF.2006.211p