Copernicus SST L3S Observation Converter#1109
Conversation
Added a new satellite SST converter that converts a netCDF-based blended L3S SST product from Copernicus into DART's obs sequence files. Temperature in the obs seq file is in C. The obs error sd are read from the incoming data. An SD lower bound is set to avoid over-fitting. Observations with bad incoming QC are discarded.
Added an rst readme file. Also, included the new page in the toctree.
hkershaw-brown
left a comment
There was a problem hiding this comment.
Hi Moha,
Looking good.
Definite changes:
- call check_namelist_read to catch invalid namelist entries
- readme.rst add the info about the temp being in C in the obs seq to the docs
- add this new converter to the list of available converters in the docs.
Then a few suggestions:
There's a couple of typos that I spotted.
A couple of your loops are looping in the slow order. I am not sure if you want to keep this order to have the obs in a certain order in obs_seq.
Defaults matching in the docs/example input.nml and the code
Cheers,
Helen
| do ilon = 1, nlon | ||
| do ilat = 1, nlat |
There was a problem hiding this comment.
Same comment on loop ordering, typically you want the inner loop to be the fastest varying.
Note I'm not sure if the observation sequence order is better for science reasons to loop along latitude first.
| - Observation time and location | ||
|
|
||
| and writes them into a DART observation sequence file using the | ||
| ``SATELLITE_BLENDED_SST`` observation kind. |
There was a problem hiding this comment.
Trying to use type for observation (I know kind is all over the place).
| ``SATELLITE_BLENDED_SST`` observation kind. | |
| ``SATELLITE_BLENDED_SST`` observation type. |
| Only observations with quality level greater than a specified threshold | ||
| (default: > 3) are retained. |
There was a problem hiding this comment.
This is a nitpick, but from the docs it seems like you can change the specified threshold, whereas it is a hardcoded parameter.
integer, parameter :: OBS_QC_LOW_QUALITY = 3
| Only observations with quality level greater than a specified threshold | |
| (default: > 3) are retained. | |
| Only observations with quality level greater than 3 are retained. |
|
Also @todo add cmems_sst_to_obs executable to .gitignore |
|
Helen, thanks for the helpful comments. I believe I addressed all of them. |
Description:
A new observation converter for a blended satellite Sea Surface Temperature product from the Copernicus Marine Environment Monitoring Service. The converter reads CMEMS NetCDF files and generates DART
obs_seqfiles containing bias-corrected blended SST observations (in K) suitable for assimilation into ocean models such as ROMS. The PR also adds a readme file containing brief documentation for the new converter.Types of changes
Documentation changes needed?
Tests
I've tested the converter with several CMEMS observation files. These can be found here:
/glade/work/gharamti/inacawo/data_snippets/SSTChecklist for merging
Add cmems_sst_to_obs to .gitignore
Updated changelog entry
Documentation updated
Update conf.py
Checklist for release
Testing Datasets