****************************** -------------------------------------------- * PALM 6.0 Rev: 4071 * ocean - run without 1D - prerun ****************************** -------------------------------------------- Date: 2019-09-02 Run: oceanml Time: 15:36:39 Run-No.: 00 Run on host: default Number of PEs: 4 Processor grid (x,y): ( 2, 2) calculated ------------------------------------------------------------------------------ Numerical Schemes: ----------------- --> Use the Moeng_Wyngaard turbulence closure (LES mode). --> Use the boussinesq approximation for the model equations. --> Solve perturbation pressure via FFT using temperton-algorithm routines perturbation pressure is calculated at every Runge-Kutta step --> Momentum advection via Wicker-Skamarock-Scheme 5th order --> Scalar advection via Wicker-Skamarock-Scheme 5th order --> Loop optimization method: cache --> Time differencing scheme: runge-kutta-3 ------------------------------------------------------------------------------ Run time and time step information: ---------------------------------- Timestep: variable maximum value: 20.000 s CFL-factor: 0.90 Start time: 0.000 s End time: 120.000 s Computational grid and domain size: ---------------------------------- Grid length: dx = 1.250 m dy = 1.250 m dz(1) = 1.250 m Domain size: x = 80.000 m y = 80.000 m z(u) = -80.000 m Number of gridpoints (x,y,z): (0: 63, 0: 63, 0: 65) Subdomain size (x,y,z): ( 32, 32, 66) Characteristic levels of the geo. wind component ug: Height: 0.0 0.0 m ug: 0.00 0.00 m/s Gradient: ------ 0.00 1/100s Gridpoint: 0 65 Characteristic levels of the geo. wind component vg: Height: 0.0 0.0 m vg: 0.00 0.00 m/s Gradient: ------ 0.00 1/100s Gridpoint: 0 65 Topography information: ---------------------- Topography: flat Boundary conditions: ------------------- p uv pt B. bound.: p(0) = p(1) | uv(0) = uv(1) | pt(0) = pt(1) T. bound.: p(nzt+1) = p(nzt) | uv(nzt+1) = uv(nzt) | pt(nzt+1) = pt(nzt) e B. bound.: e(0) = e(1) T. bound.: e(nzt+1) = e(nzt) = e(nzt-1) sa B. bound.: sa(0) = sa(1) T. bound.: sa(nzt+1) = sa(nzt) Top surface fluxes are used in diffusion terms at k=nzt Predefined constant momentumflux: u: -0.000100 m**2/s**2 v: 0.000000 m**2/s**2 Predefined constant heatflux: 0.000102 K m/s Predefined constant salinityflux: 0.000000 psu m/s Lateral boundaries: left/right: cyclic north/south: cyclic Initial profiles: ---------------- Characteristic levels of the initial temperature profile: Height: 0.0 0.0 m Temperature: 300.00 300.00 K Gradient: ------ 0.00 K/100m Gridpoint: 0 65 Characteristic levels of the initial salinity profile: Height: 0.0 0.0 m Salinity: 35.00 35.00 psu Gradient: ------ 0.00 psu/100m Gridpoint: 0 65 List output: ----------- 1D-Profiles: Output every 9999999.90 s Data output: ----------- Time series: Output format: netCDF 64bit offset Output every ******** s ------------------------------------------------------------------------------ Physical quantities: ------------------- Geograph. latitude : latitude = 55.0 degr Geograph. longitude : longitude = 0.0 degr Angular velocity : omega = 0.729E-04 rad/s Coriolis parameter : f = 0.000119 1/s f* = 0.000084 1/s Day of the year at model start : day_init = 172 UTC time at model start : time_utc_init = 43200.0 s Gravity : g = 9.8 m/s**2 Reference state used in buoyancy terms: initial_profile Cloud physics quantities / methods: ---------------------------------- LES / Turbulence quantities: --------------------------- Mixing length is limited to 1.80 * z Actions during the simulation: ----------------------------- Disturbance impulse (u,v) every : ****** s Disturbance amplitude : 0.00 m/s Lower disturbance level : -28.12 m (GP 42) Upper disturbance level : -4.38 m (GP 61) Disturbances cease as soon as the disturbance energy exceeds 0.010 m**2/s**2 Random number generator used : random-parallel *** dynamic module disabled Ocean settings: ------------------------------------------ --> Craik-Leibovich vortex force and Stokes drift switched on waveheight: 1.0 m wavelength: 40.0 m --> surface heat flux is switched off after 900.0 s ------------------------------------------------------------------------------ Run-control output: ------------------ RUN ITER. HH:MM:SS.SS DT(E) UMAX VMAX WMAX U* W* THETA* Z_I ENERG. DISTENERG DIVOLD DIVNEW UMAX(KJI) VMAX(KJI) WMAX(KJI) ADVECX ADVECY MGCYC --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 0 0 00:00:00.00 20.0000A 0.0002D 0.0002D 0.0001 0.000 0.00 0.000E+00 0. 0.122E-08 0.122E-08 0.270E-04 0.387E-14 53 1 6 58 26 28 47 9 8 0.000 0.000 0 0 1 00:00:20.00 20.0000A 0.0016 -0.0002 0.0001 0.000 0.00 0.000E+00 0. 0.401E-07 0.122E-08 0.380E-06 0.423E-16 65 55 48 60 21 57 47 9 8 0.000 0.000 0 0 2 00:00:40.00 20.0000A 0.0032 -0.0002 0.0001 0.000 0.00 0.000E+00 0. 0.156E-06 0.122E-08 0.380E-06 0.422E-16 65 56 48 61 22 63 47 9 8 0.000 0.000 0 0 3 00:01:00.00 20.0000A 0.0048 -0.0003 0.0001 0.000 0.00 0.000E+00 0. 0.347E-06 0.122E-08 0.380E-06 0.421E-16 65 56 48 64 54 31 47 9 8 0.000 0.000 0 0 4 00:01:20.00 20.0000A 0.0063 -0.0004 0.0001 0.000 0.00 0.000E+00 0. 0.611E-06 0.122E-08 0.380E-06 0.419E-16 65 55 49 64 54 31 47 9 8 0.000 0.000 0 0 5 00:01:40.00 20.0000A 0.0079 -0.0005 0.0001 0.000 0.00 0.000E+00 0. 0.945E-06 0.122E-08 0.381E-06 0.417E-16 64 55 50 65 54 31 58 16 15 0.000 0.000 0 0 6 00:02:00.00 20.0000A 0.0094 -0.0007 0.0001 0.000 0.00 0.000E+00 0. 0.135E-05 0.122E-08 0.381E-06 0.415E-16 65 55 50 64 63 26 58 16 15 0.000 0.000 0