1

2
3 UNIVERSIDAD NACIONAL DEL ALTIPLANO - PUNO

4 FACULTAD DE INGENIERÍA CIVIL Y ARQUITECTURA

5 ESCUELA PROFESIONAL DE CIENCIAS FÍSICO MATEMÁTICA

9 USO DEL MODELO ROMS PARA LA SIMULACIÓN DE LA

10 HIDRODINÁMICA DEL LAGO TITICACA – PUNO
11

12 BORRADOR DE TESIS

13 PRESENTADO POR:

14 Bach. Grimaldo Apaza Chino

15 PARA OPTAR EL TÍTULO PROFESIONAL DE:

16 0

17
18 PUNO – PERÚ
19
20 201_
21

1 1
2
22
23
24
25 UNIVERSIDAD NACIONAL DEL ALTIPLANO - PUNO
26 FACULTAD DE INGENIERÍA CIVIL Y ARQUITECTURA
27 ESCUELA PROFESIONAL DE CIENCIAS FÍSICO MATEMÁTICAS
28
29 BORRADOR DE TESIS
30 USO DEL MODELO ROMS PARA LA SIMULACIÓN DE LA HIDRODINÁMICA

31 DEL LAGO TITICACA

32 PRESENTADO POR:
33 Bach. Grimaldo Apaza Chino

34 PARA OPTAR EL TITULO PROFESIONAL DE:

35 LICENCIADO EN CIENCIAS FÍSICO MATEMÁTICAS

36APROBADO POR:
37
38PRESIDENTE: ____________________________________
39 (Abrev. de Grado) Presidente de Jurado
40
41PRIMER MIEMBRO: ____________________________________
42 (Abrev. de Grado) Primer miembro de Jurado
43
44SEGUNDO MIEMBRO: ____________________________________
45 (Abrev. de Grado) Segundo miembro de Jurado
46
47DIRECTOR / ASESOR: ____________________________________
48 (Abrev. de Grado) Director/Asesor de Tesis
49
50
51
52Área : Física aplicada
53Tema : Hidrodinámica

3 2
4
54ÍNDICE GENERAL

55
56ÍNDICE DE FIGURAS...................................................................................................3
57ÍNDICE DE TABLAS.....................................................................................................3
58ÍNDICE DE ACRÓNIMOS............................................................................................3
59RESUMEN.......................................................................................................................4
60ABSTRACT.....................................................................................................................4
61I. INTRODUCCIÓN.......................................................................................................4
62II. REVISIÓN DE LITERATURA................................................................................4
63III. MATERIALES Y MÉTODOS................................................................................4
64IV. RESULTADOS Y DISCUSIÓN..............................................................................5
65V. CONCLUSIONES......................................................................................................5
66VI. RECOMENDACIONES...........................................................................................5
67VII. REFERENCIAS......................................................................................................5
68ANEXOS..........................................................................................................................5
69
70

5 3
6
71ÍNDICE DE FIGURAS

72Figura N°1 DESCRIPCIÓN….………………………………………………Pag. N° #

73Figura N°2 : DESCRIPCIÓN ……………………..…………………………Pag. N° #

74Figura N°3 : DESCRIPCIÓN …………..……………………………………Pag. N° #

75

76Observaciones: Cada Figura debe tener un título el cual debe describir la ilustración
77presentada con la menor cantidad de palabras posible el cual debe estar ubicado en la
78parte inferior de la misma. Se debe evitar la redundancia entre el contenido de la Figura
79y su título. Si la Figura tiene símbolos o abreviaciones, éstas deben ser explicadas en el
80respectivo título. Si existen símbolos estos deben anotarse en la parte inferior (pie de
81página). Fotografías, esquemas, gráficos, micrografías, mapas, carta, diagramas deben
82ser consideradas como Figuras y numeradas arábicamente y consecutivamente
83conforme aparezcan en el texto que las cita y deben ubicarse lo más próximo posible del
84mismo. Si las Figuras no son de propiedad del autor debe considerarse la respectiva
85referencia.

86Si las Figuras tienen derechos de autor es necesario el permiso del autor(es).

87El formato de la Figura debe ser de acuerdo a la especialidad del trabajo de

88investigación.

7 4
8
 89ÍNDICE DE TABLAS

90Tabla N°1 : DESCRIPCIÓN ……………………………….…………………Pag. N°

91#

