1 | include(dom.inc)
2 |
3 | SUBROUTINE GRAY_CASE(all_k_abs,kabs_gray,ncells,nkabs,homosyst)
4 |
5 | IMPLICIT NONE
6 |
7 | ! IN
8 | CHARACTER*3 :: homosyst
9 | DOM_INT :: ncells,nkabs
10 | DOM_REAL,DIMENSION(ncells) :: kabs_gray
11 |
12 | ! OUT
13 | DOM_REAL,DIMENSION(ncells,nkabs) :: all_k_abs
14 |
15 | ! LOCAL
16 | DOM_INT :: ielt
17 |
18 | !------------------------!
19 | ! Non-homogeneous system !
20 | !------------------------!
21 |
22 | IF (homosyst=='NO') THEN
23 |
24 | DO ielt=1,ncells
25 |
26 | all_k_abs(ielt,:)=kabs_gray(ielt)
27 |
28 | ENDDO
29 |
30 | !--------------------!
31 | ! Homogeneous system !
32 | !--------------------!
33 |
34 | ELSE IF (homosyst=='YES') THEN
35 |
36 | DO ielt=1,ncells
37 |
38 | all_k_abs(ielt,:)=kabs_gray(1)
39 |
40 | ENDDO
41 |
42 | END IF
43 |
44 | END SUBROUTINE GRAY_CASE
grey_case.F could be called by: