When i run a model, I don't know Why are so much circulation?
Who can help me?
Code: Select all
#define UV_ADV /* advection terms */
#define UV_COR /* Coriolis terms */
#define UV_VIS2 /* harmonic horizontal mixing */
#define UV_LDRAG /* linear bottom friction */
#undef UV_QDRAG /* Use quadratic bottom drag law...but no waves*/
/* tracers equations */
#define TS_MPDATA /* recursive MPDATA 3D advection */
/* need a positive-definate advection scheme for sediment */
#define TS_DIF2 /* harmonic horizontal mixing */
#define SALINITY
#define NONLIN_EOS
/* pressure gradient algorithm */
#define DJ_GRADPS /* splines density Jacobian (Shchepetkin, 2000) */
/* model configuration */
#define SOLVE3D
#define CURVGRID
#define MASKING
#define SPLINES /* parabolic splines reconstruction of vertical derivatives */
#define SPHERICAL /* analytical spherical grid */
#undef AVERAGES
#undef DIAGNOSTICS_TS /* writing out tracer diagnostics */
#undef DIAGNOSTICS_UV /* writing out momentum diagnostics */
#define STATIONS /* writing out station data to file */
/* horizontal mixing of momentum*/
#define MIX_S_UV /* mixing along constant S-surfaces */
/* horizontal mixing of tracers*/
#define MIX_GEO_TS /* mixing along geopotential (constant depth) surfaces */
/* NetCDF input and output*/
#define OUT_DOUBLE /* Double precision output fields in NetCDF files*/
/* boundary conditions*/
#define NORTHERN_WALL
#define SOUTHERN_WALL
#define WESTERN_WALL
#undef EW_PERIODIC
#define ANA_M2OBC
#define ANA_FSOBC
#define ADD_FSOBC
#define ADD_M2OBC
#define EAST_M3RADIATION
#define EAST_TGRADIENT
/* define UV_TIDES to specify boundary tidal currents from data e.g. ADCIRC */
/* define UV_TIDES to compute tidal currents using reduced physics */
#define RAMP_TIDES /* ramping (over one day) tidal forcing */
#define SSH_TIDES /* read tidal SSH from file */
#define EAST_FSCHAPMAN
#define UV_TIDES /* read tidal U and V from file */
#define EAST_M2FLATHER
#undef UV_PSOURCE
#undef Q_PSOURCE
#undef TS_PSOURCE
/* sediment transport model */
#undef SEDIMENT
#ifdef SEDIMENT
# define SED_MORPH
# define SED_DENS
# define SUSPLOAD /* suspended load transport */
# undef BEDLOAD_SOULSBY
# define BEDLOAD_MPM
#endif
#ifdef SEDIMENT
#undef ANA_SEDIMENT
#endif
/*bottom boundary layer closure */
#undef MB_BBL
#undef SG_BBL
#undef SSW_BBL
/* vertical mixing of momentum and tracers */
#define MY25_MIXING
#ifdef MY25_MIXING
# define KANTHA_CLAYSON
# undef CANUTO_A
# define N2S2_HORAVG
#endif
/* */
#define ANA_SMFLUX
#define ANA_STFLUX
#define ANA_SPFLUX
#define ANA_SSFLUX /* analytical surface salt flux defaults to zero */
#define ANA_BPFLUX /*bottom passive tracers fluxes */
#define ANA_BSFLUX /* bottom temp flux...defaults to zero */
Code: Select all
Process Information:
Thread # 0 (pid= 0) is active.
Model Input Parameters: ROMS/TOMS version 3.4
Tuesday - January 18, 2011 - 10:44:20 AM
-----------------------------------------------------------------------------
ZHANJIANG Model
Operating system : CYGWIN
CPU/hardware : i686
Compiler system : gfortran
Compiler command : /usr/bin/gfortran-4
Compiler flags : -frepack-arrays -O3 -ffast-math -ffree-form -ffree-line-length-none
SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : 510:527M
Local Root : /cygdrive/e/ROMS/CASE/NEW/trunk
Header Dir : /cygdrive/e/ROMS/CASE/NEW/trunk/ROMS/Include
Header file : zhanjiang_1225.h
Analytical Dir: /cygdrive/e/ROMS/CASE/NEW/trunk/ROMS/Functionals
Resolution, Grid 01: 0218x0238x016, Parallel Threads: 1, Tiling: 001x001
Physical Parameters, Grid: 01
=============================
28800 ntimes Number of timesteps for 3-D equations.
5.000 dt Timestep size (s) for 3-D equations.
10 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
1 ERstr Starting ensemble/perturbation run number.
1 ERend Ending ensemble/perturbation run number.
0 nrrec Number of restart records to read from disk.
T LcycleRST Switch to recycle time-records in restart file.
360 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
360 nSTA Number of timesteps between the writing of data
the stations file.
T ldefout Switch to create a new output NetCDF file(s).
360 nHIS Number of timesteps between the writing fields
into history file.
2.0000E+01 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
2.0000E+01 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
1.0000E+01 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
1.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
1.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
1.0000E-05 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
1.0000E-05 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
5.0000E+00 theta_s S-coordinate surface control parameter.
4.0000E-01 theta_b S-coordinate bottom control parameter.
1.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.000 dstart Time-stamp assigned to model initialization (days).
0.000 tide_start Reference time origin for tidal forcing (days).
20101105.00 time_ref Reference time for units attribute (yyyymmdd.dd)
2.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
2.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
2.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
2.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
2.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
2.0000E+01 obcfac Factor between passive and active
open boundary conditions.
12.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
T Hout(idFsur) Write out free-surface.
T Hout(idUbar) Write out 2D U-momentum component.
T Hout(idVbar) Write out 2D V-momentum component.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idWvel) Write out W-momentum component.
T Hout(idOvel) Write out omega vertical velocity.
T Hout(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt
Output/Input Files:
Output Restart File: ocean_rst.nc
Output History File: ocean_his.nc
Output Stations File: ocean_sta.nc
Input Grid File: Data/ROMS/Grid/zhanjiang_grd.nc
Input Nonlinear Initial File: Data/ROMS/Initial/zhanjiang_ini.nc
Station positions File: ROMS/External/stations_zhanjiang.in
Tile partition information for Grid 01: 0218x0238x0016 tiling: 001x001
tile Istr Iend Jstr Jend Npts
0 1 218 1 238 830144
Tile minimum and maximum fractional grid coordinates:
(interior points only)
tile Xmin Xmax Ymin Ymax grid
0 -0.50 220.50 -0.50 240.50 RHO-points
0 -0.50 220.50 -0.50 240.50 U-points
0 -0.50 220.50 -0.50 240.50 V-points
Stations Parameters, Grid: 01
=============================
7 Nstation Number of stations to write out into stations file.
T Sout(idFsur) Write out free-surface.
T Sout(idUbar) Write out 2D U-momentum component.
T Sout(idVbar) Write out 2D V-momentum component.
T Sout(idUvel) Write out 3D U-momentum component.
T Sout(idVvel) Write out 3D V-momentum component.
T Sout(idWvel) Write out W-momentum component.
T Sout(idOvel) Write out omega vertical velocity.
T Sout(idTvar) Write out tracer 01: temp
T Sout(idTvar) Write out tracer 02: salt
Flag and positions for station 0001: 0 78.0000 143.0000
Flag and positions for station 0002: 0 67.0000 133.0000
Flag and positions for station 0003: 0 71.0000 144.0000
Flag and positions for station 0004: 0 81.0000 135.0000
Flag and positions for station 0005: 0 84.0000 145.0000
Flag and positions for station 0006: 0 104.0000 145.0000
Flag and positions for station 0007: 0 105.0000 145.0000
Activated C-preprocessing Options:
ZHANJIANG_1225 ZHANJIANG Model
ADD_FSOBC Adding tidal elevation to proccesed OBC data.
ADD_M2OBC Adding tidal currents to proccesed OBC data.
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_FSOBC Analytical free-surface boundary conditions.
ANA_M2OBC Analytical 2D momentum boundary conditions.
ANA_SMFLUX Analytical kinematic surface momentum flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ANA_STFLUX Analytical kinematic surface temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EAST_FSCHAPMAN Eastern edge, free-surface, Chapman condition.
EAST_M2FLATHER Eastern edge, 2D momentum, Flather condition.
EAST_M3RADIATION Eastern edge, 3D momentum, radiation condition.
EAST_TGRADIENT Eastern edge, tracers, gradient condition.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MY25_MIXING Mellor/Yamada Level-2.5 mixing closure.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
NORTHERN_WALL Wall boundary at Northern edge.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
OUT_DOUBLE Double precision output fields in NetCDF files.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RAMP_TIDES Ramping tidal forcing for one day.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SOUTHERN_WALL Wall boundary at Southern edge.
SPLINES Conservative parabolic spline reconstruction.
SPHERICAL Spherical grid configuration.
SSH_TIDES Add tidal elevation to SSH climatology.
STATIONS Writing out station data.
TS_MPDATA Recursive flux corrected MPDATA 3D advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_LDRAG Linear bottom stress.
UV_TIDES Add tidal currents to 2D momentum climatologies.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WESTERN_WALL Wall boundary at Western edge.
INITIAL: Configuring and initializing forward nonlinear model ...
Vertical S-coordinate System:
level S-coord Cs-curve at_hmin over_slope at_hmax
16 0.0000000 -0.0000000 0.000 -0.000 -0.000
15 -0.0625000 -0.0048948 -0.063 -0.119 -0.175
14 -0.1250000 -0.0119914 -0.125 -0.263 -0.401
13 -0.1875000 -0.0230908 -0.188 -0.453 -0.719
12 -0.2500000 -0.0409944 -0.250 -0.722 -1.194
11 -0.3125000 -0.0696342 -0.313 -1.114 -1.916
10 -0.3750000 -0.1133185 -0.375 -1.680 -2.984
9 -0.4375000 -0.1742175 -0.438 -2.443 -4.449
8 -0.5000000 -0.2489214 -0.500 -3.366 -6.232
7 -0.5625000 -0.3284417 -0.563 -4.344 -8.125
6 -0.6250000 -0.4042644 -0.625 -5.279 -9.933
5 -0.6875000 -0.4744488 -0.688 -6.150 -11.612
4 -0.7500000 -0.5437741 -0.750 -7.010 -13.271
3 -0.8125000 -0.6205544 -0.813 -7.957 -15.101
2 -0.8750000 -0.7145184 -0.875 -9.101 -17.327
1 -0.9375000 -0.8366254 -0.938 -10.569 -20.201
0 -1.0000000 -1.0000000 -1.000 -12.513 -24.025
Time Splitting Weights: ndtfast = 10 nfast = 15
Primary Secondary Accumulated to Current Step
1 0.0014850648765542 0.1000000000000000 0.0014850648765542 0.1000000000000000
2 0.0059386031399402 0.0998514935123446 0.0074236680164944 0.1998514935123445
3 0.0133457074936657 0.0992576331983506 0.0207693755101601 0.2991091267106951
4 0.0236484051180389 0.0979230624489840 0.0444177806281990 0.3970321891596791
5 0.0366959666818608 0.0955582219371801 0.0811137473100598 0.4925904110968592
6 0.0521753389587947 0.0918886252689940 0.1332890862688545 0.5844790363658532
7 0.0695217010484137 0.0866710913731145 0.2028107873172682 0.6711501277389678
8 0.0878091442019233 0.0797189212682732 0.2906199315191915 0.7508690490072409
9 0.1056214752525633 0.0709380068480808 0.3962414067717548 0.8218070558553218
10 0.1209031436506851 0.0603758593228245 0.5171445504224399 0.8821829151781463
11 0.1307902921035067 0.0482855449577560 0.6479348425259466 0.9304684601359023
12 0.1314219308195445 0.0352065157474053 0.7793567733454911 0.9656749758833076
13 0.1177312353577226 0.0220643226654509 0.8970880087032138 0.9877392985487585
14 0.0832169680811588 0.0102911991296786 0.9803049767843727 0.9980304976784371
15 0.0196950232156274 0.0019695023215627 1.0000000000000000 0.9999999999999999
ndtfast, nfast = 10 15 nfast/ndtfast = 1.50000
Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):
1.000000000000 1.083803072626 0.541901536313 1.000000000000 1.000000000000
Power filter parameters, Fgamma, gamma = 0.28400 0.00000
Minimum X-grid spacing, DXmin = 1.93639455E-01 km
Maximum X-grid spacing, DXmax = 6.04239487E-01 km
Minimum Y-grid spacing, DYmin = 2.68942210E-01 km
Maximum Y-grid spacing, DYmax = 7.28061065E-01 km
Minimum Z-grid spacing, DZmin = 6.25000000E-02 m
Maximum Z-grid spacing, DZmax = 3.82426846E+00 m
Minimum barotropic Courant Number = 1.22495929E-02
Maximum barotropic Courant Number = 4.45839358E-02
Maximum Coriolis Courant Number = 5.32083094E-04
NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(File: zhanjiang_ini.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 2.24573922E+01 Max = 2.27319504E+01)
- salinity
(Min = 3.41531521E+01 Max = 3.41673488E+01)
GET_NGFLD - tidal period
(Min = 4.30819210E+04 Max = 9.67262421E+04)
GET_2DFLD - tidal elevation amplitude
(Min = 5.73345254E-02 Max = 5.85054670E-01)
GET_2DFLD - tidal elevation phase angle
(Min = 1.36518644E+00 Max = 5.21264685E+00)
GET_2DFLD - tidal current inclination angle
(Min = 1.11310007E-05 Max = 6.28318263E+00)
GET_2DFLD - tidal current phase angle
(Min = 1.71641227E-03 Max = 6.28239915E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 9.59182533E-03 Max = 1.28520087E-01)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -1.56127638E-02 Max = 7.50047868E-03)
Maximum grid stiffness ratios: rx0 = 2.088082E-01 (Beckmann and Haidvogel)
rx1 = 3.468923E+00 (Haney)
Initial basin volumes: TotVolume = 3.4728667676E+10 m3
MinVolume = 3.8897822161E+03 m3
MaxVolume = 1.6160127023E+06 m3
Max/Min = 4.1545068915E+02
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00028800)
STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
0 0 00:00:00 0.000000E+00 5.585238E+01 5.585238E+01 3.817329E+10
DEF_HIS - creating history file: ocean_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_STATION - creating stations file: ocean_sta.nc
1 0 00:00:10 1.090756E-16 5.585238E+01 5.585238E+01 3.817329E+10
2 0 00:00:20 5.301220E-11 5.585239E+01 5.585239E+01 3.817329E+10
3 0 00:00:30 1.836297E-10 5.585240E+01 5.585240E+01 3.817330E+10
4 0 00:00:40 3.694276E-10 5.585241E+01 5.585241E+01 3.817330E+10
5 0 00:00:50 6.179113E-10 5.585243E+01 5.585243E+01 3.817331E+10
6 0 00:01:00 9.592319E-10 5.585244E+01 5.585244E+01 3.817331E+10
7 0 00:01:10 1.431474E-09 5.585247E+01 5.585247E+01 3.817332E+10
8 0 00:01:20 2.067357E-09 5.585249E+01 5.585249E+01 3.817333E+10
9 0 00:01:30 2.888218E-09 5.585252E+01 5.585252E+01 3.817334E+10
10 0 00:01:40 3.905421E-09 5.585255E+01 5.585255E+01 3.817335E+10
11 0 00:01:50 5.125385E-09 5.585258E+01 5.585258E+01 3.817336E+10
12 0 00:02:00 6.554275E-09 5.585261E+01 5.585261E+01 3.817337E+10
13 0 00:02:10 8.200629E-09 5.585265E+01 5.585265E+01 3.817338E+10
14 0 00:02:20 1.007634E-08 5.585269E+01 5.585269E+01 3.817340E+10
15 0 00:02:30 1.219682E-08 5.585273E+01 5.585273E+01 3.817341E+10
16 0 00:02:40 1.458062E-08 5.585277E+01 5.585277E+01 3.817343E+10
17 0 00:02:50 1.724820E-08 5.585282E+01 5.585282E+01 3.817345E+10
18 0 00:03:00 2.022000E-08 5.585287E+01 5.585287E+01 3.817346E+10
19 0 00:03:10 2.351452E-08 5.585292E+01 5.585292E+01 3.817348E+10
.......
ROMS/TOMS: DONE... Tuesday - January 18, 2011 - 8:56:44 PM
- topo.png (22.64 KiB) Viewed 2027 times
- zeta.png (22.99 KiB) Viewed 2027 times
- ubar.png (31.38 KiB) Viewed 2027 times
- vabr.png (31.6 KiB) Viewed 2027 times
