       DOUBLE PRECISION FUNCTION   FSAT (AL,BL,CL,DL,EL,FL,TIME)
       IMPLICIT REAL *8 (A-J,L-M,O-Z), INTEGER *4(K,N)
       PW=2512.
       PMAF=-0.00608
       PLFF= 0.085
       PMFKD=.0078
       FLUX=PMAF*AL*DCOS(PW*TIME)+PMAF*BL*DCOS(PW*TIME -2.0944) +
      *     PMAF*CL*DCOS(PW*TIME+2.0944)+PLFF*DL+PMFKD*EL
       FSAT=DABS(DCOS(PW*TIME -FL)*FLUX)
       IF (FL.EQ .1.0) FSAT=DABS(FLUX)
       RETURN
       END
