Code vapor
unknown
plain_text
2 years ago
1.7 kB
10
Indexable
""" Extracción de una serie de tiempo con WRF-python """ import numpy as np import xarray as xa from matplotlib import pyplot as plt from glob import glob from netCDF4 import Dataset import wrf import pandas as pd from datetime import datetime import sys from datetime import datetime, timedelta df = pd.read_csv('vismet_centrosur_pp_upto_20230627.csv', index_col="Fecha", parse_dates=True) files = [ Dataset(x) for x in sorted(glob('wrfout_d04_2023-04-26_00:00:00')) ] pr = wrf.getvar(files,'RAINNC', timeidx=wrf.ALL_TIMES) pc = wrf.getvar(files,'RAINC' , timeidx=wrf.ALL_TIMES) pos = wrf.ll_to_xy(files[0], -37.91, -71.45).values lon = pos[0] lat = pos[1] pr = pr.interp( west_east=lon, south_north=lat ) pc = pc.interp( west_east=lon, south_north=lat ) fig, ax = plt.subplots(2,1) plt.sca(ax[0]) plt.plot(pr.Time, pr+pc, label="Precipitación total") plt.plot(pr.Time, pc, label="Pprecipitación convectiva") plt.grid() plt.ylabel('mm') plt.legend() plt.sca(ax[1]) plt.plot(pr.Time[1:], np.diff(pr+pc), label="Tasa precipitación total") plt.plot(pr.Time[1:], np.diff(pc), label="Tasa precipitación convectiva") plt.ylabel('mm/hora') plt.grid() plt.legend() df = pd.read_csv('vismet_centrosur_pp_upto_20230627.csv', index_col="Fecha", parse_dates=True) df = df[datetime(2023,4,26):datetime(2023,4,29)] df = df["08313002-4"] df[np.isnan(df)] = 0 df.index = df.index - timedelta(hours=-3) plt.sca(ax[0]) plt.plot(df.index, df.cumsum(),'ko-',markersize=3, label="Precipitacion observada") plt.legend() plt.sca(ax[1]) plt.plot(df.index, df,'ko-',markersize=3, label="Tasa para la observacion") plt.legend() plt.show()
Editor is loading...
Leave a Comment