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

調(diào)劑小程序

登錄注冊

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

24小時熱門版塊排行榜

>論壇更新日志 (4259)
>考研 (801)
>導師招生 (483)
>蟲友互識 (163)
>基金申請 (125)
>文獻求助 (113)
>休閑灌水 (102)
>招聘信息布告欄 (52)
>考博 (50)
>博后之家 (36)
>碩博家園 (29)
>論文投稿 (21)
>論文道賀祈福 (20)
>公派出國 (13)
>教師之家 (9)
>找工作 (9)

返回列表

zouji1983

金蟲 (正式寫手)

應助: 6 (幼兒園)
金幣: 567.4
散金: 2
紅花: 3
帖子: 323
在線: 550.6小時
蟲號: 445773
注冊: 2007-10-28
性別: GG
專業(yè): 材料科學

[求助] fortran error求助

當我在debug mode下debug我的fortran code時，出現(xiàn)erro‘user breakpoint called from code 07c90120e’，同時在黑屏里出現(xiàn)'run-time error M6201:MATH-floor

OMAIN error'. 網(wǎng)上查了一些資源，也沒有找到解，請高手指點，多謝！

回復此樓

» 猜你喜歡

初試 317 已經(jīng)有3人回復
279分求調(diào)劑一志愿211 已經(jīng)有15人回復
286求調(diào)劑已經(jīng)有4人回復
279求調(diào)劑已經(jīng)有4人回復
求調(diào)劑已經(jīng)有3人回復
求調(diào)劑已經(jīng)有3人回復
299求調(diào)劑已經(jīng)有3人回復
330求調(diào)劑0854 已經(jīng)有3人回復
299求調(diào)劑已經(jīng)有6人回復
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有5人回復

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

急求fortran運行錯誤原因，在線等已經(jīng)有7人回復
請教 fortran 運行錯誤的原因已經(jīng)有13人回復
【求助】Fortran語言賦值問題？已經(jīng)有3人回復
【求助】求助：安裝abinit-6.2.3出現(xiàn)錯誤已經(jīng)有12人回復
【求助】vasp安裝問題已經(jīng)有27人回復
【求助】有很多錯誤的程序，找不出來了已經(jīng)有18人回復
【求助】Fortran 如何讀取文件中特定行的內(nèi)容已經(jīng)有16人回復
【求助】VASP 編譯出錯已經(jīng)有24人回復
【求助】運行錯誤，大家看看問題在哪？已經(jīng)有7人回復
【交流】Fortran語言答疑專帖已經(jīng)有157人回復

1樓 2012-09-12 22:47: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)化學

【答案】應助回帖

感謝參與，應助指數(shù) +1

貼代碼吧

出錯信息，不同的編譯器是不同的。

贊一下

回復此樓

2樓2012-09-13 15:00:20

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

zouji1983

金蟲 (正式寫手)

應助: 6 (幼兒園)
金幣: 567.4
散金: 2
紅花: 3
帖子: 323
在線: 550.6小時
蟲號: 445773
注冊: 2007-10-28
性別: GG
專業(yè): 材料科學

嗯。謝謝！代碼如下：
module thermodynamics

use global
use numlib

implicit none

contains

!*******************************************************************************************
! SUBROUTINE CONVERT_COMPOSITION
! This subroutine calculates the local phase concentrations at point (ix,iy) based on the
! equality of diffusion potential of that component in all phases and a mass balance.
!
! STATUS : ONLY FOR UP TO THREE COMPONENTS
!*******************************************************************************************
subroutine convert_composition(ix,iy)
integer :: ix, iy, i
real(8) :: r
real(8), dimension(1:np*(nc-1)) :: sol, RT
real(8), dimension(1:np*(nc-1),1:np*(nc-1)) :: coeff
integer, dimension(1:np*(nc-1)) :: IPIV

coeff = 0.d0

coeff(1,1) = dmudx(1,1,1,ix,iy)
coeff(1,2) = -dmudx(1,1,2,ix,iy)
coeff(1,3) = dmudx(1,2,1,ix,iy)
coeff(2,1) = dmudx(1,2,1,ix,iy)
coeff(2,3) = dmudx(2,2,1,ix,iy)
coeff(2,4) = -dmudx(2,2,2,ix,iy)
coeff(3,1) = phi(1,ix,iy)
coeff(3,2) = phi(2,ix,iy)
coeff(4,3) = phi(1,ix,iy)
coeff(4,4) = phi(2,ix,iy)

