Sediment for beginner))

Discussion of how to use ROMS on different regional and basin scale applications.

Moderators: arango, robertson

Post Reply
Message
Author
fancer
Posts: 45
Joined: Fri Feb 06, 2009 5:20 pm
Location: NNSTU, The Nizhniy Novgorod State Technical University named after R.E.Alekseev, Russia

Sediment for beginner))

#1 Post by fancer » Wed Mar 25, 2009 8:08 pm

Hello, I attend sediment transport to a problem. I have already learnt to collect ROMS for the necessary problem, whether it be sediment transport or not, and to start processing. You could not help?
1) how to set a step on time and time?
2) what should be source files and that each of them means.
3) where I can find detailed description shoreface for ROMS?

Many thanks for the help!))))
PS Probably, it not the last my questions))))

jcwarner
Posts: 816
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Sediment for beginner))

#2 Post by jcwarner » Fri Mar 27, 2009 12:44 pm

A lot of information you seek is on the wiki site:

https://www.myroms.org/wiki/index.php/D ... ion_Portal

fancer
Posts: 45
Joined: Fri Feb 06, 2009 5:20 pm
Location: NNSTU, The Nizhniy Novgorod State Technical University named after R.E.Alekseev, Russia

Re: Sediment for beginner))

#3 Post by fancer » Fri Mar 27, 2009 2:37 pm

I'm so sorry for my english)))) I still study)))

I have read all in wiki-site, but I have not found something concrete about NetCDF input files and what should be in this files. For example, in ROMS/External/ocean_sed_test1.in:
...
! Input NetCDF file names, [1:Ngrids].

GRDNAME == ocean_grd.nc
ININAME == ocean_ini.nc
ITLNAME == ocean_itl.nc
IRPNAME == ocean_irp.nc
IADNAME == ocean_iad.nc
CLMNAME == ocean_clm.nc
BRYNAME == ocean_bry.nc
FWDNAME == ocean_fwd.nc
ADSNAME == ocean_ads.nc
...
But how ROMS calculate without this files?((( What ROMS use as defaults files instead files in this list?

And how I can assign bathymetry?

ps A step on time assign in ROMS/External/ocean_sed_test1.in)))

jcwarner
Posts: 816
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Sediment for beginner))

#4 Post by jcwarner » Fri Mar 27, 2009 2:41 pm

well, roms does not need all these files for all applications. It depends on what you are doing. So it typically will need an grid file (GRDNAME) and an init file (ININAME).
Grid files can be difficult to create, but there is info on the wiki to make these.

fancer
Posts: 45
Joined: Fri Feb 06, 2009 5:20 pm
Location: NNSTU, The Nizhniy Novgorod State Technical University named after R.E.Alekseev, Russia

Re: Sediment for beginner))

#5 Post by fancer » Mon Mar 30, 2009 11:51 am

Thank You very much for help)))
I run SHOREFACE application, but it crash with segmentation fault error. Step by step:

