新人初学Fortran，有一段代码好像第一步语法就出错了。 - Fortran Coder

[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')