premsf.f

Go to the documentation of this file.

      REAL*8 FUNCTION premsf(M,T,RR)        
      REAL*8 rzams,m,t,rr
      REAL*8 rzamsf
      EXTERNAL rzamsf
*      
      rzams = rzamsf(m)
*       Coefficients in PREMSR
      IF(m.LE.1.d0)THEN
         pre1 = 0d0
         pre2 = 0d0
         pre3 = 7.432d-2 - 9.43d-2*m + 7.439d-2*m**2
      ENDIF
      IF(m.GT.1d0.AND.m.LT.2d0)THEN
         pre1 = -4.00772d0 + 4.00772d0*m
         pre2 = 8.5656d0 - 8.5656d0*m
         pre3 = -4.50678d0 + 4.56118d0*m
      ENDIF
      IF(m.GE.2d0)THEN
         pre1 = 1.60324d0 + 2.20401d0*m - 0.60433d0*m**2 
     &        + 5.172d-2*m**3
         pre2 = -4.56878d0 - 4.05305d0*m + 1.24575*m**2 
     &        - 0.10922d0*m**3
         pre3 = 3.01153 + 1.85745*m - 0.64290d0*m**2 
     &        + 5.759d-2*m**3
      END IF
      premsf = rzams*10**((pre1*t**3 + pre2*t**4 
     &       + pre3*t**5)/(1.05 - t)) - rr
      RETURN
      END