!=============================================================================
!   PARAMETER FILE FOR PARALLEL SPECTRAL MODEL.  =============================
!=============================================================================
      runlabel='default_IW'
      nx=16           !! collocation pts in x direction
      ny=1            !! collocation pts in y direction
      nz=16           !! collocation* pts in z direction
      ntrajectories=0
      numprocs=01

!     *For slip bcs, we actually use, say, 17 pts for cos transforms in z
!     but still specify nz=16. The sin transforms use only 15 pts but all
!     variables u,v,w etc, regardless of whether they are odd or even
!     functions of z, would have 17 valid data pts in the z direction.
!=============================================================================
!    Time step
!=============================================================================
      dt=600.       !! integration time step [s]
      t_start=0.  !! starting time         [s]
!!!   t_end=3000000.    !! ending time           [s]
      t_end=6000.

!=============================================================================
!     Declare how the two scalars relate to density
!     R (rho)  T (temp) S (salinity)   P (passive)
!     Legal Settings: 'TS' 'TP' 'SP' 'RP'  'PP'
!=============================================================================
!!! xlv      scalars = 'RP'  !! scalar1=rho'  scalar2=passive
      scalars = 'RP'
      s1_scale= 3.e-2  
      s2_scale= 1  
 

!=============================================================================
!    Type of z boundary conditions.
!=============================================================================
      bc_flag = 'zslip'     !! 'zperiodic'  or 'zslip' allowable
                            !!  always periodic in x,y

!=============================================================================
!    Size of computational domain.
!    x in [0,Lx)   , y in [0,Ly)   , z in [0,Lz]  zslip
!                                         [0,Lz)  zperiodic
!=============================================================================
      Lx = 2000.d3      !! x domain size [m]
      Ly = 2000.d3       !! y domain size [m]
      Lz = 4200.d0       !! z domain size [m]



!=============================================================================
!     Flag specifying whether flow is forced via user_defined_forcing.f
!=============================================================================
      force_flag='yes'
      nstep_storm=72
      istart_storm=1

!=============================================================================
!     Physical parameters
!=============================================================================
      g=9.81
      xlat=47.5
      xOmega=7.29e-5
      Rayon=6400.e3
      plan_trad=1.        ! =0 pas de composante verticale du paramètre de coriolis
      plan_beta=0.        ! =0 pas d'effet beta
      plan_nontrad=0.     ! =0 pas de composante horizontale du paramètre de coriolis
      plan_beta_tild=0.   ! =0 pas d'effet beta tild

!=============================================================================
!     Fluid properties
!=============================================================================
      rho_0=1027.             ! characteristic density [kg/m3]
      nu = 1.e-6              ! viscosity [m2/s]
      kappa_1 = 1.e-6         ! diffusivity of scalar 1 [m2/s]
      kappa_2 = 1.e-6         ! diffusivity of scalar 2 [m2/s]


!=============================================================================
!     Additional user specified scales required to nondimensionalize
!     the equations of motion. (Lz,nu,kappa)
!=============================================================================
      DGRAD=3.e-4             ! char. scale of dens gradient [kg/m4]
      U0=0.6                  ! char. velocity scale [m/s]