def test(self):
cmd = stat0A(earth_model_stat=join(self.share_dir, 'earth.model'),
stat0_out=self.out_file,
l_min=1,
l_max=10,
fault_bottom_depth=60.,
fault_top_depth=3.,
obs_dep=0.,
if_skip_on_existing_output=True,
stdout=None,
stderr=None,
cwd=join(self.outs_dir, 'cwd'),
if_keep_cwd=True)
cmd()
def test(self):
cmd = stat0A(
earth_model_stat = join(self.share_dir,'earth.model'),
stat0_out = self.out_file,
l_min = 1,
l_max = 10,
fault_bottom_depth = 60.,
fault_top_depth = 3.,
obs_dep = 0.,
if_skip_on_existing_output = True,
stdout = None,
stderr = None,
cwd = join(self.outs_dir,'cwd'),
if_keep_cwd = True)
cmd()
def add_model(model_str, l_max, model_num):
cmd1['%d' % model_num] = stat0A(earth_model_stat=join(
model_str, 'earth.model_' + model_str),
stat0_out=join(model_str, 'stat0.out'),
l_min=1,
l_max=15000,
fault_bottom_depth=fault_bottom_depth,
fault_top_depth=fault_top_depth,
obs_dep=0.,
if_skip_on_existing_output=True,
stdout=FNULL)
cmd2['%d' % model_num] = vj.pollitz.ComputeEarthModelVISCO1DNonGravity(
earth_file=join(model_str, 'earth.model_' + model_str),
l_max=lmax,
outputs_dir=model_str,
if_skip_on_existing_output=True,
stdout=FNULL,
stderr=FNULL,
)
def add_model(model_str, l_max, model_num):
cmd1['%d'%model_num]= stat0A(
earth_model_stat = join(model_str, 'earth.model_'+model_str),
stat0_out = join(model_str, 'stat0.out'),
l_min = 1,
l_max = 15000,
fault_bottom_depth = fault_bottom_depth,
fault_top_depth = fault_top_depth,
obs_dep = 0.,
if_skip_on_existing_output = True,
stdout = FNULL
)
cmd2['%d'%model_num] = vj.pollitz.ComputeEarthModelVISCO1DNonGravity(
earth_file = join(model_str, 'earth.model_'+model_str),
l_max = lmax,
outputs_dir = model_str,
if_skip_on_existing_output = True,
stdout = FNULL,
stderr = FNULL,
)
FNULL = open(os.devnull, 'w')
fid = FaultFileIO('../fault_model/fault_He50km.h5')
fault_bottom_depth = fid.depth_bottom
fault_top_depth = fid.depth_top
lmax = 810
######################
earth_file_dir = 'He50km_Vis1.0E19/'
cmd1_He50km_Vis1E19 = stat0A(earth_model_stat=join(
earth_file_dir, 'earth.model_He50km_Vis1.0E19'),
stat0_out=join(earth_file_dir, 'stat0.out'),
l_min=1,
l_max=15000,
fault_bottom_depth=fault_bottom_depth,
fault_top_depth=fault_top_depth,
obs_dep=0.,
if_skip_on_existing_output=True,
stdout=FNULL)
cmd2_He50km_Vis1E19 = ComputeEarthModelVISCO1D(
earth_file=join(earth_file_dir, 'earth.model_He50km_Vis1.0E19'),
l_max=lmax,
outputs_dir=earth_file_dir,
if_skip_on_existing_output=True,
stdout=FNULL,
stderr=FNULL,
)
import os
from os.path import join
import argparse
from viscojapan.pollitz.pollitz_wrapper import stat0A
from viscojapan.fault_model import FaultFileReader
import viscojapan as vj
FNULL = open(os.devnull, 'w')
fid = FaultFileReader('../fault_model/fault_bott60km.h5')
fault_bottom_depth = fid.depth_bottom
fault_top_depth = fid.depth_top
model_str = 'He63km_VisM1.0E19'
cmd = stat0A(
earth_model_stat = join(model_str, 'earth.model_'+model_str),
stat0_out = join(model_str, 'stat0.out'),
l_min = 1,
l_max = 15000,
fault_bottom_depth = fault_bottom_depth,
fault_top_depth = fault_top_depth,
obs_dep = 0.,
if_skip_on_existing_output = True,
stdout = FNULL
)
cmd.run()
fid = FaultFileIO('../fault_model/fault_He50km.h5')
fault_bottom_depth = fid.depth_bottom
fault_top_depth = fid.depth_top
lmax = 810
######################
earth_file_dir = 'He50km_Vis1.0E19/'
cmd1_He50km_Vis1E19 = stat0A(
earth_model_stat = join(earth_file_dir, 'earth.model_He50km_Vis1.0E19'),
stat0_out = join(earth_file_dir, 'stat0.out'),
l_min = 1,
l_max = 15000,
fault_bottom_depth = fault_bottom_depth,
fault_top_depth = fault_top_depth,
obs_dep = 0.,
if_skip_on_existing_output = True,
stdout = FNULL
)
cmd2_He50km_Vis1E19 = ComputeEarthModelVISCO1D(
earth_file = join(earth_file_dir, 'earth.model_He50km_Vis1.0E19'),
l_max = lmax,
outputs_dir = earth_file_dir,
if_skip_on_existing_output = True,
stdout = FNULL,
stderr = FNULL,
)
