Global map at surface and at 500 mt (averaged over a period and at one precise time) -> earth_diagnostics coupled with MapGenerator or s2dverification can do the work

Time series (globally averaged and for basins or specific coordinates) -> also in this case earth_diags + s2dv (or MG in the future)

Latitude/Longitude VS time -> in this case we need to do some development in MapGenerator

vertical section -> with MapGenerator

total mass for PO4 -> we can combine two diagnostics with earth_diags + s2dv

Depth integration for P04 -> earth_diags + MapGenerator