Eddies en Python

(00:06) dar inicio Muchas gracias Daniela bueno Buenos días buenas tardes a todos Bienvenidos a la tercera sesión de este jcat en el día de hoy tendremos eh tres charlas en python si queréis podéis ir entrando ya la plataforma entrando por python en lugar de por r eh Como todos los días si tenéis algún problema nos vais diciendo y las tres charlas de hoy serán de Charles tupá de Laura Navarro y Yay Santana eh la primera será sobre datos de corrientes en python bueno
(00:34) todas en python datos de corrientes la segunda simulaciones de trayectorias de plásticos Después tendremos un breve descanso y por último datos de clima en también en lenguaje python quería recordaros antes de que se me olvide me ha pasado me han pasado de la ecid una el formulario para que os apuntéis como vuestra asistencia de hoy están fijado en el chat principal hakon 2024 eh si no lo localizándose [Música]
(01:30) o si tenéis algún alguna incidencia para entrar a la plataforma o lo que sea nos vais diciendo podéis ir confirmando también que vais entrando y que ya estáis ya sabéis que para todo nos comunicamos por el discord por escrito eh pero siempre que queráis podéis abrir los micros y y la cámara si queréis y y decirnos lo que [Música] queráis bueno no veo comentarios ni ni dudas Así que si queréis vamos a dar comienzo ya a la al primer tutorial que nos dará Charles datos de corrientes en
(02:22) python cuando quieras Charles no sé si hay algún comentario ahora mismo ya conectada Perfecto Bueno pues chars cuando quieras tienes el microgo apagado ahora mismo Sí sí Perfecto Bueno pues buenos días a todos Soy soy Charles y Bueno antes de todo Quería presentarme un poco son Soy ingeniero en física y trabajo en escenografía Desde hace como 10 años si se oye un terremoto es que mis niños están jugando encima de de un poco poco difícil entonces trabajo en osc enog
(03:01) grafía empezar a trabajar primero con el lenguaje matlab que me gustaba muchísimo Antes luego cuando cambié de trabajo en el instituto donde trabajé me dijeron No no aquí no no hay licencias de mlab tiene que otro lenguaje empecé a trabajar en python fue eso en 2000 2013 durante casi 10 años estuve casi solamente con python y luego cambiando través de trabajo me dijo no no python no ti que pasarte al Julia que que
(03:33) también es otro lenguaje bastante potente que se usa mucho en tratamiento procesamiento de datos grafía y de hecho habíamos pensado hace unos meses ha tutoriales de Julia pero era ya bastante complicado con solo python y r y bueno entonces solo indicarles que que sí que Julia es un lenguaje mu muy bueno y que mce la pena echarlo un vistazo cuando cuando tienen tiempo hoy lo que quería hacer es hablar de corrientes corrientes en el mar Pues todo el mundo sabe que que hay muchísimo
(04:08) tiene mucha importancia y vamos a ver cómo se emiten las corrientes dónde bajar datos de de corrientes Y qué hacer con estos datos tengo que compartir mi pantalla porque no la he hecho todavía compartir al entonces Supongo que aquí podéis ver la la pantalla de datos de corrientes s y si si no es el caso pued entonces Gracias en esta figura me gusta el programa nombre muy sencillo aquí que vemos es el campo de viento partir del mod numérico la web está muy chula se ve
(05:03) se ve todo todo muy claro y dentro de eso tenemos opciones Air Ocean chemistry particul Ocean corrientes y y así y así nos da una una visión global de las corrientes pero en superficie y corrientes hay corrientes fuertes Entonces el color rojo naranja la intensidad corrientes fuertes como a nivel del Ecuador Aquí también y otras cosas que vemos los famosos remolinos aquí por ejemplo hay uno en el golfo de México hay hay otros
(05:42) por ahí un montón de de remolinos es algo que me un tema que me interesa mucho en estenografía y por eso quería dedicar un poco tiempo hablar de estos ronos Cómo se detectan porque realmente ustedes todo lo pueden detectar con sus ojos puede decir Bueno hay un remino por ahí un allá un más pequeño pero existen técnicas de detección automáticas de remolinos Porque si cogemos todo el océano remolinos hay cientos de ellos o
(06:07) o más hay muchísimos y nos gustará tener una técnica adecuada para detectarlos y también seguirlos porque un remolino no se queda en su en su sitio se desplaza interactúa con corrientes y transporta propiedades del agua del mar Entonces algo algo muy necesario end Pero antes de todo Sí esta web es una una web que que me gusta mucho por por este aspecto de de visualización no es no es muy reciente creo que tiene más de 10 años porque la presenté Ya varias veces pero para mí sigue siendo una de de las
(06:45) mejores Simplemente porque es sencilla Pero tenemos todo si quiero ver la temperatura del mar Pues aquí la tengo con las corrientes ya está un trabajo muy limpio y Y eso y cerramos la pequeña introducción entonces primera pregunta Que que hago Cómo se van a medir las corrientes la técnica más sencilla no hay que ser grafo para entenderlo sencillo tiro algo al mar y miro como la cosa que tira al mar se va a desplazar
(07:17) puedo ver la la posición si tiene el GPS tiene la posición cada cada minuto cada x tiempo Y a partir de la posición y del tiempo puedo calcular la la corriente la velocidad de la corriente eso se hace hay actualmente en el océano hay muchísimos muchísimas bollas de drivas que llaman en inglés drifting bus y nos dan información muy interesante porque van toda la zona del océano y lo único es que hay que entender que
(07:47) estas bollas tienen que tener una forma específica porque si no van a mucho el efecto del viento y no solo de la corriente Pero bueno no entramos en detalles pero primera idea algo el mar y ver cóm cómo se desplaza luego existen lo que llamamos los radar de alta frecuencia entonces son un sistema hecho por antenas que se colocan cerca de la costa por ejemplo cerca de un faro cerca de un instituto dond luz internet
(08:17) etcétera y estas antenas que van a hacer van a mandar una una onda hacia el océano según la las características del océano entonces oleaje corrientes la señal va a volver a la antena al sensor y se puede con un procesamiento bastante complejo estimar la velocidad de las corrientes pero únicamente en superficie por supuesto estas ondas no no entran dentro de del mar pero por lo menos tenemos con una resolución muy buena
(08:49) una una visión de de cómo está las las corrientes en superficie en una zona bastante amplia hay radares en todos los país No todos los países pero muchísimos países que ti que tienen Costa es una instalación que puede costar bastante pero son grafia es algo muy útil y se usa mucho luego otra cosa interesante es que los satélites nos pueden medir las corrientes primero no directamente pero a través de de la lo que llamamos anomalía del nivel del mar o anomales en
(09:22) inglés y a parte de eso de de de este nivel del mar podemos estimar las corrientes a grandes escalas temporales y espaciales que llamamos las corrientes geostrófico Bueno de eso es que tenemos una cobertura global o casi global y una cobertura temporal una frecuencia temporal bastante interesante tamb y otro instrumento es lo que llamamos el adcp el perfilador acústico de corrientes dopl un aparato que se puede o fondear
(09:56) en el mar o en unos barcos escenográficos tienen este este eso y miden cuando se desplazan miden la velocidad de la corriente en la una cierta capa de de la columna de agua entonces no es solo superficie pero es varias profundidades entonces la dificultad es combinar toda esta información que tenemos voas radares satélites y datos más in que tienen cada uno tiene su ventaja algunos de resolución espacial otros de resolución temporal otros solo en la superficie pues pu eso Entonces qué
(10:33) quería hacer hoy ilustrar unos de estos estos sensores no de estas plataformas o RM que permiten tener acceso a los datos en python va a cargar estos paquetes habituales para cargar fich sf se puede usar el paquete mf4 o x array vaer las figuras voy a emplearlo el cart p que dimos una charla en el taller intermedio me parce sobre sobre este tema entonces ahí se pueden enar todos los detalles aunque no hay no hay muchas
(11:11) cosas complicada un poco más y bueno ejecutar eso qué estoy haciendo aquí pues importando los paquetes voy a crear un objeto que va a ser la línea de Costa y se puede extraer de la base de datos G hs directamente indicando la la escala F significa que es full resolution entonces la resolución completa voy a crear una proyección geográfica para crear mis mapas eso existen muchas muchas proyecciones con el paquete cartopy bueno no voy a encender ahí pero si ponen en Google cartopy projection o
(11:53) proyecciones salen toda la lista y bueno cada uno libre hacer la proyección que que apce empezado por datos de de radar Como le decía radar muchas veces se encuentra cerca de Costa sistema de antenas Y en este caso he cogido datos de del sosp que es una una instituto que se encuentra las Islas vale ellos proporcionan sus datos a través un servid es decir No tengo que bajar el fichero para acceder a los datos pequeño problema si hay un fallo de internet o
(12:29) aquí o en en su ahí no tengos acceso a los datos en este caso Esperamos que funcione voy definir el URL de los datos y que me va a indicar el último campo de radar que han construido con Aquí he creado un pequeño try excep ver si conseguía cargar y si no consigo cargar significa que de deo bien tarda unos segundos en cargar porque hay una una llamada esta web baja los datos Y aquí tengo los datos si quiero tener más detalles Pues voy a ver lo que hay dentro de la estructura que se llama radar
(13:14) Data tard un pequeño rato y vemos que tenos una un fichero que contiene más de 1300 periodos de tiempo 30 longitud 27 30 latitudes 27 longitudes y aquí tenemos los atributos que much entonces que son datos de instrumento que se llam se encuentra guiza Y eso siempre importante los metadatos cuando creamos a m me inter me interesa representar el último campo de velocidad de este radar lo que hay variable y me indica que el
(13:56) último fue hoy a las 13 me parece que es hora utc o sea hace entiendo que hace como hace dos una hora o hace 20 minutos siempre me pierdo con las horas es decir es son datos de hoy Cómo representar datos en este caso voy a crear una figura bastante sencilla Ah vamos a verlo aquí que hago va a crear una figura voy a crear un subplot y lo importante aquí indicar el tipo de proyección que quiero aquí elegido la p que es una proyección que funciona
(14:37) bastante bien para todo que Europa pero se podía cambiar ponerlo aquí y aquí tenemos los tenemos Todos queremos una global podríamos hay otras un po más que se ven menos pero su utilidad por ejemplo un particular Y eso está bien saber que existen y cada una su ventaja voy a representar pu flechitas cada flecha va a representar las corrientes y el tamaño de la flecha indica intensidad corriente hacia el sur en la parte oeste del dominio ahí hacia suroeste y donde vemos
(15:34) punito significa que la corriente es muy muy muy débil como he representado he utilizado la función se llama quer que justamente representa vectores teno que proporcionar un campo longitud unaa longitudes una malla de latitudes Y en este caso las coordenadas de no los dos componentes radar componente o este la norte sur y eso Entonces es un ejemplo muy sencillo para que puedan ver qué tipo de cobertura espacial tiene una
(16:16) pinta tiene un campo de radar el primer índice aquí representa el tiempo entonces las coordenadas latitud longitud y el primero es es el tiempo entonces lo que quiero representar el menos un signica que son los datos los más recientes dentro del fichero lo que había indicado la fecha antes entonces 27 es no va no va a cambiar mucho Pues un poco de cambio Entonces lo interesante es poder representar este este para para
(16:54) varias para varios varios tiempos Luego si no hay pregunta justo ahora Perdona justo justo te iba a interrumpir porque hay un comentario en perfecto Dime el discord no sé si lo puedes ver pero parece que hay no les está yendo bien la línea que empieza con ser dir igual y el Paz El paz que sera el directorio claro No aquí est Okay pues la verdad Que de momento no sé pensaba que era algo común a todos los
(17:53) participantes tener este esta carpeta ah a ver que están poniendo más comentarios nos dice Bueno Emilio y Laura han puesto lo que debe ser a ver que utilicen este o sea sería quitar el re write and shared vamos a darles un minuto o dos para que utilicen este nuevo directorio a ver si va funcionando y nos puedes confirmar rf Vera que te funcionan no sé el nombre real de rc Vera Así que impecable le funciona perfecto Pues nada alguien más Muchas gracias Rodrigo
(18:53) alguien más todo el mundo le ha ido bien Hasta aquí si no nos vais poniendo también por el discord Pues nada Gracias puedes continuar bueno como les comenté muchos sensores capaces de proporcionarnos datos de corrientes pero mucha un inconveniente que tenemos en muchos de ellos es que no nos dan información bajo la superficie por ejemplo el satélite lo que no es que nos dan es corrientes superficiales lo mismo con los drifters y por eso es interesante trabajar con modelos
(19:37) numéricos ya hemos visto esta semana bastante veces lo de copernicus no es que trabajamos por ellos no trabajamos por copernicus pero que al final prof datos en Europa pero también productos globales modelos globales creado por por instituto Mercator en Francia y de hecho creo que son ellos los uno de los primeros que cor modelo oceanográfico global poro tanta fama aquí en esta pantalla explico más o menos como extraer parte de un producto es decir un producto global si Baja un
(20:17) fichero es bastante gordo porque tiene tiene muchas variables tiene una malla muy grande tiene do nivel de profundidad Entonces es interesante seleccionar una zona más pequeña tanto al nivel espacial que a nivel temporal en esta captura pantalla he trabajado con el producto que veis ven el código aquí global análisis Forecast phy 00124 este lo que he hecho he creado una pequeña caja aquí Esos son los datos que me gustaría visualizar pero claramente
(20:55) se puede elegir cualquier zona aquí y el tiempo hab elegido el 24 de noviembre Cuáles son las variables que me gustan dos componentes de la velocidad es hacia Este no sé eso probar vamos abrir sio aquí exp si vamos descripción de este producto características se puede visualizar dentro de maion Pro que es un otro visor de datos no el mismo anes pero Funciona muy bien y aquí tenemos el acceso a los datos aquí quizás algo más complicado porque vemos
(21:45) un montón de productos y no entendemos siempre lo que es Pero por eso nos ayuda la pequeña símbolo de información que nos dice elado de cada uno el grupo etca y nos permite entender esada es cada producto en mi caso quería trabajar con corrientes que puesto es el mod glo bajado he bajado un ficho de corriente a ver si está si no me equivo Ok se supone que que está para todo mundo Entonces es lo que extraí el otro día
(22:40) copernicus lectura del fichero voy a trabajar otra vez con x aunque se puede trajar con paete hdf4 queremos yente lo que hay Dentro de este fichero entonces tengo un paso de tiempo una profundidad 481 latitudes y 840 longitudes Eso es porque he cogido la profundidad más cerca de la superficie pero claramente se puede profundidad más baja si nos interesa es verdad que muchas aplicaciones como la que van a ver con laora después nos interesa la
(23:20) superficie a través de eso vemos todo lo que hay dentro el fichero y sobre todo los dos variables los dos componentes voy a hacer con eso voy a el primer paso primer paso de tiempo y la primera profundidad y voy a calcular la magnitud de de la corriente es dicir la norma del vector de corriente porque necesito eso porque las voy a utilizar donde la donde las las gráficas entonces hac muy rápidamente Hay muchísimas maneras de representar corrientes en el
(23:57) océano y es lo que quiero presentar ahora la primera manera es separar por corrientes voy a representar por un lado la el componente según Norte según oeste este que vemos aquí y aquí norte y sur es interesante porque aquí ya vemos lo que pasa c del Ecuador que tenemos corrientes que van en dirección contraria y aquí pues la corriente el componente Sur Norte No no es tan visible lo que está pasando pero siempre está bien para una prima visualización separ los componentes y ver por ejemplo
(24:39) la amplitud y ver que estamos con corriente bastante intensas del orden de 1 m por segundo en este caso 1,5 m por segundo en la parte más occidental del dominio entonces son corrientes intensas sin embargo con esta visualización la temida de lo que está pasando porque no es difícil en nuestra nuestra mente combinar dos componentes para ver la dirección Pero lo que vamos a hacer ahora vamos a trabajar con vectores que
(25:13) es lo que lo mismo que hemos hecho con la radar lo que indico aquí que la longitud de cada flecha va a depender de la velocidad de la corriente es decir pequeñas flechas o puntito significa que casi no hay corriente y trabajar con y luego flechas largas cor corrientes intensas en la la función quiver que usamos antes hay parámetros que se llaman regrid shape y también este parámetro del ancho de la flecha pues hay que a veces cambiar estos valores para obtener una gráfica que tiene una buena pinta es
(25:53) decir que visualmente podemos visualizarlo los vectores porque si por ejemplo vectores son demasiado si hay demasiados vectores en una zona por la figura sale negra porque cada uno buen todos los vectores están juntos en la figura y lo mismo con la flecha si tengo una flecha demasiado grande pues no se ve nada o pequeña en este caso había probado con 40 para shape y ancho 0,001 y creo que se ve bastante bien la corriente aquí impensar cerca del
(26:27) Ecuador vemos incluso me apce dosos y claro corrientes mucho menos mucho más zon voy a probar con otro valor a ver si la diferencia con 20 por ejemplo sieme ratito en calcular laolu eso es algo que que pueden probar también en su su sesión python pero un esta l el es el objetivo es represent para sir sirve para representar una flecha de referencia significa que quiero representar una fcha de referencia de 1 m por segundo lo que está indicando ahí
(27:31) Y esta flecha la voy a colocar a coordenadas relativas a 0 0,95 es decir la parte derecha figura y 0,5 es decir en medio de la figura lo hice intentó para justamente poner esta flecha donde tenemos el continente porque si no estaría mejo una esquina y podría molestar un poco ya vemos que cambiando a much una densidad de vectores bastante menor y es lo mejor para la visualización también se puede probar con los valores por defecto es decir quitamos los valores no me acuerdo si da
(28:14) un resultado satis problem buo 100 y esperamos y mientras tanto voy a explicar la otra manera que podemos representar es poner un color según la intensidad entonces elegido este este mapa este esta paleta de color que se llama c Ocean Speed no sé si ya conocen c Ocean pero lo voy a poner aquí es un grupo de de paletas de color creado específicamente para la sabemos que siempre que trabajar con temperatura han
(29:03) creado color que se llama termal que se supone pero cada uno tiene su gusto Se supone que el termal está muy bien para presentar comp de temperatura existe el que se llama paraid solar no me acuerdo para qué eseno aquí tenemos una una vista general las paletas que nos ofrecen yo en este caso dije va a representar velocidad pues va a la que se llama Speed a ver que si si me gusta Ah vale eso en la figura anterior con el rap 100 y bueno vemos que es bastante
(29:41) denso pero no está mal vemos muchos detalles me pasa aquí que incluso vemos Islas que que parece que tien efecto sobre el campo de de velocidad Pero bueno eso es tema pero interesante ver que que esta resolución de de datos del modelo copernicus buena resolu espacial entonces aquí creando una figura de nuevo con la función quier Cuál es la diferencia con código anterior Es que le puse coordenadas longitud latitud primer componente segundo
(30:19) componente segundo componente y luego la añadí como argumento la magnitud lo que había calculado antes aquí magnitud que es en realidad la norma del campo de velocidad usando este argumento la función me va a cambiar el color según la intensidad otra manera de representar y bueno no voy a decir que mejor que el anterior puo es diferente porque como tenemos colores menos intensos pues se ven menos las fechas corrientes débiles de nuevo y podemos poner en evidencia corrientes
(30:53) fuertes que haí de un75 m por segundo que sí que son corrientes muy intensas en esta zona Y por último una técnica que también bastante visual se llama streamline líneas de corrientes eh Ya voy a correr este código porque tard un rato en en hacerse Qué está haciendo eso pues realmente integra el campo de velocidad para generar lo que se llama líneas de corrientes e es una visualización como si vemos trayectorias de de de partículas que seguían líneas
(31:28) entonces lo que está haciendo ahora python es calculando estas líneas de corrientes a partir de los datos dos componentes de de de la corriente dos componentes u y v y bueno vamos a ver cómo nos representa este campo tenemos un parámetro importante se llama density y density va por supuesto controlar qué densidad de líneas tenemos Ah en la la figura puse 10 y significa que si aumentamos eso el cálculo cuesta más y aquí es el es el
(32:10) tipo de de figura que tenemos que al final representamos el mismo campo de varias maneras distintas porque la corriente que decíamos intensa aquí la vemos el el Ecuador losmolinos salen salen muy claros Aquí voy a decir es más cuestión de de gusto como quiero representar y si si es a veces se busca más un lado artístico en las figuras otro más algo para publicar y si por ejemplo pongo un valor más bajo el cálculo de la integración para calcular las de corrientes va más
(32:44) rápido se supone pero el resultado también se va se van a impactar porque tenemos menos menos densidad de corrientes durante el tiempo este cálculo si hay preguntas antes de empezar la última parte bu no s vi que en el chat había unem nos puso la enlace para los map de colores Muchas gracias Y por último la pregunta que hacer Bueno hay una pregunta que les voy a hacer pero no la contesto ahora mismo pero por solo hablamos de componentes
(33:30) horizontales de la corriente he hablado de componente norte sur y componente oeste Este pero nunca hablamos de la del componente vertical porque al final dentro del mar sí que existe una corriente vertical y la pregun es porque no no hablamos de eso dejo tiempo para pensarlo y a lo mejor pondré unos elementos de de respuesta dentro del chat qué podemos hacer con estas corrientes esto después de lo que presentó Laura les va a presentar algo sobre plásticos pues una primera
(34:04) aplicación yo la aplicación que quería centrarme sobre sobre era detección seguimiento lo puedes vol un poco el texto que puesto ahí un remolino que en inglés un Ed puede definir como una corriente circular de agua una figura aquí claro lo que pasa es que estos remolinos existen en todos los todos los océanos y tienen propiedades interesantes mucha gente nuf les estudia por satélite conos con datos de biología ver si hay más pl dentro
(34:46) remolino Pues un tópico muy muy interesante hay dos tipos de remolinos buen hay muchos tipos Pero hay dos principales lo que ll anticiclónicos y los ciclónicos simplemente es que giran de una manera de un sentido distinto los anticiclónicos se definen con los remolinos que tienen agua menos densa en el interior del remolino y los ciclónicos contrario agua más densa enor dificultad quizs es que los anticiclónicos
(35:22) yend que corregir aquí entonces el hemisfer hemisferio sur gira los anticiclónicos en sentido y los ciclónicos sentido contrario por no por nos interesan porque nos los remolinos pu ya lo dije al principio en la introducción es que transportar el agua con sus características Entonces si pueden transportar calor puede transportar zona de agua ricas en nutrientes y al final pued afectar a varias escalas espaciales
(35:59) y temporales entonces un rolino es un elemento que está presente en todos los ocanos los mares creado por varios procesos y que que tiene una una vida que puede puede ser más de un año entonces la herramienta Ed tracking consiste en detectar los trinos a partir de datos de anomalías de nivel del mar de corrientes estos datos pueden venir de satélite ah como satélites de altimetría pero también a partir Modelo Si tengo un
(36:31) modelo numérico bastante bueno con buena resolución pues puedo aplicar tracking estos estos datos estas salidas de artículo es de ya 2014 entonces no es un método muy reciente pero aquí ha sido utilizado decenas de veces para estudiar remos en muchas zonas aquí tiene la referencia mams 2014 y Bueno le indiqué aquí datos de satélite y modelos numéricos ahora entramos en la zona conflictiva por qué conflictiva Porque este este paquete de
(37:09) p Tracker Se puede instalar aquí pero no va a funcionar porque hay conflictos en las versiones de python las versiones de las versiones de ma etcétera entonces si pueden correr el código qué va a pasar vamos a ver si funciona bueno los detalles Qué pasa empezamos hemos definido un fichero que se llama Cent altimet que contiene los datos de altimetría y vamos a cargarlo con este Comando H es igual Simplemente porque la malla es regular indico las coordenadas longitud
(37:50) latitud me carga los datos Mu la segunda etapa que no es siempre necesaria es que vamos aar un filtro estos datos ad es la variable que contiene el nivel del mar en detalles absolute dynamic topography no entro en detalles Porque todo lo que es altimetría much muchas variables que que permiten detectar losmolinos y una de ellas la que llama absol topografía y la otra es pero no entro en detalles Aquí vamos a decir okay trabajamos esto
(38:27) se supone que estos valores son los que suelen utilizar en detección de remolinos y explico hasta ahora bien porque aplico este filtro pues para atar componentes de baja frecuencia Todo bien No hay hay ningún problema es ahí donde vamos a de El problema pero en teoría sí etapa final pues detectar los remolinos entonces a c son las salidas de esta identificación a va a contener una lista de Rem anticiclónicos y va contener una lista
(39:02) deinos ciclónicos algoritmo lo llamamos a través de método identification que funciona sobre lo que hemos cargado antes el H que usamos usamos los variables a topogra usamos el componente de velocidad que se llama dentro el fichero se llama u ugos ugs v gos para el otro componente el tiempo dentro para identificación que realmente eso creo que no importa eso es el paso entre dos isolíneas de detección entonces también
(39:47) es un valor que que se suele utilizar Pero podemos a veces cambiarla pero si mejora detección de remolino para este caso necesario otro parámetro Pixel liit es decir queremos que los hermonos que detectan tenga al menos cinco píxeles de de tamaño porque si no a lo mejor va a detectar armonios muy muy pequeños que realmente no existen y 500 también un límite superior Pero eso se puede aumentar y el shape error es un error
(40:13) que va a descartar los remolinos que no tienen una una forma redonda entre comillas es decir si ponemos un Ah un valor muy bajo el algoritmo va aceptar remolinos que tien una forma Ah si está muy deformado que no se parece un círculo Pues sí pues si las corrientes están cerradas aunque no tiene forma de círculo Pues si lo acceptamos como como remolino bu es el contrario valor pequeño es muy estricto y valor alto es un es poco poco estricto y
(40:50) acceptamos mucho muchos el Val 55 también no Val por defecto que que usamos Qué pasa aquí pues no encontramos ninguno he puesto len a significa la longitud del vector a que con losos pues no funciona eso es lo que comenté ayero antes es que hemos intentado instalar este paquete dentro del entorno y no por motivos de versiones de python y de otros paquetes de otros módulos no se puede correr aquí es un problema que
(41:27) hemos dedicado bastante tiempo y a lo mejor Emilio tiene una una solución que podemos implementar pero F fue complicado pero también una una conclusión que a veces comentamos cuando comparamos lenguajes que python a pesar de tener muchos años de de vida es un lenguaje que en unos casos nos puede dar muchos problemas por temas de versiones porque simplemente porque un paquete este de ped tracker que funcionaba perfectamente en mi máquina o en todas las máquinas un año o 2 años simplemente
(41:59) por el hecho de haber cambiado la versión de python y pasar a 3,2 nos da prueba Qué puedo hacer ahora pues voy a compartir mi instalación de de python y de Jupiter es decir el mismo Notebook pero corrido en mi máquina entonces todo lo mismo todo aquí es igual voy a poner De hecho voy a poner all lo único que cambi es que estoy trabajando con un entorno virtual y más seguro que todos los paquetes funcionan y que tengo la versión par 3,10 de python y no la 3
(42:40) 3,12 Entonces vamos a hacer exactamente lo mismo Ah cargando la misma el mismo fichero y lo que lo que nos va a dar aquí como como resultado y aprovecho este momento para preguntar si hay si hay dudas o preguntas porque dentro poco gastado mis 45 minutos no veo que haya preguntas pero Emilio veo que está escribiendo algo igual quiere aprovechar para comentarlo en voz em si quiere hablar es bueno el resumen Es que
(43:20) desgraciadamente llegué a a problemas de nuevo al tratar anoche Casi casi lo instalo en el en la plataforma en el Jupiter Hub Pero al final no pude Eh bueno y solo voy a Añadir que como dijo char realmente eh Comparado a r sobre todo este tipo de problemas surge más frecuentemente en python e pero más específicamente surge con más frecuencia cuando son paquetes especializados como este Que obviamente la comunidad de usuarios de python que utiliza que que y
(43:55) no me refiero a escenógrafo más sino que en el universo de usuarios de python e la comunidad que que se interesa en extraer remolino es bastante pequeñita e y eso e si alguien si los desarrolladores de este paquete eh digamos otras distracciones en sus vida hacen que cambien de tema etcétera hacen que no le dediquen el mismo esfuerzo que le dedicaron cuando lo desarrollaron entonces e no se mantienen al tanto de las versiones no solo de python sino que toda esta secuela de dependencia de
(44:34) bibliotecas que lo usan como Non y otra gran serie de cosas e Esto es lo que ocurrió con este paquete que no no lo han mantenido completamente al tanto y y sobre todo cuando Tratamos de integrar un paquete bien especializado con dependencia muy específica dentro de un entorno que es bastante enorme que es el que tenemos en el en el en el Júpiter Hub en cambio Si uno trata de crear un entorno bien reducido que se concentre nada más en esto e y unos cuantos
(45:05) paquetes generales como Panda sex eh No es muy difícil crear el entorno Charles lo ha logrado yo lo he logrado también el problema es cuando uno trata de integrar de crear entorno bien bien amplio como lo que tenemos en Jupiter H porque es necesario para que cubra todo entonces desgraciadamente eso es una realidad de python pero también una realidad de lo que ocurre probablemente ocurre hasta cierta medida en r también
(45:29) cuando eh un paquete tiene utiliza una gran serie de biblioteca e y no lo mantienen al tanto porque es una comunidad pequeña porque todo el mundo está ocupado y tiene muchas responsabilidades y al mismo tiempo no trata de como hacemos en este evento de crear un entorno universal digamos entonces entre para los que están trabajando en este proyecto Charles y yo segu estos tratando de encontrarle soluciones pero está difícil
(46:01) en la en el Jupiter Hub sí hay que comentar que estamos con ese problema No desde ayer ni desde antes de ayer algo que empezamos a mirar Hace algunos días y los dos probando varias cosas Y eso corriendo este código en mi máquina qué pasa pues en este caso obtengo 80 remolinos anti ciclónicos y cuos ciclónicos porque no que es el mismo número 88 para representar sencillo Porque los que han creado la el paquete ha creado una función display
(46:45) que permite representar el mapa y aquí es el mapa para el dí este día concreto de los remolinos que hemos encontrado como habían dicho antes hay montones de en todas partes tambos diferentes formas distintas alguno poco raro pero son remolinos intensidad distinta y lo que podemos ver Es características tipo el Rayo del del remolino a ver si es por ejemplo radio efectivo aquí estamos en metros Entonces tenemos remos que van hasta 100 km de Rayo otros más pequeños
(47:22) etcétera la amplitud que ser la deviación desde nivel del mar medio Bueno simplific un poco la las definiciones pero la amplitud es algo vertical en este caso y la etapa que propongo dentro del proyecto que que estamos haciendo el proyecto número tres es el seguimiento Porque si aplicamos este hasta detección de remolinos a varios momentos como tenemos esta información cada día si Durante un año aplico la detección de de remol no puedo aplicar luego un seguimiento es decir he
(47:55) detectado remolino en un día el día siguiente parece que se ha desplazado puedes decir rre es el mismo remolino que se ha desplazado en unos kilómetros y el día siguiente también entonces hay un algoritmo que está también publicado para no solo detectar sino seguir remolinos la ventaja es que con eso tenemos acceso a otras variables no solo el tamaño radio tenemos acceso a la zonas donde donde nacen remolinos donde se crean por ejemplo se pueden crear por
(48:27) efecto de unas Islas y donde mueren y y también la Life la vida la duración de la vida de entonces lo que proponemos hacer dentro del proyecto este proyecto número tres sobre remos y aquí acabo la la presentación Muchas gracias Charles vamos unos minutitos si queréis preguntar algo o algún comentario nos podéis comentar también si os ha corrido bien todo el todo el Script salvo Esta última parte Claro que ya hemos visto todos los problemas que hay alguien quiere comentar alguna cosa
(49:13) dudas parece que no hay comentarios bien pues muchas gracias muchas gracias al segundo tutorial con que nos va a hablar de la distribución de plásticos en el océano adelante Laura cuando quieras Hola buenos días se ve tu pantalla Laura se ve bien está configurando yo aquí también para veros vale Ok veis mi pantalla y me ois bien s Vale pues sí gracias Charles por la introducción que ayuda un poco a retomar luego este
(50:15) tutorial entonces bueno como bien decía Charles una de las aplicaciones que podemos tener para usar estos datos de corrientes es hacer simulación de trayectorias podemos decir que estas lo que vamos a simular son plásticos Pero bueno podéis elegir vuestra partícula contaminante favorita y para nuestro caso básico que vamos a ver aquí funciona igualmente entonces Eh sí como decía poder simular la trayectoria de una partícula en el océano es importante por muchas razones
(50:51) y una de ellas es por todos los temas de contaminación Marina por plásticos como ya dijimos pero también por otros contaminantes en el océano como pueden ser derrames de petróleo o también e para el seguimiento de de otras cosas naturales como podría ser el algas sargaso pero que en altas cantidades está dando muchos problemas sobre todo el el golfo de de México y bueno Y también tiene aplicaciones en en temas
(51:20) de Ecología Marina para entender la conectividad de por ejemplo larvas entre entre diferentes Islas o es marinas o sea es una herramienta que se puede usar para para muchas cosas y cada vez se está usando más y bueno existen muchas herramientas para facilitarnos un poco la vida no porque nosotros e vamos a usar estos datos de de velocidades pero necesitamos nuestras partículas virtuales que vamos a usar
(51:48) necesitamos moverlas no con estas velocidades y hay todo un código detrás para poder eh convertir nuestros datos de de corrientes eulerian en un punto de tiempo a una trayectoria y una de estas herramientas de software libre se llama Ocean parcels que vemos aquí y y bueno es bastante sencilla y hay bastantes materiales de de tutoriales y em Y sí Y bueno si es una una herramienta de software libre que tenéis bastantes tutoriales y lo único que necesitamos para poder correr estas simulaciones son
(52:33) los datos de velocidad que ya nos presentó Charles alguno algunos de los que podemos usar y por ejemplo uno de los datos el que vamos a usar en este caso Son productos de satélite en el que usan diferentes datos de de distintos satélites y los combinan para tener un unos mapas de velocidad y bueno si queréis saber más sobre Ocean parcel os he dejado aquí el el enlace de del software y y para ver un pequeño ejemplo de de algo que se puede hacer Lu
(53:09) con estas simulaciones se hice este proyecto que se llama plástica drift más para temas de divulgación a ver si me carga bien okay Y que esto nos muestra aquí un mapa global y podemos seleccionar un punto de este mapa y ver en el curso del tiempo donde puede ser que acabe en este caso nuestro pequeño patito de boma virtual y os muestra aquí al después de diferente tiempo transcurrido donde la las zonas donde es
(53:40) más probable que esté nuestro nuestro pato virtual en rojo las zonas que es más probable donde se encuentre y va disminuyendo los colores y por ejemplo esto lo hicieron para unas campañas de divulgación y todo esto es usando el la herramienta que os voy a mostrar ahora de de oce Passos y eso podemos elegir diferentes zonas y ver cómo se mueve al ir pasando los dos años y vemos cómo se distribuyen eh Por zonas muy extensivas bueno esto Solo para un
(54:16) pequeño ejemplo de cosas que se pueden llegar a hacer con esta herramienta también para temas divulgativos okay Bueno pasando ahora más al código pueso vamos a intentar usar estos campos de corrientes que nos enseñó Charles eh Para hacer una simulación de partículas e vamos a hacer una simulación e muy cortita aquí puse inicialmente 30 días pero veremos después que la haremos eh más corta para que no tarde tanto el código y luego podemos elegir el número
(54:50) de partículas que lanzamos eh os enseñaré dos dos ejemplos e sí eh aquí veis que tenéis dos notebooks Perdón esto no os lo mostré estamos en este del hay el particle set pset a y PS B es solo para mostros dos opciones de cóm podéis decidir lanzar partículas y nosotros estamos aquí en el en el pset a e antes de pasar y a correr todo el código e si alguien tiene algún problema o no ha encontrado el YouTube Notebook o algo que lo diga por el chat o lo podéis
(55:30) abrir los los micrófonos voy a dar un momentito para asegurarnos que estamos todos y todas al mismo punto bueno y así no dudéis en interrumpir cuando estemos presentando si alguna duda si surge algún problema aprovechar que ahora en vivo y en directo para entender todos los los pasos Bueno vamos a ir empezando aquí simplemente cargamos alguno de los paquetes que vamos a usar estos para hacer gráficas y estos los paquetes así
(56:13) también más básicos este de exar Rey que ya también es el que os ha enseñado Charles que estado usando y vale Y luego Aquí ya hemos cargado parcels aquí lo vuelvo a cargar para que que veáis la versión que estamos usando como comentaban Emilio Y Charles estas versiones van cambiando mucho en python entonces a veces es importante ver cuál es la versión que estamos usando ya que puede haber cambios de una versión a otra vale como
(56:43) os decía vamos a usar unos datos de corrientes eh que también fueron descargados de la plataforma copernicus y ya los tenéis en la carpeta esta de shared así que que no debería haber problema y lo que vamos a hacer es primero cargar un un fichero solo de ejemplo para ver cómo son estos datos eh En este caso cambia por cada producto pero tenemos e son datos horarios Ay horarios perdon diarios de corrientes Y tenemos un fichero por cada día entonces aunque luego usemos más ficheros para poder
(57:21) hacer nuestra simulación vamos a visualizar uno okay Lo vamos con y aquí al de Ese es nuestro dataset que hemos creado abierto con xarray podemos ver las características y todas las variables que que tenemos y nosotros vamos a usar como os he dicho necesitamos nuestras variables de velocidad para detectar nuestras partículas y vamos a usar hugos y v que es lo que os mostraba también Charles antes y vale vamos a ver cómo son estos
(57:58) campos de velocidad entonces una de las ventajas de este paquete de exar Rey es que podemos visualizar de manera muy rápida nuestros nuestros datos Y podemos simplemente elegir aquí nuestra variable hugos Y luego vemos que aquí hugos depende del tiempo latitud y longitud Entonces vamos a especificar que solo nos muestra un paso de tiempo y luego hacemos punto plot y ya por defecto sin tener que dar más detalles podemos tener esta visualización rápida de nuestros campos
(58:33) de de velocidad zona recuerdo aquí que tenemos como mostraba Charles tenemos nuestras dos variables de de corriente u de manera zonal oest este y v que es nuestra componente meridional norte sur vale Y bueno como Aquí vemos que es un campo muy muy grande son datos globales vamos a seleccionar un subzona entonces eh vamos a elegir definimos aquí nuestra longitud y latitud mínima y máxima y conexar también nos permite e subselect nuestros datos para reducir
(59:12) nuestro dataset con este temp y vamos a volver a plotear nuestro campo de velocidad ubos y vemos aquí los datos aquí en esta parte del específico y vamos a usar esta zona para para este ejemplo de de simulación entonces Bueno ya hemos visto un poco los datos de velocidad que vamos a usar eh siempre recomiendo sobre todo si son unos datos nuevos que no hemos usado antes Eh pues hacer este check rápido para ver un poco Cómo es el campo
(59:45) y estar seguros que realmente estos son los campos que nos interesan luego cómo decidir qué campo de velocidad es mejor eso se puede hacer todo un tutorial ap parte o una discusión Pero bueno lo primero que tenemos que hacer para usar pasel es definir nuestro fset nuestro campo de velocidad y el primer paso es decir qué ficheros vamos a usar para eh las velocidades entonces en nuestro caso tanto u como V están en el mismo fichero
(1:00:15) así que simplemente lo tenemos que definir dos veces el filed el directorio donde están los datos y el nombre de del fichero Entonces como veis aquí Bueno no por extenderme mucho pero por explicar un poco en vez de poner toda una lista de Todos Los ficheros e que quiero usar he usado esta anotación de aquí eh e Gracias Julia por recordarnos con seción en español ayer en nuestro proyecto se llaman expresiones regulares y nos ayuda a poder decir en plan bueno queremos
(1:00:51) todos los ficheros que del 2018 septiembre que luego seguidamente vayan de cer a uno y luego cualquier cosa entonces una manera de ver de entender un poco que estamos haciendo aquí es usar esta expresión ls de los comandos bash y poner nuestro tenemos que poner nuestro nuestro file nuestro directorio y luego copiado simplemente esto aquí y como veis nos muestra todos los ficheros perdonad que hay en este directo de septiembre de 2018 y que luego vayan
(1:01:30) el siguiente sea de cer a un y luego nos da igual cualquier cualquier dígito cualquier cosa Entonces ya le estamos diciendo Estos son todos los ficheros que vamos a usar para para nuestra simulación y bueno si queréis saber un poco más de esta anotación hay está página web que también os recomendaba Julia Ok Vale y lo siguiente es que tenemos que decir cómo se llaman nuestras no entonces las variables que
(1:01:58) hay que definir son u y v y en este en estos ficheros estas variables se llaman uos y vg y Lu hay que definir Cómo se llaman las dimensiones de estas variables No como ya vimos las velocidades dependen del tiempo de la latitud y la longitud entonces aquí los nombres también son muy sencillos la latitud latitud longitud Time Time y aquí también la importancia de lo que os comentaba antes de visualizar uno de Los
(1:02:27) ficheros que vamos a usar para ver que hemos puesto bien los nombres de de todas las variables que necesitamos okay Y y ya podemos usar este Comando donde simplemente cargamos los datos de velocidad con los nombres de las variables de las eh Y dimensiones y ya tenemos nuestro fset o sea aquí Simplemente ya tenemos listo el campo de velocidades que vamos a usar eem parece que sean muchas cosas y es cierto pero es verdad que son códigos así muy sencillos tanto con este tutorial como con los otros tutoriales
(1:03:10) que podéis usar ver en la página de parcels que están hechos para que bueno no haya que saberse todos estos comandos de memoria sino eh se puede reciclar copiar el código y adaptarlo a los a los datos que que tenemos hasta aquí hay alguna duda o cosa que no que no esté claro por el discord no tenemos nada apuntado si alguien quiere abrir el micrófono y hacer alguna pregunta pero parece que todo claro de momento
(1:03:46) Okay Vale pues que ahora qué tenemos que hacer No pues definir nuestro hemos definido El Fil el campo de velocidades que vamos a usar y ahora particle set que son las partículas virtuales que vamos a lanzar en nuestro en nuestros datos entonces hay como comentaba antes hay varias maneras de de hacer esto varios métodos Aquí voy a mostrar uno en que le decimos a qué longitud y latitud queremos lanzar partículas Entonces yo elegí aquí pues
(1:04:20) una serie una longitud men 90 y lanzar partículas en esta en esta longitud en una línea entre 4 y 6 gros de de latitud norte entonces simplemente le digo pues defino particle set que se van a afectar estas partículas con el campo de velocidad que definimos esto simplemente por el tipo de advección no hace falta que entremos en estos detalles y luego Aquí digo mi en las posiciones de las partículas inicial donde las quiero lanzar
(1:04:57) [Música] os he ido hablando y no está ocurriendo todo un momento e y luego Aquí podemos ver las partículas que que hemos lanzado que hemos definido o sea como os decía no a 90 gr oeste queremos lanzar esta línea de partículas en4 y 6 gr Norte solo para visualizar esto un poco y bueno Y aquí recordaros Que bueno para hacer este plot usamos el parquete cartop que si alguna persona os acordáis de cómo funcionaba esto tenéis aquí el
(1:05:48) enlace al a un tutorial que hizo Charles en la en la edición anterior Entonces ya tenemos casi todas las piezas necesarias para para nuestra simulación lo último que nos falta es definir Eh cuántos días queremos afectar nuestras partículas como os dije en este caso vamos a hacer solo 10 días y luego tenemos el nombre del fichero de salid yo lo he llamado de la gr simulación La Grana el particle set el ejemplo a y bueno un primer test porque Bueno A veces hay problemitas Así
(1:06:23) que empecé por testun y e la última versión de parcels antes guardaba Los ficheros por defecto en formato Ne cdf que vimos que es como tenemos los datos de velocidad y ahora por defecto los guarda en otro formato que se llama punzar Ah no me quiero meter en muchos detalles de esto pero solo dec aquí tenéis alguna información más pero funciona de la se puede cargar con exar rey de la misma manera con un cdf y es un formato de
(1:06:54) datos que cada vez está usando más sobre todo por Cuando hacemos cosas en la nube como como este hacat por ejemplo porque bueno pesa menos y ayuda un poco a lo al acceso a a los datos sobre todo cuando son multidimensionales okay Y ahora esto es la última parte cuando finalmente e hacemos la simulación como veis con Solo dos líneas de código Bueno aquí lo he separado en varias líneas hacemos nuestras simulación de de partículas
(1:07:28) entonces e aquí definimos un poco el el nuestro fichero de salida entonces este es el nombre de nuestro fichero de salida y aquí le digo Cada cuánto quiero e que me diga la posición de las partículas entonces hacemos una estamos haciendo una simulación y yo le digo que me dé la posición de las partículas en este caso cada 24 horas O sea cada día y luego Aquí ya es cuando realmente ejecutamos nuestra simulación en este caso como estamos haciendo una simulación muy simple de
(1:08:04) simplemente es una partícula pasiva no tenemos ningún comportamiento de Por ejemplo si fuera ya queremos hacer una simulación de plásticos más realista de por ejemplo la flotabilidad de de la partícula eem pondríamos otros comportamientos que los llamamos como el kerner no de queé otro componente de velocidad estamos poniendo o por ejemplo si también quisiéramos poner que la partícula se está moviendo por las corrientes y por el viento Pues añadiríamos otra cosa
(1:08:37) pero bueno aquí es el caso más simple entonces simplemente adamos nuestras partículas con nuestros campos de corrientes esto Solo significa que lo adamos con un mro que es el runje cuta 4 Son detalles más técnicos no hace falta que entremos en eso aquí Cuántos días queremos que que adar las partículas los 10 días esto también son temas más de El dt que es otro parámetro que hay que poner aquí por la advección e son
(1:09:10) parámetros más técnicos si alguien tiene dudas podemos entrar en ello y por último el fichero output que queremos entonces con estos con solo cuatro parámetros aquí ya podemos ejecutar la simulación ahora perdon se me había olvidado ejecutar esta línea y Bueno está bien que no sale aquí que va avanzando como veis va bastante rápido y Ya está nos dice out output Fil ya generados el nombre del fichero y veis la extensión pun Star y
(1:09:51) os debería aparecer aquí este fichero a todo el mundo le ha funcionado Alguien ha tenido algún problema e yo por Bueno aquí en este caso en el Cloud podemos ver directamente que aparece la carpeta yo si no lo que a veces hago en el Jupiter notb usar estos comandos B solo añadiendo aquí el signo de exclamación y simplemente me muestra Cuáles son los últimos ficheros que se han generado en nuestro directorio y
(1:10:29) aquí me dice que a las 3:15 hora local Se nos generó nuestro output file alguna duda algún problema técnico parece que nada no hay de comentarios okay Bueno est sí yo tengo una duda Pero tiene que ver más con el con el respaldo de los archivos Cómo yo puedo sacar o respaldar o copiar el archivo de salida sar a mi directorio del computador por ejemplo Bueno desde el desde el Cloud una manera sencilla si haces clic derecho te da la opción puser ajusto no
(1:11:23) da la opción de Download ex Bueno aquí para este punto para este Notebook si me da la opción de Download para esto no s si al ser un directorio ya no me ves ya igual te pid Ah no deja Descargar directorios Si fichero por por fichero puedo Añadir un comentario permite Descargar directorio o carpeta pero tal vez una manera y no sé si quieres tratar de demostrarlo es comprimir el directorio no sé si tenés que ir a la terminal o si se puede hacer desde ahí s
(1:12:07) De seguro hay que ir a la terminal e no sé si vale la pena hacerlo ahorita pues tal vez lo podemos mandar las instrucciones por el discord Sí después si no tiene mucho que con Sí luego no sé a ver cuánto ocupa exactamente este chentito ya aquí también al poner el H es para que nos diga esto el tamaño del fichero en en ya en tamaño que los humanos podamos Ender no no sol las máquinas y son solo Sí 627 kilb así que es bastante ligero Entonces igual sí que es algo si yo por ejemplo pensé que era más si
(1:12:56) no los ni los subí al github pero es una otra opción es via el git para ficheros ligeros también os lo podríais Descargar Okay alguna otra duda Bueno pues si no la última parte que os quería Mostrar es luego como podemos un poco eh echar un vistazo a estos datos que hemos generado pues podemos usar otra vez xarray para cargar los datos se usa igual que lo hemos hecho antes con un ncdf Simplemente hay que Añadir esta parte de decir que usamos que son datos t para
(1:13:48) cargarlos Y ya y el formato de estos datos un poco particular no como son datos lagrangianos que tamban en el tiempo em tenemos e 20 trayectorias o sea 20 partículas esto a veces siempre es un poco confuso son 20 trayectorias que o sea 20 partículas y observaciones refiere a a los 10 pasos de tiempo no Entonces tenemos 20 partículas Y tenemos su trayectoria durante 10 días y qué variables tenemos Pues tenemos la la posición a la que se
(1:14:24) encuentran la longitud luego la la profundidad que en nuestro caso es constante porque solo hemos usado las corrientes de superficie para albitar las partículas y el tiempo entonces Bueno podemos visualizar ahora nuestras nuestras trayectorias con este plot donde aquí ploteo la posición inicial de las partículas y aquí la la trayectoria vale se va un poco pequeño como veis aquí tengos todos los puntos de la de las partículas iniciales y aquí vemos Las trayectorias que siguen
(1:15:07) durante estos 10 días y bueno para entender para entender un poco por cómo Por qué se han movido así porque han seguido estas trayectorias podemos repetir esta gráfica pero con el campo de de velocidad de de fondo eh Aquí voy a poner el campo de velocidad de la velocidad zonal de fondo ya es el mismo plot que antes Entonces sí Aquí hacemos graficamos nuestro campo de velocidad zonal y aquí tenemos la la la posición
(1:15:48) de las partículas final en rojo y la posición de las partículas eh inicial en negro entonces podemos ver no el campo rojo significa que las partícul es positivo Entonces está yendo hacia el este Entonces está todas estas partículas que empezaron aquí nos aparecen a la derecha y en cambio aquí nuestro campo de velocidad es azul es negativo están yendo hacia el oeste y todas estas partículas de aquí nos
(1:16:16) aparecen más hacia hacia lo esto y si queréis también podemos hacer el mismo plot eh con nuestro campo de velocidad meridional Entonces copiamos esto de aquí y simplemente podemos cambiar en vez de la variable ugos nuestra velocidad zonal podemos poner vgs nuestra velocidad meridional y aquí tenemos lo lo equivalente entonces vemos to toda esta zona es es azul y todo significa que nuestras partículas se han desplazado hacia el sur y eso nos sale
(1:16:55) que en verdad si vemos con respecto a la posición inicial que las partículas en negro todas aparecen más hacia abajo pero bueno aquí de fondo solo mostramos un campo de velocidad en el día el primer día o sea el día el día 6 de de septiembre pero si cargan podemos cargar luego otro de de nuestros datos de otro día y ver en relación a ese campo que tenemos o otra cosa que podéis probar hacer es de calcular la la velocidad absoluta que os mostraba ch la
(1:17:35) Norma y tener ese campo de de fondo también laora tienes una pregunta en el chat que te leo se puede poner una gran cantidad de partículas en un solo punto estoy pensando en un derrame de petróleo Gracias de Harold sí se puede hacer eso lo único clo que hay que tener en cuenta si realmente te interesa hacer eso O no no Claro si tú pones todas las partículas en el mismo punto se van a si no añades algun otro tipo de de parámetro esas
(1:18:08) partículas igual por ruido numérico se separan pero pero si no no O sea lo que solemos hacer hemos hecho también pruebas con derrames de petróleo si por ejemplo en este punto rojo que vemos aquí ha habido un derrame en vez de solo poner partículas ahí ponerla poner muchas pero al en ese punto y alrededor No para que así tengas también un poco Esa esa variabilidad pero sí se pueden yo ahora lo he hecho con muy
(1:18:38) pocas partículas Pero podemos poner muchas más el límite en el número de partículas va a depender de la de la potencia del ordenador o o o la máquina que estés usando para tu simulación pero general es bastante potente sobre todo con con lo que comentaba ahora de de que Los ficheros sean punto zar en vez de punto NC eh alguna otra duda si no esta pregunta creo que hace como enlazar bien a otra vale Okay está res el problema de pero bueno
(1:19:24) Gracias Harold esto esta pregunta hace enlazar bien al otro ejemplo que os he puesto aquí en el Hub que es con el particle set B que es todo lo mismo la única diferencia es de cómo vamos a lanzar nuestras nuestras partículas no como os dije antes ahí pusimos como partículas en en unas coordenadas que que nos interesaban Pero hay otra manera en la que podemos hacer esto Vale ahora voy a volver a ver todo
(1:19:55) eh si vais corriendo también este Notebook a ver si os surge alguna otra eh duda cualquier cosa no no hay preguntas cuentas usamos los mismos campos de velocidad todo igual hasta aquí en este punto no antes pues otra vez antes lanzamos las partículas en en unas coordenadas y ahora lo que vamos a hacer en vez de en una longitud y latitud precisa vamos a decir pues en en en cada punto de malla queremos lanzar una partícula no y vamos a hacer como una malla de partículas y vamos a
(1:20:42) decir en cada punto de donde tengamos datos eh de nuestro producto de velocidad vamos a lanzar una partícula entonces aquí tenemos el fils en nuestro campo de de velocidades y estos valores aquí son valores de de los índices pues una subselect de datos voy a mostraros esto a ver si queda un poco más claro entonces aquí tenemos esta malla de partículas donde en cada punto donde tenemos una partícula tenemos un punto de datos de velocidad porque eso podría ser que
(1:21:20) lancemos e varias partículas en el mismo punto de malla de nuestros datos pero en verdad todas esas partículas van a estar moviéndose con la misma velocidad entonces simplemente aquí definimos Tenemos aquí un vector de longitud un vector de latitud y con este Comando creamos la malla Y tenemos ya aquí definido nuestro nuestro nuevo set de partículas que nos parece así no sé si queda muy claro este punto voy a intentar enseñaros un ejemplo con los
(1:22:01) datos aquí por ejemplo si usamos la el campo de velocidad zonal si usamos estos índices para ahora plotear nuestro campo de velocidad aquí son todos los datos que tenemos de nuestra de nuestra subzona Y esto es si usamos los los índices que hemos puesto para lanzar las partículas entonces aquí veos de longitud y aquí latitud vale Sí otra vez esto es nuestro campo inicial perd que estamos viendo la
(1:23:11) subregión pero aquí estos índices corresponden a a los datos globales entono Estos son todos los datos que tenemos le decimos vale de nuestro campo de velocidad que es este de aquí elegimos estas posiciones en longitud y en latitud que si las visualizamos ahora con con esta gráfica de aquí tenemos las posiciones en latitud las posiciones en longitud y estos son los datos Entonces esta gráfica de aquí lo que les ha en
(1:23:45) cada punto aquí de nuestros datos en cada cuadradito vamos a lanzar una partícula y entonces estas son nuestras partículas iniciales más o menos queda claro Hay alguna duda Okay pues bueno el resto del código es todo igual vamos a inyectar estas partículas durante 10 días definimos el nombre del fichero de salida aquí nos tenemos que acordar de de cambiarlo para que no se sobreescribi el otro fichero que hemos
(1:24:27) generado y como veis antes teníamos 10 partículas ahora tenemos unas cuantas más pero aún así va un poco más lento pero aún así el código es bastante rápido y ya tenemos generado aquí nuestro nuestro nuevo fichero punto punto zar y otra vez como os comenté podemos usar este Comando y vemos Cuál ha sido el último fichero lo que se ha generado en este caso se nos ha generado correctamente nuestro fichero de de
(1:25:03) salida y por último otra vez vamos a visualizar un poco los los datos los cargamos con xarray vamos a ver pues mira aquí si antes teníamos perdonad no 10 pero 20 partículas Ahora aquí tenemos 100 cinco veces más y nada ha tardado igual dos o TR segundos más vamos a visualizarlos las trayectorias aquí veis otra vez las posiciones inerciales con los puntitos negros y las trayectorias que sigue cada partícula y otra vez vamos a ver con el campo de velocidad zonal de fondo si
(1:25:49) entendemos Cómo se mueven y como veis aquí Tenemos también una velocidad zonal muy fuerte hacia el este entonces todas nuestras partículas iniciales aquí en negro se desplazan hacia hacia el este y tenemos estas posiciones finales en en rojo Okay y bueno esto era lo que os quería Mostrar para hoy es una simulación muy básica eh Como dije lo llamamos simulaciones de plásticos pero en este punto podría ser de eh cualquier cosa que esté en la la superficie del mar ya que
(1:26:25) no hemos añadido ninguna otra cosa que ningun otra comportamiento a nuestras partículas aparte del efecto de las de las corrientes pero desde esto sería como la parte inicial de cómo crear esta simulación Y desde ahí las podemos ir haciendo más y más complejas e hay alguna duda o también puede ser una curiosidad de cosas que se podrían llegar a hacer como como comentaba el compañero para el caso de de rames de
(1:27:03) petróleo no veo ninguna pregunta por el chat si hay alguien que quiera abrir micrófono y comentar alguna cosa también si queréis parece que hay una pregunta iba por justamente lo lo último que mencionas no O sea hay diferentes aplicaciones para todo este campo de de todo este campo entre ellos está lo de contaminación y en mi caso a m me interesa lo que es huevos huevos y larvas de peces no cómo seleccionar si esto es bien lo básico o
(1:27:47) sea cuál sería el siguiente paso para seleccionar las características importantes para para las partículas O sea que en este caso sea son huevos no Sí ahí sería ver yo eso Sé que hay hay compañeras que han estado trabajando más en simulaciones de huevos y y larvas no es mi caso pero sería ver qué qué afecta el movimiento de de de los huevos o larvas particulares que que quieres ver tú son son por ejemplo huevos que se
(1:28:22) mantienen la superficie de del océano o son Realmente algo que estén luego en la en la columna de agua si por ejemplo si se desplazan en la en la columna de agua habría que Añadir aquí aparte de nuestras velocidad eh u y v habría que Añadir en W no desplazamiento vertical y entonces usar otro otro producto de datos que nos de estas corrientes en en tres dimensiones y y y por ejemplo si la temperatura si el campo de temperatura es algo que afecta
(1:29:01) a las larvas y a su desplazamiento también es algo que se puede Añadir o sea sería un poco más ver e una manera en la que a veces definimos esta herramienta de parel como si fueran unos bloques de legos no Entonces tenemos todos estas piezas que podemos ir juntando ahora tenemos esta pieza base y luego se haría Yo sí recomendaría un poco identificar Cuáles son las otras piezas que que te faltan eso para definir el el movimiento de de los
(1:29:34) huevos o larbas que quieres estudiar e si con un compañero Comenta me preguntaba también el otro día sobre Por qué no definimos el tamaño de las partículas o Cuánto pesa Y eso Y es que primero antes de definir las características de la de la partícula hay que definir el impacto no de esto si ahora aquí decimos que la partícula mide 5 cm no nos va a cambiar nada porque no hemos e implementado ese impacto No si por ejemplo decimos bueno las partículas
(1:30:11) que miden 5 cm me lo invento pues van más rápido más lentas con la corriente e Pues habría que Añadir e esta parte de de código como otra pieza de de legos Aquí más o menos no sé si queda muy claro Sí gracias [Música] nada Gracias Laura hay alguna otra pregunta pregunta o curiosidad como queráis bueno mientras si como me quedan igual 5 minutos Marian Eh bueno IB vamos a hacer justo ahora un descansito pero si tienes algo rápido que quieras enseñar Ah sí sí pores si le surgí alguna duda si no eso
(1:31:08) que tanto a mí como char también nos pueden poner si se les ocurre una duda más tarde como siempre en el Canal de discord eso solo quería Mostrar por último aquí en este enlace que os puse al principio os lleva la página de parcels y como comentaba tenéis aquí mucho material de tutoriales y documentación donde hay muchos ejemplos e esto es un ejemplo de cómo empezar e Cómo usar diferentes datos de de modelos eh Cómo tener en cuenta la la profundidad hay
(1:31:43) muchos ejemplos diferentes que podéis explorar y ver qué se puede o no se puede hacer con esta herramienta y ya está estás muteada Perdón Así decía que no sé si has copiado este enlace en el discord O querías o lo vas a poner y pueden ir echando un vistac Yo muy bien pues ahora aprovechamos a hacer un descansito en mi ordenador pone y 38 si queréis paramos 10 minutitos vale Mirad la hora que tenéis en vuestro
(1:32:23) vuestro ordenador porque cambia en cada sitio así que dentro de 10 minutitos nos vemos otra vez hasta ahora muchas gracias muchas gracias mariam eh Bueno Buenos días y buenas tardes a todos primero pedirles un poquito de disculpas porque no estoy pudiendo asistir a a las sesiones en grupo me ha surgido un viaje y la verdad que mi mi agenda está un poco un poco complicada Pero pero bueno Perdonen espero al menos
(1:32:56) estar el el viernes para las presentaciones de del trabajo de los grupos y nada Me gustaría aprovechar dos minutitos para presentarme porque entiendo que la mayoría no me conocen me llamo Jay Santana Falcón eh soy escenógrafo y hasta hace poquito estaba trabajando Me supongo que será mi afiliación donde la que les aparece en mete France pero hace nada me mudé a a las palmas aquí en España Así que estoy trabajando ahí en la universidad ahora
(1:33:21) mismo y nada dedico a el tema del cambio climático a pesar de que soy anrao pues estudio un poco el cambio climático Y utilizando modelos de sistema tierra o System models esm que es lo que vamos a ver hoy un poco para contarles qué son modelos esm son modelos que intentan simular todos aquellos procesos y todas aquellas interacciones entre los diferentes sistemas de la tierra entre las atmósfera la superficie continental junto con su río su química su física su
(1:34:00) su biota eh También la interfaz entre la atmósfera y el océano y en el caso de que nos ataña ahora mismo el océano tanto la parte física la química como también toda la biología el c y eh Y de nuevo todas las interacciones entre ellos particularmente yo estoy estudiando desde hace unos 3 cu años cómo eh cambios producidos en el océano pueden afectar al a los organismos Por qué con estas herramientas porque estas
(1:34:34) herramientas como ya hemos visto son complejas como ya les digo Son complejas y nos permiten ver el pasado o el estado del pasado el presente y también simular o proyectar diferentes futuros para para la Tierra así que en mi caso lo que hago es ver cómo esos cambios en el futuro podrían afectar a los organismos y ya pasando al a lo que vamos a ver hoy eh primero repasaremos un poco lo que ya vimos con con Julia el lunes que sería la descarga de los datos de de estos
(1:35:08) modelos eh concretamente de los que están dentro del proyecto semip 6 los proyectos de semip son proyectos de intercomparación vale los mip los model intercomparison project y se porque es la sexta versión de este de estos proyectos No es el único mip que existe este sería el cip c viene de cle porque son modelos acoplados entre los diferentes sistemas de la de la tierra pero también tendríamos otros otros mips tendríamos el highest mip que se encarga
(1:35:40) de modelos de alta resolución o el escenario mip que también lo veremos que son eh proyectos de intercomparación en modelos simulando eh diferentes narrativas del futuro diferentes escenarios y y otros otro otros semip eh antes de de comenzar si tienen alguna duda si tienen alguna pedirme alguna aclaración o o simplemente algún comentario eh Por favor interrumpan en cualquier cualquier momento estaré intentaré estar atento al al chat
(1:36:15) y y si no pues directamente abran el micrófono y y me y me interrumpen vale Bueno pues comenzamos con como lo que digo con la descarga primero de estos de estos modelos eh Como ya le les diría Julia lo primero que vamos a hacer es cargar las librerías básicas que yo siempre cargo para este tipo de de Pues de datos que son las de xarray porque vamos a manejar datos de tres o más dimensiones eh pandas nampi que son
(1:36:50) las dos básicas y m lip en este caso para ver eh las figuras de una manera así rápida y básica este Api es una aplicación eh disponible eh públicamente No no es una aplicación que yo haya hecho está en pangeo que es una una página en la que además de bueno de muchísimas otras cosas de servir como repositorio tienen también varios eh varias herramientas varios tutoriales y de uno de ellos Pues he conseguido esta Api que funciona hay varias apis por ahí
(1:37:25) circulando en internet y la verdad que esta eh Por lo menos a mí me ha funcionado esperemos que ahora en vivo y en directo Pues también funcione lo que hace esta Api pues directamente Busca en los servidores de isgf que es los servidores en los que todos los eh participantes en estos proyectos cuelgan los resultados de sus modelos Aquí tengo abierto la página del isgf como vemos aquí les dije que tenemos
(1:37:56) diferentes diferentes mips vamos al caso que nos daña cip y podríamos aquí Buscar directamente por el por el nombre del del modelo por ejemplo vamos al modelo británico que tiene bastante resultado y como Ven a medida que nosotros seleccionamos el el modelo el mip al que pertenece eh o el experimento tenemos por ejemplo la histórica que ahora Les explico van eh inmediatamente cargándose los resultados que tenemos disponibles para esos para esas simulaciones pues lo que hace
(1:38:40) el el Api disponible en pangeo como les digo y lo hemos copiado aquí es buscar en esos en esos servidores en los servidores de isgf eh los resultados que nosotros definimos Cómo definimos Pues por actividad como hemos visto en la en la página la variable que nosotros queremos el experimento al cual pertenece la simulación que nosotros estamos buscando etcétera y también aquí vemos el el el miembro el miembro no es más que el tipo de configuración que
(1:39:19) eh eh a la que pertenece esa simulación normalmente cada modelo eh lanza una serie de eh simulaciones dentro de un ensemble lo que significa es que se inician de manera diferente y así podemos tener de una una misma simulación varios resultados entonces la r la R1 es la simulación 1 y a partir de ahí pues se van lanzando diferentes simulaciones la i sería el inicio qué inicio se ha eh tenido en cuenta eh Y luego las otras dos serían la física si está incluida y el forzamiento qué forzamiento se está utilizando un
(1:40:01) forzamiento de CO2 y un forzamiento de otro tipo entonces con todos ellas con todas estas eh filtros que vamos aplicando a nuestra búsqueda Pues el Api va a la página que les comenté y nos da los resultados que obtiene de este filtro aquí y como nos ha dado bastante resultado Simplemente nos centramos en uno de ellos y lo abrimos para ver el fichero que nos ha descargado de la de la página de isgf este sería el fichero es un fichero que lo he abierto en en
(1:40:41) xarray que está para toda la para todo el globo Aquí vemos que en latitud y longitud tenemos 180 gr y 360 gr Porque en estos ulciones vamos a tener normalmente resoluciones en torno a a un grado de la de latitud y longitud y tenemos 1980 datos eh de tiempo y estos datos temporales están en meses de hecho aquí en el Api hemos aclarado que lo que queríamos era unos ficheros mensuales pertenecientes esta o mayúscula al océano
(1:41:24) podemos ver qué pinta tiene la variable que en este caso es la temperatura del océano y simplemente pues pintar un poco para ver qué pinta tiene nuestra variable seleccionando el el primer año el primer perdón el primer registro que en este caso era 1960 creo que estoy equivocado creo que es50 Pero bueno nos vamos a 1960 enero de 1960 y vemos tiene eh Perdón nuestro nuestro fichero del mismo modo podemos descargarnos otros ficheros en este caso
(1:42:08) lo que lo que indica aquí es el fichero del área para la zona del para para el océano ya que solamente vamos a tener en cuenta como vemos aquí el grid oceánico mientras que son todos vamos a tener Val o ceros en el en el Continente pues este fichero lo que nos va a a permitir es calcular como luego veremos eh las medias eh para el océano global Y de nuevo Aunque está hecho ploteos el el área en este caso en el caso del modelo
(1:42:42) que hemos seleccionado como vemos la los grid no tienen exactamente el mismo el mismo área sino que eh aquí en las zonas eh subtropicales van a tener una un área un poco diferente porque no es un grid eh regular sino que dependiendo de la latitud va a variar eh su su resolución es bastante típico en lo en los modelos de este de este de cip y vamos a tener que que tenerlo en cuenta he ido un poquito rápido porque
(1:43:17) esto lo vimos el el lunes pero si hay alguna duda hasta aquí se ha ido muy rápido todo bien vale de momento no hay comentarios por escrito Vale pues entonces podemos seguir una vez tenemos estos ficheros descargados lo que quería mostrarles es algunos cálculos muy básicos que se pueden hacer con con estos y que Y explicarles un poquito las diferentes simulaciones que tenemos en en semip Entonces si podemos abrir la el fichero de El Notebook de
(1:44:08) visualizar he incluido un poco de una un resumen unas explicaciones muy muy sencillas de lo que sería el tanto el cip como el el como les dije el el entee una comparativa de todos aquellos modelos o de todos aquellos sí modelos de todos aquellas instituciones que desarrollan modelos acoplados del sistema tierra la sexta versión comenzó en 2014 y hoy en día supone el grueso de los modelos climáticos usados supone el grueso también de los estudios de los que se
(1:44:49) basan los informes del ipcc los informes del clima ya que son modelos que muchas instituciones eh desarrollan pueden ser eh comparados entre ellos ya que ya que el el semip el proyecto se encarga de lanzar una serie de de estandarizaciones de protocolos que tienen que seguir todas las instituciones y gracias a ellos pueden ser comparados y pueden llevarse a cabo estudios con una con bastante robustez digámoslo así eh estadística
(1:45:24) así que de nuevo para ya sumergirnos en el en el uso de estos modelos vamos a cargar estas librerías básicas que que les comenté antes xarray pandas nampi y la para hacer figuras básicas y lo que vamos a hacer Es irnos a la al semip y también utilizar el escenario mip por un lado utilizaremos la la simulación histórica que es aquella a la que o la que se usa para estudiar el pasado reciente del del del del sistema tierra
(1:46:04) se considera que comienza a partir de El de la era industrial y eso se considera estandarizado que el año 1850 comienza el momento en el que el hombre o la actividad humana eh comienza a alterar el clima terrestre y ya que el semip 6 se estableció en 2014 se ha establecido que las simulaciones del historical llegan hasta 2014 y como Pongo aquí pues se caracteriza porque siguen las concentraciones históricas de gases de
(1:46:40) efecto invernadero de los ghg en inglés a partir de ahí a partir del año 2014 comenzaría lo que denominamos escenarios que son una narrativa son inventados una narrativa de la evolución futura del planeta se parte del año 2014 No termina en diciembre del 2014 Y a partir de 2015 empiezan estas estas narrativas obviamente estas narrativas tienen unos consensos y están basadas en datos no reales pero en en datos
(1:47:16) posibles y se van a establecer en diferentes eh subgrupos en diferentes sub escenarios hay escenarios en los que se presupone que la humanidad va a tomar un camino basado en la sostenibilidad perdona parece que hay también una pregunta con el Paz como ha habido antes in igual a si puedes no sé si puedes ver el el discord el Ah en el discord perdón bueno es el lo que has llamado inpath equal y pone users Santana Falcón downloads si está correcto supongo que
(1:47:50) no no eh dónde está ese impaz entiendo que aquí igual si quieres decir Dónde está localizado Esto sí sí sí está en el archivo visualizar dos pero la presentación que tiene no es el mismo link por eso de acuerdo sí puede ser que me haya que haya subido la la versión uno en lugar de esta versión no nos pu copiar el enlace en el discord puedes copiarlo tú sí sí sí sí sí sí por supuesto Ah bueno pero lo mismo que antes en vez de
(1:48:38) share rewrite creo que va a ser share sin y ya está es exactamente ahí Verdad sh que está así correcto sí sí exact pues comprobar Rodrigo que te va bien Ahora si pensaba que había subido en los tutoriales esta versión seguramente subí la primera versión que tenía lo cuando trabajaba local bueno Supongo que así ya le va funcionando si no nos comentas por disc Rodrigo vale si en cualquier momento hay otra eh otro otro fallo del del Notebook
(1:49:27) díganmelo y lo y lo y lo corrijo subo de nuevo esta versión de todas maneras cuando termine ahora de de dar la explicación Eh bueno Pues seguimos con los escenarios lo que les decía hay diferentes subgrupos dentro de los escenarios en el que se supone que la humanidad va a seguir un camino de s sostenibilidad nos vamos a quedar a medio camino no el middle of the road que es un poco cerca de la evolución que estamos en la que los cambios no son muy
(1:49:59) significativos con respecto a lo al histórico un mundo de regionalismo un mundo donde donde se favorece la defensa en detrimento de de otr de otras cuestiones un mundo de disparidades entre entre países y un mundo en el que nos desarrollamos eh basándonos principalmente en combustible fósiles y por tanto Pues un mundo un poco más pesimista no un mundo en el que contaminamos eh bastante nosotros vamos a centrar en el en en dos escenarios uno en el que seguimos a la senda de que que
(1:50:38) está marcada por el por los datos históricos y un escenario pesimista el 585 si ven eh artículos en el que se utilicen System model verán seguramente o les será era Familiar El 585 ya que hasta hace muy poco era lo que venía a llamarse el Cómo era el business as usual no el el todo sigue igual el 585 realmente es un es un escenario bastante pesimista pero en los estudios se suele se suele utilizar ahora veremos realmente Por qué se utiliza bastante este este
(1:51:15) escenario Igualmente vamos a cargar el el archivo are celo el archivo que nos va a dar el la área de del océano para luego utilizarlo y lo primero que hacemos Es calcular las eh series temporales globales lo que vamos a hacer es por un lado tanto los datos históricos como los de los escenarios la temperatura del de superficie y vamos a multiplicarla por el área vamos a sumar latitud y longitud y dividirla por el área total de esta manera calcularemos
(1:51:47) una serie temporal de la media ponderada ponderada eh Para ponderar todos los el valor de todos los de todos los puntos de todo el grid eh Ya que si el grid fuera eh regular pues seguramente las diferenci no sería muchas pero como vimos antes el gren no es regular Así que vamos a tener que tener en cuenta la el área de todas y cada uno de los píxeles de cada uno de los de los píxeles del grid y nada vamos
(1:52:20) a observar simplemente el la media durante el periodo histórico que empezaba en 1850 como les dije y llegamos hasta finales del 2014 vamos a hacer un plot bastante sencillo y además vamos a a incorporarle eh un Rolling min el Rolling min eh lo vimos en los talleres intermedios no estoy seguro si todos ustedes estuvieron allí o no pero eh lo que básicamente hace es suavizar la la serie temporal realmente aquí la vemos es vamos a a
(1:52:51) darle como ven tenemos todas las en azul datos mensuales con toda su variabilidad Cuando hacemos un Rolling me lo que hacemos Es suavizar es decir que cada dato es la media de los eh de los de 12 meses de 6 meses hacia atrás 6 meses hacia delante y con eso conseguimos un suavizado muy parecido si hubiésemos hecho un valores eh anuales Así que esta es la evolución de la eh temperatura de oceánica superficial desde que empezó la acción
(1:53:26) del hombre como vemos más o menos en los años se mantiene más o menos constante hasta los años 1980 1970 80 en el que ya comienza a verse que hay una subida bastante grande de la temperatura Un aumento de la lo que no el Ocean warming el aumento de la temperatura del océano y vamos a a ver a cuantificar este este Ocean warming o por lo menos a a a visualizarlo normalmente los niveles de calentamiento global tanto en el océano como en el en el planeta tierra
(1:54:05) pues se hacen en relación a la temperatura global que había antes de la actividad antropogénica de la actividad del hombre es ese momento como vimos antes por conceso por consenso se considera que es anterior al año 1850 lo que se suele hacer en los estudio es mirar la temperatura de una simulación como conocida como pi control Perdón este de aquí este pi control realmente lo que hace es intentar estabilizar el el los modelos Entonces
(1:54:42) se se introduce el forzamiento que en este caso son los flujos de de carbono atmosférico y se deja eh Por protocolo unos 500 años de pi control a partir de esa simulación de ese pi control van iniciándose eso que les decía antes del ensemble esos diferentes eh simulaciones esas simulaciones históricas entonces lo que se hace normalmente es se cogen en una cantidad grande de o o el pi control entero los 500 años o se coge una
(1:55:16) selección grande de años dentro de pi control se hace la media y a partir de ahí pues se calcula el calentamiento global como nosotros no tenemos capacidad de de preferí no cargar un un pi control en el en el Júpiter Hub vamos a aproximar a la temperatura media del año 1850 es una buena aproximación Así que simplemente calculamos la temperatura del océano para el año 1850 hacemos la media de todo lo de todo el año y la plotea y vemos Cuál es el el
(1:55:54) estado del del océano Sería más o menos este para el el momento en el que comenzó la actividad humana hacemos Igual que antes calculamos la media ponderada en este caso es solamente un momento un instante una media del año 1850 y lo que tenemos es un valor para el océano global el océano global de media estaba Este modelo a unos 18 com 26 18,3 gr Cent de media y ahora lo que vamos a hacer es construir una figura en la que vamos a
(1:56:38) eh calcular lo el calentamiento global a partir del de este valor cargamos las librerías en este caso voy a utilizar la librería cbor y bueno grid Spec Es simplemente para para para seleccionar en lugar de hacer un plot hacer varios subplot pero se puede hacer Igualmente con con plot me gusta sorn para este tipo de de figuras para estas series temporales que es lo que vamos a hacer aquí porque nos permite pues darle un poco de de estilo
(1:57:11) darle un poco de de una visualización bonita al a nuestra figura pero se puede hacer Igualmente con o simplemente usando las herramientas de de xarray definimos la figura y definimos una serie de de parámetros del estilo del de la figura que vamos a a plotear ponemos su título y demás y aquí calculamos eh con esta con esta parte de aquí calculamos el calentamiento global simplemente restamos a nuestra serie
(1:57:47) temporal en este caso histórica el valor que acabamos de calcular el valor preindustrial y lo mismo vamos a hacer para el la simulación 245 y para la simulación 585 Además de eso les he dejado aquí no entiendo que no hace falta que lo que lo explique simplemente es una un un código que me parece interesante que es bueno de utilizar de vez en cuando o para algunas de estas figuras porque de nuevo mejora la visualización de la figura que es indicar dentro de la figura
(1:58:27) que qué estamos ploteando en este caso vamos a plotear la serie histórica el 245 y el 585 Así que lo indicamos tarda un poquito aquí simplemente al final he establecido los límites de la de la figura tanto de los valores de calentamiento global como de los años en los que quiero que se represente y esta sería la eh figura Aquí vemos que nuestra figura parte del cero ya que hemos eliminado el valor preindustrial y como vimos antes la serie histórica que sería el azul empieza a calentarse el planeta
(1:59:07) alrededor del de los años 80 y si seguimos un escenario 245 que como vimos alineado con la serie histórica eh nos encontraríamos que a final de siglo en el año 2100 estaríamos En unos 2 com 2 gros Cent por encima del valor preindustrial y no sé si ustedes están familiarizados pero seguramente lo habrán oído que el acuerdo de París que fue reafirmado hace pocos años en en Glasgow establecía que tenemos que intentar que el planeta no se calentaste
(1:59:42) por encima de 1,5 gr cados con lo cual incluso siguiendo un escenario cercano al al actual estaríamos muy mu por encima de ese por lo menos en el océano muy por encima de esos eh de esos límites si como les dije siguemos el escenario más pesimista el que han llamado business as usual estaríamos hablando que el océano se llegaría a calentar hasta 4 gr en relación a lo que estaba eh antes de que la acción del ser
(2:00:16) humano comenzase en 1850 alguna duda hasta aquí vale efectivamente Emilio aquí lo pone exactamente que que sí que gran parte de los estudios del último semip están o están basados en en resultados de o en estudios realizados con resultados de estos de estos modelos históricos ahora lo que podríamos hacer es ya que sabemos en la serie temporal del océano global Cuál será el de media la el aumento de temperatura a final de siglo en relación al al inicial podríamos verlo en en dos dimensiones
(2:01:15) Así que lo que vamos a hacer es calcular las anomalías desde final de siglo al momento preindustrial cargo de nuevo librería sem es para los colores y cartop lo hemos visto en otras en otros tutoriales que es para plotear mapas Y en este caso esta variable que les puse Aquí es simplemente un fichero para pintar la superficie continental en el mapa que vamos a crear las anomalías normalmente necesitan igual que las climatología rangos temporales de 30 años o más para que sean representativos de un periodo
(2:01:58) Así que lo que vamos a hacer es quitarle o sustraerle al final de siglo para cada uno de los escenarios eh 30 años de final de siglo menos los primeros 30 años de la sección histórica y así calculamos nuestras anomalías tanto para la el escenario 585 como para el escenario 245 Espero que no tarde mucho y mientras voy explicando que aquí de nuevo vamos a definir la figura elegimos una proyección en este caso Robinson que es
(2:02:31) una proyección bastante usada en en estos en los estudios que utilizan cip seguramente que si ven algunos eh artículos verán que esta esta proyección est es bastante usada y simplemente pues en este caso eh utilizando un Loop entre nuestros dos anomalías eh vamos a pintar los dos escenarios aquí lo que tendríamos es la anomalía la diferencia entre el mundo siguiendo un escenario eh o siguiendo nuestra actividad tal
(2:03:15) cual la tenemos hoy en día y el mundo antes de que comenzase nuestra actividad y por otro lado un mundo más pesimista en relación a Cuando comenzó nuestra actividad aquí lo ven como ven La anomalía básicamente es positiva Siempre casi todo rojo aunque como vemos el rojo es más intenso cuando la eh contaminación por parte o el uso de los combustibles fósiles por parte del de la humanidad es mucho mayor pero si se dan
(2:03:44) cuenta los patrones son muy similares es decir que independientemente de nuestra acción el océan va a calentarse de la misma manera la única lo que más va a variar es la intensidad de ese de ese calentamiento también vemos una característica muy pues muy característica valga la redundancia que es esta parte de aquí entre Islandia y Groenlandia que es lo que se llama el el el warming H que es una zona en la
(2:04:17) que existe una anomalía y que en lugar de calentarse se está eh enfriando especialmente si seguimos un escenario de de eh de cierta digamos de cierta mitigación de de una contaminación no muy extrema como el 585 es una eh característica que aún está en estudio que tiene que ver con la con la amoc con la con la circulación meridional pero que todavía su o hasta donde ellos sus orígenes no están del todo eh
(2:04:56) determinadas Y de la misma manera podríamos ver algo quizás más interesante que sería cómo será el mundo a final de siglo en relación al mundo eh de hoy en día aquí les dejo la manera de de de calcularlo Aunque ahora en las veces que lo he intentado hacer con Jupiter H no no me no me ha funcionado porque eh concatenar que es lo que teníamos que hacer dos pues dos variables en este caso dos dos dos array la histórica y el
(2:05:30) escenario Pues lleva mucho tiempo es verdad que se podría directamente seleccionar los años Sería mucho más sencillo pero bueno lo vamos a aproximar de nuevo a cómo sería el mundo en relación al final de la serie histórica al final de la semip miro por si hay alguna alguna duda vale veo aquí que han compartido el el problema que vamos a tener con con con el calentamiento solo una una cosa No no he leído este este
(2:06:07) esta esta noticia pero comentarles que a pesar de que el año 2024 sea un año caluroso y por encima del 1,5 gr centí Eso no significa que no podamos cumplir con lo establecido en el en el en en el consenso de París con eso de 1,5 gr a final de siglo Porque eem primero estas estimaciones se basan en en periodos de 30 años como hemos visto y eh también porque eh podemos dejar de contaminar o existen una serie de tecnologías para intentar eh paliar parte del calentamiento que
(2:06:46) son los cdr y bueno no no no me lío por este lado pero vamos que que no significa que hay que ser optimistas y no significa que que este año estemos por encima del del 1,5 que no podamos conseguirlo Pero bueno vamos a nuestro a nuestras anomalías vamos a calcularla eh Para no irnos muy lejos de los datos actuales vamos a calcularla con 15 años como visto no sería lo lo lo mejor lo mejor es hacerlo eh con 30 años o más
(2:07:17) pero bueno es una pequeña aproximación que tomamos una pequeña licencia que nos tomamos para este cálculo calculamos nuestras anomalías Y de nuevo eh pintamos la figura Exactamente igual que lo teníamos antes y lo que veríamos ya para eh finalizar Ya les dejo en paz cómo será el mundo en relación a cómo era eh a principios de de este siglo como vemos si seguimos eh contaminando a ritmo actual el mundo será más cálido o el
(2:07:47) océano Perdón será más cálido pero si eh nuestros gobierno y nosotros mismos no tomamos conciencia y eh seguimos contaminando a un ritmo aún mayor utilizando únicamente combustibles fósiles y sin poner ningún eh remedio nuestro mundo es un mundo Muy cálido un mundo más peligroso y un mundo en el que seguramente no querremos vivir y nada más por mi parte si tienen alguna duda algún comentario pues estoy encantado de de intentar
(2:08:32) resolver muy bien muchas gracias Jay alguien tiene alguna pregunta comentario todos habéis sido capaz todos habéis sido capaces de correr bien el código verdad aparte de algún problemilla con con los directorios parece que todo ha ido bien muy bien pues muchas gracias Geral parece que no hay preguntas para ti nada Ustedes aquí terminamos la sesión hoy de tutoriales quería antes de que paséis ya a la parte de eh de los proyectos quería comentaros bueno a ver también que
(2:09:10) participen Los tutores de proyectos que a ver qué pensáis la jornada de mañana en principio la tenemos dedicada a trabajo colaborativo también vuestro Pero tenemos la opción de de ponernos en bueno de de hacer una pequeña sesión común de media hora quizás o lo que haga falta pero para poner en común Pues eso dudas que hayáis ido teniendo eh A veces también eso lo de trabajar en grupos pequeños va bien para trabajar así más
(2:09:35) ágilmente pero a veces va bien poner también en común pues esto los problemas que habéis tenido cada uno de los grupos y tal si son cosas muy cortitas que quizás lleve poco tiempo lo podemos comentar ahora o si veis que es interesante que nos reunamos mañana pues hacemos así me comentáis Los tutores como lo veis y el resto de organizadores si habéis bueno comentarnos por ejemplo empezamos con el grupo de pesca creo que estás tú Marina
(2:10:06) No sé si está Héctor ya por aquí Marina si nos quieres comentar problemas que hayáis tenido si hay algo que os vendría bien lo tienes apagado el micrófono ahora mismo Hola buenas s nosotros nosotros ayer le dedicamos bastante tiempo a a ver cómo trabajar en local y luego subirlo a github y al final nos dimos cuenta de que no iba que no iba a ser posible que todos trabajás a la vez y que se actualizas en github que iba a ser mejor que cada uno trabajase Bueno
(2:10:39) estábamos acabando de decidir Cómo iban a hacerlo cada uno empezó a trabajar en su r lo que que hicimos fue con zoom que yo lo tengo ilimitado crear salas de grupo para que pudiesen agruparse en distintos grupos de trabajo porque nosotros somos siete y pensab que era mucha gente y y luego retomar en la sala central para compartir lo que se había hecho y y y bueno y ahora leyendo los comentarios que teníamos en la sala de profesorado Eh sí que era buena idea crear un
(2:11:14) dropbox o que que ellos creasen un dropbox o un Google Drive para compartir documentación o que nosotros viéramos a al sher de giub estábamos y y empezaron a ver ya decidieron que iban a trabajar con datos de Pacífico de tunidos eh de zona eh del ciat creo que se llamaba que es parecido al icat pero para la zona de latinoamérica inter América y y bueno y estaban empezando a trabajar con con variables ambientales revisando eh tutoriales de los tutoriales intermedios
(2:11:53) porque no ninguno había participado en los tutoriales intermedios Entonces les hemos pasado todos los vídeos y estaban tratando de ponerse al día y con la sesión de Jer de hoy pues seguro que les ha sido super útil también y el resto claro muy bien pues pasamos al segundo grupo contaminación Julia o Laura o ambas si queréis comentar có lo ves Hola Bueno no sé si Laura está la por aquí el micro apagado Laura Pero bueno hablas
(2:12:24) la que queráis vale Pues habíamos decidido que trabajaran en grupo igual una persona compartiendo y que fueran comentando a la vez por los temas de subir al repositorio no como ha dicho Marina que no sab muy bieno organizarlo y estaban escogiendo la zona para para hacer la simulación de de las trayectorias de plásticos estamos un poco en ese punto y a ver H dónde podemos avanzar avanzar hoy no sé si tienen otras dudas que quieran comentar luego
(2:12:58) con los otros equipos o no sé si alguien de nuestro grupo quiere decir algo en torno a eso o o es decisión nuestra esto marián Bueno un poco de todo Si hay gente que le parezca bien pues eso tener una pequeña sesión mañana o lo quiere comentar ahora es un poco la opinión de todos Laura tú querías Añadir algo más igual Eh bueno por lo que que comentan por el canal que también fue creo que un problema que nosr F el año pasado todo el tema de identificación en github con los tokens Y eso eh que parece que igual
(2:13:35) lo solucionaron Pero bueno no sé si lo quiero indicar por aquí o por el canal del jaat Si siguen teniendo si es un error común no hay un error un problema común esto de de al hacer Push al al github este problema con el password o el token y toda esta identificación que según el Software que usas puedes ya meterlo en en un apartado para que no te lo pregunte cada vez Pero supongo que cambiará un poco de unos a otros Sí yo creo que es más la configuración inicial luego ya si se
(2:14:12) pueden hacer cosas para que no se tea que poner cada vez pero creo que había problemas a la hora de de conseguir este este password toen Laura o alguien me podrían recordar a través del canal de profesores creo que leí algo al respecto de cómo facilitar el uso de get dentro de la plataforma que Leí algo en la la compañía que maneja este sistema para que sea más fácil que se congeló Laura Hola ahora Sí sí dices que si Bueno creo que pasaron por aquí por el canal alguna de
(2:14:57) las posibles soluciones lo reenviamos por el principal y si no si yo tenía pendiente Buscar que el año pasado creo que también nos pasó algunos grupos y hicimos algo a ver si recuerdo que den dio instrucciones más detalladas vale la pena buscarlas de nuevo Lo que descubrí hace un par de semanas que parece que la el grupo la compañía que maneja esta la plataforma en sí tiene instrucciones sobre cómo eh facilitarlo Para evitar tener que lidar con Tok y cosas así así es que yo diría
(2:15:33) que alub sí eh el uso de que ha en el si es que hagamos las dos cosas diría veamos qué hicimos el año pasado y también yo puedo investigar esto eh específico para la plataforma el problema que yo que yo tuve que parece que varios lo tuvieron es de que al momento de hacer un comit el passw de hth Hub no lo reconocía sino que pedía un token Hay que entrar a la configuración de Hub e crear un token que es básicamente otra contraseña más larga y
(2:16:10) que guardarla en algún archivo y después cuando uno quiera hacer el comit volver y poner esa esa clave y eso es todo Sí reconozco ese dolor de cabeza y me parece que esto que estoy mencionando lo lo reduce y es doble contraseña al final básicamente pero más o menos lo habéis ido resolviendo o sur o es algo que os surge cada así recurente Rodrigo no una vez que ya creé el el token lo guardé ahora cada vez que hago un comit lo ingreso en vez de la contraseña y ahí no
(2:16:43) hay problema muy bien Este era un problema que También tenían otros grupos solo pu por creo que la situación sería igual para todos Sí seguramente alguien más le le ha tenido que surgir el problema sal en tu como Ah bueno Ah vale Marina pensaba que era alguien más no si querías comentar yí en nuestro grupo también a mí me salía por ejemplo y no logré no logré arreglarlo eh Ahora completo tu ordenador Sí sí pero bueno
(2:17:20) porque r me estaba dando problemas pero sí me daba el problema de la password y del token y ahí ya decidí dejar a hctor que Héctor No le estaba dando problemas e trabajar en su local Pues sí sería bueno ponerlo en algún apartado y guardarlo para otros años porque es algo que seguirá surgiendo muy bien Charles nos querías comentar algo de tu grupo o cualquier yo puse la página yo puse la página del token ahí
(2:17:49) en el chat del del discord del grupo contaminación busc el chat ahí debe estar la página para configurar el token individual por supuesto cada uno Configura su token Ah ya lo veo pues Si quieres copiarlo tú mismo al hakat 2024 al Canal general y así ya lo tiene todo el mundo Gracias Rodrigo lo compartí el enlace que menciona Rodrigo al Canal este ya está ahí muy bien más cosas Charles no hoy nada que comentar porque
(2:18:32) no Hemos llegado al punto de subir subir ficheros alit Entonces no hemos encontrado los problemas estamos viendo un poco como compartir Los ficheros de datos de copernic porque al final son muchos gigas de datos entonces vo a ver cómo solucionar esto muy bien Ah Diana tienes la mano levantada Sí gracias Igual para comentar lo del token también a mí me ayudó mucho que la clase de Jorge Cornejo de markdown y gi ahí viene super bien explicado en la en la parte inicial
(2:19:11) Entonces yo seguí los pasos eso me ayudó mucho perfecto alen más quiere comentar algún problema de los que hemos comentado o algú otro que haya tenido siempre va bien ponerlo en común y Qué es mejor la el general New token el Classic o el beta del enlace que habéis mandado Rodrigo yo usé el clásico por lo menos las vers bet siempre arrojan errores Vale gracias vale Diana ya nos ha puesto el enlace a lo que comentaba de de Jorge Así que para los que
(2:20:00) usáis ahí tenéis la [Música] información muy bien alguna otra algún otro comentario que queráis hacer al respecto de los proyectos o de los tutoriales o cualquier cosa que queráis comentar pues parece que no hay más comentarios yo no sé al final creo que hemos hecho un poco lo que estaba planeado para mañana ahora mismo si os parece Entonces mañana haríais reunión también o sea hoy ahora terminamos esta sesión veis con vuestro grupo Igual que ayer si hacéis
(2:20:45) un descanso más largo más corto y hacéis la sesión de proyecto y mañana sería lo mismo hemos puesto una hora orientativa pero creo que cada grupo se puede organizar como quiera a ver en concreto teníamos trabajo en grupo de proyectos a partir de las 10 hora Guatemala pero bueno eso lo habláis también entre vosotros y organizáis las horas como queráis trabajar en grupo si queréis y el viernes tendríamos ya la presentación de de los grupos como
(2:21:17) dijimos nos bueno estaría bien que todo el mundo hiciera una parte de la presentación en cada grupo para que todo el mundo tenga un poco tenga un poco la oportunidad de presentar su trabajo Bueno una parte de este proyecto colaborativo Así que si os parece Pues eso os preparáis una presentación compartida el viernes tenemos hora de comienzo de sesión a las 8 cu de Guatemala tenéis alguna duda comentario
(2:21:53) Los tutores organizadores queréis al se me olvida alguna cosa que queráis comentar Sí si alguien más tiene dudas del grupo de pesca que yo he mencionado esa pero no sé si alguien más tiene más dudas es un buen momento para mencionarla sobre todo si mañana no va a haber Ronda de dudas puede ser buen un buen momento lo mismo no tenéis eh pero por si acaso De todas formas por el discord podéis seguir escribiendo también
(2:22:29) comentarios o lo que sea y y os vamos ayudando entre todos Sí yo solo comentar una cosa que también lo comentamos en nuestro grupo de que no duden de si tienen dudas No solo de los tutoriales de esta vezo si es también de los tutoriales intermedios o iniciales que tamb aprovechen que esta semana estamos disponibles para para ponernos dudas por el disco Muy bien pues no sé si Emilio quieres comentar alguna cosa más de
(2:23:08) parte por la parte de la plataforma o no principio nada más perfecto pues lo dejaríamos entonces aquí Muchas gracias a todos y como digo seguimos en contacto por el discord y ya en reiniciamos sesión el viernes a las 88:15 hora de Guatemala Muchas gracias a Daniela también Gracias hasta el viernes seguimos en contacto por escrito Hasta luego ahora mandaremos enl mandaremos enlace para la reunión