版塊導航: 正在加載中...

調(diào)劑小程序

登錄注冊

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

24小時熱門版塊排行榜

>論壇更新日志 (1444)
>考研 (153)
>導師招生 (132)
>論文投稿 (40)
>蟲友互識 (34)
>博后之家 (30)
>公派出國 (23)
>基金申請 (22)
>教師之家 (22)
>考博 (22)
>碩博家園 (11)
>找工作 (10)
>論文道賀祈福 (9)
>材料工程 (3)
>簽證指南 (3)
>材料綜合 (2)

返回列表

曾經(jīng)電工

金蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 2090.4
散金: 68
帖子: 119
在線: 224.2小時
蟲號: 582687
注冊: 2008-07-25
性別: GG
專業(yè): 燃燒學

[交流] 【求助】FORTRAN中增加計算區(qū)域后出錯，已有2人參與

開始我計算點為50×100，程序能夠正常編譯和運行，結(jié)果也還是可以。但是當計算點增加到500×100時，可以運行，而出來的結(jié)果完全不正確。看了一下數(shù)據(jù)，發(fā)現(xiàn)是在計算濃度場中，增加計算點后其計算值出現(xiàn)負數(shù)。不知怎么回事？另外還有一個問題就是計算過程中本來值為1.000000000000000，但是顯示結(jié)果卻為1.00000000000001，再繼續(xù)運行下去的話，會出現(xiàn)累計現(xiàn)象，從而影響結(jié)果的正確性。參數(shù)的精度我都定義為雙精度的，這是怎么一個情況？

回復此樓

» 猜你喜歡

284求調(diào)劑已經(jīng)有10人回復
一志愿山東大學藥學學碩求調(diào)劑已經(jīng)有4人回復
07化學280分求調(diào)劑已經(jīng)有4人回復
298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑已經(jīng)有12人回復
求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復
335求調(diào)劑已經(jīng)有5人回復
求調(diào)劑已經(jīng)有7人回復
一志愿吉大化學322求調(diào)劑已經(jīng)有4人回復
環(huán)境學碩288求調(diào)劑已經(jīng)有8人回復
341求調(diào)劑(一志愿湖南大學070300) 已經(jīng)有6人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

急求fortran運行錯誤原因，在線等已經(jīng)有7人回復
請教 fortran 運行錯誤的原因已經(jīng)有13人回復
fortran轉(zhuǎn)換CHGCAR出錯已經(jīng)有4人回復
FORTRAN新手求助主程序循環(huán)問題已經(jīng)有10人回復
我編的Simpson積分法fortran程序給不出結(jié)果，大俠們看看哪里出了問題？已經(jīng)有4人回復
請教一個fortran小程序編譯出錯的問題，謝謝已經(jīng)有9人回復
【求助】Fortran語言賦值問題？已經(jīng)有3人回復
Fortran的格式化輸入輸出問題已經(jīng)有14人回復
有關fortran的一次而問題，希望大家能幫幫忙，謝謝已經(jīng)有4人回復
【求助】用fortran解矩陣問題【已解決】已經(jīng)有5人回復

1樓 2010-09-14 09:23:25

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學

★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
余澤成(金幣+1):謝謝參與應助！ 2010-09-14 16:02:06
曾經(jīng)電工(金幣+5):謝謝您的幫助！ 2010-09-14 16:07:21
曾經(jīng)電工: 回帖置頂 2011-12-06 09:50:19

第一個問題，無源碼無真相。

第二個問題，基本上不用擔心，所有的浮點計算都有一定的誤差，但這點誤差基本上可以忽略不計，我想

贊一下(2人)

回復此樓

2樓2010-09-14 14:20:36

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

曾經(jīng)電工

金蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 2090.4
散金: 68
帖子: 119
在線: 224.2小時
蟲號: 582687
注冊: 2008-07-25
性別: GG
專業(yè): 燃燒學

引用回帖:

Originally posted by snoopyzhao at 2010-09-14 14:20:36:
第一個問題，無源碼無真相。

第二個問題，基本上不用擔心，所有的浮點計算都有一定的誤差，但這點誤差基本上可以忽略不計，我想

