160 likes | 229 Views
TOMS Inaugural Meeting, NCAR, Aug. 23, 2001. Model Comparison and Evaluation Studies. Tal Ezer (with H. Arango & A. Shchepetkin). Goal: to evaluate and test new numerical elements in terrain-following ocean models.
E N D
TOMS Inaugural Meeting, NCAR, Aug. 23, 2001 Model Comparison and Evaluation Studies Tal Ezer (with H. Arango & A. Shchepetkin)
Goal: to evaluate and test new numerical elements in terrain-following ocean models. • Why: to help developers in improving codes and users in selecting options and parameterizations. • How: by comparing different schemes and parameterizations in POM and ROMS/TOMS. • Criteria: numerical errors, accuracy, stability, compatibility between schemes, computational cost.
Tested Elements • Advection schemes • Time stepping algorithms • Pressure gradient schemes • Others
Seamount configuration test Very steep case h=4050m, w=50km s=0.36, r=14.2 moderately steep h=2700m, w=100km s=0.07, r=2.7 S=max(DH/2H) r=max(sDH/HDs) grid=(64x64x20), Dx=8km
Zonal flow -50cm 0 +50cm Topography Sea surface height
Effect of advection scheme on model: Surface elevation anomaly ADV2- 2nd ord. cnt. POM ROMS ADV4- 4th ord. cnt. ADV3- 3rd ord. upst. ROMS ROMS
Time-stepping schemes(split mode: baroclinic/internal and barotropic/external)
Coupling of barotropic (external) and baroclinic (internal) modes in ROMS weights Un+1 = SamUm Un+½ = S bm’Um’ DTE 1<m<N, N=DTI/DTE DTI
Sensitivity to internal (DTI) & external (DTE) time steps ROMS UNSTABLE STABLE TDI/DTE POM CFL=13s
The adjustment process in POM and ROMS: forced case (zonal flow)
Roms- sensitivity of adjustment procesto time step choices DTE=12s DTE=24s DTI=360s DTI=720s
Structure ofV (cm/s) in ROMS for different PG schemes (medium seamount case) R-DJ (Vmax=3.7) R-PJQ (Vmax=0.03) R-FPJ (Vmax=30) R-DJC (Vmax=0.06) R-WDJ (Vmax=0.3)
More model elements and schemes need to be tested • Boundary conditions • Data assimilation • Nesting and air-sea coupling • Horizontal and vertical mixing schemes • Generalized coordinate system: comparisons between sigma, z-level & isopycnal grids.