wolfhece.mar.Interface_MAR_WOLF_objet

Created on Tue Jan 3 16:31:47 2023

@author: jbrajkovic

Module Contents

class wolfhece.mar.Interface_MAR_WOLF_objet.MAR_input_data(xsummits=np.zeros(0), ysummits=np.zeros(0), date_debut=datetime.datetime(2020, 7, 11, 5), date_fin=datetime.datetime(2020, 7, 11, 5), directory='~/BUP_srv7/', directory_hist_sim='~/BUP_srv7/', model_name='MIROC6', var='MBRR', var_unb='E', UnborNot=0, syu=1981, eyu=2010, mod_ydays=1, generate_quantiles=1)[source]

directory_hist_sim = '~/BUP_srv7/'[source]

UnborNot = 0[source]

var_unb = 'E'[source]

var = 'MBRR'[source]

xsummits[source]

ysummits[source]

date_debut[source]

date_fin[source]

directory = '~/BUP_srv7/'[source]

mod_ydays = 1[source]

fn[source]

ds[source]

lons[source]

lats[source]

Lb72[source]

mask[source]

vec_data[source]

directory_unbiasing = '/srv7_tmp1/jbrajkovic/These/IRM/'[source]

fn_quant_ev = '/srv7_tmp1/jbrajkovic/These/Unbiasing/evapotranspiration_quantiles_1981_2010.nc'[source]

fn_quant_pr = '/srv7_tmp1/jbrajkovic/These/Unbiasing/precipitation_quantiles_1981_2010.nc'[source]

syu = 1981[source]

eyu = 2010[source]

generate_quantiles = 1[source]

model_name = 'MIROC6'[source]

time[source]

mask_rectangles()[source]: Creates the rectangular mask so MAR values can be extracted only for the precised zone

plot_mask()[source]

select_MARdata()[source]

Inputvar:nom de la variable hydro MAR (string): date_debut:date initiale (vecteur[heure,jour,mois,année] date_fin:idem pour date finale directory:répertoire avec simus MAR (en fonction du GCM/scénario) mask: masque spatiale(matrice de 0 et 1 de la zone d’intéret)

Description : Sélectionne la variable hydro MAR, pour les pixels du masque. Retourne une matrice 2D avec toutes les valeurs MAR pour tous les pas de temps exemple: 5 pas de temps et 100 pixels , output = matrice de dimensions(100,5)

MAR_unbiasing()[source]

find_timestep()[source]: Routine qui trouve le time step de MAR en heures

make_time()[source]: formatte une matrice avec la date pour chaque pas de temps en heure,jour,mois,année à redévelopper si pas de temps inférieurs à l’heure

MAR_summits()[source]: utilise les longitudes et latitudes des centres des pixels MAR pour calculer les coordonnées des sommets des pixels en Lambert 72 outputs: deux matrices contenant pour chaque pixels les 4 coordonnées des 4 sommets

MAR_shapefile(name, dirout1)[source]: cette routine sort les pixels MAR au format shapefile le nom donné dans le sous-dossier GRID

MAR_TextOutputs(dirout1, extension='.rain.txt')[source]: sortie au format texte 1 fichier par polygone nom du fichier = ID du polygone.rain

MAR_BinaryOutputs(dirout1, extension='.rain.dat')[source]: sortie au format texte 1 fichier par polygone nom du fichier = ID du polygone.rain

MAR_BinaryOutputs_V2(dirout1, extension='.rain.dat')[source]: Sortie au format binaire 1 fichier par polygone nom du fichier = ID du polygone + extension

MAR_BinaryOutputs_temperatures(dirout1, extension='.rain.dat')[source]: Sortie au format binaire 1 fichier par polygone nom du fichier = ID du polygone + extension

MAR_BinaryOutputs_temperatures_daily(dirout1, extension='.rain.dat')[source]: Sortie au format binaire 1 fichier par polygone nom du fichier = ID du polygone + extension

MAR_BinaryOutputs_temperatures_daily_text(dirout1, extension='.rain.dat')[source]: Sortie au format binaire 1 fichier par polygone nom du fichier = ID du polygone + extension

mask_subbasin(bassin_shp='', buff=100)[source]

cumulated_volume(fp='/climato_tmp1/jbrajkovic/BUP_srv7/Drainage_basin_tabreux.shp')[source]: Area of the shapefile

precip_time(time_ch=datetime.datetime(2000, 1, 1, 0), maille_num=5)[source]

test_prepro(xor1, yor1, dimx, dimy, res=100)[source]: Test si les volumes sont les mêmes que ceux générés par wolf sur une grill 100x100

test_prepro_mat(xor1, yor1, dimx, dimy, res=100, sdate=datetime.datetime(2020, 1, 1, 0), edate=datetime.datetime(2020, 1, 1, 0), out_fn='', pref_out='')[source]: Returns the interpolated matrix at 100 m between the two dates

test_vect()[source]

compute_coverage_and_coords(fn1, fn2, epsg=31370)[source]

Compute coverage fraction and return coordinates of intersected polygons.

Returns:

fractions (np.ndarray) – Fraction (0–1) of each polygon in fn2 covered by fn1
clipped_coords (list) – List of (x, y) coordinate arrays for each intersected polygon (None if no intersection)

check_input_volume(fn1, fn2, epsg=31370, dir_out=None)[source]

wolfhece.mar.Interface_MAR_WOLF_objet.dir_ds = '/phypc11_tmp3/MARv3.14/MARv3.14-EUi-MIROC6-5km-ssp585/'[source]