1) in makefie:
ROMS_APPLICATION ?= SHOREFACE
2) in cppdefs.h no one word about SHOREFACE (((
3) and use ROMS/External/ocean_shoreface.in:
oceanS < ROMS/External/ocean_shoreface.in >& ocean_shoreface.log &

Here log-file (Where is "attach file"-button??):

Process Information:

Thread # 0 (pid= 31509) is active.

Model Input Parameters: ROMS/TOMS version 3.2
Monday - March 30, 2009 - 4:01:16 PM
-----------------------------------------------------------------------------

Shore Face Planar Beach Test Case

Operating system : Linux
CPU/hardware : i686
Compiler system : gfortran
Compiler command : /usr/bin/gfortran
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 : 333M

Local Root : /home/fancer/ROMS/build/roms
Header Dir : /home/fancer/ROMS/build/roms/ROMS/Include
Header file : shoreface.h
Analytical Dir: /home/fancer/ROMS/build/roms/ROMS/Functionals

Resolution, Grid 01: 0059x0007x020, Parallel Threads: 1, Tiling: 001x001


Physical Parameters, Grid: 01
=============================

4320 ntimes Number of timesteps for 3-D equations.
5.000 dt Timestep size (s) for 3-D equations.
20 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.
720 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).
360 nHIS Number of timesteps between the writing fields
into history file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
720 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
1 ntsDIA Starting timestep for the accumulation of output
time-averaged diagnostics data.
720 nDIA Number of timesteps between the writing of
time-averaged data into diagnostics file.
0.0000E+00 visc2 Horizontal, harmonic mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
5.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
5.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.000 gls_p GLS stability exponent.
1.500 gls_m GLS turbulent kinetic energy exponent.
-1.000 gls_n GLS turbulent length scale exponent.
7.6000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
1.0000E-12 gls_Pmin GLS minimum value of dissipation.
5.4770E-01 gls_cmu0 GLS stability coefficient.
1.4400E+00 gls_c1 GLS shear production coefficient.
1.9200E+00 gls_c2 GLS dissipation coefficient.
-4.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
1.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
1.3000E+00 gls_sigp GLS constant Schmidt number for PSI.
1400.000 charnok_alpha Charnok factor for Zos calculation.
0.500 zos_hsig_alpha Factor for Zos calculation using Hsig(Awave).
0.250 sz_alpha Factor for Wave dissipation surface tke flux .
100.000 crgban_cw Factor for Craig/Banner surface tke flux.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
1.5000E-03 rdrg2 Quadratic bottom drag coefficient.
1.5000E-02 Zob Bottom roughness (m).
2.0000E-02 Zos Surface roughness (m).
1.0000E-01 Dcrit Minimum depth for wetting and drying (m).
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
0.0000E+00 theta_s S-coordinate surface control parameter.
0.0000E+00 theta_b S-coordinate bottom control parameter.
0.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.00 time_ref Reference time for units attribute (yyyymmdd.dd)
0.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
0.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
0.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
0.0000E+00 obcfac Factor between passive and active
open boundary conditions.
10.000 T0 Background potential temperature (C) constant.
30.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).
7.6000E-04 Scoef Saline contraction coefficient (1/PSU).
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
T Hout(idUbms) Write out bottom U-momentum stress.
T Hout(idVbms) Write out bottom V-momentum stress.
T Hout(idW2xx) Write out 2D radiation stress, Sxx.
T Hout(idW2xy) Write out 2D radiation stress, Sxy.
T Hout(idW2yy) Write out 2D radiation stress, Syy.
T Hout(idU2rs) Write out total 2D u-radiation stress.
T Hout(idV2rs) Write out total 2D v-radiation stress.
T Hout(idU2Sd) Write out 2D u-momentum stokes velocity.
T Hout(idV2Sd) Write out 2D v-momentum stokes velocity.
T Hout(idW3xx) Write out 3D horizonrtal radiation stress, Sxx.
T Hout(idW3xy) Write out 3D horizonrtal radiation stress, Sxy.
T Hout(idW3yy) Write out 3D horizonrtal radiation stress, Syy.
T Hout(idW3zx) Write out 3D vertical radiation stress, Szx.
T Hout(idW3zy) Write out 3D vertical radiation stress, Szy.
T Hout(idU3rs) Write out total 3D u-radiation stress.
T Hout(idV3rs) Write out total 3D v-radiation stress.
T Hout(idU3Sd) Write out 3D u-momentum stokes velocity.
T Hout(idV3Sd) Write out 3D v-momentum stokes velocity.
T Hout(idWamp) Write out wave height.
T Hout(idWlen) Write out wave length.
T Hout(idWdir) Write out wave direction.
T Hout(idBott) Write out bottom property 01: grain_diameter
T Hout(idBott) Write out bottom property 02: grain_density
T Hout(idBott) Write out bottom property 03: settling_vel
T Hout(idBott) Write out bottom property 04: erosion_stress
T Hout(idBott) Write out bottom property 05: ripple_length
T Hout(idBott) Write out bottom property 06: ripple_height
T Hout(idBott) Write out bottom property 07: bed_wave_amp
T Hout(idBott) Write out bottom property 08: Zo_def
T Hout(idBott) Write out bottom property 09: Zo_app
T Hout(idBott) Write out bottom property 15: active_layer_thickness
T Hout(idVvis) Write out vertical viscosity coefficient.
T Hout(idMtke) Write out turbulent kinetic energy.
T Hout(idMtls) Write out turbulent generic length-scale.

