[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)

      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
          darklat  = sunlat + 90.0

      print *, 'Dark latitude: ', darklat

      CALL exit(1)

Attachment: darklat.ps
Description: PostScript document

reply via email to

[Prev in Thread] Current Thread [Next in Thread]