Package 'mesoda'

Title: Mesoscale Data Assimilation - PhD Thesis
Description: En el futuro: mi tesis de doctorado. Por ahora será un compendio de análisis ordenados.
Authors: Paola Corrales [aut, cre]
Maintainer: Paola Corrales <[email protected]>
License: MIT + file LICENSE
Version: 0.0.0.9000
Built: 2024-11-04 05:05:51 UTC
Source: https://github.com/paocorrales/mesoda

Help Index

Calculate areal coverage

Description

Calculate areal coverage

Usage

calculate_areal_coverage(
  exp,
  run,
  ini_date,
  ciclos,
  acumulado = 1,
  q = c(1, 5, 10, 25, 50),
  lon_band = c(-66.5, -61.5),
  lat_band = c(-35.5, -29),
  wrf_path = "/home/paola.corrales/datosmunin3/EXP/",
  out_path = "/home/paola.corrales/datosmunin3/EXP/derived_data/ppacum/"
)

Arguments

exp

name of the experiment
run

if "ana" or "fcst
ini_date

init date and time
ciclos

number of hourly cycles/fcst
acumulado

time of accumulated
q

pp threshold
lon_band

longitud limits
lat_band

latitude limits
wrf_path

where the simulations are located
out_path

where the output file will be saved

Wrap para FSS

Description

Usa la función fss del paquete verification pero previamente requiere que las variables estén en matrices. También puede iterar para distintos q (valor de pp) y w (tamaño de la caja = w2+1)

Usage

calculate_fss(fcst, obs, q, w, binary = TRUE)

Arguments

fcst

matriz con precipitación pronosticada
obs

matriz con precipitación observada (mismas dimensiones que fcst)
q

número o vector numérico con umbrales a ser calculados
w

número o vector numérico con escalas, correspondiente al ancho de la caja
binary

bolean, TRUE = ambos campos son binarios.

Calculate the FSS for non-binary fields

Description

Adapted from the SpatialVx. Calculates the value of Fraction Skill Score (FSS) for multiple neighborhood sizes.

Usage

calculate_FSSvector_from_nonbinary_fields(fbin1, fbin2, nvector)

Arguments

nvector

A numeric vector containing neighborhood sizes for which the FSS values are to be calculated. Only positive odd values are allowed in the vector. A square neighborhood shape is assumed and the specified value represents the length of square side.
fbind1

A numeric matrix representing the first field. Any values are allowed in the matrix.
fbind2

A numeric matrix representing the first field. Any values are allowed in the matrix. The matrix needs to have the same dimensions as fbin1.

Calculate probability matching

Description

Calculate probability matching

Usage

calculate_probability_marching_pp(
  exp,
  run,
  ini_date,
  ciclos,
  acumulado = 1,
  wrf_path = "/home/paola.corrales/datosmunin3/EXP/",
  out_path = "/home/paola.corrales/datosmunin3/EXP/derived_data/ppacum/"
)

exp

Arguments

exp

name of the experiment
run

if "ana" or "fcst
ini_date

init date and time
ciclos

number of hourly cycles/fcst
acumulado

time of accumulated
wrf_path

where the simulations are located
out_path

where the output file will be saved

Format

An object of class character of length 1.

Radiance to BT

Description

Radiance to BT

Usage

calculate_rad(rad, ncfile)

Arguments

rad

radiance
ncfile

NetCDF file

Discretiza una variable en cajas y devuelve labes lindos

Description

Discretiza una variable en cajas y devuelve labes lindos

Usage

cut_round(x, breaks)

Arguments

x

variable a discretizar
breaks

vector con breaks

Calculate Reduced Centered Random Variable (RCRV)

Description

Calculate Reduced Centered Random Variable (RCRV)

Usage

get_RCRV(dt, tipo.var = "superficie", calculo = "temporal")

Arguments

dt

data.table with the diag file
tipo.var

character, superficie == variables de superficie, perfil = variables en altura, radianzas = randianzas.
calculo

