[問題] 解暫態方程式
我要使用差分法解一個波動方程式
假設是二維的
這樣要解的是u
所以宣告u(100,100,100)
前兩個分別是X方向和Y方向的網格數,後面那個是t
也就是u(x,y,t)
假設我給dt是0.01
總時間長T給定1
這樣就有了100格,一切都很順利
我想改的問題來了~
今天我想寫一個副程式sub來解上面的波動方程式
但是時間方面我想改成以下這樣
do it=1, 1000
call sub
end do
it每跑一步就是跑了一個dt
我想讓他跑1000步,也就是不會被限制宣告那邊就要給他設定t的網格數
所以就不用設定時間總長
請問我該如何做呢?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.117.58.27
→
04/26 21:43, , 1F
04/26 21:43, 1F
→
04/26 21:44, , 2F
04/26 21:44, 2F
→
04/26 21:45, , 3F
04/26 21:45, 3F
→
04/26 21:46, , 4F
04/26 21:46, 4F
→
04/26 21:46, , 5F
04/26 21:46, 5F
→
04/26 21:49, , 6F
04/26 21:49, 6F
→
04/26 21:49, , 7F
04/26 21:49, 7F
→
04/26 21:51, , 8F
04/26 21:51, 8F
→
04/26 21:54, , 9F
04/26 21:54, 9F
D大你說到我的問題點了,第一種的話 太容易超出過記憶體
第二種時間run太長也不行,照我目前的需求,時間可能至少需要10萬格
兩種方法我都有試,但是都無法順利解決
推
04/26 22:16, , 10F
04/26 22:16, 10F
→
04/26 22:16, , 11F
04/26 22:16, 11F
→
04/26 23:18, , 12F
04/26 23:18, 12F
t大說到我的需求了,我並不想把所有的t網格都存下來
我想改成只看需要的那幾個dt就夠了
ex.設定跑1000步,每100步輸出資料
推
04/27 09:44, , 13F
04/27 09:44, 13F
※ 編輯: lprsdmn 來自: 140.117.58.27 (04/27 10:39)
推
04/27 11:04, , 14F
04/27 11:04, 14F
→
04/27 11:05, , 15F
04/27 11:05, 15F
→
04/27 11:06, , 16F
04/27 11:06, 16F
→
04/27 11:06, , 17F
04/27 11:06, 17F
t大可以幫我看看我程式該如何修改嗎
你真的了解
program wave
implicit none
real(8) u(0:101,0:101,0:10001)
real(8) lamda,vel,c,dx,dy,dt,x,y,z,tfinal
integer i,j,k,nxstep,nystep,ntime,L
real, parameter :: pi=3.14159
c=2.9d8
L=1
tfinal=2.1d-9
dx=0.01*L
dy=0.01*L
dt=2.1d-13
nxstep=(L/dx)
nystep=(L/dy)
ntime=int(tfinal/dt)+1
!b.c
do k=0,ntime
do j=0,nystep
u(0,j,k)=0.
u(nxstep,j,k)=0.
end do
end do
do k=0,ntime
do i=0,nxstep
u(i,0,k)=0
u(i,nystep,k)=0
end do
end do
!i.c
do i=1,nxstep-1
do j=1,nystep-1
x=dx*i
y=dy*j
u(i,j,0)=sin(pi*x)*sin(pi*y)
u(i,j,1)=sin(pi*x)*sin(pi*y)
end do
end do
!==============================
do i=1,nxstep-1
do j=1,nystep-1
do k=1,ntime
u(i,j,k+1)=2*u(i,j,k)-u(i,j,k-1)+c*c*((dt/dx)**2*(u(i+1,j,k)-2*u(i,j,k)&
+u(i-1,j,k))+(dt/dy)**2*(u(i,j+1,k)-2*u(i,j,k)+u(i,j-1,k)))+100*dt*dt
end do
end do
end do
!=========================================================
open (7, file='ya1.plt')
open (8, file='ya2.plt')
open (9, file='ya3.plt')
do i=0,nxstep
write(7,*) i , u(i,50,5)
write(8,*) i , u(i,50,20)
write(9,*) i , u(i,50,100)
end do
stop
end
※ 編輯: lprsdmn 來自: 140.117.58.27 (04/27 14:54)
討論串 (同標題文章)
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
