sam/sam.h at master · tslund/sam · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
      parameter( mLx=2050, mLy=2050, mLz=2050 )
      parameter( Ls=54, L_params=100 )
      parameter( maxp=1024)

      real    lambda_z_fs, Nx_inv, Ny_inv, Nz_inv, NxNy_inv, NxNyNz_inv,
     &        l_cbc, m_cbc, k_truncate, k_truncate_sq, lapse0, lapse,
     &        lambda_x, lambda_z, k_w, m_w, N_sq
      complex iunit
      integer fh
      character(12) labels
      logical l_root, required, fixed, found, squak

      common/static/ xL, yL, zL, c_smag, dt0, cfl0, k_truncate, z0, to,
     & lapse0, grav, vis, pr, amplitude, lambda_x, Gam, omega, Uo,
     & shear, shear_ratio, flct_u, flct_t, lambda_z_fs, Ri_fs,
     & Nx, Ny, Nz, i_restart, nt_restart, i_strat, i_les, n_steps,
     & n_skip_h, n_skip_p, n_skip_v, n_skip_s, k_xy_plane(9),
     & j_xz_plane(9), i_yz_plane(9), i_stat, i_cfl, nrk_max, n_dealias,
     & iubc_z, i_prob, i_force, k_force, i_fs

      common/dynamic/ time, dt, t_stat,
     & nt, n_frame_p, n_hist

      common/consts/ pi, two_pi, iunit

      common/inputs/   labels(L_params), values(L_params),
     &               required(L_params),  fixed(L_params),
     &                  found(L_params),  squak(L_params)

      common/strat/ buoy_fac_x, buoy_fac_z, Pr_inv, To_inv

      common/gravity_wave/ N_sq, k_w, m_w, omega_i, lambda_z,
     &                     sin_theta, cos_theta, lapse,
     &                     i_gw_type

      common/dimensions/ Lu, Lr, Lw, Lm, Nze, ipad_r,
     &                   Nx_min, Ny_min, Nz_min,
     &                   n_inputs_r, n_inputs_i, n_inputs,
     &                   n_dynpar_r, n_dynpar_i, n_dynpar,
     &                   n_params

      common/mesh/ dx, dy, dz, dx_inv, dy_inv, dz_inv,
     &             xL_inv, yL_inv, zL_inv

      common /mpi_stuff/ myid, numprocs, l_root, i_root,
     &                   ixs, ixe, izs, ize,
     &                   nxp, nzp, fh(100),
     &                   ix_s( 0:maxp), ix_e( 0:maxp), nx_p(0:maxp),
     &                   iz_s( 0:maxp), iz_e( 0:maxp), nz_p(0:maxp)

      common/tstep/ gamma(4), zeta(4), beta(4), t_start,
     &              adv_cfl_limit, vis_cfl_limit,
     &              cfl_x, cfl_y, cfl_z, cfl_vis,
     &              vis_fac_x(mLx), vis_fac_y(mLy), vis_fac_z(mLz),
     &              nrk, nt_start

      common/les/ E_cl, a_cl, b_cl, c_cl, k_cl

      common/cbc/ l_cbc, m_cbc, t0_cbc, u_inf_cbc, u_0_cbc, tfact

      common/de_alias/ z_shift, k_truncate_sq

      common/fft/ trigx( 2*mLx+15), trigy( 4*mLy+15), trigz( 4*mLz+15),
     &            trigzs(3*mLz+15), trigzc(3*mLz+15), trigzr(2*mLz+15),
     &            trigyr(2*mLy+15),
     &            wave_x, wave_y, wave_z,
     &            k_x(mLx), k_y(mLy), k_z(mLz),
     &            i_symm(4), Nx_inv, Ny_inv, Nz_inv,
     &            NxNy_inv, NxNyNz_inv,
     &            kx_max, ky_max, kz_max,
     &            Nx2, NxNy, NxNyNz

      common/planes/ n_xy_planes, n_xz_planes, n_yz_planes

      common/statistics/ t_stat0,
     &                   u_avg(4), u_var(4), u_max(4), vis_max, energy,
     &                   stat(Ls), stat_tmp(Ls)