Parent Directory | Revision Log | Revision Graph | Patch
revision 1.6 by stephd, Thu Apr 1 22:58:07 2004 UTC | revision 1.7 by stephd, Thu Aug 25 22:32:44 2005 UTC | |
---|---|---|
# | Line 68 C$TAF INIT bio_export = static, 10 | Line 68 C$TAF INIT bio_export = static, 10 |
68 | lit=sfac(j) | lit=sfac(j) |
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)*hFacC(i,j,k,bi,bj)*.5) |
72 | if (k.gt.1) atten=atten+(k0*drF(k-1)*.5) | if (k.gt.1) atten=atten+(k0*drF(k-1) |
73 | & *hFacC(i,j,k-1,bi,bj)*.5) | |
74 | lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj)) | lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj)) |
75 | if (lit.lt.0.d0.or.lit.gt.150) then | if (lit.lt.0.d0.or.lit.gt.150) then |
76 | print*,'QQ lit', lit | print*,'QQ lit', lit |
|
ViewVC Help | |
Powered by ViewVC 1.1.22 |