Output/Input Files:

Output Restart File: ocean_rst.nc
Output History File: ocean_his.nc
Output Averages File: ocean_avg.nc
Output Diagnostics File: ocean_dia.nc
Input Forcing File 01: Data/ROMS/Forcing/swan_shoreface_angle_forc.nc

Tile partition information for Grid 01: 0059x0007x0020 tiling: 001x001

tile Istr Iend Jstr Jend Npts

0 1 59 1 7 8260

Tile minimum and maximum fractional grid coordinates:
(interior points only)

tile Xmin Xmax Ymin Ymax grid

0 -0.50 60.50 -3.50 10.50 RHO-points

0 -0.50 60.50 -3.50 10.50 U-points

0 -0.50 60.50 -3.50 10.50 V-points


Sediment Parameters, Grid: 01
=============================


Size Sd50 Csed Srho Wsed Erate poros
Class (mm) (kg/m3) (kg/m3) (mm/s) (kg/m2/s) (nondim)

1 1.0000E+00 0.0000E+00 2.6500E+03 1.0000E+00 1.0000E-05 5.0000E-01

tau_ce tau_cd tnu2 tnu4 Akt_bak Tnudg
(N/m2) (N/m2) (m2/s) (m4/s) (m2/s) (day)

1 7.0000E-02 2.0000E-02 0.0000E+00 0.0000E+00 5.0000E-06 0.0000E+00

morph_fac
(nondim)

1 1.0000E+00

New bed layer formed when deposition exceeds 0.10000E-01 (m).
Two first layers are combined when 2nd layer smaller than 0.00000E+00 (m).
Rate coefficient for bed load transport = 0.50000E-01

T Hout(idTvar) Write out sediment01: sand_01
T Hout(idfrac) Write out bed fraction, sediment 01: sandfrac_01
T Hout(idfrac) Write out mass, sediment 01: sandmass_01
T Hout(idUbld) Write out bed load at U-points, sediment 01: bedload_Usand_01
T Hout(idVbld) Write out bed load at V-points, sediment 01: bedload_Vsand_01
T Hout(idSbed) Write out BED property 01: bed_thickness
T Hout(idSbed) Write out BED property 02: bed_age
T Hout(idSbed) Write out BED property 03: bed_porosity

Activated C-preprocessing Options:

SHOREFACE Shore Face Planar Beach Test Case
ANA_BPFLUX Analytical bottom passive tracers fluxes.
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_FSOBC Analytical free-surface boundary conditions.
ANA_GRID Analytical grid set-up.
ANA_INITIAL Analytical initial conditions.
ANA_M2OBC Analytical 2D momentum boundary conditions.
ANA_MASK Analytical Land/Sea Masking.
ANA_SEDIMENT Analytical sediment initial conditions.
ANA_SMFLUX Analytical kinematic surface momentum flux.
ANA_SPFLUX Analytical surface passive tracer fluxes.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ANA_STFLUX Analytical kinematic surface temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged fields.
AVERAGES_BELOAD Writing out time-averaged bed load.
AVERAGES_NEARSHORE Writing out time-averaged nearshore radiation fields.
BEDLOAD_MPM Activate bed load sediment transport Meyer-Peter-Mueller.
DIAGNOSTICS_UV Computing and writing momentum diagnostic terms.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EASTERN_WALL Wall boundary at Eastern edge.
GLS_MIXING Generic Length-Scale turbulence closure.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MIX_S_UV Mixing of momentum along constant S-surfaces.
NEARSHORE_MELLOR Nearshore RAdiation Stress Terms.
NONLINEAR Nonlinear Model.
!NONLIN_EOS Linear Equation of State for seawater.
NS_PERIODIC North-South periodic boundaries.
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.
!RST_SINGLE Double precision fields in restart NetCDF file.
SEDIMENT Cohesive and noncohesive sediments.
SUSPLOAD Activate suspended sediment transport.
SOLVE3D Solving 3D Primitive Equations.
SPLINES Conservative parabolic spline reconstruction.
THREE_GHOST Using three Ghost Points in halo regions.
TS_MPDATA Recursive flux corrected MPDATA 3D advection of tracers.
UV_ADV Advection of momentum.
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_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WEST_FSCHAPMAN Western edge, free-surface, Chapman condition.
WEST_M2FLATHER Western edge, 2D momentum, Flather condition.
WEST_M3GRADIENT Western edge, 3D momentum, gradient condition.
WET_DRY Wetting and drying activated.

