Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

Question about how iron_flux is converted into dissolved Fe tendencies (which module handles it?)

yujiezhang

yujie zhang
New Member
Hi everyone,


I am currently working on modifying iron inputs in CESM2’s ocean biogeochemistry (POP2 + MARBL), and I would like to understand exactly how iron_flux is converted into tracer tendencies for dissolved iron (Fe) inside the model(e.g., mmol Fe m⁻³ s⁻¹).
 

mlevy

Michael Levy
CSEG and Liaisons
Staff member
The surface fluxes for all tracers (which includes those returned from MARBL, which includes iron) are used by the ocean GCM in the vertical diffusion and mixing processes.

For POP2, I believe all the fluxes are converted to tendencies in the vdiff() subroutine, specifically by the multiplication of dzr(k) (1/layer thickness, which has units of m⁻¹). Note that it's easy to get lost in the variable names, but STF refers to surface tracer fluxes and VTF = STF in the top layer of ocean cells (there is a merge() statement setting VTF earlier in the routine)

Code:
         VDTK(:,:,n) = merge((VTF(:,:,n,bid) - VTFB)*dzr(k), &
                             c0, k <= KMT(:,:,bid))
 
Top