My question now is about light attenuation, which is given by an exponential: exp[-(kw + kchl . CHL) z],
where CHL is the integrated chlorophyll concentration in the column from surface to level z.
In BioFasham.in, kchl (AttChl) is given in units of 1/(mgChl m^2). In the code, the computation of light attenuation is
Code: Select all
Att=EXP(-0.5_r8*(AttSW(ng)+AttChl(ng)*Bio(i,k,iChlo))* & & (z_w(i,j,k)-z_w(i,j,k-1)))
Thanks a lot for the input,