character, temporal == calculo en el tiempo, box = por cajas, nada = por variable

Goes Proyection

Description

Goes Proyection

Usage

goes_projection(x, y, ncfile)

Arguments

x

regular grid
y

regular grid
ncfile

NetCDF file

Input conventional observation errors from errtable to diag files

Description

Input conventional observation errors from errtable to diag files

Usage

input_obs_error(variable, type, level, path_to_errtable = "errtable.csv")

Arguments

variable

character
type

numeric, BUFR code
level

numeric
path_to_errtable

path to errtable file

Radiance to BT

Description

Radiance to BT

Usage

rad_to_tb(rad, ncfile, h = 6.629e-34, c = 299800000, kb = 1.381e-23)

Arguments

rad

radiance
ncfile

NetCDF file
h

constant
c

constant
kb

constant

Lectura de archivos diag para observaciones convencionales

Description

Internamente define la fecha. El formato de nombre de archivo debería ser "asim_conv_date.mem"

Usage

read_diag_conv(
  file_list,
  exp,
  member = "000",
  variable = c("uv", "p", "t", "q")
)

Arguments

file_list

lista de archivos
exp

nombre del experimento asociado
member

miembro del ensamble en caracteres. Si es la media correspode "000"
variable

vector con las variables que se quieren obtener.

Lectura de archivos diag para radianzas

Description

Internamente define el sensor, plataforma, la fecha y el miembro del ensamble asociado. El formato de nombre de archivo debería ser "asim_sensorplatdate.mem"

Usage

read_diag_rad(file_list, this_exp)

Arguments

file_list

lista de archivos
exp

nombre del experimento asociado

Lee archivo "satbias" de coeficientes

Description

Esta función parsea el archivo satbias y lo ordena en formato tidy.

Usage

read_satbias(file)

Arguments

file

archivo a leer

Humedad relativa

Description

Humedad relativa

Usage

rh(QVAPOR, P, T)

Arguments

QVAPOR

vector con la humedad especifica en kg/kg
P

vector con la presión en Pa
T

vector con la temperatura en kelvin

Color scale for cloud tops

Description

Color scale for cloud tops

Usage

scale_color_topes(
  colours = hex,
  limits = c(-90, 50),
  breaks = seq(-90, 50, 20),
  guide = guide_colorbar(barheight = 15),
  ...
)

Arguments

colours

hex format
limits

vector with limits
breaks

vector with breaks
...

others

Transformación de escalas a symlog

Description

Permite transformar una escala (colores, ejes, etc.) a symlog.

Usage

symlog_trans(base = 10, thr = 1, scale = 1)

References

Un buen samaritano en https://stackoverflow.com/a/14674703

Dew Point

Description

Dew Point

Usage

td(QVAPOR, P)

Arguments

QVAPOR

vector con la humedad especifica en kg/kg
P

vector con la presión en Pa

Temperatura

Description

Temperatura

Usage

tk(T, P, T_BASE = 300)

Arguments

T

vector con la temperatura potencial en kelvin
P

vector con la presión en Pa
T_BASE

valor numérico con la temperatura de base del WRF, usualmente 300 K

Rota los vientos e interpola a la retícula de masa

Description

La función calcula la rotación correcta e interpola a la reticula de masa. Solo funciona para Labert

Usage

uvmet(ncfile)

Arguments

ncfile

path al archivo netCDF donde están los datos

Proyección de coordenadas

Description

Proyecta y reproyecta coordenadas utilizando el paquete proj4

Usage

wrf_project(
  lon,
  lat,
  inverse = FALSE,
  round = c(-3, -4),
  map_proj = norargentina_lambert
)

Arguments

lon, lat

longitud y latitud o x e y para calcular la inversa
inverse

si es FALSE la función devuelve las coordenadas cartográficas (distancia al centro del dominio). Si es TRUE devuelve lon, lat.
round

igual a TRUE resuelve problemas de redondeo
map_proj

recibe un string con la proyección correspondiente en el formato aceptado por proj4