2 #include saprc99_mosaic_4bin_vbs2.spc
3 #include saprc99_mosaic_4bin_vbs2.eqn
8 !__________________________________________________
10 REAL(KIND
=dp
) FUNCTION
Keff ( A0
,B0
,C0
, TEMP
,X1
,X2
,y1
,y2
)
11 REAL(KIND
=dp
),INTENT(IN) :: X1
,X2
,y1
,y2
12 REAL(KIND
=dp
),INTENT(IN) :: TEMP
13 REAL(KIND
=dp
),INTENT(IN):: A0
,B0
,C0
14 Keff
= A0
* EXP(- B0
/TEMP
) &
15 *(TEMP
/300._dp
)**C0
*(y1
*X1
/(X1
+ X2
+ 1.0e-35) &
16 +y2
*(1-X1
/(X1
+ X2
+ 1.0e-35)))
18 !__________________________________________________
20 REAL(KIND
=dp
) FUNCTION
Keff2 ( C0
,X1
,X2
,y1
,y2
)
21 REAL(KIND
=dp
),INTENT(IN) :: X1
,X2
,y1
,y2
22 REAL(KIND
=dp
),INTENT(IN):: C0
23 Keff2
= C0
*(y1
*X1
/(X1
+ X2
+ 1.0e-35) &
24 +y2
*(1-X1
/(X1
+ X2
+ 1.0e-35 )))
26 !__________________________________________________