92Tabla N°2 : DESCRIPCIÓN ………………………………….………………Pag. N°

93#

94Tabla N°3 : DESCRIPCIÓN …………………………………………………Pag. N°

95#

96Cada Tabla debe tener un título el cual debe describir la ilustración presentada con la
97menor cantidad de palabras posible el cual debe estar ubicado en la parte superior de la
98misma. Se debe evitar la redundancia entre el contenido de la Tabla y su título. Las
99Tablas deben ser numeradas arábiga y consecutivamente de acuerdo al orden que
100aparezcan en el texto que las cita y deben ubicarse lo más próximo posible del mismo.
101Si existen símbolos estos deben anotarse en la parte inferior (pie de página). Si las
102Tablas no son de propiedad del autor debe considerarse la respectiva referencia.

103El formato de la Tabla debe ser de acuerdo a la especialidad del trabajo de

104investigación.

105Si las Tablas tienen derechos de autor es necesario el permiso del autor(es).

106Verificar si Figuras y Tablas tienen derechos de autor.

107

9 5
10
108ÍNDICE DE ACRÓNIMOS

109Observación: Abreviaciones y acrónimos deben ser especificados únicamente la

110primera vez que se use en el texto entre paréntesis y a continuación del texto que
111describa su significado
112ROMS: Regional Ocean Modeling System
113IRD: Institut de recherche pour le développement
114INRIA: Institut national de recherche en informatique et en automatique
115AGRIF: Adaptative grid refinement in Fortran
116CROCO: Coastal and Regional Ocean Community model
117NetCDF: Network Common Data Form
118FORTRAN: The IBM Mathematical Formula Translating System
119MATLAB: Matrix Laboratory

11 6
12
120RESUMEN

121El trabajo de tesis debe tener una descripción del problema lugar y fecha concisa del
122trabajo considerando la justificación, los objetivos, los principales métodos, los
123resultados más relevantes y la conclusión del trabajo utilizando para ello 350 palabras
124como máximo, incluyendo los conectores. En el resumen se debe evitar el uso de
125acrónimos y abreviaciones que no sean de uso común en la especialidad, así como las
126referencias bibliográficas. El resumen debe ser escrito en un solo párrafo continuado,
127sin el uso de puntos aparte.

128Inmediatamente después del resumen y en un párrafo aparte se deben citar como

129máximo cinco (05) palabras clave (keywords) que describan el trabajo de investigación
130realizado, en el orden de importancia del trabajo.

131Ejemplo:

132Palabras Clave : tesis, formato, investigación.

13 7
14
133ABSTRACT

134Resumen del trabajo en idioma inglés considerando también al final y en párrafo aparte
135las palabras claves (keywords)

136

137Keywords: Thesis, Format, Research.

15 8
16
138

17 9
18
139I. INTRODUCCIÓN

140Sección del trabajo de tesis donde en un máximo de dos páginas se justifica la

141realización de la investigación. Se pueden Utilizar referencias y se debe considerar en
142los párrafos finales el objetivo general y los objetivos específicos de la investigación.

143En este estudio se trabaja con la variante de ROMS llamado Roms_Agrif, el cual es un
144modelo hidrodinámico ampliamente utilizado para este tipo de estudios, particularmente
145para aplicaciones de alta resolución. Enmarcado en el objetivo de este trabajo de
146investigación se apunta a calificar su capacidad para representar la circulación en el
147Lago Titicaca, a niveles espaciales de submesoescala, es decir a resoluciones espaciales
148de xxx kilómetros. Las pruebas realizadas con el modelo de las variables temperatura,
149salinidad y velocidad de corrientes, sugieren que es capaz de representar adecuadamente
150la variabilidad temporal y espacial de las mismas.

151Por su versatilidad, la modelación numérica es una herramienta que suma utilidad tanto
152para diagnósticos como pronósticos de procesos y eventos extremos.