RT(1) = dmudx(1,1,1,ix,iy)*xp(1,1,ix,iy) + dmudx(1,2,1,ix,iy)*xp(2,1,ix,iy) - mu(1,1,ix,iy) - dmudx(1,1,2,ix,iy)*xp(1,2,ix,iy) - dmudx(1,2,2,ix,iy)*xp(2,2,ix,iy) + mu(1,2,ix,iy)
RT(2) = dmudx(2,2,1,ix,iy)*xp(2,1,ix,iy) + dmudx(1,2,1,ix,iy)*xp(1,1,ix,iy) - mu(2,1,ix,iy) - dmudx(2,2,2,ix,iy)*xp(2,2,ix,iy) - dmudx(1,2,2,ix,iy)*xp(1,2,ix,iy) + mu(2,2,ix,iy)
RT(3) = x(1,ix,iy,old)
RT(4) = x(2,ix,iy,old)

sol = RT(

call ludcmp(coeff,np*(nc-1),np*(nc-1),IPIV,r)
call lubksb(coeff,np*(nc-1),np*(nc-1),IPIV,sol)

xp(1,1:2,ix,iy) = sol(1:2)
xp(2,1:2,ix,iy) = sol(3:4)

end subroutine convert_composition

!*******************************************************************************************
! FUNCTION DMUDX
! This function retrieves the derivative of the diffusion potential mu_{iu}^{ip} with
! respect to the concentration variable x_{iv}^{ip}
!
! STATUS : FINAL FOR 3-COMPONENT SYSTEM --> only i,j
!*******************************************************************************************
function dmudx(iu,iv,ip,ix,iy)
integer :: iu, iv, ip, ix, iy, i, j
real(8) :: dmudx
i = floor(xp(1,ip,ix,iy)/thermo_dx+0.5d0)+1 ! x1
j = floor(xp(2,ip,ix,iy)/thermo_dx+0.5d0)+1 ! x2

if (ip == 2) then
! Wollastonite @ 1320C
if (iu == iv) then
dmudx = 2.d11
else
dmudx = 0.d0
end if
else
if (iu == iv) then
dmudx = (((A(iu,thermo_p(ip),i+1,j+1)-A(iu,thermo_p(ip),i,j+1))/(thermo_dx)*(xp(1,ip,ix,iy)-(i-1)*thermo_dx)+A(iu,thermo_p(ip),i,j+1))-((A(iu,thermo_p(ip),i+1,j)-A(iu,thermo_p(ip),i,j))/(thermo_dx)*(xp(1,ip,ix,iy)-(i-1)*thermo_dx)+A(iu,thermo_p(ip),i,j)))/thermo_dx*(xp(2,ip,ix,iy)-(j-1)*thermo_dx)+((A(iu,thermo_p(ip),i+1,j)-A(iu,thermo_p(ip),i,j))/(thermo_dx)*(xp(1,ip,ix,iy)-(i-1)*thermo_dx)+A(iu,thermo_p(ip),i,j))
else
dmudx = (((B(iu,iv,thermo_p(ip),i+1,j+1)-B(iu,iv,thermo_p(ip),i,j+1))/(thermo_dx)*(xp(1,ip,ix,iy)-(i-1)*thermo_dx)+B(iu,iv,thermo_p(ip),i,j+1))-((B(iu,iv,thermo_p(ip),i+1,j)-B(iu,iv,thermo_p(ip),i,j))/(thermo_dx)*(xp(1,ip,ix,iy)-(i-1)*thermo_dx)+B(iu,iv,thermo_p(ip),i,j)))/thermo_dx*(xp(2,ip,ix,iy)-(j-1)*thermo_dx)+((B(iu,iv,thermo_p(ip),i+1,j)-B(iu,iv,thermo_p(ip),i,j))/(thermo_dx)*(xp(1,ip,ix,iy)-(i-1)*thermo_dx)+B(iu,iv,thermo_p(ip),i,j))
end if
end if

end function dmudx

!*******************************************************************************************
! FUNCTION MU
! This function retrieves the diffusion potential for component k in phase i from linear
! interpolation from the descretized tables for the chemical potential
! (CHPOT in the program and input/MU_.txt files)
!
! STATUS : FINAL FOR 3-COMPONENT SYSTEM --> only i,j
!*******************************************************************************************
function mu(iu,ip,ix,iy)
integer :: iu, ip, ix, iy, i, j
real(8) :: mu
i = floor(xp(1,1,ix,iy)/thermo_dx+0.5d0)+1 ! x1
j = floor(xp(2,1,ix,iy)/thermo_dx+0.5d0)+1 ! x2

if (ip == 2) then
! Wollastonite @ 1320C
if (iu == 1) then
mu = dmudx(1,1,2,ix,iy)*xp(1,ip,ix,iy) - .8915873511e11
else
mu = dmudx(2,2,2,ix,iy)*xp(2,ip,ix,iy) - .3652507619e11
end if
else
mu = (((CHPOT(iu,thermo_p(1),i+1,j+1)-CHPOT(iu,thermo_p(1),i,j+1))/(thermo_dx)*(xp(1,1,ix,iy)-(i-1)*thermo_dx)+CHPOT(iu,thermo_p(1),i,j+1))-((CHPOT(iu,thermo_p(1),i+1,j)-CHPOT(iu,thermo_p(1),i,j))/(thermo_dx)*(xp(1,1,ix,iy)-(i-1)*thermo_dx)+CHPOT(iu,thermo_p(1),i,j)))/thermo_dx*(xp(2,1,ix,iy)-(j-1)*thermo_dx)+((CHPOT(iu,thermo_p(1),i+1,j)-CHPOT(iu,thermo_p(1),i,j))/(thermo_dx)*(xp(1,1,ix,iy)-(i-1)*thermo_dx)+CHPOT(iu,thermo_p(1),i,j))
end if

end function mu

!*******************************************************************************************
! FUNCTION F
! This function retrieves the gibbs energy for phase i from linear interpolation from
! the descretized tables for the gibbs free energy
!
! STATUS : FINAL FOR 3-COMPONENT SYSTEM --> only i,j
!*******************************************************************************************
function f(ip,ix,iy)
integer :: ip, ix, iy, i, j
real(8) :: f
i = floor(xp(1,ip,ix,iy)/thermo_dx+0.5d0)+1 ! x1
j = floor(xp(2,ip,ix,iy)/thermo_dx+0.5d0)+1 ! x2

if (ip == 2) then
! Wollastonite @ 1320C
f = dmudx(1,1,2,ix,iy)/2.d0*xp(1,ip,ix,iy)**2+dmudx(2,2,2,ix,iy)/2.d0*xp(2,ip,ix,iy)**2-.8915873511e11*xp(1,ip,ix,iy)-.3652507619e11*xp(2,ip,ix,iy) - .2093397165e11
else
f = (((GIBBS(thermo_p(ip),i+1,j+1)-GIBBS(thermo_p(ip),i,j+1))/(thermo_dx)*(xp(1,ip,ix,iy)-(i-1)*thermo_dx)+GIBBS(thermo_p(ip),i,j+1))-((GIBBS(thermo_p(ip),i+1,j)-GIBBS(thermo_p(ip),i,j))/(thermo_dx)*(xp(1,ip,ix,iy)-(i-1)*thermo_dx)+GIBBS(thermo_p(ip),i,j)))/thermo_dx*(xp(2,ip,ix,iy)-(j-1)*thermo_dx)+((GIBBS(thermo_p(ip),i+1,j)-GIBBS(thermo_p(ip),i,j))/(thermo_dx)*(xp(1,ip,ix,iy)-(i-1)*thermo_dx)+GIBBS(thermo_p(ip),i,j))
end if

end function f

end module thermodynamics

贊一下

回復此樓

3樓2012-09-13 15:37:41

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

zouji1983

金蟲 (正式寫手)

應助: 6 (幼兒園)
金幣: 567.4
散金: 2
紅花: 3
帖子: 323
在線: 550.6小時
蟲號: 445773
注冊: 2007-10-28
性別: GG
專業(yè): 材料科學

代碼比較多，謝謝！

回復此樓

4樓2012-09-13 15:38:45

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

相關版塊跳轉(zhuǎn) 我要訂閱樓主 zouji1983 的主題更新

返回列表

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

[考研] 286求調(diào)劑 +4	Faune 2026-03-21	4/200	2026-03-21 13:22 by zhukairuo
[考研] 299求調(diào)劑 +3	某某某某位 2026-03-21	3/150	2026-03-21 13:11 by zhukairuo
[考研] 一志愿華中科技大學071000，求調(diào)劑 +3	沿岸有貝殼6 2026-03-21	3/150	2026-03-21 10:35 by 暮云清寒
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 藥學383 求調(diào)劑 +3	藥學chy 2026-03-15	5/250	2026-03-20 22:11 by 云游重陽
[考研] 329求調(diào)劑 +9	想上學吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 本人考085602 化學工程專碩 +19	不知道叫什么！ 2026-03-15	21/1050	2026-03-20 20:48 by zhukairuo
[考研] 0817 化學工程 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 學員8dgXkO
[考研] 一志愿南京航空航天大學大學，080500材料科學與工程學碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 0703化學調(diào)劑，六級已過，有科研經(jīng)歷 +13	曦熙兮 2026-03-15	13/650	2026-03-20 19:35 by Dream007008
[考研] 08工學調(diào)劑 +5	用戶573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 梁成偉老師課題組歡迎你的加入 +9	一鴨鴨喲 2026-03-14	11/550	2026-03-19 17:22 by ！本暗一次！
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 288求調(diào)劑 +4	奇點0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky

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

24小時熱門版塊排行榜

zouji1983

[求助] fortran error求助

» 猜你喜歡

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

snoopyzhao

【答案】應助回帖

zouji1983

zouji1983

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