| 1 | C***************************************************************************** |
|---|
| 2 | c PARAMETER FILE FOR PARALLEL SPECTRAL MODEL. ***************************** |
|---|
| 3 | c***************************************************************************** |
|---|
| 4 | pi=4.*atan(1.) |
|---|
| 5 | |
|---|
| 6 | c Initial conditions: either name a file or specify none |
|---|
| 7 | c xlv netcdf_file='/home/alpha0/kraig/NSF_IWS/input/ics/GM_cdf_2145.000' |
|---|
| 8 | c netcdf_file='/workdir/rech/dyc/rdyc711/cfd_2D_000500' |
|---|
| 9 | netcdf_file='none' ! none means ICS defined in initialize.f |
|---|
| 10 | |
|---|
| 11 | c Flag specifying whether flow is forced via user_defined_forcing.f |
|---|
| 12 | c force_flag='yes' |
|---|
| 13 | force_flag='yes' |
|---|
| 14 | |
|---|
| 15 | |
|---|
| 16 | c***************************************************************************** |
|---|
| 17 | c Type of z boundary conditions. |
|---|
| 18 | c***************************************************************************** |
|---|
| 19 | bc_flag = 'zslip' ! 'zperiodic' or 'zslip' allowable |
|---|
| 20 | ! always periodic in x,y |
|---|
| 21 | c============================================================================= |
|---|
| 22 | c Declare how the two scalars relate to density |
|---|
| 23 | c R (rho) T (temp) S (salinity) P (passive) |
|---|
| 24 | c Legal Settings: 'TS' 'TP' 'SP' 'RP' 'PP' |
|---|
| 25 | c============================================================================= |
|---|
| 26 | scalars = 'TS' ! scalar1=rho' scalar2=passive |
|---|
| 27 | s1_scale= 1.0 |
|---|
| 28 | s2_scale= 1.0 |
|---|
| 29 | |
|---|
| 30 | c***************************************************************************** |
|---|
| 31 | c Size of computational domain. |
|---|
| 32 | c x in [0,Lx) , y in [0,Ly) , z in [0,Lz] zslip |
|---|
| 33 | c [0,Lz] zperiodic |
|---|
| 34 | c***************************************************************************** |
|---|
| 35 | Lx = 150000 ! x domain size [m] |
|---|
| 36 | Ly = Lx ! y domain size [m] (unused for 2d see below) |
|---|
| 37 | Lz = 4000. ! z domain size [m] |
|---|
| 38 | dz = Lz/float(nz) ! vertical grid spacing [m] (*$*) |
|---|
| 39 | c***************************************************************************** |
|---|
| 40 | c***************************************************************************** |
|---|
| 41 | |
|---|
| 42 | |
|---|
| 43 | c***************************************************************************** |
|---|
| 44 | c Physical parameters |
|---|
| 45 | c***************************************************************************** |
|---|
| 46 | g=9.81 ! gravity [m/s2] |
|---|
| 47 | c xlv xlat=pi/4 ! latitude in radians |
|---|
| 48 | xlatit=47.5*pi/180 |
|---|
| 49 | xOmega=2.*pi/(24.*3600.) ! earth rotation frequency [1/s] |
|---|
| 50 | f=2*xOmega*sin(xlatit) ! Coriolis parameter [1/s] |
|---|
| 51 | c beta=(1/R)2(Omega)cos(xlatit); R~6400km ! beta effect pba June 2005 |
|---|
| 52 | c xlv beta=(1./(6400*1.e3))*2*xOmega*cos(xlatit) ! [1/ms] |
|---|
| 53 | beta=0. |
|---|
| 54 | |
|---|
| 55 | c***************************************************************************** |
|---|
| 56 | c***************************************************************************** |
|---|
| 57 | |
|---|
| 58 | |
|---|
| 59 | |
|---|
| 60 | c***************************************************************************** |
|---|
| 61 | c Fluid properties |
|---|
| 62 | c***************************************************************************** |
|---|
| 63 | rho_0=1000. ! characteristic density [kg/m3] |
|---|
| 64 | nu = 1.e-6 ! viscosity [m2/s] (3),(4) |
|---|
| 65 | kappa_1 = nu/1. ! diffusivity for scalar 1 [m2/s] |
|---|
| 66 | kappa_2 = nu/1. ! diffusivity for scalar 2 [m2/s] |
|---|
| 67 | c***************************************************************************** |
|---|
| 68 | c dissipation scheme |
|---|
| 69 | c for DNS, i.e. nu grad^2 u, kappa grad^2 rho |
|---|
| 70 | c choose 'isotropic' and p=2, higher orders of p |
|---|
| 71 | c and anisotropic are undocumented currently |
|---|
| 72 | diss_flag='isotropic' ! isotropic or anisotropic |
|---|
| 73 | p=2. ! order of (hyper)viscosity operator |
|---|
| 74 | T_diss=3000 ! [s] decay time scale at maximim wavenumber |
|---|
| 75 | ! ignored if diss_flag=isotropic and p=2 i.e. for DNS |
|---|
| 76 | c***************************************************************************** |
|---|
| 77 | |
|---|
| 78 | c***************************************************************************** |
|---|
| 79 | c Additional user specified scales required to nondimensionalize |
|---|
| 80 | c the equations of motion. |
|---|
| 81 | c***************************************************************************** |
|---|
| 82 | DGRAD=3.e-4 ! char. scale of dens gradient [kg/m4] |
|---|
| 83 | c xlv s1_scale= 1.0 ! char. scale of scalar 1 e.g. would be deg C if s1->Temp. |
|---|
| 84 | c xlv s2_scale= 1. ! char. scale of scalar 2 e.g. would be psu if s2->Salinity |
|---|
| 85 | U0=0.01 ! char. velocity scale [m/s] |
|---|
| 86 | bfreq = sqrt((g/rho_0)*DGRAD) ! char. buoyancy frequency [1/s] (*$*) |
|---|
| 87 | c***************************************************************************** |
|---|
| 88 | c***************************************************************************** |
|---|
| 89 | |
|---|
| 90 | c***************************************************************************** |
|---|
| 91 | c Time step and execution control parameters |
|---|
| 92 | c***************************************************************************** |
|---|
| 93 | dt=600 ! integration time step [s] |
|---|
| 94 | t_start=0*dt ! starting time [s] |
|---|
| 95 | t_end=100*dt ! ending time [s] (test wave per = 4.89 hrs --> 3 pers.) |
|---|
| 96 | vort_flag='no' ! write out vorticity fields as well as velocity/density |
|---|
| 97 | t_stat=20*dt ! time increment for calls to energetics [s] |
|---|
| 98 | c***************************************************************************** |
|---|
| 99 | c***************************************************************************** |
|---|
| 100 | |
|---|
| 101 | c****************************************************************************** |
|---|
| 102 | c Lagrangian "float" parameters |
|---|
| 103 | c****************************************************************************** |
|---|
| 104 | t_floats_on =t_start ! float insertion time |
|---|
| 105 | z_close=0.5*dz/Lz ! min. distance to upper/lower bdries |
|---|
| 106 | ! that floats can attain [d'less]=[1] |
|---|
| 107 | z_offset=-0.5*dz/Lz ! up offset for second float dens. measurement |
|---|
| 108 | delta_w=0.0/U0 ! float velocity due to residual buoyancy [1] |
|---|
| 109 | c****************************************************************************** |
|---|
| 110 | c****************************************************************************** |
|---|
| 111 | |
|---|
| 112 | c****************************************************************************** |
|---|
| 113 | c Perturbed, "noisy" ICS parameter |
|---|
| 114 | c****************************************************************************** |
|---|
| 115 | efactor=0*1.e-4 |
|---|
| 116 | c if efactor is nonzero, the (linear) vortical mode field is perturbed |
|---|
| 117 | c with a spectrum of "noise" (decays linearly with horiz. wavenumber) such |
|---|
| 118 | c that the total kinetic energy of the noise is a specified fraction of |
|---|
| 119 | c the horizontal kinetic energy of the deterministically specified ICs |
|---|
| 120 | c efactor=0.01 --> noise field has 1% of the HKE in the deterministic fields |
|---|
| 121 | c****************************************************************************** |
|---|
