Index: trunk/doc/SDToolBox - Test Plan/testAntolinez/testAntolinez.py =================================================================== diff -u -r104 -r105 --- trunk/doc/SDToolBox - Test Plan/testAntolinez/testAntolinez.py (.../testAntolinez.py) (revision 104) +++ trunk/doc/SDToolBox - Test Plan/testAntolinez/testAntolinez.py (.../testAntolinez.py) (revision 105) @@ -39,9 +39,24 @@ lon = ds['lon'][:] #era5 = pd.read_pickle(os.path.join(pth.data,'NS_WAM_ERA5.p')) plt.plot(lat,lon,'or');plt.axis('equal') -#%% Extract data +#%% Extract data non-gridded # 1st initialize the class inda = sd_eda.InputData(input_coordinates={'LON':lon,'LAT':lat}, input_variables=['swh'],input_scenarios=['era5'], input_years=[1986]) ds = sd_eda.ExtractData.get_era_5(directory_path=os.path.join(pth.era5),input_data=inda) +dsxr = ds.generate_netcdf(dir_path=pth.out,base_name='test',dataset_code='era5') + +#%% Extract data gridded +lonmin=-5;lonmax=11;dlon=0.5 +latmin=50;latmax=62;dlat=0.5 +lon = np.arange(latmin,latmax+dlat,dlat) +lat = np.arange(lonmin,lonmax+dlon,dlon) +#lon,lat = np.meshgrid(lon,lat) +#plt.plot(lat,lon,'or');plt.axis('equal') +del(inda) +inda = sd_eda.InputData(input_coordinates={'LON':lon,'LAT':lat}, + input_variables=['msl_p'],input_scenarios=['era5'], + input_years=[1986],is_gridded=True) + +ds = sd_eda.ExtractData.get_era_5(directory_path=os.path.join(pth.era5),input_data=inda)