Opened 17 years ago
Closed 17 years ago
#115 closed upgrade (Done)
Detided time-averaged fields
Reported by: | arango | Owned by: | arango |
---|---|---|---|
Priority: | major | Milestone: | Release ROMS/TOMS 3.1 |
Component: | Nonlinear | Version: | 3.1 |
Keywords: | Cc: |
Description
A new CPP option (AVERAGES_DETIDE) has been coded to write detided free-surface and 2D/3D momentum. These new fields (zeta_detided, ubar_detided, vbar_detided, u_detided, v_detided) are written to output time-averaged NetCDF (AVGNAME).
The detiding is somewhat elaborated. It solves the following least-squares problem:
F = Favg + SUM[A(k)*SIN(omega(k)*t)] + SUM[B(k)*COS(omega(k)*t)], k=1:Ntides
where Favg is the time-averaged field, omega(k) are the tidal frequencies, A(k) and B(k) are the fitting coefficients, and Ntides is the number of tidal harmonics.
The detiding improves with time. A restart mechanism is provided by writing time-accumulated harmonics into input tidal forcing NetCDF file. More information about detinding will be provided later in wikiROMS.
Many thanks to John Wilkin and Gordon Zhang for discussions about how to implement this option.