153Es importante conocer los alcances y limitaciones de los modelos, como las
154incertidumbres tanto en casos históricos como en tiempo real. Desde luego un modelo
155oceánico de alta resolución proporciona una mejor representación de las características
156locales y la variabilidad de la temperatura superficial del cuerpo de agua. El
157acoplamiento entre la atmósfera y el océano da una mejor representación de la
158interacción atmósfera-océano gracias a una mejor simulación de los flujos de calor,
159masa y momento.
1601.1. Objetivo General
161Analizar las ecuaciones que describen la hidrodinámica para determinar los patrones de
162circulación del Lago Titicaca verificando la intervención de forzantes atmosféricos
163mediante la implementación y calibración del modelo ROMS.
1641.2. Objetivos Específicos
165Desarrollar la simulación idealizada preliminar para el análisis de las características
166hidrodinámicas del Lago Titicaca, basándose en la batimetría idealizada y real, en el
167estrés superficial uniforme del viento durante un cierto periodo.
168Analizar los movimientos del agua a pequeña escala, en un periodo espacio temporal.
169Evaluar e interpretar la circulación de agua (corrientes), ya sea en periodos temporales
170en el Lago Titicaca, basándose en la batimetría real, el estrés del viento, los flujos de
171calor de superficie, la radiación solar y el flujo de afluencia (salida) de ríos.

172II. REVISIÓN DE LITERATURA

173Sección del trabajo de tesis donde se citan las referencias teóricas y antecedentes que
174sustentan el trabajo de investigación para cada uno de los objetivos propuestos de una
175forma crítica, explicando su importancia y relevancia con el trabajo realizado. Evitar la
176información irrelevante que no contribuya al logro de los objetivos planteados.

19 10
20
1772.1. FUNDAMENTOS BÁSICOS DE LA DINÁMICA OCEÁNICA
1782.2. MODELADO NUMÉRICO
1792.3. DEFINICIONES Y CONCEPTOS IMPORTANTES
180Aguirre, C. A. (2009) en su tesis de investigación “Simulación numérica de la
181hidrodinámica del Lago Titicaca, frente a las costas de la ciudad de Puno en el Perú”
182cuyo objetivo principal fue utilizar el modelo numérico Estuary, Lake and Coastal
183Ocean Model (ELCOM), como una herramienta de análisis para comprender mejor la
184hidrodinámica del Lago Titicaca, logrando realizar simulaciones con viento,
185condiciones iniciales de temperatura y el efecto Coriolis. Concluye que el gradiente de
186temperatura es el principal factor que influye en las corrientes que se presentan en el
187lago. Así mismo, de que por el tamaño del lago el efecto de la fuerza de Coriolis, es un
188factor importante que afecta la dinámica y el efecto del viento, fue notable solo en la
189superficie libre. Se debe considerar que tales resultados fueron productos de realizar
190simulaciones entre el primer y quinto día del mes de Enero del 2005.
191Estudios referidos al patrón de corrientes superficiales fueron realizados por Guzmán,
192E. J. (2006) en su tesis “Modelamiento de corrientes superficiales inducidas por viento
193y marea en la bahía Ferrol - Chimbote” en cuya investigación los resultados muestran
194la influencia del viento y la marea en la generación de corrientes, concluyendo que
195debido a la configuración de la bahía, la circulación que se presenta es muy lenta y que
196la marea ejerce mayor influencia de circulación dentro de la bahía Ferrol, mientras que
197la influencia del viento es notoria fuera de la bahía. No obstante, el modelo está sujeto a
198una revalidación de los resultados, debido a que los valores medidos en el campo
199corresponden a periodos muy cortos de medición.
200En cuanto al uso del modelo ROMS; Escobar, M. (2012) en su tesis de investigación de
201título “Validación del modelo numérico ROMS en la región ecuatorial comprendida
202entre la costa ecuatoriana y las Islas Galápagos” concluye que el modelo tiene la
203capacidad de pronosticar la salinidad del mar con menos error que la temperatura del
204mar, como resultado de usar métricas estadísticas entre los valores pronosticados y las
205mediciones in situ de la temperatura y la salinidad del mar para tres cruceros
206oceanográficos; y advierte que mientras se incrementa la profundidad mayor error en el
207modelo.

208III. MATERIALES Y MÉTODOS

209Sección de la tesis donde se describe con detalle el material experimental utilizado en la

210investigación. Así mismo, los métodos materiales y técnicas empleados para cada uno
211de los objetivos propuestos (evitando repeticiones). En el caso de material de
212laboratorio y equipos se debe indicar entre paréntesis las especificaciones técnicas
213(modelo, marca, número de serie y procedencia de los mismos). En el caso de reactivos
214la marca, lote y fecha de vencimiento de los mismos. De igual manera, se debe señalar
215la metodología experimental empleada y el análisis estadístico utilizado en la
216interpretación de los datos de la investigación.
2173.1. TIPO Y DISEÑO DE LA INVESTIGACIÓN
2183.2. VARIABLES
2193.2.1. VARIABLE INDEPENDIENTE:
220Modelo ROMS (Regional Ocean Modeling System)

