OceanBioME.jl icon indicating copy to clipboard operation
OceanBioME.jl copied to clipboard

PAR units

Open jagoosw opened this issue 1 year ago • 2 comments

Hi all,

I think we need to clarify the units for PAR. LOBSTER (and other BGC models) seem to use W/m $^2$ while e.g. the kelp model use einsteim/m $^2$ /day.

Since PAR is defined as an integral between 400 and 700nm (by MODIS at least) then we can take the very rough approximaiton at the central wavelength for putting it into the BGC models, but I don't know if this is the best solution.

At the moment we're implicitly assuming W/m $^2$ and converting for the kelp model:

https://github.com/OceanBioME/OceanBioME.jl/blob/70ff47f3326cbe668eb92fdf0f6222a1ae96cdd2/src/Models/Individuals/SLatissima.jl#L536

(This is, as far as I can tell, also being done wrong). I'm not exactly sure when this got added in but we need to a) agree on units for the defaults, and b) agree how to convert for the different models.

I think the above should read:

...* day / (3.99e-10 * (299792458 / 550e-9))

(the second factor is ~545e12 so we've lost a factor somewhere in the past)

jagoosw avatar Mar 23 '23 18:03 jagoosw