def testCW(numiter):
params_matrix = get_params_matrix()
output = rice_model.cw(rice_model.mai1(params_matrix.get('mai')[numiter], params_matrix.get('a')[numiter]), params_matrix.get('dw')[numiter], params_matrix.get('dsed')[numiter], params_matrix.get('osed')[numiter], params_matrix.get('pb')[numiter], params_matrix.get('kd')[numiter])
print "Test of function name: %s expected: %i != calculated: %i" % ("cw", params_matrix.get('cw_out')[numiter], output)
assert round(output, 3) == round(params_matrix['cw_out'][numiter], 3)
def html_table(row_inp,iter):
dsed_temp=float(row_inp[0])
dsed.append(dsed_temp)
a_temp=float(row_inp[1])
a.append(a_temp)
pb_temp=float(row_inp[2])
pb.append(pb_temp)
dw_temp=float(row_inp[3])
dw.append(dw_temp)
osed_temp=float(row_inp[4])
osed.append(osed_temp)
mai_temp=float(row_inp[5])
mai.append(mai_temp)
kd_temp=float(row_inp[6])
kd.append(kd_temp)
mai1_temp=rice_model.mai1(mai_temp, a_temp)
mai1_out.append(mai1_temp)
cw_temp=rice_model.cw(mai1_temp,dw_temp,dsed_temp,osed_temp,pb_temp,kd_temp)
cw_out.append(cw_temp)
Input_header="""<table border="1">
<tr><H3>Batch Calculation of Iteration %s</H3></tr><br>
<tr>
<td><b>Input Name</b></td>
<td><b>Input value</b></td>
<td><b>Unit</b></td>
</tr>"""%(iter)
Input_mai="""<tr>
<td>Mass of Applied Ingredient Applied to Paddy</td>
<td>%s</td>
<td>kg</td>
</tr>""" %(mai_temp)
Input_dsed="""<tr>
<td>Sediment Depth</td>
<td>%s</td>
<td>m</td>
</tr>""" %(dsed_temp)
Input_a="""<tr>
<td>Area of the Rice Paddy</td>
<td>%s</td>
<td>m<sup>2</sup></td>
</tr>""" %(a_temp)
Input_pb="""<tr>
<td>Bulk Density of Sediment</td>
<td>%s</td>
<td>kg/m<sup>3</sup></td>
</tr>""" %(pb_temp)
Input_dw="""<tr>
<td>Water Column Depth</td>
<td>%s</td>
<td>m</td>
</tr>""" %(dw_temp)
Input_osed="""<tr>
<td>Porosity of Sediment</td>
<td>%s</td>
<td>-</td>
</tr>""" %(osed_temp)
Input_kd="""<tr>
<td>Water-Sediment Partitioning Coefficient</td>
<td>%s</td>
<td>L/kg</td>
</tr></table><br>""" %(kd_temp)
Output_header="""<table border="1">
<tr>
<td><b>Output Name</b></td>
<td><b>Output value</b></td>
<td><b>Unit</b></td>
</tr>"""
Output_mai1="""<tr>
<td>Application Rate</td>
<td>%0.2E</td>
<td>kg a.i./A</td>
</tr>""" %(mai1_temp)
Output_cw="""<tr>
<td>Peak & Chronic EEC</td>
<td>%0.2E</td>
<td>µg/L</td>
</tr></table><br>""" %(cw_temp)
Inout_table = Input_header+Input_mai+Input_dsed+Input_a+Input_pb+Input_dw+Input_osed+Input_kd+Output_header+Output_mai1+Output_cw
return Inout_table
def testMAI1(numiter):
params_matrix = get_params_matrix()
output = rice_model.mai1(params_matrix.get('mai')[numiter], params_matrix.get('a')[numiter])
print "Test of function name: %s expected: %i != calculated: %i" % ("mai1", params_matrix.get('mai1_out')[numiter], output)
assert round(output, 3) == round(params_matrix['mai1_out'][numiter], 3)
