Add SFe_climatology.py

[dougiesquire]

This PR adds a very simple script to calculate the monthly climatology of soluble Iron from CESM-MIMI_1980-2015_CAM4-6MEAN_MonthlyDep_Hamiltonetal2020.nc (https://doi.org/10.7298/xqqj-qk90).

Example output from the script can be found at: /g/data/ik11/inputs/WOMBAT/CESM-MIMI_1980-2015_CAM4-6MEAN_MonthlyDep_Hamiltonetal2020_clim.nc

Inputs at different model resolutions will be created from this climatology using the regrid_forcing.py script, e.g. see here: /g/data/tm70/ds0092/model/inputs/wombat/SFe_Hamiltonetal2020_monthly_clim.nc

I've successfully tested input files created using this workflow in an ACCESS-OM3+WOMBAT configuration

[anton-seaice]

Looks good

Climatology_bounds and cell_methods could be included - https://cfconventions.org/Data/cf-conventions/cf-conventions-1.12/cf-conventions.html#climatological-statistics but as we don't / probably can't achieve cf-compliance its fairly moot

[dougiesquire]

Thanks. I think I will add these.

[dougiesquire]

I've added cell_methods and climatology_bounds in cd2e695 - updated file is here: /g/data/ik11/inputs/WOMBAT/CESM-MIMI_1980-2015_CAM4-6MEAN_MonthlyDep_Hamiltonetal2020_clim.nc

@anton-seaice, could you please re-review/approve?

[anton-seaice]

We've made it this far...

--------------------------------------------------------------------------------
                         IOOS Compliance Checker Report                         
                                 Version 5.3.0                                  
                     Report generated 2025-09-03T05:50:45Z                      
                                    cf:1.10                                     
http://cfconventions.org/Data/cf-conventions/cf-conventions-1.10/cf-conventions.html
--------------------------------------------------------------------------------
                               Corrective Actions                               
CESM-MIMI_1980-2015_CAM4-6MEAN_MonthlyDep_Hamiltonetal2020_clim.nc has 7 potential issues


                                     Errors                                     
--------------------------------------------------------------------------------
§2.2 Data Types
* The variable climatology_bounds failed because the datatype is int64

§3.3 Standard Name
* Attribute long_name or/and standard_name is highly recommended for variable SFe
* Attribute long_name or/and standard_name is highly recommended for variable time
* Attribute long_name or/and standard_name is highly recommended for variable lon
* Attribute long_name or/and standard_name is highly recommended for variable lat


                                    Warnings                                    
--------------------------------------------------------------------------------
§2.3 Naming Conventions
* global attribute Dust+Comb Dep should begin with a letter and be composed of letters, digits, and underscores
* global attribute Fire+Anth Dep should begin with a letter and be composed of letters, digits, and underscores

§2.6 Attributes
* §2.6.2 global attribute title should exist and be a non-empty string
* §2.6.1 Conventions field is not present

§4.1 Latitude Coordinate
* latitude variable 'lat' should define standard_name='latitude' or axis='Y'

§4.2 Longitude Coordinate
* longitude variable 'lon' should define standard_name='longitude' or axis='X'

§7.4 Climatological Statistics
* Climatology variable coordinates are in improper order: ('nv', 'time'). Bounds-specific dimensions should be last
* Climatology variable coordinates are in improper order: ('nv', 'time'). Bounds-specific dimensions should be last
* The "time: method within years/days over years/days" format is not correct in variable SFe.
* The "time: method within years/days over years/days" format is not correct in variable SFe.
WARNING: The following exceptions occurred during the cf checker (possibly indicate compliance checker issues):
cf.check_time_coordinate: 'NoneType' object has no attribute 'group'

[dougiesquire]

Thanks for checking. I'm gonna call that good enough for now