謝謝您的回答，我的濃度計算程序如下，麻煩您幫我看一下是否是程序中的問題呢？我認為這種浮點計算的誤差是在積累，最后都影響到結(jié)果的正確性了。
module const
implicit none
integer,parameter :: N=100
integer,parameter :: L=50
integer,parameter :: tmax=900000
double precision,parameter :: Dis=2.0d-3
double precision,parameter :: dx=4.0d-6
double precision,parameter :: Ds=1.68d-9
double precision,parameter :: Cs0=10.8d0
double precision,parameter :: Cxm=20.0d0
double precision,parameter :: Yxs=0.61d0
double precision,parameter :: Ks=5.204d0
double precision,parameter :: dts=1.0d0
end module const

!=====================================
!             主程序
!=====================================
   program main
   use const
   implicit none
   integer :: it,x,z,ci(1:N,1:L)
   double precision :: Sd(1:N,1:L),Cd(1:N,1:L)

   call init_Field(Sd,Cd,ci)

   do it=1,tmax

      call substrate(Sd,Cd,ci)

   end do

stop
end

   subroutine init_Field(Sd,Cd,ci)
   use const
   implicit none
   double precision :: a,Sd(1:N,1:L),Cd(1:N,1:L)
   integer :: ci(1:N,1:L)
   integer :: x,z

   do x=1,N
      do z=1,L
      ci(x,z)=0
         Sd(x,z)=1.0d0
         Cd(x,z)=0.0d0
      end do
   end do

   call random_seed()
   do x=1,N
      z=1
   call random_number(a)
      Cd(x,z)=a
      if (Cd(x,z)<=0.5d0) then
         ci(x,z)=1
      else
      Cd(x,z)=0.0d0
   end if
   end do

   return
   end
!=========================================================================================
! 濃度場迭代計算: 在每個坐標方向上分別采用追趕法進行求解三對角方程組，
! 其他方向則按照顯示處理的交替方向隱式方法ADI(alternating direction implicit)
!=========================================================================================
   subroutine substrate(Sd,Cd,ci)
   use const
   implicit none
   integer :: x,z,ci(1:N,1:L),Zm,Zb,k(1:L)
   double precision :: Sd(1:N,1:L),Cd(1:N,1:L),Ud(1:N,1:L),Md(1:N,1:L),a1(1:N,1:L),a2(1:N,1:L),a3(1:N,1:L),aa1(1:N,1:L),aa2(1:N,1:L),aa3(1:N,1:L)  ! a1,a2,a3,b1為三對角陣方程系數(shù)(x方向)；ps為底物消耗量
   double precision :: ps2(1:N,1:L),f(1:N,1:L),fe(1:N,1:L),y(1:N,1:L),m(1:N,1:L),me(1:N,1:L),r(1:N,1:L),tt,kk,I,b1(1:N,1:L),b2(1:N,1:L),ps1(1:N,1:L)
   double precision :: um,ms

   ms=0.562137d0
   um=0.25986d0
   Zm=0
   k(1:L)=0
   do z=1,L
      do x=1,N
      k(z)=k(z)+ci(x,z)
   end do
   if (k(z)==0) then
Zm=z-1
goto 300
   else
      continue
   end if
   end do

300 Zb=Zm+10

   do x=1,N
      do z=Zb+1,L
      Sd(x,z)=1.0d0
   end do
   end do

   do x=1,N
      do z=1,L
      Ud(x,z)=Sd(x,z)
   end do
   end do


   !!! X-direction
   tt=um
   kk=ms
   a1(1:N,1:Zb)=-1.0d0
   a3(1:N,1:Zb)=-1.0d0
   a2(1:N,1:Zb)=2.0d0+(2.0d0*(Dis**2)*(dx**2))/(dts*Ds)
   do z=2,Zb
      do x=2,N-1
      b1(x,z)=Sd(x,z+1)+((2.0d0*(Dis**2)*(dx**2))/(dts*Ds)-2.0d0)*Sd(x,z)+Sd(x,z-1)-ps1(x,z)
