************************************************************************ * Recombination Fitting * ************************************************************************ real*8 zeta,nu,j *----- Read Parameter ------------------------- open ( unit = 80, file = 'init.par' , status = 'unknown' & , access = 'sequential', form = 'formatted' ) read(80,2000) zeta read(80,2000) nu close ( 80 ) 2000 format( 10x,f8.0) *----- Header --------------------------------- open ( unit = 80, file = 'recomb.dat' , status = 'unknown' & , access = 'sequential', form = 'formatted' ) write (80,1100) 1100 format(/,3x,' zeta ',' nu ', & ' J ( in equation (6.12) ) ') *----- Calculation --------------------------- call f(zeta,nu,j) write(80,1000) zeta,nu,j 1000 format(3x,1p3e17.8) close (80) *----- Main Routine End ----------------------- end *----- Subroutine ----------------------------------------------------- subroutine f(zeta,nu,j) real*8 zeta,nu,j *----- Set Parameter -------------------------- if ( nu.ge.-20.d0 .and. nu.lt. 0.d0) then include 'm20nu0.txt' elseif ( nu.ge. 0.d0 .and. nu.lt.10.d0) then include '0nu10.txt' else write(0,*) 'Bat nu Input! (-20