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