Hello,
me again. I just saw on in the browser that you can't see the attached
files. So here they are inline:
*************************************************************************
*** main.f *************************************************************
program main
include 'CASECOM'
include 'COMMON'
character*24 hedvis,hedden,hedx,hedy,hedz
real*4 dtime, etime
real*4 darray(2),dd,tarr(3)
parameter (ittelaps=10)
real*4 ttstepi,ttstep,ttstept,ttsteps,ttelaps(ittelaps)
dimension ustart(jt)
character*8 ttnames(ittelaps)
data ttnames/' calcu ',' calcv ',' calcw ',' conv '
&,' calcpe ',' calcte ',' vist ',' calced',' calcph'
&,' SUM '/
c------define integer pointers used in the vectors phi, urf, nsweep etc.
u=1
v=2
w=3
p=4
pp=5
te=6
ed=7
t=8
f1=9
f2=10
end
*************************************************************************
*** COMMON *************************************************************
c parameter(it=164,jt=52,kt=100,nphit=8,nphito=8)
parameter(it=478,jt=303,kt=215,nphit=8,nphito=8)
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/pres/ipref,jpref,kpref
1/prin1/plane
1/prin2/indmon,imon,jmon,kmon
1/vecti/nsweep(nphit)
1/vectl/solve(nphit)
1/vectc/head(nphit)
1/vectr/phi(it,jt,kt,nphit),phio(it,jt,kt,nphito),prt(nphit),
2 resor(nphit),reref(nphit),urf(nphit),ustarv(it,jt,kt),
3 phic(it,jt,kt,6),cs(it,jt,kt),tracel(it,jt,kt),
4 csk(it,jt,kt),
6 conve(it,jt,kt),convn(it,jt,kt),convh(it,jt,kt),
7 smp(it,jt,kt),dtfals(nphit)
1/alli/ni,nj,nk,nim1,njm1,nkm1,iter,maxit,nphmax,itstep,ntstep,
2 iturb,itime,jst1
1/allr/great,small,sormax,betap,time,dt(100000),tscoef,acrank
1/allc/scheme,schtur,name,linsol(nphit)
1/alll/save,restrt,cycli,cyclk,ortogo,steady,echo,les,expl,dyn,
2 dynk,multig
common
1/geom/xc(it),yc(jt),zc(kt),vol(it,jt,kt),
2 areaex(jt,kt),areany(it,kt),areahz(it,jt),
5 fx(it),fy(jt),fz(kt)
1/flupr/urfvis,viscos,vismin,densit,vis(it,jt,kt)
1/turb/cmucd,cmu,cd,c1,c2,cappa,elog,pfun,prandl
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)
1/point/u,v,w,p,pp,m,te,ed,t,f1,f2,f3
integer u,v,w,p,pp,m,te,ed,t,f1,f2,f3
logical solve,save,restrt,cycli,cyclk,steady,ortogo,echo,les
logical dynk,multig,expl,dyn
character plane*2,head*24,scheme*2,schtur*2,name(nphit)*2,linsol*2
Greetings!
Fabian
> > There's something missing in this example, as this doesn't compile.
> > Please send us a complete, compilable example.
> >
>
> Yes, the example is attached.
>
> Thanks!
> Fabian