ps1(x,z)=((Dis**2)*(dx**2)/Ds)*((tt*Cxm*Cd(x,z)*Sd(x,z))/(Yxs*Cs0*(Ks/Cs0+Sd(x,z)))+Cxm*Cd(x,z)*kk/Cs0)
   end do
   end do


do z=2,Zb
   f(2,z)=a3(2,z)/a2(2,z)
   fe(2,z)=a2(2,z)
      do x=3,N-2
   fe(x,z)=a2(x,z)-(a1(x,z)*f(x-1,z))
   f(x,z)=a3(x,z)/fe(x,z)
      end do
fe(N-1,z)=a2(N-1,z)-(a1(N-1,z)*f(N-2,z))
!----------------------------------------------------------
      y(2,z)=(b1(2,z)-a1(2,z)*Sd(1,z))/a2(2,z)
do x=3,N-2
   y(x,z)=(b1(x,z)-a1(x,z)*y(x-1,z))/fe(x,z)
end do
      y(N-1,z)=(b1(N-1,z)-a3(N-1,z)*Sd(N,z)-a1(N-1,z)*y(N-2,z))/fe(N-1,z)
!----------------------------------------------------------
Ud(N-1,z)=y(N-1,z)
do x=N-2,2,-1
   Ud(x,z)=y(x,z)-f(x,z)*Ud(x+1,z)
      end do
!-----------------------------------------------------------
Ud(N,z)=(4.0d0*Ud(2,z)+4.0d0*Ud(N-1,z)-Ud(N-2,z)-Ud(3,z))/6.0d0
Ud(1,z)=Ud(N,z)
end do

do x=1,N
   Ud(x,1)=(4.0d0*Ud(x,2)-Ud(x,3))/3.0d0
   Ud(1,1)=Ud(N,1)
end do

do x=1,N
   do z=1,L
      Md(x,z)=Ud(x,z)
   end do
end do



   !!! Z-direction
   aa1(1:N,1:Zb)=-1.0d0
   aa3(1:N,1:Zb)=-1.0d0
   aa2(1:N,1:Zb)=2.0d0+(2.0d0*(Dis**2)*(dx**2))/(dts*Ds)
   do x=2,N-1
      do z=2,Zb
      b2(x,z)=Md(x+1,z)+((2.0d0*(Dis**2)*(dx**2))/(dts*Ds)-2.0d0)*Md(x,z)+Md(x-1,z)-ps2(x,z)
ps2(x,z)=((Dis**2)*(dx**2)/Ds)*((tt*Cxm*Cd(x,z)*Md(x,z))/(Yxs*Cs0*(Ks/Cs0+Md(x,z)))+Cxm*Cd(x,z)*kk/Cs0)
      end do
   end do

      do x=2,N-1
         m(x,2)=aa3(x,2)/aa2(x,2)
         me(x,2)=aa2(x,2)
         do z=3,Zb-1
         me(x,z)=aa2(x,z)-(aa1(x,z)*m(x,z-1))
         m(x,z)=aa3(x,z)/me(x,z)
         end do
me(x,Zb)=aa2(x,Zb)-(aa1(x,Zb)*m(x,Zb-1))
!------------------------------------------------------------------------------------
      r(x,2)=(b2(x,2)-aa1(x,2)*Md(x,1))/aa2(x,2)
do z=3,Zb-1
   r(x,z)=(b2(x,z)-aa1(x,z)*r(x,z-1))/me(x,z)
end do
         r(x,Zb)=(b2(x,Zb)-aa3(x,Zb)*Md(x,Zb+1)-aa1(x,Zb)*r(x,Zb-1))/me(x,Zb)
!-----------------------------------------------------------------------------------
Sd(x,Zb)=r(x,Zb)
do z=Zb-1,2,-1
         sd(x,z)=r(x,z)-m(x,z)*Sd(x,z+1)
         end do
!-----------------------------------------------------------------------------------
Sd(x,1)=(4.0d0*Sd(x,2)-Sd(x,3))/3.0d0
      end do
   do z=1,Zb
      Sd(N,z)=(4.0d0*Sd(N-1,z)+4.0d0*Sd(2,z)-Sd(3,z)-Sd(N-2,z))/6.0d0