INITIAL: Configuring and initializing forward nonlinear model ...


Vertical S-coordinate System:

level S-coord Cs-curve at_hmin over_slope at_hmax

20 0.0000000 0.0000000 0.000 0.000 0.000
19 -0.0500000 -0.0500000 0.150 -0.219 -0.588
18 -0.1000000 -0.1000000 0.300 -0.438 -1.175
17 -0.1500000 -0.1500000 0.450 -0.656 -1.763
16 -0.2000000 -0.2000000 0.600 -0.875 -2.350
15 -0.2500000 -0.2500000 0.750 -1.094 -2.938
14 -0.3000000 -0.3000000 0.900 -1.313 -3.525
13 -0.3500000 -0.3500000 1.050 -1.531 -4.113
12 -0.4000000 -0.4000000 1.200 -1.750 -4.700
11 -0.4500000 -0.4500000 1.350 -1.969 -5.288
10 -0.5000000 -0.5000000 1.500 -2.188 -5.875
9 -0.5500000 -0.5500000 1.650 -2.406 -6.463
8 -0.6000000 -0.6000000 1.800 -2.625 -7.050
7 -0.6500000 -0.6500000 1.950 -2.844 -7.638
6 -0.7000000 -0.7000000 2.100 -3.063 -8.225
5 -0.7500000 -0.7500000 2.250 -3.281 -8.813
4 -0.8000000 -0.8000000 2.400 -3.500 -9.400
3 -0.8500000 -0.8500000 2.550 -3.719 -9.988
2 -0.9000000 -0.9000000 2.700 -3.938 -10.575
1 -0.9500000 -0.9500000 2.850 -4.156 -11.163
0 -1.0000000 -1.0000000 3.000 -4.375 -11.750

Time Splitting Weights: ndtfast = 20 nfast = 29

Primary Secondary Accumulated to Current Step

1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
3-0.0011514592651645 0.0501156998674241-0.0034654566136461 0.1501639558342180
4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649003
5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792859
6 0.0029141799764308 0.0501428056118168 0.0000580677400948 0.3006719858911027
7 0.0054132615310267 0.0499970966129953 0.0054713292711215 0.3506690825040980
8 0.0084687837865132 0.0497264335364439 0.0139401130576347 0.4003955160405419
9 0.0120633394191050 0.0493029943471183 0.0260034524767397 0.4496985103876602
10 0.0161716623600090 0.0486998273761630 0.0421751148367487 0.4983983377638232
11 0.0207585511322367 0.0478912442581626 0.0629336659689854 0.5462895820219857
12 0.0257765478740990 0.0468533167015507 0.0887102138430844 0.5931428987235364
13 0.0311633730493853 0.0455644893078458 0.1198735868924698 0.6387073880313822
14 0.0368391158442262 0.0440063206553765 0.1567127027366960 0.6827137086867586
15 0.0427031802506397 0.0421643648631652 0.1994158829873357 0.7248780735499238
16 0.0486309868367616 0.0400292058506332 0.2480468698240973 0.7649072794005570
17 0.0544704302037592 0.0375976565087951 0.3025173000278565 0.8025049359093521
18 0.0600380921294286 0.0348741349986072 0.3625553921572851 0.8373790709079593
19 0.0651152103984763 0.0318722303921357 0.4276706025557614 0.8692513013000950
20 0.0694434033194840 0.0286164698722119 0.4971140058752453 0.8978677711723069
21 0.0727201499285570 0.0251442997062377 0.5698341558038024 0.9230120708785446
22 0.0745940258796570 0.0215082922098099 0.6444281816834594 0.9445203630883545
23 0.0746596950216180 0.0177785909158270 0.7190878767050775 0.9622989540041815
24 0.0724526566618460 0.0140456061647461 0.7915405333669234 0.9763445601689277
25 0.0674437485167025 0.0104229733316538 0.8589842818836260 0.9867675335005816
26 0.0590334053485720 0.0070507859058187 0.9180176872321980 0.9938183194064003
27 0.0465456732896125 0.0040991156383901 0.9645633605218105 0.9979174350447904
28 0.0292219798521905 0.0017718319739095 0.9937853403740009 0.9996892670186999
29 0.0062146596259991 0.0003107329813000 1.0000000000000000 0.9999999999999998

