应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 求fortran95编写的数值程序
51/1返回列表
查看: 995  |  回复: 4
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xmch2011

铁虫 (小有名气)

[求助] 求fortran95编写的数值程序

想学学用fortran95写的数值程序，那位同学有没有这样的程序，学习一下，谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
数学，力学，物理
1楼2012-04-01 22:01:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wll778824

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wll778824

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
小雨萌萌: 金币+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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wll778824

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
小雨萌萌: 金币+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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmch2011

铁虫 (小有名气)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小雨萌萌: 金币-10, 谢人家也不把悬赏的金币发一下呀？帮你发放金币！ 2012-04-05 16:33:33
引用回帖:
4楼: Originally posted by wll778824 at 2012-04-02 12:13:01:
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楼2012-04-03 11:35:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xmch2011 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