diff options
Diffstat (limited to 'ast/ast_tester/testskyframe.f')
-rw-r--r-- | ast/ast_tester/testskyframe.f | 89 |
1 files changed, 0 insertions, 89 deletions
diff --git a/ast/ast_tester/testskyframe.f b/ast/ast_tester/testskyframe.f deleted file mode 100644 index 6359d63..0000000 --- a/ast/ast_tester/testskyframe.f +++ /dev/null @@ -1,89 +0,0 @@ - program testskyframe - implicit none - include 'SAE_PAR' - include 'AST_PAR' - - integer status, sf1, sf2, fs - double precision vals(5) - - status = sai__ok - - sf1 = ast_skyframe( 'system=fk5,epoch=2015.0', status ) - sf2 = ast_skyframe( 'system=fk5,epoch=2015.1', status ) - fs = ast_convert( sf1, sf2, 'SKY', status ) - if( fs .eq. AST__NULL ) then - call stopit( status, 'Error 0' ) - end if - - if( .not. ast_isaunitmap( ast_getmapping( fs, ast__base, - : ast__current, status ), - : status ) ) then - call stopit( status, 'Error 1' ) - end if - - if( ast_getd( sf1, 'SkyTol', status ) .ne. 0.001D0 ) then - call stopit( status, 'Error 2' ) - end if - - call ast_setd( sf2, 'SkyTol', 1.0D-6, status ) - fs = ast_convert( sf1, sf2, 'SKY', status ) - - if( ast_isaunitmap( ast_getmapping( fs, ast__base, - : ast__current, status ), - : status ) ) then - call stopit( status, 'Error 3' ) - end if - - sf2 = ast_skyframe( 'system=fk5,epoch=2016.6', status ) - fs = ast_convert( sf1, sf2, 'SKY', status ) - if( ast_isaunitmap( ast_getmapping( fs, ast__base, - : ast__current, status ), - : status ) ) then - call stopit( status, 'Error 4' ) - end if - - - vals(1) = 6.1D0 - vals(2) = 6.15D0 - vals(3) = 6.2D0 - vals(4) = 6.25D0 - vals(5) = 6.3D0 - call ast_axnorm( sf1, 1, 0, 5, vals, status ) - if( vals(1) .ne. 6.1D0 .or. - : vals(2) .ne. 6.15D0 .or. - : vals(3) .ne. 6.2D0 .or. - : vals(4) .ne. 6.25D0 .or. - : vals(5) .ne. 6.3D0 - 2*AST__DPI ) then - call stopit( status, 'Error 5' ) - end if - - call ast_axnorm( sf1, 1, 1, 5, vals, status ) - if( vals(1) .ne. 6.1D0 - 2*AST__DPI .or. - : vals(2) .ne. 6.15D0 - 2*AST__DPI .or. - : vals(3) .ne. 6.2D0 - 2*AST__DPI .or. - : vals(4) .ne. 6.25D0 - 2*AST__DPI .or. - : vals(5) .ne. 6.3D0 - 2*AST__DPI ) then - call stopit( status, 'Error 6' ) - end if - - if( status .eq. sai__ok ) then - write(*,*) 'All SkyFrame tests passed' - else - write(*,*) 'SkyFrame tests failed' - end if - - end - - subroutine stopit( status, text ) - implicit none - include 'SAE_PAR' - integer status - character text*(*) - - if( status .ne. sai__ok ) return - status = sai__error - write(*,*) text - - end - - |