43 integer(i4b) :: n, n_data_kei_max
45 real(dp) :: r_val, d_dummy
48 n_data_kei_max = 10000
54 open(unit=11, iostat=ios, &
55 file=inpath//
'/general/kei.dat', &
57 if (ios /= 0) stop
' read_kei: Error when opening the kei file!'
59 read(unit=11, fmt=
'(a)') ch_dummy
60 read(unit=11, fmt=
'(15x,f5.0)') kei_r_max
61 read(unit=11, fmt=
'(15x,f5.0)') kei_r_incr
62 read(unit=11, fmt=
'(a)') ch_dummy
64 n_data_kei = nint(kei_r_max/kei_r_incr)
66 if (n_data_kei > n_data_kei_max) stop
' read_kei: Array kei too small!'
68 read(unit=11, fmt=
'(a)') ch_dummy
69 read(unit=11, fmt=
'(a)') ch_dummy
70 read(unit=11, fmt=
'(a)') ch_dummy
72 do n=-n_data_kei, n_data_kei
73 read(unit=11, fmt=*) d_dummy, kei(n)
76 close(unit=11, status=
'keep')
Declarations of kind types for SICOPOLIS.
Declarations of global variables for SICOPOLIS (for the ANT domain).
subroutine read_kei()
Reading of the tabulated kei function.
Declarations of global variables for SICOPOLIS.