IDL Files
The file jkmoore_idl.tar was missing the file: oceanarea_x3p.dat
The file x3p_routines.pro was updated to include the "get_x3p_coords" routine.
The tar file jkmoore_idl.tar has been updated to reflect the above changes.
Input Files
The file jkmoore_inputs.tar was missing the input file: woa98_sst_surf_m_x3p.dat
There were three NaN values in the two files: bio_slope_x3p_globalV5.dat and hmix_levpd_min25_x3p.dat. Both of these files were corrected and replaced in the tar file: jkmoore_inputs.tar
The tar file jkmoore_inputs.tar has been updated to reflect the above changes.
Final References
DSR pubs are now out, please cite these papers if you use this code:
2002. Moore, J.K., Doney, S.C.,Kleypas, J.A., Glover,D.M., Fung, I.Y., An intermediate complexity marine ecosystem model for the global domain, Deep-Sea Res. II, 49: 403-462.
2002. Moore, J.K., Doney, S.C., Glover, D.M., Fung, I.Y., Iron cycling
and nutrient limitation patterns in surface waters of the world ocean,
Deep-Sea Res. II, 49: 463-508.
Updates to Mixed Layer Model Since Publication
Variable thetaNmax
There is evidence that the maximum chl/C (chl/N) ratio can be higher
in diatoms than other phytoplankton groups (Geider et al., 1997). I suggest
using a thetaNmax of 4.0 for the diatoms and 3.0 for the other phytoplankton
groups. This improves global chlorophyll fields in terms of winter season
gyre values (too high around borders of the gyres in Moore et al. 2002).
In depth-resolved applications this also improves the behavior of diatom
biomass at the deep chlorophyll max.
New lambda function - for all phytoplankton groups
replace: lambda = lambda_no3 * max(frat_sp, 0.5) ! used
in DSR papers
with: lambda = lambda_no3
* (frat_sp * 0.5 + 0.5)
This is a smoother function more true to the data presented by Geider (1992).