不好意思,前面有的,我没注意到,就在第一行:
INTEGER,PARAMETER::iwp=SELECTED_REAL_KIND(15)
INTEGER::fixed_freedoms,i,iel,iters,k,limit,nband,ndim=2,nels,neq,nip=4, &
nn,nod=8,nodf=4,nodof=3,nr,ntot=20,nxe,nye
REAL(iwp)::det,one=1.0_iwp,penalty=1.e20_iwp,pt5=0.5_iwp,rho,tol,ubar, &
vbar,visc,zero=0.0_iwp |