21 11
22
2213.2.2. VARIABLE DEPENDIENTE:
222Simulación de la Hidrodinámica del Lago Titicaca
2233.3. METODOLOGÍA
224En la presente sección se describe los procedimientos considerados para el desarrollo de
225la tesis, que incluyen la ubicación geográfica del estudio (sub-sección 3.3.1), la
226descripción del modelo ROMS (sub-sección 3.3.2), la verificación del modelo (sub-
227sección 3.3.3), la adaptación del modelo para el dominio de estudio (sub-sección 3.3.4),
228el forzamiento del modelo para el Lago Titicaca (sub-sección 3.3.5), las características
229de la simulación del Lago (sub-sección 3.3.6).
2303.3.1. UBICACIÓN GEOGRÁFICA DEL ESTUDIO
231El Lago Titicaca se encuentra localizado al sur del país de Perú, geográficamente entre:
23215°13’19’’-16°35’37’’S; 68°33’36’’-70°02’13’’W. Para el presente proyecto se
233considera como dominio de estudio todo el espejo de agua que constituye la cuenca del
234Lago Titicaca.

235
236 Figura 1. Ubicación geográfica del Lago Titicaca delimitado con coordenadas en grados, minutos y
237 segundos.

238El Lago Titicaca está dividido en dos cuencas lacustres separadas por el estrecho de
239Tiquina (800m de ancho). Al norte está el Lago Mayor o Chucuito que tiene una
240superficie de 6 400km2 y una profundidad máxima de alrededor 283m. Al Oeste del
241Lago Mayor, en el sector peruano, está la bahía de Puno, bordeada por las penínsulas de
242Capachica al Norte y Copaya al Sur, en cuyo interior se asienta la ciudad de Puno. Al

23 12
24
243sur está el Lago Menor o Huiñaimarca que tiene una superficie de 2 100km2 y una
244profundidad máxima de alrededor de 40m.
2453.3.2. DESCRIPCIÓN DEL MODELO
246El modelo ROMS (Regional Ocean Modeling System), como bien lo describen sus
247creadores Shchepetkin y McWilliams (2005), es un modelo de circulación oceánica de
248nueva generación que se diseñó para simular de manera precisa sistemas oceánicos
249regionales (Por ejemplo, Blanke et al., 2002; Di Lorenzo et al., 2003;Haidvogel et al.,
2502000; MacCready et al., 2002; Marchesiello et al., 2003; Penven et al., 2001). Salinas y
251Maya (2011) estructuran de manera general los modelos numéricos de este tipo de la
252siguiente manera:
253 Tabla 1
254 Estructuración de los modelos numéricos oceánicos
Modelos numéricos oceánicos
El modelo
Información de numérico Resolución Salida de
entrada propiamente numérica resultados
dicho
Condiciones Ecuaciones: Esquematización: Salidas numéricas
iniciales y de - Continuidad - Redes o Salidas gráficas:
frontera: - Navier-Stokes mallado - Tendencia del
- Batimetría - Movimiento patrón de
- Corrientes - Dispersión Discretización e circulación.
- Mareas integración: - Campo de
- Vientos Parámetros y - Diferencias temperatura.
constantes. finitas - Concentración de
- Elementos sustancias.
Finitos
Calibración (comparación con
Validación (mediante observación y
mediciones, ajuste, modificaciones y
datos históricos)
predicciones)
255 Nota: Adaptado de “SERVICIOS MENORES PARA LA RED DE DESASTRES
256 HIDROMETEOROLÓGICOS Y CLIMÁTICOS: PROYECTO TH1115.6”, por Salinas, J. A.,
257 Maya, M. E., 2011, Instituto Mexicano de Tecnología del Agua, p. 4.

258En referencia a la tabla X, ROMS es un modelo numérico oceánico basado en

