      parameter(it=4,jt=200,kt=102,nphit=10,nphito=10)
c it    : i-dimension of arrays (must be >= ni)
c jt    : j-dimension of arrays (must be >= nj)
c kt    : k-dimension of arrays (must be >= nk)
c nphit : dimension of last index in phi-vector (must be >= number of variables)
c nphito: dimension of last index in phio-vector (if steady=.false., it must
c         be >= number of variables; if steady=.true. it can be set to 1 in
c         order to save memory)
c
      common
     1/prin2/imon,jmon,kmon
     1/vectc/head(nphit)
     1/vecti/nsweep(nphit)
     1/vectl/solve(nphit),isol(7)
     1/vectr/phi(it,jt,kt,nphit),phio(it,jt,kt,nphito),prt(nphit),
     2        resor(nphit),reref(nphit),
     3        smp(it,jt,kt),delta(it,jt,kt),prt_lam(nphit),
     4        fk(it,jt,kt),fe(it,jt,kt),cdiss(it,jt,kt)
     1/vectr1/conve(it,jt,kt),convn(it,jt,kt),convh(it,jt,kt)
     1/alli/ni,nj,nk,nim1,njm1,nkm1,iter,maxit,nphmax,itstep,ntstep,
     2      iprint_start,nfiles_restart,minit,
     3      jl0,jformax,jhub
     1/for_r/cdfor,fmax_leaf_area_dens,formaxdens,
     2      forheight,delta_leaf,f_lad(it,jt,kt),ustarv(it,jt,kt),
     3      forceu(it,jt,kt),forcev(it,jt,kt),forcew(it,jt,kt),
     4      qfluxh,extinct,heat_forest,dqdy(jt),grav,volextcoef,
     5      fcori,u_geo,w_geo,flowangle,hWind,
     6      alpha_geo,rrl(it,jt,kt)
     1/allr/great,small,sormax,betap,time,dt(100000),acrank,
     2      acrank_conv,blend,c_omega1,c_omega2,ynoll,tnoll
     1/allc/scheme,schtur,name,schemet
     1/alll/save,restrt,cycli,cyclk,steady,echo,les,wale,pans,kom,
     2      les_abl,stability,keq_abl,keq,des,iddes
      common
     1/geom/xc(it,jt,kt),yc(it,jt,kt),zc(it,jt,kt),vol(it,jt,kt),
     2      areaex(it,jt,kt),areaey(it,jt,kt),areaez(it,jt,kt),
     3      areanx(it,jt,kt),areany(it,jt,kt),areanz(it,jt,kt),
     4      areahx(it,jt,kt),areahy(it,jt,kt),areahz(it,jt,kt),
     5      fx(it,jt,kt),fy(it,jt,kt),fz(it,jt,kt),
     6      xp(it,jt,kt),yp(it,jt,kt),zp(it,jt,kt),dist2d(it,jt)
     1/flupr/urfvis,viscos,densit,vis(it,jt,kt)
     1/turb/cmucd,cmu,cd,c1,c2,cappa,elog,pfun
     1/coef/ap(it,jt,kt),an(it,jt,kt),as(it,jt,kt),ae(it,jt,kt),
     2      aw(it,jt,kt),ah(it,jt,kt),al(it,jt,kt),sp(it,jt,kt),
     3      su(it,jt,kt),suv(it,jt,kt),suw(it,jt,kt),spvw(it,jt,kt),
     4      gente(it,jt,kt)
     1/point/u,v,w,p,pp,m,te,ed,t,f1,f2,f3,om,conc
      integer u,v,w,p,pp,m,te,ed,t,f1,f2,f3,om,conc
      logical solve,save,restrt,cycli,cyclk,steady,echo,les,kom
      logical pans,wale,les_abl,stability,keq_abl,keq,des,iddes
      character plane*2,head*24,scheme*2,schtur*2,name(nphit)*2,
     1schemet*2
