8 implicit real*8 (a-h,o-z)
31 call geom(crftd1(1),crftd1(2),crftd1(3),crr1,crftd1(4))
32 call geom(crftu1(1),crftu1(2),crftu1(3),crr1,crftu1(4))
37 if(ftmp(k,i,j).le.1.0d-15)
then
38 fiolvl1(k,i,j)=ftmp(k,i,j)
40 call geom(ftmpb(k,i,j),ftmpa(k,i,j),ftmp(k,i,j),
41 1 crftrd1,fiolvl1(k,i,j))
50 if(ftmp(k,i,j).le.1.0d-15)
then
51 fiolvl1(k,i,j)=ftmp(k,i,j)
53 call geom(ftmpb(k,i,j),ftmpa(k,i,j),ftmp(k,i,j),
54 1 crftru1,fiolvl1(k,i,j))
72 call geom(crftd2(1),crftd2(2),crftd2(3),crr2,crftd2(4))
73 call geom(crftu2(1),crftu2(2),crftu2(3),crr2,crftu2(4))
78 if(ftmp(k,i,j).le.1.0d-15)
then
79 fiolvl2(k,i,j)=ftmp(k,i,j)
81 call geom(ftmpb(k,i,j),ftmpa(k,i,j),ftmp(k,i,j),
82 1 crftrd2,fiolvl2(k,i,j))
91 if(ftmp(k,i,j).le.1.0d-15)
then
92 fiolvl2(k,i,j)=ftmp(k,i,j)
94 call geom(ftmpb(k,i,j),ftmpa(k,i,j),ftmp(k,i,j),
95 1 crftru2,fiolvl2(k,i,j))
103 call xntpln(pcrft,pl(lvlcrft),pl(lvlcrft+1),crftd1(4),
105 call xntpln(pcrft,pl(lvlcrft),pl(lvlcrft+1),crftu1(4),
111 call xntpln(pcrft,pl(lvlcrft),pl(lvlcrft+1),
112 1 fiolvl1(k,i,j),fiolvl2(k,i,j),fiocrft(k,i,j))
121 cfdown(ksza)=crftfd*
pi
125 crftzd(ksza,i,j)=fiocrft(1,i,j)+fiocrft(2,i,j)
131 crftzu(ksza,m,j)=fiocrft(1,i,j)+fiocrft(2,i,j)
139 ef=amuo*efactb(nolyr+1)
143 cttdn(i,j)=(fiocrft(1,i,j)+fiocrft(2,i,j))
149 cttup(m,j)=(fiocrft(1,i,j)+fiocrft(2,i,j))
154 crfttup(ksza,i,j)=cttup(i,j)*ftot
155 crfttdn(ksza,i,j)=cttdn(i,j)*ftot
161 ef=amuo*efactb(nolyr+1)
165 crfttup(ksza,i,j)=cttup(i,j)*ftot
166 crfttdn(ksza,i,j)=cttdn(i,j)*ftot
