Arts2 {
VectorCreate( data_f_grid )
VectorCreate( data_t_grid )
VectorSet( data_f_grid, [ 230e9, 240e9 ] )
VectorSet( data_t_grid, [ 220, 250, 270] )
complex_refr_indexIceMatzler06(
        data_f_grid = data_f_grid,
        data_T_grid = data_t_grid )
#WriteXML("ascii", complex_refr_index)
GriddedField3Create(ref)
ReadXML(ref, "TestRefice.complex_refr_indexREFERENCE.xml")
Compare(complex_refr_index, ref, 1e-3)
}