c............................................................................
c This program solves the barotropic vorticity equation in non-dimensional
c form using finite differences.
c
c The model has incorporated the "partial" slipping boundary conditions.
c............................................................................
im=202 % number of grid points in the zonal direction
jm=102 % number of grid points in the meridional direction
ds=0.05 % grid step
dt=0.05 % time step
Ro=0.0 % Rossby number (measures non-linearity of the flow)
eps=0.3 % non-dimensional coefficient representing bottom friction
Ah=0.0 % non-dimensional coeff. of horizontal Laplacian mixing
Bh=0.0 % non-dimensional coeff. of horizontal bi-harmonic mixing
gamma=0.0 % coeff. of "intermediate slipping" used as boundary cond.
nst=1 % start time step number
nend=2000 % end time step number
nlpt=100 % frequency (time steps) for saving output
MCF=0 % matrix (0) or column (1) output
ncrit=4000 % number of steps allowed to do the relaxation (sub. helm)
pcrit=0.1 % criterium to stop the relaxation
BFP=1 % Beta (BFP=1) or F plane (BFP=0)
GYR=1 % Simple Gyre (GYR=1) or Double Gyre (GYR=2)
HEM=-1 % North Hemisphere Gyre (HEM=1) or South Hemisphere Gyre (HEM=-1)