source: anuga_work/production/dampier_2006/check_data.py @ 6498

from os import sep, getcwd, access, F_OK, mkdir, getenv
from anuga.abstract_2d_finite_volumes.util import get_data_from_file
import os
from Numeric import zeros, array, allclose

#makes the csv files from the evolved model

home = getenv('INUNDATIONHOME') #Sandpit's parent dir   
#user = get_user_name()
data = 'data'
state = 'western_australia'
scenario_name = 'dampier'
scenario = 'dampier_tsunami_scenario_2006'
#scenario = 'test_dampier'
an = 'anuga'
bo = 'boundaries'

#run_time = '20070305_055256_run'
#run_time = '20070305_134855_run'
run_time = '20070321_164505_run'

topo = 'topographies'
out = 'outputs'
urs = 'urs'
gridded = '1_10000'

gauge_boundary_filename = 'gauges_time_series_first.csv'
gauge_evolved_filename = 'gauges_time_series_first.csv'

boundary_dir_filename = os.path.join(home,data,state,scenario,an,bo,gauge_boundary_filename)

evolved_dir_filename= os.path.join(home,data,state,scenario,an,out,run_time,gauge_evolved_filename)


print'evolved_dir_filename',evolved_dir_filename
e_time, e_stage, e_momentum, e_speed, e_elevation = get_data_from_file(evolved_dir_filename)

print'boundary_dir_filename',boundary_dir_filename
b_time, b_stage, b_momentum, b_speed, b_elevation = get_data_from_file(boundary_dir_filename)

 
# compares the 2 models
j=0
b_sample = []
for i in range(len(b_time)):
    if j<(len(e_time)) and b_time[i] == e_time[j]:
        b_sample.append(float(b_stage[i]))
        print 'time e equal b:', b_time[i], b_stage[i],i, j, b_sample[j], e_stage[j],(len(e_time)-1)
        j = j +1


assert allclose (b_sample, e_stage, 2, 0.1) 
print "test successful"