ndtfast, nfast = 20 29 nfast/ndtfast = 1.45000

Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):

1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.14200

Minimum X-grid spacing, DXmin = 2.00000000E-02 km
Maximum X-grid spacing, DXmax = 2.00000000E-02 km
Minimum Y-grid spacing, DYmin = 2.00000000E-02 km
Maximum Y-grid spacing, DYmax = 2.00000000E-02 km
Minimum Z-grid spacing, DZmin = -1.50000000E-01 m
Maximum Z-grid spacing, DZmax = 5.87500000E-01 m

Minimum barotropic Courant Number = 1.23806805E-02
Maximum barotropic Courant Number = 1.89792238E-01
Maximum Coriolis Courant Number = 0.00000000E+00

GET_2DFLD - wind-induced wave direction, t = 0 00:00:00
(Rec=0001, Index=1, File: swan_shoreface_angle_forc.nc)
(Tmin= 0.0000 Tmax= 0.0000)
(Min = 3.14159265E+00 Max = 4.73490388E+00)
GET_2DFLD - wind-induced significant wave height, t = 0 00:00:00
(Rec=0001, Index=1, File: swan_shoreface_angle_forc.nc)
(Tmin= 0.0000 Tmax= 0.0000)
(Min = 0.00000000E+00 Max = 2.20400000E+00)
GET_2DFLD - wind-induced mean wavelength, t = 0 00:00:00
(Rec=0001, Index=1, File: swan_shoreface_angle_forc.nc)
(Tmin= 0.0000 Tmax= 0.0000)
(Min = 1.00000000E+01 Max = 7.05299988E+01)

Maximum grid stiffness ratios: rx0 = 1.685714E+00 (Beckmann and Haidvogel)
rx1 = 6.574286E+01 (Haney)


Initial basin volumes: TotVolume = 7.22750000000000E+05 m3
MinVolume = -5.50833333333334E+01 m3
MaxVolume = 2.30083333333334E+02 m3
Max/Min = -4.17700453857791E+00

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00004320)

STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME

0 0 00:00:00 0.000000E+00 3.856303E+01 3.856303E+01 7.732200E+05

jcwarner
Posts: 816
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Sediment for beginner))

#6 Post by jcwarner » Mon Mar 30, 2009 3:20 pm

why did it crash?? you say it crashes with seg fault error, but i do not see that error in the log file.

fancer
Posts: 45
Joined: Fri Feb 06, 2009 5:20 pm
Location: NNSTU, The Nizhniy Novgorod State Technical University named after R.E.Alekseev, Russia

Re: Sediment for beginner))

#7 Post by fancer » Mon Mar 30, 2009 7:07 pm

After last line:
0 0 00:00:00 0.000000E+00 3.856303E+01 3.856303E+01 7.732200E+05
In terminal next message:
segmentation fault

