42 integer(i4b) :: n, n_data_kei_max
44 real(dp) :: r_val, d_dummy
47 n_data_kei_max = 10000
53 open(unit=11, iostat=ios, &
54 file=inpath//
'/general/kei.dat', &
56 if (ios /= 0) stop
' read_kei: Error when opening the kei file!'
58 read(unit=11, fmt=
'(a)') ch_dummy
59 read(unit=11, fmt=
'(15x,f5.0)') kei_r_max
60 read(unit=11, fmt=
'(15x,f5.0)') kei_r_incr
61 read(unit=11, fmt=
'(a)') ch_dummy
63 n_data_kei = nint(kei_r_max/kei_r_incr)
65 if (n_data_kei > n_data_kei_max) stop
' read_kei: Array kei too small!'
67 read(unit=11, fmt=
'(a)') ch_dummy
68 read(unit=11, fmt=
'(a)') ch_dummy
69 read(unit=11, fmt=
'(a)') ch_dummy
71 do n=-n_data_kei, n_data_kei
72 read(unit=11, fmt=*) d_dummy, kei(n)
75 close(unit=11, status=
'keep')