Arts2 {

# First run basic tests of implementation
#####
Tensor4Create( out_rad )
RT4Test( out_rad=out_rad )
#WriteXML( in=out_rad, output_file_format="ascii" )

Tensor4Create( ref_rad )
ReadXML( out=ref_rad, filename="artscomponents/polradtran/refRT4" )
Compare( out_rad,ref_rad,5e-12 )

}