<commit_before>from __future__ import division

import numpy as np
import joblib

from scipy.interpolate import interp1d
from scipy.interpolate import UnivariateSpline as USpline

import sim
import nuclide
import na
import production

constraints = {
    'n_samples':      10,
    'h_surface':      0.0,   # surface elevation, assumed constant, m
    'lat':            65.0,
    'rho':            2.67,  # g/cm2
    # define max allowable parameter values to normalize to
    'max_dz_m':       10,    # m
    'min_dz_m':       0.01,  # m
    't_postgl':       12000, #yr
    'nuclide':        nuclide.Be10Qtz(),
    'n_gl':           30,     # number of glaciations
    'alt':            0,
}

con = constraints

# evenly spaces samples in m
con['bottom_depth_m'] = 10
con['sample_depths_m'] = np.linspace(0, con['bottom_depth_m'], con['n_samples'])
con['sample_depths'] = con['sample_depths_m'] * 100 * con['rho']
con['bottom_depth'] = con['sample_depths'][-1] 
con['max_dz'] = con['max_dz_m'] * 100 * con['rho']
con['min_dz'] = con['min_dz_m'] * 100 * con['rho']
  
# Define our "true" model
# 25 m per 20 kyr glaciation, 80kyr interglacials
con['t_gl'] = 20000
con['t_int'] = 80000
t_gl = con['t_gl']
t_int = con['t_int']
con['dz_true_m']  = 2.5 # m / glaciation
dz_true = np.ones(con['n_gl']) * (con['dz_true_m'] * con['rho'] * 100)
assert con['max_dz'] > dz_true[0]

mtrue = dz_true / con['max_dz']
<commit_msg>Experiment with lowering the max limit on depths removed.<commit_after>from __future__ import division

import numpy as np
import joblib

from scipy.interpolate import interp1d
from scipy.interpolate import UnivariateSpline as USpline

import sim
import nuclide
import na
import production

constraints = {
    'n_samples':      10,
    'h_surface':      0.0,   # surface elevation, assumed constant, m
    'lat':            65.0,
    'rho':            2.67,  # g/cm2
    # define max allowable parameter values to normalize to
    'max_dz_m':       3,    # m, setting a lower max to see if we still get
                            # high old erosion rate
    'min_dz_m':       0.01,  # m
    't_postgl':       12000, #yr
    'nuclide':        nuclide.Be10Qtz(),
    'n_gl':           30,     # number of glaciations
    'alt':            0,
}

con = constraints

# evenly spaces samples in m
con['bottom_depth_m'] = 10
con['sample_depths_m'] = np.linspace(0, con['bottom_depth_m'], con['n_samples'])
con['sample_depths'] = con['sample_depths_m'] * 100 * con['rho']
con['bottom_depth'] = con['sample_depths'][-1] 
con['max_dz'] = con['max_dz_m'] * 100 * con['rho']
con['min_dz'] = con['min_dz_m'] * 100 * con['rho']
  
# Define our "true" model
# 25 m per 20 kyr glaciation, 80kyr interglacials
con['t_gl'] = 20000
con['t_int'] = 80000
t_gl = con['t_gl']
t_int = con['t_int']
con['dz_true_m']  = 2.5 # m / glaciation
dz_true = np.ones(con['n_gl']) * (con['dz_true_m'] * con['rho'] * 100)
assert con['max_dz'] > dz_true[0]

mtrue = dz_true / con['max_dz']

{'fraction_non_alphanumeric': 0.145069000336587, 'fraction_numerical': 0.041400201952204646, 'mean_word_length': 3.180028129395218, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'dir': 'git-commits-cleaned', 'id': '1955732', 'n_tokens_mistral': 1198, 'n_tokens_neox': 1098, 'n_words': 284}