-
Notifications
You must be signed in to change notification settings - Fork 16
Description
u_raw = xr.open_dataset('G:/data/ERA5/ERA5_2024.grib')['u']
v_raw = xr.open_dataset('G:/data/ERA5/ERA5_2024.grib')['v']
u_stda = metdig.utl.xrda_to_gridstda(
u_raw[0,0],lon_dim='longitude',lat_dim='latitude',level=[850])
v_stda = metdig.utl.xrda_to_gridstda(
v_raw[0,0],lon_dim='longitude',lat_dim='latitude',level=[850])
metdig.io.custom.split_stda_to_cache_sfc(u_stda,var_name='u',is_overwrite=False,data_name='cz')
metdig.io.custom.split_stda_to_cache_sfc(v_stda,var_name='v',is_overwrite=False,data_name='cz')
init_time = datetime.datetime(2024,6,1,0)
u = get_model_grid(init_time=init_time,fhour=None,data_name='custom',var_name='u',level='850')
v = get_model_grid(init_time=init_time,fhour=None,data_name='custom',var_name='v',level='850')
caldata = mdgcal.shear(u, v)