I am running the tides only model and keeping constant temp. But I have got the surface temperature as follow. I also found the temperature of all boundaries in the bottom layers are not constant. What's the reason for this?
Please find the attachment for my model's configure.
how to keep constant temp when running tide only modal?
-
xiaozhu557
- Posts: 62
- Joined: Fri Sep 11, 2009 1:48 pm
- Location: nmefc
how to keep constant temp when running tide only modal?
- Attachments
-
- Temperature on the bottom layer
-
- SST
-
xiaozhu557
- Posts: 62
- Joined: Fri Sep 11, 2009 1:48 pm
- Location: nmefc
Re: how to keep constant temp when running tide only modal?
Though I have not defined TS_FIXED, I have set the initial conditions, OBC to T0(ng), and set the surface and bottom fluxes to 0, by ana_initial, ana_tobc.h, ana_stflux.h and ana_btflux.h.
So I really can not understand why the temperature can be changed still during the LOOP?
So I really can not understand why the temperature can be changed still during the LOOP?
Re: how to keep constant temp when running tide only modal?
You don't need to define TS_FIXED for tide-only modeling. Just constant tempe/salt as you specified.
It looks to me that your problem may be something to do with open boundaries.
What's your cppdefs and ocean.in?
Also, why did you edit ana_tobc.h? I don't think I did that for tide-only modeling.
Cheers,
DJ@TAMU
It looks to me that your problem may be something to do with open boundaries.
What's your cppdefs and ocean.in?
Also, why did you edit ana_tobc.h? I don't think I did that for tide-only modeling.
Cheers,
DJ@TAMU
Re: how to keep constant temp when running tide only modal?
For a really simple tides-only problem you can turn off SOLVE3D. Then you won't have any temperature to be evolving. The one thing you would have to change in ocean.in is to set DT to your fast DT and NDTFAST to 1.
-
xiaozhu557
- Posts: 62
- Joined: Fri Sep 11, 2009 1:48 pm
- Location: nmefc
Re: how to keep constant temp when running tide only modal?
My cppdefs are including but not limited as follow,d.kobashi wrote:You don't need to define TS_FIXED for tide-only modeling. Just constant tempe/salt as you specified.
It looks to me that your problem may be something to do with open boundaries.
What's your cppdefs and ocean.in?
Also, why did you edit ana_tobc.h? I don't think I did that for tide-only modeling.
Cheers,
DJ@TAMU
UV_* ANA_BSFLUX ANA_BTFLUX VISC_GRID MIX_S_UV ANA_FSOBC ANA_M2OBC ANA_M3OBC ANA_TOBC ANA_INITIAL ANA_STFLUX ANA_SMFLUX, RADIATION_2D SPLINES SOLVE3D CURVGRID MASKING DJ_GRADPS.
For other cppdefs, you can check the attached file "batch-41443.out", as follow,
Code: Select all
xcn069
Model Input Parameters: ROMS/TOMS version 3.7
Monday - February 9, 2015 - 4:17:44 PM
-----------------------------------------------------------------------------
The new SCS Operational Forecasting System, 1/30 degree resolution
Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /export/home/zhuxm/mpich-3.1/bin/mpif90
Compiler flags : -heap-arrays -fp-model precise -ip -O3 -free -free -free
Input Script :
SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : exported
Local Root : /mnt/swgfs/zhuxm/roms/ROMS_trunk_20141023
Header Dir : /export/home/zhuxm/swgfs/roms/ofs_scs
Header file : scs.h
Analytical Dir: /export/home/zhuxm/swgfs/roms/ofs_scs
Resolution, Grid 01: 0985x0793x030, Parallel Nodes: 128, Tiling: 016x008
Physical Parameters, Grid: 01
=============================
10000000 ntimes Number of timesteps for 3-D equations.
90.000 dt Timestep size (s) for 3-D equations.
30 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.
480 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
40 nHIS Number of timesteps between the writing fields
into history file.
172800 ndefHIS Number of timesteps between creation of new
history files.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
0 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
172800 ndefAVG Number of timesteps between creation of new
time-averaged file.
1.0000E+01 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
T LuvSponge Turning ON sponge on horizontal momentum.
F LtracerSponge(01) Turning OFF sponge on tracer 01: temp
F LtracerSponge(02) Turning OFF sponge on tracer 02: salt
1.0000E-05 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-05 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.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.0000E-03 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
1.0000E-03 Zob Bottom roughness (m).
2 Vtransform S-coordinate transformation equation.
4 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.
30.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1027.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).
20000101.00 time_ref Reference time for units attribute (yyyymmdd.dd)
3.6000E+02 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
3.6000E+02 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
3.6000E+02 Znudg Nudging/relaxation time scale (days)
for free-surface.
3.6000E+02 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
3.6000E+02 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.2000E+02 obcfac Factor between passive and active
open boundary conditions.
T VolCons(1) NLM western edge boundary volume conservation.
T VolCons(2) NLM southern edge boundary volume conservation.
T VolCons(3) NLM eastern edge boundary volume conservation.
T VolCons(4) NLM northern edge boundary volume conservation.
14.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1027.000 R0 Background density (kg/m3) used in linear Equation
of State.
1.7000E-04 Tcoef Thermal expansion coefficient (1/Celsius).
0.0000E+00 Scoef Saline contraction coefficient (1/PSU).
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
F LtracerSrc(01) Turning OFF point Sources/Sinks on tracer 01: temp
F LtracerSrc(02) Turning OFF point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
F LtracerCLM(01) Turning OFF processing of climatology tracer 01: temp
F LtracerCLM(02) Turning OFF processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
F LnudgeTCLM(01) Turning OFF nudging of climatology tracer 01: temp
F LnudgeTCLM(02) Turning OFF nudging of climatology tracer 02: salt
T Hout(idFsur) Write out free-surface.
T Aout(idFsur) Write out averaged free-surface.
T Aout(idUbar) Write out averaged 2D U-momentum component.
T Aout(idVbar) Write out averaged 2D V-momentum component.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idWvel) Write out averaged W-momentum component.
T Aout(idTvar) Write out averaged tracer 01: temp
T Aout(idTvar) Write out averaged tracer 02: salt
Output/Input Files:
Output Restart File: scs_rst.nc
Prefix for History Files: scs_his0
Prefix for Averages Files: scs_avg
Input Grid File: scs_grd.nc
Input Forcing File 01: scs_tide_20000101.nc
Tile partition information for Grid 01: 0985x0793x0030 tiling: 016x008
tile Istr Iend Jstr Jend Npts
Tile minimum and maximum fractional coordinates for Grid 01:
(interior points only)
tile Xmin Xmax Ymin Ymax grid
RHO-points
U-points
V-points
Maximum halo size in XI and ETA directions:
HaloSizeI(1) = 148
HaloSizeJ(1) = 224
TileSide(1) = 106
TileSize(1) = 7208
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Chapman Imp Chapman Imp Chapman Imp Chapman Imp
ubar 1 Flather Flather Flather Flather
vbar 1 Flather Flather Flather Flather
u 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud
v 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud
temp 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud
salt 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud
tke 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud
Activated C-preprocessing Options:
SCS The new SCS Operational Forecasting System, 1/30 degree resolution
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_FSOBC Analytical free-surface boundary conditions.
ANA_INITIAL Analytical initial conditions.
ANA_M2OBC Analytical 2D momentum boundary conditions.
ANA_M3OBC Analytical 3D momentum boundary conditions.
ANA_SMFLUX Analytical kinematic surface momentum flux.
ANA_STFLUX Analytical kinematic surface temperature flux.
ANA_TOBC Analytical tracers boundary conditions.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
INLINE_2DIO Processing 3D IO level by level to reduce memory needs.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
MY25_MIXING Mellor/Yamada Level-2.5 mixing closure.
NONLINEAR Nonlinear Model.
!NONLIN_EOS Linear Equation of State for seawater.
NO_WRITE_GRID Not Writing grid arrays into NetCDF ouput files.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RADIATION_2D Use tangential phase speed in radiation conditions.
RAMP_TIDES Ramping tidal forcing for one day.
!RST_SINGLE Double precision fields in restart NetCDF file.
SOLVE3D Solving 3D Primitive Equations.
SPLINES Conservative parabolic spline reconstruction.
SSH_TIDES Add tidal elevation to SSH climatology.
TIDES_ASTRO Compute astronomical (long-period tide) corrections.
TS_C4HADVECTION Fourth-order centered horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection 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_QDRAG Quadratic bottom stress.
UV_TIDES Add tidal currents to 2D momentum climatologies.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
VISC_GRID Horizontal viscosity coefficient scaled by grid size.
Process Information:
INITIAL: Configuring and initializing forward nonlinear model ...
*******
Vertical S-coordinate System, Grid 01:
level S-coord Cs-curve Z at hmin at hc half way at hmax
30 0.0000000 0.0000000 0.000 0.000 0.000 0.000
29 -0.0333333 -0.0002307 -0.143 -0.503 -1.676 -2.372
28 -0.0666667 -0.0009291 -0.286 -1.014 -4.742 -7.537
27 -0.1000000 -0.0021144 -0.430 -1.532 -9.256 -15.608
26 -0.1333333 -0.0038189 -0.574 -2.057 -15.313 -26.779
25 -0.1666667 -0.0060892 -0.719 -2.591 -23.052 -41.328
24 -0.2000000 -0.0089871 -0.864 -3.135 -32.657 -59.624
23 -0.2333333 -0.0125913 -1.009 -3.689 -44.362 -82.137
22 -0.2666667 -0.0169996 -1.155 -4.255 -58.457 -109.450
21 -0.3000000 -0.0223307 -1.302 -4.835 -75.296 -142.273
20 -0.3333333 -0.0287278 -1.449 -5.431 -95.303 -181.459
19 -0.3666667 -0.0363614 -1.597 -6.045 -118.986 -228.028
18 -0.4000000 -0.0454338 -1.747 -6.682 -146.947 -283.187
17 -0.4333333 -0.0561830 -1.897 -7.343 -179.892 -348.356
16 -0.4666667 -0.0688882 -2.049 -8.033 -218.652 -425.203
15 -0.5000000 -0.0838752 -2.203 -8.758 -264.196 -515.673
14 -0.5333333 -0.1015226 -2.358 -9.523 -317.648 -622.026
13 -0.5666667 -0.1222681 -2.516 -10.334 -380.310 -746.874
12 -0.6000000 -0.1466152 -2.676 -11.199 -453.679 -893.225
11 -0.6333333 -0.1751391 -2.839 -12.127 -539.465 -1064.512
10 -0.6666667 -0.2084930 -3.006 -13.127 -639.610 -1264.635
9 -0.7000000 -0.2474113 -3.177 -14.211 -756.296 -1497.978
8 -0.7333333 -0.2927113 -3.352 -15.391 -891.953 -1769.421
7 -0.7666667 -0.3452901 -3.532 -16.679 -1049.250 -2084.319
6 -0.8000000 -0.4061157 -3.719 -18.092 -1231.062 -2448.452
5 -0.8333333 -0.4762095 -3.912 -19.643 -1440.427 -2867.917
4 -0.8666667 -0.5566161 -4.112 -21.349 -1680.449 -3348.952
3 -0.9000000 -0.6483588 -4.320 -23.225 -1954.172 -3897.665
2 -0.9333333 -0.7523753 -4.537 -25.286 -2264.382 -4519.654
1 -0.9666667 -0.8694298 -4.764 -27.541 -2613.351 -5219.481
0 -1.0000000 -1.0000000 -5.000 -30.000 -3002.500 -6000.000
Time Splitting Weights for Grid 01: ndtfast = 30 nfast = 42
==================================
Primary Secondary Accumulated to Current Step
1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
3-0.0017877524645903 0.0334071600137066-0.0040025528757869 0.1001008081383191
4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463536
7-0.0011849025289723 0.0336508794757796-0.0107112868023600 0.2343461315221331
8-0.0005032751608631 0.0336903762267453-0.0112145619632232 0.2680365077488784
9 0.0003887272597151 0.0337071520654408-0.0108258347035081 0.3017436598143192
10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
12 0.0043042707117221 0.0335513677379153-0.0022367614257356 0.4026337758325831
13 0.0060106451121704 0.0334078920475245 0.0037738836864348 0.4360416678801076
14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
15 0.0099910761708920 0.0329439123123590 0.0216737068001213 0.5021931174029188
16 0.0122483446563884 0.0326108764399960 0.0339220514565097 0.5348039938429148
17 0.0146692120341107 0.0322025982847830 0.0485912634906204 0.5670065921276978
18 0.0172400033810439 0.0317136245503127 0.0658312668716643 0.5987202166780105
19 0.0199444086685725 0.0311389577709445 0.0857756755402368 0.6298591744489550
20 0.0227631639997064 0.0304741441486588 0.1085388395399432 0.6603333185976138
21 0.0256737146312911 0.0297153720153352 0.1342125541712342 0.6900486906129490
22 0.0286498597812016 0.0288595815276255 0.1628624139524359 0.7189082721405746
23 0.0316613792205220 0.0279045862015855 0.1945237931729578 0.7468128583421600
24 0.0346736416507075 0.0268492068942347 0.2291974348236653 0.7736620652363948
25 0.0376471948657328 0.0256934188392112 0.2668446296893981 0.7993554840756060
26 0.0405373376992232 0.0244385123436867 0.3073819673886213 0.8237939964192927
27 0.0432936737565710 0.0230872677537126 0.3506756411451924 0.8468812641730054
28 0.0458596469320356 0.0216441452951603 0.3965352880772280 0.8685254094681656
29 0.0481720587108284 0.0201154903974257 0.4447073467880565 0.8886408998655914
30 0.0501605672561820 0.0185097551070648 0.4948679140442384 0.9071506549726561
31 0.0517471682814030 0.0168377361985254 0.5466150823256415 0.9239883911711815
32 0.0528456577069106 0.0151128305891453 0.5994607400325521 0.9391012217603267
33 0.0533610761022577 0.0133513086655816 0.6528218161348098 0.9524525304259084
34 0.0531891349131379 0.0115726061288397 0.7060109510479478 0.9640251365547481
35 0.0522156244733761 0.0097996349650684 0.7582265755213239 0.9738247715198165
36 0.0503158038019030 0.0080591141492892 0.8085423793232269 0.9818838856691057
37 0.0473537721847153 0.0063819206892258 0.8558961515079423 0.9882658063583315
38 0.0431818225418188 0.0048034616164019 0.8990779740497611 0.9930692679747334
39 0.0376397765791564 0.0033640675316746 0.9367177506289175 0.9964333355064080
40 0.0305543017255206 0.0021094083123694 0.9672720523544381 0.9985427438187774
41 0.0217382098544504 0.0010909315881854 0.9890102622088885 0.9996336754069628
42 0.0109897377911118 0.0003663245930371 1.0000000000000004 0.9999999999999999
ndtfast, nfast = 30 42 nfast/ndtfast = 1.40000
Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):
1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000
Power filter parameters, Fgamma, gamma = 0.28400 0.18933
Metrics information for Grid 01:
===============================
Minimum X-grid spacing, DXmin = 3.23632803E+00 km
Maximum X-grid spacing, DXmax = 9.17057348E+00 km
Minimum Y-grid spacing, DYmin = 3.59466836E+00 km
Maximum Y-grid spacing, DYmax = 8.96506546E+00 km
Minimum Z-grid spacing, DZmin = 1.43021931E-01 m
Maximum Z-grid spacing, DZmax = 7.80518607E+02 m
Minimum barotropic Courant Number = 3.73923787E-03
Maximum barotropic Courant Number = 2.81036039E-01
Maximum Coriolis Courant Number = 6.21345144E-03
Horizontal mixing scaled by grid size, GRDMAX = 9.05374848E+00 km
GET_NGFLD - tidal period
(Min = 4.30819210E+04 Max = 9.67262421E+04)
GET_2DFLD - tidal elevation amplitude
(Min = 2.52492159E-05 Max = 1.81496844E+01)
GET_2DFLD - tidal elevation phase angle
(Min = 7.93344988E-06 Max = 6.28318465E+00)
GET_2DFLD - tidal current inclination angle
(Min = 3.70293851E-07 Max = 6.28318519E+00)
GET_2DFLD - tidal current phase angle
(Min = 1.48087661E-06 Max = 6.28318352E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 1.46277260E-06 Max = 3.88200810E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -4.84188479E-01 Max = 4.11129338E-01)
Basin information for Grid 01:
Maximum grid stiffness ratios: rx0 = 4.152431E-01 (Beckmann and Haidvogel)
rx1 = 8.935458E+00 (Haney)
Initial basin volumes: TotVolume = 3.7659566873E+16 m3
MinVolume = 1.7915665593E+06 m3
MaxVolume = 5.3941176727E+10 m3
Max/Min = 3.0108385561E+04
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 10000000)
STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed
0 0 00:00:00 0.000000E+00 2.142222E+04 2.142222E+04 3.772388E+16
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: scs_his0_0001.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
1 0 00:01:30 1.658685E-30 2.142222E+04 2.142222E+04 3.772388E+16
(666,679,01) 1.379357E-16 1.050798E-16 0.000000E+00 1.323787E-14
Thanks for your suggestions.
- Attachments
-
ana_tobc.h- (11.45 KiB) Downloaded 298 times
-
xiaozhu557
- Posts: 62
- Joined: Fri Sep 11, 2009 1:48 pm
- Location: nmefc
Re: how to keep constant temp when running tide only modal?
I just want to run 3D model, not 2D, though I know 2D is enough for tide only case.kate wrote:For a really simple tides-only problem you can turn off SOLVE3D. Then you won't have any temperature to be evolving. The one thing you would have to change in ocean.in is to set DT to your fast DT and NDTFAST to 1.
Thanks. Any other suggestions still?