版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

xmch2011

鐵蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 184
散金: 10
帖子: 75
在線: 45.2小時(shí)
蟲(chóng)號(hào): 1264781
注冊(cè): 2011-04-13
性別: GG
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

[求助] 求fortran95編寫(xiě)的數(shù)值程序

想學(xué)學(xué)用fortran95寫(xiě)的數(shù)值程序，那位同學(xué)有沒(méi)有這樣的程序，學(xué)習(xí)一下，謝謝！

回復(fù)此樓

» 猜你喜歡

303求調(diào)劑已經(jīng)有6人回復(fù)
306求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有6人回復(fù)
275求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
考研調(diào)劑已經(jīng)有4人回復(fù)
285求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有7人回復(fù)
318求調(diào)劑已經(jīng)有6人回復(fù)
305分求調(diào)劑（食品工程）已經(jīng)有5人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助計(jì)算機(jī)高人，編寫(xiě)個(gè)Fortran程序已經(jīng)有13人回復(fù)
【求助】求助用Fortran編寫(xiě)模擬中的接觸角已經(jīng)有4人回復(fù)
【求助】如何寫(xiě)FORTRAN程序?qū)崿F(xiàn)求平均最近鄰距離已經(jīng)有34人回復(fù)
【求助】用fortran怎么描述下面這樣的情況，我是新手，請(qǐng)大家?guī)兔?/a> 已經(jīng)有8人回復(fù)

【求助】有關(guān)fortran進(jìn)行數(shù)值模擬已經(jīng)有10人回復(fù)
【求助】如何實(shí)現(xiàn)編寫(xiě)的fortran程序讓別人使用，但不讓別人知道源代碼已經(jīng)有6人回復(fù)

數(shù)學(xué)，力學(xué)，物理

1樓 2012-04-01 22:01:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wll778824

銀蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1013.4
散金: 242
紅花: 1
帖子: 273
在線: 150.9小時(shí)
蟲(chóng)號(hào): 917314
注冊(cè): 2009-12-01
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

【答案】應(yīng)助回帖

★ ★ ★
小雨萌萌: 金幣+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

贊一下

回復(fù)此樓

3樓2012-04-02 12:12:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 5 個(gè)回答

wll778824

銀蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1013.4
散金: 242
紅花: 1
帖子: 273
在線: 150.9小時(shí)
蟲(chóng)號(hào): 917314
注冊(cè): 2009-12-01
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
小雨萌萌: 金幣+5, 幫發(fā)放金幣咯~ 2012-04-05 16:34:13

給你一個(gè)有限體積算法求解Burgers  方程的fortran程序，源代碼用固定格式寫(xiě)的
!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

以上為主程序，以下為子程序

贊一下

回復(fù)此樓

2樓2012-04-02 12:12:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wll778824

銀蟲(chóng) (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1013.4
散金: 242
紅花: 1
帖子: 273
在線: 150.9小時(shí)
蟲(chóng)號(hào): 917314
注冊(cè): 2009-12-01
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

【答案】應(yīng)助回帖

★ ★ ★
小雨萌萌: 金幣+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

贊一下

回復(fù)此樓

4樓2012-04-02 12:13:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xmch2011

鐵蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 184
散金: 10
帖子: 75
在線: 45.2小時(shí)
蟲(chóng)號(hào): 1264781
注冊(cè): 2011-04-13
性別: GG
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小雨萌萌: 金幣-10, 謝人家也不把懸賞的金幣發(fā)一下呀？幫你發(fā)放金幣！ 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) ...

謝謝你啊，學(xué)習(xí)一下

贊一下

回復(fù)此樓

數(shù)學(xué)，力學(xué)，物理

5樓2012-04-03 11:35:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 5 個(gè)回答

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +3	酥酥魚(yú).. 2026-03-21	3/150	2026-03-22 11:31 by 楊楊楊紫
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	4/200	2026-03-22 11:04 by 搏擊518
[考研] 0856材料專碩353求調(diào)劑 +4	NIFFFfff 2026-03-20	4/200	2026-03-22 09:49 by 2026paper
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 二本跨考鄭大材料306英一數(shù)二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 22408 344分求調(diào)劑一志愿華電計(jì)算機(jī)技術(shù) +4	solanXXX 2026-03-20	4/200	2026-03-20 23:49 by alg094825
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 317求調(diào)劑 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	3/150	2026-03-20 20:35 by 學(xué)員8dgXkO
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無(wú)懈可擊111
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu