[求助] 新人初学Fortran，有一段代码好像第一步语法就出错了。

Last change: N     17 April 2003   10.00 am
! Program to Calculate water surface profileof ice covered
! channel and non-equilibrium ice coverthickness
! Nimal C. Jayasundara and Prof. Hung TaoShen
! Department of Civil and EnvironmentalEngineering
! Clarkson University
      program ne_jam
    parameter(kr=65,kn=50, kt=100)
    commonnhecu(kr),ywidu(kr,kn),zbu(kr,kn),
    $      nhecd(kr),ywidd(kr,kn),zbd(kr,kn)
    realbman(kr),iman(kr),xl(kr),lengthx(kr), x(kr)
    realheud, h0(kr), hu1, hu2
    realr(kr), nc(kr), hh(kr), hw(kr)
    realb(kr),p(kr),area(kr),dd(kr), v(kr),k(kr)
    realsf(kr), error(kr), eed
    realhini, q, tini, tsave(kr)
    realvmax, rhow, rhoi, rhorat, coef, alpha, g, pi, phi, e
    realk0,k1,xkk,kx,gama_e
    realti(kr),  t1,t21,t22, tavg
    realaa,bb, cc, top, ravg
    realric(kr), sw, dtdx
    realhw1(kr), ze(kr)
    realerrorv, tolv, tolh, tolice
    realhtop(kr),hbotom(kr)
    realni0, ti0, ni(kr), nimax, tti(kr)
    realti_del(kr), area_req
    realt_step(kt), Qice(kt)
    realicevol(kr), cover_vol(kr)
    realvol_diff, ice_supply
    integercount, itern, istep, vstep, ni_option
    open(unit = 2, file = 'thames-g.dat')
    open(unit = 3, file = 'thames-p.dat')
    open(unit = 4, file = 'thames-ini.dat')
    open(unit = 7, file = 'thames_icev.dat')