Update src/seismic_hazard_forecasting.py

This commit is contained in:
ftong 2025-06-30 18:21:22 +02:00
parent 4aef02a31e
commit d904101ba1

View File

@ -1,12 +1,7 @@
# -*- coding: utf-8 -*-
from datetime import timedelta, datetime
import numpy as np
from igfash.rate import apply_beast, datenum_to_datetime
def datenum_to_datetime(datenum):
# MATLAB epoch starts at year 0, Python's datetime starts at year 1
matlab_epoch = datetime(1, 1, 1)
days = datenum - 367 # Subtracting the offset for MATLAB's epoch BUT THIS COULD BE PROBLEMATIC
return matlab_epoch + timedelta(days=days)
def calc_bins2(dates, unit, bin_dur, dates_calc, rate_forecast, rate_unc_high, rate_unc_low, multiplicator, filename = "activity_rate.png", **kwargs):