Estaciones meteorológicas

Observado, ERA5-Land crudo y la serie corregida. La escala cambia la serie, la dispersión, la diferencia y la distribución empírica entre diaria, mensual y anual. El método parte en el ganador; puedes ver otro. En dirección del viento también aparece la rosa.

import { fmt } from "../_js/portal.js"
import { VARS, nombreVar, nombreEst, parseEsc, panelMeteo, aniosDisponibles, filtrarVentana, metodosConDato, aplicarMetodo } from "../_js/escalamiento.js"

kge = FileAttachment("../datos/series/esc_kge.csv").csv({typed: true})
elegido = FileAttachment("../datos/series/esc_elegido.csv").csv({typed: true})
finales = FileAttachment("../datos/series/esc_final_metodo.csv").csv({typed: true})
catalogo = FileAttachment("../datos/series/esc/catalogo.csv").csv({typed: true})
estaciones = catalogo.filter(d => d.tipo === "meteo").map(d => d.estacion).sort((a, b) => nombreEst(a).localeCompare(nombreEst(b), "es"))

viewof estacion = Inputs.select(estaciones, {
  label: "Estación",
  format: nombreEst,
  value: estaciones.includes("Villarica") ? "Villarica" : estaciones[0]
})

filasTodas = {
  const res = await fetch(`../datos/series/esc/meteo/${encodeURIComponent(estacion)}.csv`);
  if (!res.ok) throw new Error("No está el CSV de " + estacion);
  return parseEsc(d3, await res.text());
}

variables = [...new Set(filasTodas.map(d => d.variable))].sort((a, b) => nombreVar(a).localeCompare(nombreVar(b), "es"))
viewof variable = Inputs.select(variables, {
  label: "Variable",
  format: nombreVar,
  value: variables.includes("t2m") ? "t2m" : variables[0]
})
filasBase = filasTodas.filter(d => d.variable === variable)
metodos = metodosConDato(filasBase)
opt = elegido.find(d => d.estacion === estacion && d.variable === variable)
escrito = finales.find(d => d.estacion === estacion && d.variable === variable)
ganadorKge = kge.find(d => d.estacion === estacion && d.variable === variable && (d.ganador === true || d.ganador === "TRUE"))
ganador = escrito?.metodo_correccion_script15 || ganadorKge?.metodo || opt?.metodo_optimo || metodos[0]
viewof metodo = Inputs.select(metodos, {
  label: "Método",
  format: (m) => m === ganador ? `${m} (ganador)` : m,
  value: metodos.includes(ganador) ? ganador : metodos[0]
})
viewof escala = Inputs.radio(["Diaria", "Mensual", "Anual"], {
  label: "Escala de la serie",
  value: "Diaria"
})
anios = aniosDisponibles(filasBase)
opciones = ["Últimos 3 años", "Periodo común", ...anios.map(String)]
viewof ventana = Inputs.select(opciones, { label: "Ventana de la serie diaria", value: "Últimos 3 años" })
filas = aplicarMetodo(filasBase, metodo)

meta = VARS[variable] || { nombre: nombreVar(variable), unidad: "", agg: "mean" }
filasWs = variable === "wd" ? aplicarMetodo(filasTodas.filter(d => d.variable === "ws"), metodo) : null
html`<p>${nombreEst(estacion)}, ${meta.nombre}. ${metodo === ganador ? `La línea verde es <strong>${metodo}</strong>, el método ganador.` : `La línea verde es <strong>${metodo}</strong>. El ganador de esta variable es <strong>${ganador}</strong>.`} La climatología, el año hidrológico, el diagrama de Taylor y las tablas usan el periodo con observación, no la ventana diaria.</p>`

filasKge = kge.filter(d => d.estacion === estacion && d.variable === variable)
panelMeteo(Plot, {
  filas,
  filasVentana: filtrarVentana(filas, ventana),
  meta,
  lugar: nombreEst(estacion),
  filasKge,
  escala,
  metodo,
  ganador,
  filasWs,
  d3
})