259diferencias finitas de superficie libre que resuelve las ecuaciones primitivas con la
260aproximación de Boussinesq con respecto al balance de momento vertical hidrostático.
261El modelo utiliza coordenadas verticales de seguimiento del terreno que permiten un
262estiramiento diferencial, dentro de coordenadas horizontales esféricas, curvilíneas y
263ortogonales. Las ecuaciones se resuelven con pasos de tiempo cortos para la dinámica
264barotrópica y pasos de tiempo mucho mayores para la dinámica baroclínica en un
265esquema explícito dividido. Utiliza tensores rotados para reducir la mezcla diapicnal y
266mejorar el cálculo del gradiente de presión horizontal. que sigue modelo de cierre
267turbulento KPP para la superficie y fondo, obteniendo así un esquema de alta advección.
268Utiliza la paralelización mediante particiones de subdominios bidimensionales, lo que
269proporciona una paralelización compartida y distribuida (Nikolaidis et al., 2019).
270Entre sus principales características tenemos:

25 13
26
271  Ecuaciones primitivas con temperatura potencial, salinidad y una ecuación de
272 estado
273  Aproximación hidrostática y de Boussinesq
274  Coordenadas ortogonales-curvilineas
275  Esquema Arakawa C
276  Cuenca cerrada periódica prescritos en radiación
277  Condiciones de frontera abierta
278  Coordenadas sigma ajustadas al terreno en la vertical y curvilíneas ortogonales
279 en la horizontal.
280  Diversos esquemas de mezcla
281  Superficie libre
282  Solución tri-diagonal con un tratamiento implícito de la viscosidad vertical y
283 difusión.
284Se remite al lector los creadores mencionados en el párrafo anterior la descripción en
285detalle de los algoritmos no lineales y el núcleo computacional; en cuanto a las
286ecuaciones lineal tangencial, adjunto y plataformas están descritas en Moore et al.
287(2004). Incluye algoritmos de física precisos y eficientes, varios modelos acoplados de
288biogeoquímicos (Powell et al., 2006; Fennel et al., 2006), bio-óptica (Bissett et al.,
2892005), de sedimentos (Warner et al., 2008), y aplicaciones del hielo marino (Budgell,
2902005). También están incluidos varios esquemas de la mezcla vertical, las redes de
291múltiples niveles de anidamiento y compuesto.
292Desde su creación ROMS se fue modificando con los aportes de diversos
293investigadores, resultando en diferentes versiones, la versión utilizada en el presente
294trabajo es la implementada por el IRD/INRIA (Institut de recherche pour le
295développement / Institut national de recherche en informatique et en automatique)
296llamada Roms_Agrif v3.1.1 la cual es la última versión estable actualizada con soporte
297hasta el año 2014. Posteriormente cambió el nombre a CROCO (Coastal and Regional
298Ocean Community model) y que se mantiene hasta la actualidad.
299Roms_Agrif es el modelo ROMS con la corrección implementada por AGRIF
300(Adaptative grid refinement in Fortran) el cual es un paquete escrito en Fortran90 para
301la integración de características de refinamiento de malla adaptativa AMR (Adaptative
302mesh refinement) dentro de los códigos de diferencias finitas existentes (Debreu et al.,
3032008); para ver con mas detalle las correcciones realizadas por AGRIF se sugiere
304revisar a Debreu y Patoum (2016).
305El código fuente del modelo (programa principal y subrutinas) está escrito en el
306lenguaje de programación FORTRAN 90/95. Utiliza preprocesamiento C para activar
307las diversas opciones físicas y numéricas. Se puede trabajar en el sistema operativo
308GNU/LINUX, e inclusive WINDOWS (utilizando una máquina virtual), en tanto se
309cuente con librerías NetCDF (Network Common Data Form), que es un conjunto de
310bibliotecas de software y formatos de datos autodescriptivos e independientes de la
311máquina que admiten la creación, el acceso y el intercambio de datos científicos
312orientados a matrices; las cuales son utilizadas extensivamente por ROMS para manejar
313los archivos de entrada y salida. Para el presente trabajo con el objetivo de realizar
314simulaciones numéricas se utilizó la variante UBUNTU MATE 16.04 LTS Xenial
315Xerus de LINUX.

