anchor.diff:  C
         SZAD = SZA * 180. / PI
         DIPD = DIP * 180. / PI
  -      
  +
   C      write (6,444) IDATE, UT, GLAT, GLONG, F107, F107A, AP
   C  444 FORMAT (' Date=',i5,' UT=',f6.0,' Lat=',f5.1,' Lon=',f6.1,
   C     >        ' F107=',f4.0,' F107A=',f4.0,' Ap=',f4.0)
  
positive.code: C   Python connection
        PyZ = Z
  C     Hack alert, didn't get fancier due to newer GLOW version coming soon enough
        idate=Pyidate; ut=Pyut; glat=Pyglat; glong=Pyglong; 
        f107a=Pyf107a; f107=Pyf107; f107p=Pyf107p
        ap=Pyap; ef=Pyef; ec=Pyec
  C
  
negative.code: C  445 FORMAT (' SZA=',F5.1,' LST=',F5.2,' Dip=',F5.1,
  C     >        ' Ec=',F6.3,' Ie=',I1)
  C
  C Output photoionization, electron impact ionization,
  C electron density, and ion densities:
  C
  C     write (6,690)
  