title:
        SED TOY Mixed bed (erosion and sediment resuspension) TEST MODEL
time_stepping: NTIMES   dt[sec]  NDTFAST  NINFO
                216000      2.      20      40

S-coord: THETA_S,   THETA_B,    Hc (m)
           1.0         0.8       0.

start_date:
2019-01-01 00:00:00
end_date:
2019-01-05 00:00:00
output_time_steps: DT_HIS(H), DT_AVG(H), DT_RST(H)
                     1         1         24

initial: NRREC  filename
          0 
                          fileini

restart:          NRST, NRPFRST / filename
                   216000    -1
                          sed_toy_resusp_rst.nc

                                                    
history: LDEFHIS, NWRT, NRPFHIS / filename
            T      3600     0
                          sed_toy_resusp_his.nc

averages: NTSAVG, NAVG, NRPFAVG / filename
            1      3600     0
                          sed_toy_resusp_avg.nc

primary_history_fields: zeta UBAR VBAR  U  V   wrtT(1:NT)
                          T    T   T   T  T    30*T
auxiliary_history_fields:   rho Omega  W  Akv  Akt  Aks Visc3d Diff3d  HBL HBBL Bostr Wstr Ustr Vstr Shfl Swfl rsw rlw lat sen  HEL
                             50*T              
gls_history_fields:   TKE  GLS  Lscale
                      T T T

rho0:
      1025.d0

lin_EOS_cff:  R0 [86kg/m3], T0 [Celsius], S0 [PSU], TCOEF [1/Celsius], SCOEF [1/PSU]
              27            10.            35.       1.7d-4               7.6d-4 

lateral_visc:   VISC2,    VISC4    [m^2/sec for all]
                 0.       0.

tracer_diff2: TNU2(1:NT)           [m^2/sec for all]
               30*0.d0 

tracer_diff4: TNU4(1:NT)           [m^4/sec for all]
               30*0.d11


vertical_mixing: Akv_bak, Akt_bak [m^2/sec]
                    1.d-05   1.d-06

bottom_drag:     RDRG [m/s],  RDRG2,  Zob [m],  Cdb_min, Cdb_max
                  3.0d-04      3.0d-03      0.0005     1.d-6  5.d-1 

gamma2:
                 1.d0

sediments: input file
               TEST_CASES/sediment_sed_toy_resusp.in

sediment_history_fields: bed_thick bed_poros bed_fra(sand,silt)
                            20000*T