27 14
28
316El manejo de Roms_Agrif se hace con el conjunto de herramientas de nombre
317ROMSTOOLS, dichas herramientas están integradas en scripts (documento que
318contiene instrucciones, escritas en códigos de programación) para MATLAB (Matrix
319Laboratory), desarrolladas para generar la malla del dominio a simular, forzamiento
320superficial, condiciones iniciales, condiciones de borde abierto, simulaciones
321interanuales y climatológicas de las variables del océano (Escobar, 2012). Cabe
322mencionar también que el software MATLAB (versión 2015a) es requisito del modelo
323sobre el cual se programan, ejecutan y visualizan, las entradas y salidas del modelo.
324ROMSTOOLS es un software gratuito, puede redistribuirlo y/o modificarlo según los
325términos de la Licencia pública general GNU publicada por la Free Software
326Foundation; se sugiere consultar la Licencia pública general GNU para obtener más
327detalles. En ese sentido el modelo es de dominio público y tanto como su código fuente
328como ROMSTOOLS se encuentra disponible en la siguiente dirección web:
329https://www.croco-ocean.org/download/roms_agrif-project/
330Hoy en día el modelo es usado por cerca de
331Ecuaciones del modelo
332Estructura del modelo
3333.3.3. VERIFICACIÓN DEL MODELO
334Para la verificación tanto de la instalación como del correcto funcionamiento del
335modelo se editaron los archivos relevantes para las tres etapas: pre-procesamiento,
336procesamiento y post-procesamiento según la configuración indicada en la guía del
337usuario de Roms_Agrif desarrollado por Penven et al. (2010), y que se detalla a
338continuación:
339Pre-procesamiento
340En esta etapa, la cual también es llamada procesamiento de los datos de entrada, se
341ejecutan los archivos necesarios para la siguiente etapa respetando el siguiente orden:
342make_grid.m, make_forcing.m, make_clim.m. El archivo donde se determina la
343configuración de los parámetros necesarios para los archivos del pre-procesamiento de
344un dominio dado es el romstools_param.m, este archivo genera las entradas en NetCDF.
345En ese entender se editó el archivo según se muestra en los cuadros 1 y 2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% 1 - Configuration parameters
% used by make_grid.m (and others..)
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% ROMS title names and directories
%
ROMS_title = 'Prueba Benguela';
ROMS_config = 'Benguela_LR';
%
% Grid dimensions:
%
lonmin = 8; % Minimum longitude [degree east]
lonmax = 22; % Maximum longitude [degree east]
latmin = -38; % Minimum latitudeF [degree north]

29 15
30
 latmax = -26; % Maximum latitude [degree north]
%
346 Cuadro 1. Configuración del nombre y del dominio a simular en coordenadas
347 geográficas (decimales).

348El nombre que se le dio a la simulación fue “Prueba Benguela” que cubre un área mayor
349a la del estudio, delimitado por las coordenadas (38°S–26°S, 8°E–22°E), por lo que se
350evita errores de estabilización del modelo en las fronteras estén incluidos en el dominio
351de interés (ver Cuadro I). Si bien la verificación no se hizo sobre una cuenca, se debe a
352que el modelo no cuenta con datos predeterminados de lagos, está en los objetivos del
353presente trabajo configurar el modelo para simular un escenario como ese.
% Grid resolution [degree]
%
dl = 1/3;
%
% Number of vertical Levels (! should be the same in param.h !)
%
N = 32;
%
% Vertical grid parameters (! should be the same in roms.in !)
%
theta_s = 6.;
theta_b = 0.;
hc = 10.;
vtransform = 1.; % s-coordinate type (1: old- ; 2: new- coordinates)
%
% Minimum depth at the shore [m] (depends on the resolution,
% rule of thumb: dl=1, hmin=300, dl=1/4, hmin=150, ...)
% This affect the filtering since it works on grad(h)/h.
%
hmin = 75;
%
% Maximum depth at the shore [m] (to prevent the generation
% of too big walls along the coast)
%
hmax_coast = 500;
%
% Maximum depth [m] (cut the topography to prevent
% extrapolations below WOA data)
%
hmax = 5000;
%
% Slope parameter (r=grad(h)/h) maximum value for topography
smoothing
%
rtarget = 0.25;
354 Cuadro 2. Configuración de la resolución de la grilla, número de niveles verticales, parámetros verticales
355 de la grilla, profundidad mínima en la costa, profundidad máxima en la costa y el parametro de pendiente.

356La resolución de la longitud en grados de la configuración es de 1/3°36.7km para la

357construcción de las coordenadas horizontales, y de 32 niveles de profundidad para las
358coordenadas verticales sigma con estrechamiento. Los parámetros de estiramiento tanto
359superior como inferior de la coordenada-S vertical fueron theta_s = 6 y theta_b = 0
360respectivamente, es decir con mayor resolución en los niveles cercanos a la superficie y

