Parent Directory
|
Revision Log
|
Revision Graph
|
Patch
revision 1.3 by edhill, Thu Oct 9 04:19:19 2003 UTC | revision 1.4 by jmc, Tue Oct 28 22:19:53 2003 UTC | |
---|---|---|
# | Line 55 C$TAF INIT bio_export = static, 10 | Line 55 C$TAF INIT bio_export = static, 10 |
55 | lit=sfac(j) | lit=sfac(j) |
56 | DO k=1,nlev | DO k=1,nlev |
57 | C$TAF STORE lit = bio_export | C$TAF STORE lit = bio_export |
58 | atten=(k0*delz(k)*.5) | atten=(k0*drF(k)*.5) |
59 | if (k.gt.1) atten=(k0*delz(k-1)*.5) | if (k.gt.1) atten=(k0*drF(k-1)*.5) |
60 | lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj)) | lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj)) |
61 | if (lit.lt.0.d0.or.lit.gt.150) then | if (lit.lt.0.d0.or.lit.gt.150) then |
62 | print*,'QQ lit', lit | print*,'QQ lit', lit |
|
ViewVC Help | |
Powered by ViewVC 1.1.22 |