def test_scalar(self):
"""Should return a two-tuple of scalars for scalar input"""
SA = 35.0
CT = 10.0
output = gsw.entropy_first_derivatives(SA, CT)
self.assertEqual(type(output), tuple)
self.assertEqual(len(output), 2)
self.assertTrue(np.isscalar(output[0]))
def test_standard_values(self):
"""Test some standard values"""
SA = [34.7118, 34.8915, 35.0256, 34.8472, 34.7366, 34.7324]
CT = [28.8099, 28.4392, 22.7862, 10.2262, 6.8272, 4.3236]
eta_SA_standard = np.array((-0.263286800711655,
-0.263977276574528,
-0.255367497912925,
-0.238066586439561,
-0.234438260606436,
-0.232820684341694))
eta_CT_standard = np.array((13.221031210083824,
13.236911191313675,
13.489004628681361,
14.086599016583795,
14.257729576432077,
14.386429945649411))
eta_SA, eta_CT = gsw.entropy_first_derivatives(SA, CT)
self.assertTrue(np.all(abs(eta_SA - eta_SA_standard) < 1.0e-15))
self.assertTrue(np.all(abs(eta_CT - eta_CT_standard) < 1.0e-16))
#CT_SA_SA, CT_SA_pt, CT_pt_pt = gsw.CT_second_derivatives(SA_chck_cast, pt)
#test_print("CT_SA_SA")
#test_print("CT_SA_pt")
#test_print("CT_pt_pt")
#h_SA, h_CT, h_P = gsw.enthalpy_first_derivatives(SA_chck_cast, CT_chck_cast, gsw_cv.p_chck_cast)
#test_print("h_SA")
#test_print("h_CT")
#test_print("h_P")
#h_SA_SA, h_SA_CT, h_CT_CT = gsw.enthalpy_second_derivatives(SA_chck_cast, CT_chck_cast, gsw_cv.p_chck_cast)
#test_print("h_SA_SA")
#test_print("h_SA_CT")
#test_print("h_CT_CT")
eta_SA, eta_CT = gsw.entropy_first_derivatives(SA_chck_cast, CT_chck_cast)
test_print("eta_SA")
test_print("eta_CT")
eta_SA_SA, eta_SA_CT, eta_CT_CT = gsw.entropy_second_derivatives(SA_chck_cast, CT_chck_cast)
test_print("eta_SA_SA") #FIXME: diffs are not found in the original
test_print("eta_SA_CT") #FIXME: diffs are not found in the original
test_print("eta_CT_CT") #FIXME: diffs are not found in the original
#pt_SA, pt_CT = gsw.pt_first_derivatives(SA_chck_cast, CT_chck_cast)
#test_print("pt_SA")
#test_print("pt_CT")
#pt_SA_SA, pt_SA_CT, pt_CT_CT = gsw.pt_second_derivatives(SA_chck_cast, CT_chck_cast)
#test_print("pt_SA_SA")
#test_print("pt_SA_CT")
