| 查看: 1058 | 回复: 4 | ||
xmch2011铁虫 (小有名气)
|
[求助]
求fortran95编写的数值程序
|
| 想学学用fortran95写的数值程序，那位同学有没有这样的程序，学习一下，谢谢！ |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助计算机高人，编写个Fortran程序
已经有13人回复
- 【求助】求助用Fortran编写模拟中的接触角
已经有4人回复
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】有关fortran进行数值模拟
已经有10人回复
- 【求助】如何实现编写的fortran程序让别人使用，但不让别人知道源代码
已经有6人回复
wll778824
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1013.4
- 散金: 242
- 红花: 1
- 帖子: 273
- 在线: 150.9小时
- 虫号: 917314
- 注册: 2009-12-01
- 专业: 计算数学与科学工程计算
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
小雨萌萌: 金币+5, 帮发放金币咯~ 2012-04-05 16:34:13
感谢参与，应助指数 +1
小雨萌萌: 金币+5, 帮发放金币咯~ 2012-04-05 16:34:13
|
给你一个有限体积算法求解Burgers 方程的fortran程序，源代码用固定格式写的 !Solve Burgers equation u_t+(u^2/2)_x=0 using Finite Volume Method !The time is discretized by using RK3 module constant implicit double precision(a-h,o-z) parameter pi=3.1415926,dt=0.001,nw=2000 parameter n=20,h=2.0*pi/dble(n),san=0.0 end module program FVM use constant implicit double precision(a-h,o-z) common u0(-1:n+2),u1(-1:n+2),u2(-1:n+2) common/u/ u(-1:n+2),uminus(0:n),uadd(0:n) !uminus=u^-,uadd=u^+ common/flux/ flux(0:n) common/ua/ ua1(n),ua2(n),state(n) common u3(-1:n+2) do j=1,n u0(j)=1.0/3.0+2.0/(3.0*h) &*(cos((dble(j-1))*h+san)-cos(dble(j)*h+san)) enddo u0(0)=u0(n) u0(-1)=u0(n-1) u0(n+1)=u0(1) u0(n+2)=u0(2) u=u0 call RECONS do ntime=1,nw call AVERAGE u1=0 u0(0)=0 u0(n+1)=0 u(0)=0 u(n+1)=0 u1=u0+dt*u u=u1 u(0)=u(n) u(-1)=u(n-1) u(n+1)=u(1) u(n+2)=u(2) call RECONS !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! call AVERAGE u2=0 u0(0)=0 u0(n+1)=0 u1(0)=0 u1(n+1)=0 u(0)=0 u(n+1)=0 u2=0.75*u0+0.25*u1+0.25*dt*u u=u2 u(0)=u(n) u(-1)=u(n-1) u(n+1)=u(1) u(n+2)=u(2) call RECONS !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! call AVERAGE u3=0 u0(0)=0 u0(n+1)=0 u2(0)=0 u2(n+1)=0 u(0)=0 u(n+1)=0 u3=(u0+2.0*u2+2.0*dt*u)/3.0 u=u3 u(0)=u(n) u(-1)=u(n-1) u(n+1)=u(1) u(n+2)=u(2) u0=u call RECONS write(*,*)ntime if(ntime==1400)then open(11,file='1.4uadd.dat') do i=0,n write(11,'(2f15.6)')san+dble(i)*h,uadd(i) enddo open(12,file='1.4uminus.dat') do i=0,n write(12,'(2f15.6)')san+dble(i)*h,uminus(i) enddo elseif(ntime==1500)then open(13,file='1.5uadd.dat') do i=0,n write(13,'(2f15.6)')san+dble(i)*h,uadd(i) enddo open(14,file='1.5uminus.dat') do i=0,n write(14,'(2f15.6)')san+dble(i)*h,uminus(i) enddo elseif(ntime==2000)then open(15,file='2uadd.dat') do i=0,n write(15,'(2f15.6)')san+dble(i)*h,uadd(i) enddo open(16,file='2uminus.dat') do i=0,n write(16,'(2f15.6)')san+dble(i)*h,uminus(i) enddo elseif(ntime==1000)then open(17,file='1uadd.dat') do i=0,n write(17,'(2f15.6)')san+dble(i)*h,uadd(i) enddo open(18,file='1uminus.dat') do i=0,n write(18,'(2f15.6)')san+dble(i)*h,uminus(i) enddo endif enddo !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! end 以上为主程序，以下为子程序 |
2楼2012-04-02 12:12:10
wll778824
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1013.4
- 散金: 242
- 红花: 1
- 帖子: 273
- 在线: 150.9小时
- 虫号: 917314
- 注册: 2009-12-01
- 专业: 计算数学与科学工程计算
【答案】应助回帖
★ ★ ★
小雨萌萌: 金币+3, 3Q! 2012-04-05 16:33:56
小雨萌萌: 金币+3, 3Q! 2012-04-05 16:33:56
|
subroutine AVERAGE use constant implicit double precision(a-h,o-z) common/u/ u(-1:n+2),uminus(0:n),uadd(0:n) common/flux/ flux(0:n) flux=0 do i=0,n flux(i)=0.5*(0.5*uminus(i)*uminus(i) &+0.5*uadd(i)*uadd(i)-0.6*(uadd(i)-uminus(i))) enddo do j=1,n u(j)=(flux(j-1)-flux(j))/h enddo end |
3楼2012-04-02 12:12:44
wll778824
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1013.4
- 散金: 242
- 红花: 1
- 帖子: 273
- 在线: 150.9小时
- 虫号: 917314
- 注册: 2009-12-01
- 专业: 计算数学与科学工程计算
【答案】应助回帖
★ ★ ★
小雨萌萌: 金币+3, 3Q! 2012-04-05 16:33:44
小雨萌萌: 金币+3, 3Q! 2012-04-05 16:33:44
|
subroutine RECONS use constant implicit double precision(a-h,o-z) common/u/ u(-1:n+2),uminus(0:n),uadd(0:n) common/flux/ flux(0:n) uminus=0 uadd=0 do j=0,n uminus(j)=-1.0/6.0*u(j-1)+5.0/6.0*u(j)+1.0/3.0*u(j+1) enddo do j=0,n uadd(j)=1.0/3.0*u(j)+5.0/6.0*u(j+1)-1.0/6.0*u(j+2) enddo end |
4楼2012-04-02 12:13:01
xmch2011
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 184
- 散金: 10
- 帖子: 75
- 在线: 45.2小时
- 虫号: 1264781
- 注册: 2011-04-13
- 性别: GG
- 专业: 计算数学与科学工程计算
5楼2012-04-03 11:35:00
