 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1065 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
xmch2011鐵蟲 (小有名氣)
|
[求助]
求fortran95編寫的數(shù)值程序
|
||
| 想學(xué)學(xué)用fortran95寫的數(shù)值程序，那位同學(xué)有沒(méi)有這樣的程序，學(xué)習(xí)一下，謝謝！ |
回復(fù)此樓» 猜你喜歡
材料求調(diào)劑
已經(jīng)有5人回復(fù)
-
354求調(diào)劑
已經(jīng)有9人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有22人回復(fù)
-
考研調(diào)劑
已經(jīng)有3人回復(fù)
-
326求調(diào)劑
已經(jīng)有8人回復(fù)
-
333求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑
已經(jīng)有3人回復(fù)
-
考研化學(xué)學(xué)碩調(diào)劑，一志愿985
已經(jīng)有7人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有15人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級(jí)已過(guò)，有科研經(jīng)歷
已經(jīng)有14人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助計(jì)算機(jī)高人，編寫個(gè)Fortran程序
已經(jīng)有13人回復(fù)
- 【求助】求助用Fortran編寫模擬中的接觸角
已經(jīng)有4人回復(fù)
- 【求助】如何寫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)編寫的fortran程序讓別人使用，但不讓別人知道源代碼
已經(jīng)有6人回復(fù)
wll778824
銀蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1013.4
- 散金: 242
- 紅花: 1
- 帖子: 273
- 在線: 150.9小時(shí)
- 蟲號(hào): 917314
- 注冊(cè): 2009-12-01
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
【答案】應(yīng)助回帖
★ ★ ★
小雨萌萌: 金幣+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
wll778824
銀蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1013.4
- 散金: 242
- 紅花: 1
- 帖子: 273
- 在線: 150.9小時(shí)
- 蟲號(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
感謝參與，應(yīng)助指數(shù) +1
小雨萌萌: 金幣+5, 幫發(fā)放金幣咯~ 2012-04-05 16:34:13
|
給你一個(gè)有限體積算法求解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
銀蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1013.4
- 散金: 242
- 紅花: 1
- 帖子: 273
- 在線: 150.9小時(shí)
- 蟲號(hào): 917314
- 注冊(cè): 2009-12-01
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
【答案】應(yīng)助回帖
★ ★ ★
小雨萌萌: 金幣+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
xmch2011
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 184
- 散金: 10
- 帖子: 75
- 在線: 45.2小時(shí)
- 蟲號(hào): 1264781
- 注冊(cè): 2011-04-13
- 性別: GG
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
5樓2012-04-03 11:35:00
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 326求調(diào)劑 +4 | mlpqaz03 2026-03-15 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 本人考085602 化學(xué)工程 專碩 +20 | 不知道叫什么！ 2026-03-15 | 22/1100 |
|
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +6 | Liyouyumairs 2026-03-21 | 6/300 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！ +4 | lemonzzn 2026-03-17 | 7/350 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +3 | 沿岸有貝殼6 2026-03-21 | 3/150 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|