Skip to content

Latest commit

 

History

History
48 lines (40 loc) · 2.03 KB

README.md

File metadata and controls

48 lines (40 loc) · 2.03 KB
Raw

pyWOMBAT

Here we provide 0D chemostat and 1D water column models of the Australian ocean biogeochemical model WOMBAT (World Ocean Model of Biogeochemistry and Trophic-dynamics). It is written in python and intended for development purposes.

*** 0D model functional *** *** 1D model non-functional ***

0D model components

  • run_0D.py [Master file]
  • bgc_sms_0D.py [BGC routine]
  • mass_balance_0D.py [Checks for mass balance conservation]

1D model components

  • run.py [Master file]
  • bgc_sms.py [BGC routine]
  • advect_diff.py [Mixing and advection]
  • phyparams.py [Physical input parameters and z-grid]
  • bgcparams.py [Biogeochemical input parameters]
  • tra_init.py [Tracer initialisation]
  • mass_balance.py [Checks for mass balance conservation]

Structure of run_0D.py

  1. Set timestep, run length and logicals
  2. Initialise tracers
  3. Loops through model [calls bgc_sms_0D].. [calls mass_balance_0D]
  4. Final call of BGC routine [calls bgc_sms_0D]
  5. Save restart to output file
  6. Visualise the results

Structure of run.py

  1. Set timestep, run length and logicals
  2. Loads physical parameters and vertical grid [calls phyparams]
  3. Initialise tracers [calls tra_init]
  4. Loads biogeochemical parameters [calls bgcparams]
  5. Loops through model [calls advect_diff].. [calls mix_mld].. [calls bgc_sms].. [calls massb_c].. [calls massb_n].. [calls massb_p].. [calls massb_fe]
  6. Final call of BGC routine [calls bgc_sms]
  7. Save output to output file
  8. Visualise the results