| 查看: 1325 | 回复: 4 | ||
xmch2011铁虫 (小有名气)
|
[求助]
求fortran95编写的数值程序
|
| 想学学用fortran95写的数值程序，那位同学有没有这样的程序，学习一下，谢谢！ |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助计算机高人，编写个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
200