Sd(1,z)=Sd(N,z)
   end do

      return
      end

[ Last edited by 曾經(jīng)電工 on 2010-9-14 at 16:13 ]

贊一下

回復此樓

3樓2010-09-14 16:04:57

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

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學

你這個程序的輸出在哪里呢？

贊一下

回復此樓

4樓2010-09-14 19:54:31

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

曾經(jīng)電工

金蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 2090.4
散金: 68
帖子: 119
在線: 224.2小時
蟲號: 582687
注冊: 2008-07-25
性別: GG
專業(yè): 燃燒學

輸出的話，我可以在主程序里面加入的：
program main
   use const
   implicit none
   integer :: it,x,z,ci(1:N,1:L)
   double precision :: Sd(1:N,1:L),Cd(1:N,1:L)

   call init_Field(Sd,Cd,ci)

   do it=1,tmax

      call substrate(Sd,Cd,ci)
      do x=1,N
   do z=1,L
      write(*,*) x,z,Sd(x,z)
   end do
end do
   end do

stop
end

贊一下

回復此樓

5樓2010-09-15 08:35:07

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

zyj8119

木蟲 (著名寫手)

應助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
散金: 1440
紅花: 35
帖子: 2936
在線: 1329.4小時
蟲號: 664177
注冊: 2008-11-29
性別: GG
專業(yè): 理論和計算化學

★
曾經(jīng)電工(金幣+1):謝謝參與
曾經(jīng)電工(金幣+5):謝謝您的幫助！ 2010-11-10 08:21:13

引用回帖:

Originally posted by 曾經(jīng)電工 at 2010-09-14 16:04:57:

謝謝您的回答，我的濃度計算程序如下，麻煩您幫我看一下是否是程序中的問題呢？我認為這種浮點計算的誤差是在積累，最后都影響到結(jié)果的正確性了。
module const
implicit none
integer,parameter :: N=100 ...

我覺得你最好在程序里面打開一個文件，然后把所有的數(shù)據(jù)都存放在文件中，我用CVF6.6運行你的程序什么東西都沒有。

贊一下(1人)

回復此樓

好好學習，天天向上。

6樓2010-11-09 19:53:29

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

相關版塊跳轉(zhuǎn) 我要訂閱樓主曾經(jīng)電工的主題更新

返回列表

普通表情龍兔虎貓高級回復 (可上傳附件)

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

[考研] 材料調(diào)劑 +3	匹克i 2026-03-23	3/150	2026-03-23 23:18 by peike
[考研] 材料專碩英一數(shù)二306 +8	z1z2z3879 2026-03-18	8/400	2026-03-23 20:49 by baobaoye
[考研] 303求調(diào)劑 +4	元夕元 2026-03-20	4/200	2026-03-23 19:00 by macy2011
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +5	jiajunX 2026-03-22	5/250	2026-03-23 18:07 by YMU施老師
[考研] 求調(diào)劑一志愿武漢理工大學材料工程（085601） +3	WW.' 2026-03-23	5/250	2026-03-23 17:18 by 楓翼ljj
[考研] 08工學調(diào)劑 +7	用戶573181 2026-03-20	11/550	2026-03-23 15:47 by 我愛學習學習使?/a>
[考研] 0854電子信息求調(diào)劑 324 +3	Promise-jyl 2026-03-23	3/150	2026-03-23 13:43 by wangkm
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-17	10/500	2026-03-23 13:05 by 醉在風里
[考研] 070300，一志愿北航320求調(diào)劑 +3	Jerry0216 2026-03-22	5/250	2026-03-23 09:16 by 。。堂堂
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 318求調(diào)劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學員8dgXkO
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 288求調(diào)劑，一志愿華南理工大學071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

曾經(jīng)電工

[交流] 【求助】FORTRAN中增加計算區(qū)域后出錯， 已有2人參與

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

snoopyzhao

曾經(jīng)電工

snoopyzhao

曾經(jīng)電工

zyj8119

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[交流] 【求助】FORTRAN中增加計算區(qū)域后出錯，已有2人參與