31 16
32
361sobre los estratos correspondientes a hc, en donde hc es 10. de la profundidad de
362transición entre los niveles superficiales y los de fondo, siempre menor que la
363profundidad mínima hmin; estos parámetros deberán ser los mismos en el archivo
364roms_inter.in. Para las simulaciones se eligió una capa de esponja en las fronteras de
365800m ancho, dicha configuración forma parte del archivo roms_inter.in.
366Si bien las coordenadas de seguimiento de terreno permiten simular una topografía
367realista, es necesario una suavización de las pendientes del terreno para evitar errores
368del gradiente de presión, dicho parámetro rtarget se ha mantenido con la configuración
369predeterminada, teniendo cuidado que no sea mayor a 0.2, ver Cuadro 2
370El modelo se corrió para condiciones de frontera abierta dejando verificar los resultados
371de funcionalidad del modelo permitiendo evaluar las variables trazadoras o pasivas
372(temperatura y salinidad). Como consecuencia de ejecutar el primer archivo
373make_grid.m, se generó la siguiente batimetría del área de estudio (ver figura2), del cual
374los parámetros calculados para los archivos de procesamiento fueron LLm0=41,
375MMm0=42, N=32; que corresponden al número de cuadrículas de la malla en las
376direcciones x, y, z respectivamente.

377
378 Figura 2. Batimetría generada de ejecutar el archivo make_grid.m.

379El segundo archivo que se ejecutó fue el make_forcing.m, el cual crea el archivo que
380contiene los diferentes flujos superficiales y que procesa del directorio COADS05 (atlas
381mundial de datos marinos de superficie), este archivo también es necesario para realizar
382corridas climatológicas, que por defecto el modelo viene configurado con el intervalo de
383tiempo cada 15 días para 345 días, que corresponde a mitad de cada mes; se utiliza este
384tiempo para interpolar linealmente las variables forzadas en el tiempo. Las siguientes
385figuras representan los vectores de estrés de viento, el flujo neto de calor superficial, el
386flujo de agua fresca superficial, temperatura superficial del mar, salinidad superficial del

33 17
34
387mar, sensibilidad al flujo de calor neto superficial, radiación solar de onda corta en
388cuatro periodos diferentes del año.

389
390 Figura 3. Estrés de viento para cuatro periodos del año.
391

392
393 Figura 4. Flujo neto de calor superficial para cuatro periodos del año.

35 18
36
394
395 Figura 5. Flujo de agua fresca superficial para cuatro periodos del año.

396

397
398 Figura 6. Temperatura superficial del mar para cuatro periodos del año.

37 19
38
399
400 Figura 7. Salinidad superficial del mar para cuatro periodos del año.
401

402
403 Figura 8. Sensibilidad al flujo de calor neto superficial para cuatro periodos del año.

39 20
40
404
405 Figura 9. Radiación solar de onda corta para cuatro periodos del año.

406De acuerdo al orden mencionado el tercer archivo ejecutado fue el make_clim.m, el cual
407generará la información necesaria para las condiciones iniciales y limites abiertos; fue
408necesario la verificación en el archivo romstools_param.m que la configuración para los
409parámetros makeini, makeclim y makebry sea igual a uno. (ver Cuadro )
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%
% 4 - Open boundaries and initial conditions parameters
% used by make_clim.m, make_biol.m, make_bry.m
% make_OGCM.m and make_OGCM_frcst.m
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%
% Open boundaries switches (! should be consistent with
cppdefs.h !)
%
obc = [1 1 1 1]; % open boundaries (1=open , [S E N W])
%
% Level of reference for geostrophy calculation
%
zref = -1000;
%
% initial/boundary data options (1 = process)
% (used in make_clim, make_biol, make_bry,
% make_OGCM.m and make_OGCM_frcst.m)
%
makeini = 1; % initial data
makeclim = 1; % climatological data
% (for boundaries and nudging layers)
makebry = 1; % lateral boundary data

41 21
42
 makenpzd = 0; % initial and boundary data
% for NChlPZD and N2ChlPZD2 models
makebioebus= 0; % initial and boundary data for BioEBUS model
makepisces = 0; % initial and boundary data for PISCES model
%
410 Cuadro 3. Condiciones iniciales y límites abiertos configurados.

411Como resultado de muestra se obtienen algunos plots de cuatro diferentes secciones de

412temperatura para el archivo de condiciones iniciales del dominio “Prueba Benguela”,
413ver figura

