param.f 3.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119
  1. subroutine param(filen,psc,pce,pse,bb,tausc,tauce,tause,
  2. 1 fsc,fce,fse,aa)
  3. C implicit none
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. INCLUDE "param.inc"
  6. character(len=8) filen
  7. C double precision psc,pce,tausc,tauce,tause,theta,val,dRsc
  8. C double precision dRce,dRse,dTsc,dTce,dTse,fsc,fce,fse,Nsc,Nce
  9. C double precision sRsc,sRce,sRse,sTsc,sTce
  10. C double precision Nse,asct,acet,aset,ascr,acer,aser,cw,phi
  11. C double precision alpha,beta,sc,gama
  12. integer i,N
  13. character(len=30) line,ln,varname
  14. C double precision var,varvalue
  15. PRINT *, filen
  16. open(unit=1,file=filen)
  17. 1 read(1,'(a30)',end=3) LINE
  18. n=scan(line,'#',.false.)
  19. if(n.GT.0) goto 1
  20. n=scan(line,'=',.false.)
  21. if(n.eq.0) goto 1
  22. read(line(1:n-1),'(a5)') varname
  23. read(line(n+1:40),*) varvalue
  24. if(varname.eq."psc") then
  25. psc=varvalue
  26. elseif(varname.eq."pce") then
  27. pce=varvalue
  28. elseif(varname.eq."tausc") then
  29. tausc=varvalue
  30. elseif(varname.eq."tauce") then
  31. tauce=varvalue
  32. elseif(varname.eq."tause") then
  33. tause=varvalue
  34. elseif(varname.eq."theta") then
  35. theta=varvalue
  36. elseif(varname.eq."val") then
  37. val=varvalue
  38. elseif(varname.eq."dRsc") then
  39. dRsc=varvalue
  40. elseif(varname.eq."dRce") then
  41. dRce=varvalue
  42. elseif(varname.eq."dRse") then
  43. dRse=varvalue
  44. elseif(varname.eq."dTsc") then
  45. dTsc=varvalue
  46. elseif(varname.eq."dTce") then
  47. dTce=varvalue
  48. elseif(varname.eq."dTse") then
  49. dTse=varvalue
  50. elseif(varname.eq."sRsc") then
  51. sRsc=varvalue
  52. elseif(varname.eq."sRce") then
  53. sRce=varvalue
  54. elseif(varname.eq."sRse") then
  55. sRse=varvalue
  56. elseif(varname.eq."sTsc") then
  57. sTsc=varvalue
  58. elseif(varname.eq."sTce") then
  59. sTce=varvalue
  60. elseif(varname.eq."sTse") then
  61. sTse=varvalue
  62. elseif(varname.eq."fsc") then
  63. fsc=varvalue
  64. elseif(varname.eq."fce") then
  65. fce=varvalue
  66. elseif(varname.eq."fse") then
  67. fse=varvalue
  68. elseif(varname.eq."Nsc") then
  69. POPSC=varvalue
  70. elseif(varname.eq."Nce") then
  71. POPCE=varvalue
  72. elseif(varname.eq."Nse") then
  73. POPSE=varvalue
  74. elseif(varname.eq."MUsc") then
  75. MUSC=varvalue
  76. elseif(varname.eq."MUce") then
  77. MUCE=varvalue
  78. elseif(varname.eq."MUse") then
  79. MUSE=varvalue
  80. elseif(varname.eq."asct") then
  81. asct=varvalue
  82. elseif(varname.eq."acet") then
  83. acet=varvalue
  84. elseif(varname.eq."aset") then
  85. aset=varvalue
  86. elseif(varname.eq."ascr") then
  87. ascr=varvalue
  88. elseif(varname.eq."acer") then
  89. acer=varvalue
  90. elseif(varname.eq."aser") then
  91. aser=varvalue
  92. elseif(varname.eq."cw") then
  93. cw=varvalue
  94. elseif(varname.eq."phi") then
  95. phi=varvalue
  96. elseif(varname.eq."alpha") then
  97. aa=varvalue
  98. elseif(varname.eq."beta") then
  99. bb=varvalue
  100. elseif(varname.eq."sc") then
  101. sc=varvalue
  102. elseif(varname.eq."gama") then
  103. gg=varvalue
  104. elseif(varname.eq."ef") then
  105. ef=varvalue
  106. endif
  107. goto 1
  108. 3 continue
  109. return
  110. end