1 : C FAU_M2 SOURCE BECC 10/05/05 21:15:05 6674 2 : function fau_m2(M,coef) 3 : C 4 : real*8 fau_m2, M, coef 5 : C 6 : C coef = 1.0D0 for m1_plus 7 : C coef = -1.0D0 for m1_minus 8 : C 9 : fau_m2 = M + coef 10 : fau_m2 = fau_m2 * fau_m2 11 : fau_m2 = coef * 0.25D0 * fau_m2 12 : C 13 : return 14 : end 15 : 16 : 17 :
© Cast3M 2003 - All rights reserved.
Disclaimer