1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
| REAL, DIMENSION(1:2), INTENT(IN) :: SOLAD !incoming direct solar radiation (w/m2)
REAL, DIMENSION(1:2), INTENT(IN) :: SOLAI !incoming diffuse solar radiation (w/m2)
REAL, DIMENSION(1:2), INTENT(IN) :: FABD !flux abs by veg (per unit incoming direct flux)
REAL, DIMENSION(1:2), INTENT(IN) :: FABI !flux abs by veg (per unit incoming diffuse flux)
REAL, DIMENSION(1:2), INTENT(IN) :: FTDD !down dir flux below veg (per incoming dir flux)
REAL, DIMENSION(1:2), INTENT(IN) :: FTID !down dif flux below veg (per incoming dir flux)
REAL, DIMENSION(1:2), INTENT(IN) :: FTII !down dif flux below veg (per incoming dif flux)
REAL, DIMENSION(1:2), INTENT(IN) :: ALBGRD !ground albedo (direct)
REAL, DIMENSION(1:2), INTENT(IN) :: ALBGRI !ground albedo (diffuse)
REAL, DIMENSION(1:2), INTENT(IN) :: ALBD !overall surface albedo (direct)
REAL, DIMENSION(1:2), INTENT(IN) :: ALBI !overall surface albedo (diffuse)
REAL, DIMENSION(1:2), INTENT(IN) :: FREVD !overall surface albedo veg (direct)
REAL, DIMENSION(1:2), INTENT(IN) :: FREVI !overall surface albedo veg (diffuse)
REAL, DIMENSION(1:2), INTENT(IN) :: FREGD !overall surface albedo grd (direct)
REAL, DIMENSION(1:2), INTENT(IN) :: FREGI !overall surface albedo grd (diffuse)
! output #输出与吸收的辐射，为可见光波段与近红外波段之和
REAL, INTENT(OUT) :: PARSUN !average absorbed par for sunlit leaves (w/m2)
REAL, INTENT(OUT) :: PARSHA !average absorbed par for shaded leaves (w/m2)
REAL, INTENT(OUT) :: SAV !solar radiation absorbed by vegetation (w/m2)
REAL, INTENT(OUT) :: SAG !solar radiation absorbed by ground (w/m2)
REAL, INTENT(OUT) :: FSA !total absorbed solar radiation (w/m2)
REAL, INTENT(OUT) :: FSR !total reflected solar radiation (w/m2)
REAL, INTENT(OUT) :: FSRV !reflected solar radiation by vegetation
REAL, INTENT(OUT) :: FSRG !reflected solar radiation by ground
NBAND = 2
! zero summed solar fluxes
SAG = 0.
SAV = 0.
FSA = 0.
! loop over nband wavebands
DO IB = 1, NBAND
! absorbed by canopy
CAD(IB) = SOLAD(IB)*FABD(IB)
CAI(IB) = SOLAI(IB)*FABI(IB)
SAV = SAV + CAD(IB) + CAI(IB)
FSA = FSA + CAD(IB) + CAI(IB)
! transmitted solar fluxes incident on ground
TRD = SOLAD(IB)*FTDD(IB)
TRI = SOLAD(IB)*FTID(IB) + SOLAI(IB)*FTII(IB)
! solar radiation absorbed by ground surface
ABS = TRD*(1.-ALBGRD(IB)) + TRI*(1.-ALBGRI(IB))
SAG = SAG + ABS
FSA = FSA + ABS
END DO
|