[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
atan2 converted to cabs?
From: |
James Kuyper |
Subject: |
atan2 converted to cabs? |
Date: |
Mon, 13 Jan 2003 15:07:58 -0500 |
Attached is a simple Fortran program I tried to print. I've also
attached the postscript file that was created by v4.12 of a2ps. For some
reason, the 'atan2' in my Fortran code has been converted to a 'cabs' in
the postscript.
--
James Kuyper
MODIS Level 1 Lead
Science Data Support Team
(301) 352-2150
include 'PGS_CBP.f'
include 'PGS_SMF.f'
integer pgs_cbp_earth_cb_vector
character*27 asciiutc
double precision offsets
double precision sun_vector(3)
double precision rho
double precision sunlat
double precision darklat
integer retval
data offsets /0.0/
data asciiutc /'2002-03-22T12:00'/
retval = pgs_cbp_earth_cb_vector(1, asciiutc, offsets, pgsd_sun,
1 sun_vector)
IF (retval .NE. PGS_S_SUCCESS) THEN
print *, 'PGS_CBP_Earth_CB_Vector(', asciiutc, ') returned',
1 retval
C CALL exit(0)
ENDIF
rho = sqrt(sun_vector(1)**2 + sun_vector(2)**2)
sunlat = atan2(sun_vector(3), rho) * 45.0/atan(1.0)
IF(sunlat .GT. 0.0) THEN
darklat = sunlat - 90.0
ELSE
darklat = sunlat + 90.0
ENDIF
print *, 'Dark latitude: ', darklat
CALL exit(1)
END
darklat.ps
Description: PostScript document
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- atan2 converted to cabs?,
James Kuyper <=