Parent Directory | Revision Log | Revision Graph | Patch
revision 1.5 by stephd, Thu Feb 12 16:11:46 2004 UTC | revision 1.6 by stephd, Thu Apr 1 22:58:07 2004 UTC | |
---|---|---|
# | Line 69 C$TAF INIT bio_export = static, 10 | Line 69 C$TAF INIT bio_export = static, 10 |
69 | DO k=1,nlev | DO k=1,nlev |
70 | C$TAF STORE lit = bio_export | C$TAF STORE lit = bio_export |
71 | atten=(k0*drF(k)*.5) | atten=(k0*drF(k)*.5) |
72 | if (k.gt.1) atten=(k0*drF(k-1)*.5) | if (k.gt.1) atten=atten+(k0*drF(k-1)*.5) |
73 | lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj)) | lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj)) |
74 | if (lit.lt.0.d0.or.lit.gt.150) then | if (lit.lt.0.d0.or.lit.gt.150) then |
75 | print*,'QQ lit', lit | print*,'QQ lit', lit |
|
ViewVC Help | |
Powered by ViewVC 1.1.22 |