def _log_autoclean(input_file,
network_model,network_params,
station_model,station_params,
outlier_tol,
output_file):
msg = '\n'
msg += '------------- PYGEONS AUTOCLEAN RUN INFORMATION --------------\n\n'
msg += 'input file : %s\n' % input_file
msg += 'network :\n'
msg += ' model : %s\n' % ', '.join(network_model)
msg += ' parameter units : %s\n' % ', '.join(composite_units(network_model,gpnetwork.CONSTRUCTORS))
msg += ' east parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_params['east']])
msg += ' north parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_params['north']])
msg += ' vertical parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_params['vertical']])
msg += 'station :\n'
msg += ' model : %s\n' % ', '.join(station_model)
msg += ' parameter units : %s\n' % ', '.join(composite_units(station_model,gpstation.CONSTRUCTORS))
msg += ' east parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_params['east']])
msg += ' north parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_params['north']])
msg += ' vertical parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_params['vertical']])
msg += 'outlier tolerance : %s\n' % outlier_tol
msg += 'output file : %s\n\n' % output_file
msg += '--------------------------------------------------------------\n'
logger.info(msg)
return
def _log_reml_results(input_file,
network_model,network_params,network_fix,
station_model,station_params,station_fix,
likelihood,output_file):
msg = '\n'
msg += '-------------------- PYGEONS REML RESULTS --------------------\n\n'
msg += 'input file : %s\n' % input_file
msg += 'network :\n'
msg += ' model : %s\n' % ' '.join(network_model)
msg += ' parameter units : %s\n' % ', '.join(composite_units(network_model,gpnetwork.CONSTRUCTORS))
msg += ' fixed parameters : %s\n' % ', '.join(network_fix.astype(str))
msg += ' optimal east parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_params['east']])
msg += ' optimal north parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_params['north']])
msg += ' optimal vertical parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_params['vertical']])
msg += 'station :\n'
msg += ' model : %s\n' % ' '.join(station_model)
msg += ' parameter units : %s\n' % ', '.join(composite_units(station_model,gpstation.CONSTRUCTORS))
msg += ' fixed parameters : %s\n' % ', '.join(station_fix.astype(str))
msg += ' optimal east parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_params['east']])
msg += ' optimal north parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_params['north']])
msg += ' optimal vertical parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_params['vertical']])
msg += 'log likelihood :\n'
msg += ' east : %s\n' % likelihood['east']
msg += ' north : %s\n' % likelihood['north']
msg += ' vertical : %s\n' % likelihood['vertical']
msg += 'output file : %s\n\n' % output_file
msg += '--------------------------------------------------------------\n'
logger.info(msg)
return msg
def _log_reml(input_file,
network_model,network_params,network_fix,
station_model,station_params,station_fix,
output_file):
msg = '\n'
msg += '---------------- PYGEONS REML RUN INFORMATION ----------------\n\n'
msg += 'input file : %s\n' % input_file
msg += 'network :\n'
msg += ' model : %s\n' % ', '.join(network_model)
msg += ' parameter units : %s\n' % ', '.join(composite_units(network_model,gpnetwork.CONSTRUCTORS))
msg += ' fixed parameters : %s\n' % ', '.join(network_fix.astype(str))
msg += ' initial east parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_params['east']])
msg += ' initial north parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_params['north']])
msg += ' initial vertical parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_params['vertical']])
msg += 'station :\n'
msg += ' model : %s\n' % ', '.join(station_model)
msg += ' parameter units : %s\n' % ', '.join(composite_units(station_model,gpstation.CONSTRUCTORS))
msg += ' fixed parameters : %s\n' % ', '.join(station_fix.astype(str))
msg += ' initial east parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_params['east']])
msg += ' initial north parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_params['north']])
msg += ' initial vertical parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_params['vertical']])
msg += 'output file : %s\n\n' % output_file
msg += '--------------------------------------------------------------\n'
logger.info(msg)
return msg
def _log_strain(input_file,
network_prior_model,network_prior_params,
network_noise_model,network_noise_params,
station_noise_model,station_noise_params,
start_date,stop_date,output_id,rate,vertical,
covariance,output_dx_file,output_dy_file):
msg = '\n'
msg += '--------------- PYGEONS STRAIN RUN INFORMATION ---------------\n\n'
msg += 'input file : %s\n' % input_file
msg += 'network prior :\n'
msg += ' model : %s\n' % ', '.join(network_prior_model)
msg += ' parameter units : %s\n' % ', '.join(composite_units(network_prior_model,gpnetwork.CONSTRUCTORS))
msg += ' east parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_prior_params['east']])
msg += ' north parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_prior_params['north']])
msg += ' vertical parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_prior_params['vertical']])
msg += 'network noise :\n'
msg += ' model : %s\n' % ' '.join(network_noise_model)
msg += ' parameter units : %s\n' % ', '.join(composite_units(network_noise_model,gpnetwork.CONSTRUCTORS))
msg += ' east parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_noise_params['east']])
msg += ' north parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_noise_params['north']])
msg += ' vertical parameters : %s\n' % ', '.join(['%0.4e' % i for i in network_noise_params['vertical']])
msg += 'station noise :\n'
msg += ' model : %s\n' % ', '.join(station_noise_model)
msg += ' parameter units : %s\n' % ', '.join(composite_units(station_noise_model,gpstation.CONSTRUCTORS))
msg += ' east parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_noise_params['east']])
msg += ' north parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_noise_params['north']])
msg += ' vertical parameters : %s\n' % ', '.join(['%0.4e' % i for i in station_noise_params['vertical']])
msg += 'start date for output : %s\n' % start_date
msg += 'stop date for output : %s\n' % stop_date
msg += 'number of output positions : %s\n' % len(output_id)
msg += 'ignore vertical deformation : %s\n' % (not vertical)
msg += 'return strain rates : %s\n' % rate
msg += 'return covariances : %s\n' % covariance
msg += 'output east derivative file : %s\n' % output_dx_file
msg += 'output north derivative file : %s\n\n' % output_dy_file
msg += '--------------------------------------------------------------\n'
logger.info(msg)