To create a better initialization field from our Adriatic CTD data, I was thinking about taking the OA'd field, and running the model forward for a couple of days with the observed (weak) wind forcing but strongly nudging back to the initial field along the sections of initial CTD data.
The reason is that while we had dense CTD coverage in the Northern Adriatic, we only cut through the mid and south adriatic pits with one detailed section each (sort of a "feature oriented") approach to sampling. So I though if we assimilate hard to these sections, the density structure from these sections will be advected by the geostrophic flow around the basins and come up with something more reasonable than that obtained by simply using OA.
But to assimilate strongest at the data sections seems to require setting the nudging coefs with short time scales at the CTD locations, and large time scales elsewhere. I thought I could set this in a NetCDF file, but it looks like right now it's only possible to do analytically in set_nudgcof.F.
Should I hack reading some file into set_nudgcof.F, or is there a better way?
Richard P. Signell