ps I can give more info, if it necessary))

jcwarner
Posts: 816
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Sediment for beginner))

#8 Post by jcwarner » Mon Mar 30, 2009 7:08 pm

i do not use gfortran. may be try to compile with debug on, run it, and see if that gives more information.

fancer
Posts: 45
Joined: Fri Feb 06, 2009 5:20 pm
Location: NNSTU, The Nizhniy Novgorod State Technical University named after R.E.Alekseev, Russia

Re: Sediment for beginner))

#9 Post by fancer » Tue Mar 31, 2009 4:06 pm

I do not know why, but on other computer (at home) all work))) without seg fault)))
Can You answer me, please?
1) How ROMS calculate without grid, frc and others file or give me link with answer?
2) What mean variables rho and psi? This is density and flow function (current function) accordingly?

jcwarner
Posts: 816
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Sediment for beginner))

#10 Post by jcwarner » Tue Mar 31, 2009 4:58 pm

1) ROMS can not calculate without grid. That is required.
So you must be doing something different.
I strongly suggest that you look more carefully at info on the wiki. Perhpas at first, you should try to just run the upwelling case (the default case) and make sure that works for you.
Then move on to try anohter example.
2) the variable 'rho' is density perturbation.
the variable 'psi' is most likely the length scale term of the GLS turbulence closure. Did you define GLS_MIXING?

User avatar
kate
Posts: 3718
Joined: Wed Jul 02, 2003 5:29 pm
Location: IMS/UAF, USA

Re: Sediment for beginner))

#11 Post by kate » Tue Mar 31, 2009 5:46 pm

Psi is a horizontal streamfunction from a previous generation of the model. It is located at the corners of the grid boxes while rho is located at the centers of grid boxes.

fancer
Posts: 45
Joined: Fri Feb 06, 2009 5:20 pm
Location: NNSTU, The Nizhniy Novgorod State Technical University named after R.E.Alekseev, Russia

Re: Sediment for beginner))

#12 Post by fancer » Wed Apr 01, 2009 8:41 pm

2 jcwarner:
1) Thanks You, I more careful look the wiki page and see many interesting:))))
2) I already started sed_test1 and shoreface and upwelling case too and all work correctly))
2 kate:
Thanks You for your help)))

Question for all:
Where more easy generating grid-, force-, initial-files (in Linux) for beginner?)))

fancer
Posts: 45
Joined: Fri Feb 06, 2009 5:20 pm
Location: NNSTU, The Nizhniy Novgorod State Technical University named after R.E.Alekseev, Russia

Re: Sediment for beginner))

#13 Post by fancer » Fri May 01, 2009 6:37 pm

Hello all developers)))
I've any questions about SHOREFACE example:
1) I don't understand about wave direction as a whole. How explain what a wave come to beach with 270 degrees? I know what coordinates illustrate like this:
|0
|
|
-----------
270 | 90
|
|180
Or not? And why 180 degrees on a beach? And angle with what objects this degrees illustrate?
2) How can I see or what dimension illustrate changes bathymetry with time? May be some dimension h(x,y,t)???

Thanks a lot for all answers)))
ps I'm sorry for my english)))

jcwarner
Posts: 816
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Sediment for beginner))

#14 Post by jcwarner » Sat May 02, 2009 4:00 pm

SWAN is a full 3D model.We typically use SET NAUTICAL so that wave and winds are the directions that they are coming from. So a wave direction of 270 would mean waves are coming from the west.
IF you activate
#define SED_MORPH
then a variable 'bath(t,z,y,x)' can be available in the history files. You need to set:
Hout(idBath) == T ! time-dependent bathymetry
in the ocean*.in file.

fancer
Posts: 45
Joined: Fri Feb 06, 2009 5:20 pm
Location: NNSTU, The Nizhniy Novgorod State Technical University named after R.E.Alekseev, Russia

Re: Sediment for beginner))

#15 Post by fancer » Mon May 04, 2009 5:43 am

Thanks a lot. You very helped)))

Post Reply