ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c     lagr - defines Semi-lagrangin (1) or Eulerian (0 - DEFAULT) model
c ------------------------------------------------------------------------- 
c  requires consistent specification of precompiler directive SEMILAG
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (lagr=0,ieul=1-lagr)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c  n,m,l  - grid sizes in X,Y,Z respectively
c ------------------------------------------------------------------------- 
C   **** WARNING: if using multiprocessing, there exist constraint
C  ****          on values of m,n. See note below on "msg.inc" file
C  ****          and values of nprocx and nprocy.
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

c     parameter (n=128,m=64,l=21)                ! grid dimensions 
      parameter (n=128,m=64,l=51)                ! grid dimensions 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c ---  dx00, dy00, dz00 - domain grid increments/dimensions
c ---  dt00 - time step increment in seconds
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

!     parameter (dx00=10.)
!     parameter (dy00=10.)
!     parameter (dz00=5.)
!     parameter (dt00=0.1)

      parameter (dx00=10.)
      parameter (dy00=10.)
      parameter (dz00=10.)
      parameter (dt00=0.25)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c     lxyz - define meaning of dx00, dy00, dz00
c---> lxyz=0 physical grid increments, need to specify values of dx,dy,dx
c---> lxyz=1 physical domain lengths,  need to specify values of dx,dy,dx
c---> lxyz=2 normalized grid increments, do not need values of dx,dy,dx
c -------------------------------------------------------------------------
c ---  Cartesian (all units in meters and seconds)
c ---  Spherical (all units in meters and seconds)
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (lxyz=0)

ccoccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c ---     nt: number of timesteps
c ---    nth: max NT for histories
c ---  noutp: print diagnostic output/graphics plot after every "noutp" step
c ---  nplot: store abbreviated tape every "nplot" step (iowritesh)
c --- nstore: store full dataset after every "nstore" step
c---> nslice: store data for the graphics postprocessing       
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (nt=600,noutp=10,nplot=10,nstore=600,nslice=0)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c Basic physical/numerical parameters           
c
c ---   implgw = 1 advect theta perturbation (default)
c ---            0 advect full theta
c ---  isphere = 0 reference coordinates are Cartesian
c ---             (need to specify dx,dy,dz in blanelas)
c ---            1 reference coordinates are spherical
c ---              (need to specify only dz in blanelas)
c ---   icorio = 0 no coriolis accelerations
c ---            1 incorporate coriolis accelerations
c ---              (need to set icorio even if isphere=1)
c ---  icylind =0  reference coordinates are Cartesian
c ---              reference coordinates are cylindrical
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      parameter (implgw=1,isphere=0,icylind=0,icorio=0)
      parameter (intz=1-implgw)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c ---  iwrite: 1 make history tape on unit 9 (fort.9)
c --- iwrite0: 0 do not save/1 save the initial data on restart 
c ---    irst: 1 restart run from history tape (fort.10)
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (iwrite=0,iwrite0=0,irst=0)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c              CREATE RESTART/GRAPHICS/ANALYSIS PARAMETERS
c -------------------------------------------------------------------------
c ---    nfil: number of output files generated in the tape (if iwrite=1)
c ---   nfilm: parameter used to exclude initial "0" time step (if iwrite0=1)
c ---  dt_fil: timestep used in the tape run
c --- nt_fil0: time step number for the first output:
c              1) =0 initial (diagnostic) time step (if iwrite0=1) 
c              2) =nt_film0 all outputs from prognostic time steps (iwrite0=0)
c ---nt_film0: time step increments for consecutive outputs
c              nt_film0=nstore (value of "nstore" from the tape run)
c --- nfstart: first record in tape file to read
c              1) ANALIZ=0,irst=0 for initialization (fort.10)
c              2) ANALIZ>0        for analysis       (fort.10 or fort.12)
c ---  nanlfl: switch between output tapes
c              1) =0 short tape fort.12 for graphisc only 
c              2) =1 full tape fort.10 for graphics and diagnostic
c ---    ndia: flag for diagnostics during analyasis (only if nanlfl=1)
c ---    nplo: frequency of the plots, analyse every nplo'th record
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (nfil =3,             nfilm =nfil -1)
      parameter (nfilo=3,             nfilom=nfilo-1)
      parameter (dt_fil =600.e0,nt_fil0 =0,nt_film0  =720)
      parameter (dt_filo=600.e0,nt_filo0=0,nt_filom0 =720)

      parameter (nfstart=1)
      parameter ( nanlfl=1)
      parameter (   idia=1 *nanlfl)
      parameter (   nplo=1)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c     Treatment of the boundaries
c -------------------------------------------------------------------------
c ---     icyx,icyy,icyz: boundaries are open (=0) or periodic (=1)
c ---                     in x,y,z coordinates
c ---        irelx,irely: do not apply (=0) or apply damping (=1) 
c ---                     along x,y boundaries
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (icyx=1)
      parameter (icyy=1)
      parameter (icyz=0)
      parameter (irelx=0,irely=0)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c     Dry/Moist model parameters
c -------------------------------------------------------------------------
c --- moist = 0 - dry model
c --- moist = 1 - warm/cold rain moist model depending on ice options below
c                 requires setup of the precompiler MOIST directive
c --- ice   = 1 - activate simple model ice 
c --- iceab = 1 - two classe of ice A, B
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      parameter (moist=0)                         ! Dry/Moist model
      parameter (ice=0)                  ! warm rain =0/ simple ice =1
      parameter (iceab=0)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c Viscous/Inviscid model parameters
c -------------------------------------------------------------------------
c --- ivs0 = 0 - invisicd model or ILES
c --- ivs0 = 1 - viscous, requires setup of the precompiler SGS directive
c --- irid0 = 1 - computes diagnosticis of the Richardson number
c --- itke0 = 0 - Smagorinsky model
c --- itke0 = 1 - prognostic TKE
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (ivs0=1)                          ! Viscous/Inviscid model
      parameter (irid0=0)                         ! Viscous/Inviscid model
      parameter (itke0=1,itke=ivs0*itke0)         ! Smagorinsky/TKE SGS model

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c Chemistry/passive tracer on-off
c -------------------------------------------------------------------------
c     nspc - number of chemical spices/tracers
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (ichm=1,nspc=1)              ! chemical spices

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c spectral preconditioner
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter(ispcpr=0)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c immersed boundary; requires set IMRSB precompiler directive
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter(imrsb=1)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c icw ontrols lateral open boundaries: 
c icw=0 "overspecified" enforces component w=0 at the open boundaries; 
c        icw=0 makes sense only with absorbers ON
c icw=1 "correct" does not invoke components other than normal; 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (icw=1)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
control parameters and fortran statement functions for vertical stretching
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter (istr=0,SD=7000.)

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c ---  time of time integration
c--> itraj0=0 adams-bashf., itraj0=1 runge-kutta
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

      parameter(itraj=0)                          ! Ad-Bs/Rg-Kt advecting flow

ccoccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c ---    nth: max NT for histories
c ---    nthv: max NT for histories in viscous/chemistry model
ccoccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      parameter (nth=nt+1)
      parameter (nthv=nth,nthv2=2*nthv)
      parameter (nthch=nspc*nthv)