414
415 Figura 10. Cuatro secciones diferentes de temperatura para el archivo de condiciones iniciales.

416Procesamiento
417Después de generar los archivos NetCDF necesarios (roms_grd.nc, roms_frc.nc,
418roms_ini.nc, roms_clm.nc) con información de la grilla, condiciones de forzamiento,
419condiciones iniciales y de límite abierto, la siguiente etapa es el proceso de compilar el
420modelo; los parámetros que fueron configurados en el anterior archivo deben ser
421editados con las mismas características en los scripts param.h (cambiar la línea con los
422parámetros de la grilla LLm0=41, MMm0=42, N=32) y cppdesfs.h.
423Para compilar el modelo se ejecutó primeramente el archivo UNIX rcsh: jobcom; el cual
424es un programa informático interprete de órdenes. Y luego se ejecutó el archivo roms.in,
425cuyos valores se verificaron sean los mismos que en romstools_param.m en los
426parámetros (THETA_S, THETA_B, HC).

43 22
44
427Como resultado de esta etapa se generaron los archivos roms_his.nc (archivo histórico)
428y el roms_avg.nc (archivo de promedios), los cuales fueron dispuestos a verificación en
429la siguiente etapa.
430Post-procesamiento
431En esta etapa se visualizan las salidas como resultado de la etapa anterior, se utilizó la
432herramienta grafica roms_gui de MATLAB. Esta herramienta permitió visualizar dos
433tipos de variables: las variables que dispone ROMS y las variables derivadas de las
434salidas del modelo. Las figuras a continuación muestran algunas variables seleccionadas
435y de interés de verificación.

436
437 Figura 11. Plot de temperatura en °C (variable que dispone ROMS).

438

439
440 Figura 12. Sección vertical del plot de temperatura.
441

45 23
46
442
443 Figura 13. Campo de velocidades de corrientes oceánicas de corte horizontal en m/s2.

444 Con la verificación de la funcionalidad del modelo se procedió a realizar ensayos

445experimentales de configuración en los diferentes archivos que componen el modelo
446para poder adaptarlo y realizar simulaciones numéricas con el Lago Titicaca.
4473.3.4. ADAPTACIÓN DEL MODELO
448Implementación de la Batimetría en el modelo ROMS
449La batimetría se obtuvo por medio de un
45072541 datos
4513.3.5. FORZAMIENTO DEL LAGO
4523.3.6. CARACTERÍSTICAS DE LA SIMULACION DEL LAGO
4533.4. PROCEDIMIENTOS Y ANÁLISIS ESTADÍSTICO

454IV. RESULTADOS Y DISCUSIÓN

455Los resultados pueden ser presentados en forma de Tablas o Figuras según corresponda,
456(cuidando de no presentar los mismos resultados en ambas modalidades) de la forma
457más clara y concisa posible. Se deben comparar o contrastar los resultados obtenidos
458con aquellos previamente señalados en la sección de Revisión de Literatura, destacando
459la nueva información lograda a partir del trabajo de investigación desarrollado.

460V. CONCLUSIONES

461Se deben presentar tantas conclusiones como objetivos específicos tenga el trabajo de
462investigación, evitando colocar resultados en forma numérica o indicando que las

47 24
48
463diferencias fueron significativas. De preferencia, las conclusiones deben expresarse en
464forma cualitativa. Las conclusiones no deben citar ninguna referencia.

465VI. RECOMENDACIONES

466Señalar el futuro de la investigación

467
468VII. REFERENCIAS

469Las referencias deben estar en el formato de la especialidad del trabajo de investigación

470y se debe tener cuidado que absolutamente todas las referencias citadas en el texto
471aparezcan en esta sección. No se recomienda el uso de comunicaciones personales ni la
472cita de trabajos no publicados en ninguna sección de la tesis. Se acepta la cita de
473trabajos que se encuentren en proceso de impresión aceptados para publicación.

474Recomendación: Citar tesis de su escuela profesional y/o del repositorio institucional.

475(www.repositorio.unap.pe)

476ANEXOS

477Los Anexos serán identificados mediante letras mayúsculas A, B, C y lo demás que

478corresponda. Si hay ecuaciones en los Anexos serán enumeradas como Ec. A.1 ; Ec.
479B.2 o lo que corresponda, igual para Tablas y Figuras (Tabla A.1, Fig. A.1 y los demás.)
480correlativamente.

481

49 25
50