Light attenuation in Fasham

Discussion about coupled ecosystem models

Moderators: arango, robertson

Post Reply
Message
Author
LauraB
Posts: 23
Joined: Mon Jul 30, 2007 9:13 pm
Location: Dalhousie University
Contact:

Light attenuation in Fasham

#1 Unread post by LauraB »

Hello all,

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)))
Bio(i,k,iChlo) units should be mgChl*m, and should represent the integrated chlorophyll in a grid cell... Does Bio(i,k,iChlo) have these characteristics? As far as I understand fasham.h, Bio(i,k,iChlo) is chlorophyll in mg/m^3 (somewhere a comment says that "all the tracer at input have transport units: m Tunits....", and then creates Bio(i,k,indx) by dividing by Hz)

Thanks a lot for the input,
Laura

kfennel
Posts: 37
Joined: Thu Jun 12, 2003 6:55 pm
Location: Dalhousie University
Contact:

#2 Unread post by kfennel »

Hi Laura,

In the equation exp[-(kw + kchl . CHL) z] CHL represents the mean chlorophyll concentration between z and the surface (the chl integral divided by z).

The correct unit for kCHL is 1/(mg CHL m^{-2}).

Also note that the factor 1/z in front of the chl integral is missing in equation 5 in Fennel, Wilkin et al. (GBC 2006).

Hope this helps!
-Katja

longmtm
Posts: 55
Joined: Tue Dec 13, 2005 7:23 pm
Location: Univ of Maryland Center for Environmental Science

Re: Light attenuation in Fasham

#3 Unread post by longmtm »

Hi Laura,

I want to confirm that Katja is right. The Unit of AttChl in the bioFasham.in file should be[m2/(mg_Chl) rather than [1/(mg_Chl m2).This should be corrected in the file ROMS/External/bioFasham.in

and also Bio(i,k,iChlo) should have unit [milligram CHL/m3]

Wen

pmaccc
Posts: 74
Joined: Wed Oct 22, 2003 6:59 pm
Location: U. Wash., USA

Re: Light attenuation in Fasham

#4 Unread post by pmaccc »

I think that the units of AttChl are still given incorrectly in the latest version of ROMS/External/bio_Fennel.in:

! Light attenuation by chlorophyll [1/(mg_Chl m2)], {0.02486d0}.

AttChl == 0.02486d0

Should be [m2 / (mg Chl)] right?

User avatar
wilkin
Posts: 872
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: Light attenuation in Fasham

#5 Unread post by wilkin »

You are correct. The units of AttChl are m^2 / mg_Chl. This is because they are inverse meters per chlorophyll concentration.

At line 711 in fennel.h where AttChl is multiplied by the chlorophyll concentration (mg_Chl / m^3) it has to give units of inverse length ... to match units of AttSW ... and so that when Att (the sum of all attenuation factors) is multiplied by delta_z we have a nondimensional quantity to pass to the EXP( ) function.

Code: Select all

                Att=(AttSW(ng)+                                         &
     &               AttChl(ng)*Bio(i,k,iChlo)+                         &
     &               AttFac)*                                           &
     &               (z_w(i,j,k)-z_w(i,j,k-1))
                ExpAtt=EXP(-Att)
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

Post Reply