Recent Changes - Search:

Participantes

Proyectos

Cursos

Bitacoras

Support

upload Folder

edit SideBar

Rutinas MATLAB Isaac Rodriguez

Todas las rutinas de postprocesamiento, reprocesamiento, graficado de corrientes, GDOP, etc. se encuentran disponibles en el servidor OCOMEX en el directorio: /home/wera/work/matlab/hfr_process y tienen terminacion *_ir.m. Varias de estas rutinas solo son adaptaciones o modificaciones de la rutina original que no tiene terminacion *_ir.m

Listado y breve descripcion de las rutinas:

  • colorArrows_ir.m %Me permite hacer un quiver con las flechas de colores. El colorbar lo tengo que ingresar aparte y seria el que tiene matlab por default. (Se utiliza para graficar velocidades totales).
  • colorArrows2_ir.m % Me permite hacer un quiver con las flechas de colores, sin embargo el colorbar se genera en la misma rutina con colores especificos para cada rango de velocidad. (Se utiliza para graficar velocidades totales).
  • colorArrows_seismic_ir.m %Me permite hacer un quiver tomando en cuenta valores negativos y positivos (con 2 colores). El colorbar se genera en la misma rutina. (Se utiliza para graficar velocidades radiales).
  • fake_radials_ir.m %Esta rutina genera la cobertura espacial radial idealizada entre los sitios de SAU, ERE y BSQ. Para ese entonces BSQ todavia no estaba instalado y mediante esta rutina se genero una malla polar con radiales ficticios para este sitio.
  • gen_tot_RT_ir.m %Esta rutina sirve para generar archivos horarios de corrientes totales para el interior de la Bahia Todos Santos (como los que se encuentran en ocomex en ~/realtime/BTS) para un intervalo de tiempo dado. Lo unico que se debe de especificar e ingresar manualmente en el script, es la fecha inicial y final del periodo de tiempo que se quiere reprocesar las corrientes totales y el directorio en donde se desean guardar los datos reprocesados (actualmente esta en ocomex en ~/realtime/BTS/datos_reprocesados). Nota importante: Se pueden reeprocesar datos a partir del 26-may-2015 ya que es cuando los archivos *.RAD_Beam_av.mat y *.RAD_DirF_av.mat estan disponibles y generados al mismo tiempo para los 3 sitios (quinta, villarino y conalep). Aqui mismo se puede especificar si se desean reprocesar las corrientes totales tomando en cuenta Sau o si se desea tomar en cuenta la mascara que elimina los datos que estan fuera de la islas Todos santos. Si por alguna razon el script no corre, quiere decir que para la fecha inicial ingresada no hay datos (sugerencia: cambiar el dia inicial hasta que corra el programa).
  • gen_tot_RT_15_ir.m %Esta rutina carga los archivos generados cada 15 min de corrientes radiales (*.RAD_DirF.mat) de los 3 sitios (qta,vil y con; sau es opcional) y mediante la rutina rad2tot_ir.m hace la conversion de corrientes radiales a totales. Lo unico que se debe de especificar e ingresar manualmente en el script, es la fecha inicial y final del periodo de tiempo (en dias julianos) que se quiere reprocesar las corrientes totales y el directorio en donde se desean guardar los datos reprocesados (actualmente esta en ocomex en /home/wera/realtime/BTS/datos_reprocesados_15min).
  • gen_tot_CC_ir.m %Esta rutina sirve para generar archivos horarios de corrientes totales para la corriente de California (como los que se encuentran en ocomex en ~/realtime/CC) para un intervalo de tiempo dado. Lo unico que se debe de especificar e ingresar manualmente en el script, es la fecha inicial y final del periodo de tiempo que se quiere reprocesar las corrientes totales y el directorio en donde se desean guardar los datos reprocesados (actualmente esta en ocomex en ~/realtime/CC/datos_reprocesados). Nota importante: Se pueden reeprocesar datos a partir del 08-jun-2017 ya que es cuando los archivos *.RAD_DirF_av.mat estan disponibles y generados al mismo tiempo para los 3 sitios (sauzal, erendira y san quintin).
  • make_tot_CC_ir.m %Esta rutina carga el promedio horario de las corrientes radiales (*.RAD_DirF_av.mat) de los 3 sitios (sau, ere y bsq) y mediante la rutina rad2tot_CC_ir.m hace la conversion de corrientes radiales a totales. Nota importante: Mucho cuidado!!! con hacer modificaciones a esta rutina ya que es la que se ejecuta desde el crontab de ocomex mediante un link simbolico con la etiqueta make_tot_RT_xf. Esta rutina genera las corrientes totales en formato .mat .txt y una imagen en .jpg para la corriente de california.
  • make_tot_RT_ir.m %Esta rutina carga el promedio horario de las corrientes radiales (*.RAD_Beam_av.mat y *.RAD_DirF_av.mat) de los 3 sitios (qta, vil y con) y mediante la rutina rad2tot_ir.m hace la conversion de corrientes radiales a totales. Nota importante: Mucho cuidado!!! con hacer modificaciones a esta rutina ya que es la que se ejecuta desde el crontab de ocomex mediante un link simbolico con la etiqueta make_tot_RT_xf. Esta rutina genera las corrientes totales en formato .mat .txt y una imagen en .jpg para la Bahia de Todos Santos. Esta rutina no es inteligente!! si se cae por ejemplo la quinta, no se van a seguir generando corrientes totales entre vil y con. Por eso hay que especificar bien entre que estaciones se van a generar corrientes totales. En la actualidad la rutina esta configurada para obtener corrientes totales unicamente a partir de las estaciones de qta y con.
  • make_tot_RT_2_ir.m %Esta rutina es la misma que make_tot_RT_ir.m sin embargo esta hecha para poder hacer pruebas y que los cambios no sean guardados. Mediante la rutina rad2tot_2_ir.m se hace la conversion de corrientes radiales a totales. En esta rutina estan contempladas las 4 estaciones dentro de la Bahia de todos Santos (qta, vil, con y sau) y se puede configurar la rutina para que se generen corrientes totales omitiendo ciertas estaciones.
  • make_tot_RT_2_sin_sauzal_ir.m %Esta rutina es exactamente la misma que make_tot_RT_2_ir.m , solo que explicitamente esta omitida la estacion de sauzal. Es decir solo contempla las corrientes totales obtenidas utilizando las estaciones de qta, vil y con (aunque de igual manera se puede cambiar).
  • make_tot_RT_gdl_ir2.m %Esta rutina es como make_tot_RT_2_ir.m solo que esta hecha de manera generalizada de modo que en teoria esta diseņada para poder poner el numero de estaciones que sean dentro de la Bahia Todos Santos.
  • matrixUVT_ir.m %Esta rutina es una funcion que me permite generar una matriz tridimensional que contiene los siguientes datos de la Bahia de todos Santos: longitud (X), latitud (Y), corriente zonal (Um), corriente meridional (Vm) y tiempo (Tm) ingresando unicamente la fecha inicial y final del periodo de tiempo. Las corrientes totales se obtienen cada hora de ocomex en ~/realtime/BTS y se encuentran a partir del 23/01/2014
  • matrixUVT_2_ir.m %Esta rutina hace exactamente lo mismo que matrixUVT_ir.m solo que a diferencia de la anterior, no es una funcion, y se tiene que ingresar manualmente la fecha inicial y final del periodo de tiempo. Las corrientes totales se obtienen cada hora de ocomex en ~/realtime/BTS y se encuentran a partir del 23/01/2014
  • matrixUVT_3_ir.m %Esta rutina es la misma que matrixUVT_2_ir.m solo que carga los av.mat reprocesados (por la rutina gen_tot_RT_ir.m) que estan almacenados en una sola carpeta (/home/wera/realtime/BTS/datos_reprocesados)
  • matrixUVT_CC_ir.m %Esta rutina es similar a matrixUVT_3_ir.m solo que la matriz tridimensioanl generada contiene datos de la corriente de california. Actualmente esta configurada para cargar los datos reprocesados del crucero CICSE-NASA (/home/wera/realtime/CC/datos_reprocesados/crucero_cicese_nasa/).
  • plot_BTS_currents_ir.m %Rutina que grafica el ultimo archivo de corrientes radiales (para qta, vil, con y sau) y de corrientes totales dentro de la Bahia de Todos Santos.
  • plot_channels_ir.m %funcion que grafica los canales I y Q de todas las antenas a partir de un archivo comprimido .mat . Cada figura grafica NCHIRP (xaxis) vs MT (yaxis). El regreso del chirp no es considerado. Ejemplo de uso: plot_channels_ir('20162211630_vil.mat','vil').
  • plot_raw_data.m %funcion que grafica los chirps de los canales I y Q de todas las antenas a partir de un archivo comprimido .mat . Se puede graficar MT (xaxis) vs bits (yaxis) o MT (xaxis) vs miliVolts (yaxis). Ejemplo de uso: plot_raw_data('20171300130_bsq.mat','volts');
  • plot_gdop_CC_ir.m %Esta rutina grafica elipses de GDOP para los sitios que muestrean la corriente de california (sau, ere y bsq).
  • plot_gdop_ir.m %Esta rutina grafia elipses de GDOP para los sitios que muestrean el interior de la Bahia Todos Santos (qta, vil y con; sau es opcional)
  • GM_gdop_ir.m %Esta rutina grafica elipses de GDOP para los sitios que muestrean el Golfo de Mexico.
  • plot_radials_ir.m % Grafica corrientes radiales a partir de un archivo *RAD_DirF.mat o *RAD_DirF_av.mat. Las corrientes se grafican con 2 colores: Los valores positivos (color rojo) representan las corrientes radiales que se alejan respecto al radar y los valores positivos (color azul) representan las corrientes que se acercan en direccion hacia el radar. Esta rutina solo sirve para las estaciones que se encuentran de lado del Pacifico.
  • plot_radials2_ir.m % Grafica corrientes radiales a partir de un archivo *RAD_DirF.mat o *RAD_DirF_av.mat. Las corrientes se grafican solo en funcion de la magnitud (no hay valores negativos). Esta rutina solo sirve para las estaciones que se encuentran de lado del Pacifico.
  • plot_radials_gm_ir.m % Hace lo mismo que la rutina plot_radials_ir.m solo que con las estaciones que se encuentran en el Golfo de Mexico.
  • plot_totals_ir.m % Grafica las velocidades totales dentro de la Bahia Todos Santos.
  • plot_vel_mean_ir.m % Carga la matriz tridimensional de corrientes totales de CODAR y hace un grafico de corrientes totales promedio para toda la serie de tiempo.
  • rad2tot_ir.m % Esta rutina se encarga de hacer la conversion de corrientes radiales a totales para las estaciones que se encuentran dentro de la Bahia de todos Santos (qta, vil, con; sau es opcional). Esta rutina ademas emplea varios criterios para la generacion de corrientes totales: solo genera corrientes totales si entre 2 radiales existe un angulo mayor a 35 grados y menor a 145 grados (a excepcion de vil y con, donde hay un criterio extra que logra identificar si la corriente paralalela a estos 2 sitios es real), elimina corrientes con un RMSE mayor a 10. Nota importante: Mucho cuidado!!! con hacer modificaciones a esta rutina ya que es la que se ejecuta desde el crontab de ocomex mediante un link simbolico con la etiqueta make_tot_RT_xf y se encuentra dentro de la rutina make_tot_RT_ir.m
  • rad2tot_2_ir.m % Esta rutina hace exactamente lo mismo que rad2tot_ir.m solo que esta dentro de la rutina make_tot_RT_2_ir.m y esta hecha para poder hacer pruebas y graficar las corrientes totales encima de las corrientes radiales.
  • rad2tot_2_sin_sauzal_ir.m % Esta rutina es la misma que rad2tot_2_ir.m, solo que esta dentro de la rutina make_tot_RT_2_sin_sauzal_ir.m e ignora la estacion de sauzal para la generacion de corrientes totales.
  • rad2tot_gdl_ir2.m % esta rutina es similar a rad2tot_2_ir.m solo que esta hecha de manera generalizada para poder tomar en cuenta el numero de estaciones que sean dentro de la Bahia Todos Santos.
  • rad2tot_CC_ir.m % Esta rutina se encarga de hacer la conversion de corrientes radiales a totales para las estaciones que se encuentran en la corriente de california (sau, ere y bsq). Esta rutina ademas emplea varios criterios para la generacion de corrientes totales: solo genera corrientes totales si entre 2 radiales existe un angulo mayor a 40 grados y menor a 140 grados, y ademas elimina corrientes con un RMSE mayor a 10. Nota importante: Mucho cuidado!!! con hacer modificaciones a esta rutina ya que es la que se ejecuta desde el crontab de ocomex mediante un link simbolico con la etiqueta make_tot_RT_xf y se encuentra dentro de la rutina make_tot_CC_ir.m
  • rad2tot_CC_sin_criterios_ir.m % esta rutina es la misma que rad2tot_CC_ir.m solo que hace la conversion de corrientes radiales a totales se hace sin utilizar ningun criterio.
Edit - History - Print - Recent Changes - Search
Page last modified on August 14, 2017, at 05:42 PM PST