24小时热门版块排行榜    

查看: 817  |  回复: 5

欢迎你来了

新虫 (小有名气)

[求助] 运行结果不随输入变化

结果可以运行,但是结果不会随着输入变化,这是为什么呀?

运行结果不随输入变化


发自小木虫Android客户端
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

raymond137

金虫 (正式写手)

2楼2021-05-21 07:06:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

欢迎你来了

新虫 (小有名气)

module contants    implicit none    real::a1, a2    real::ed = 9.91, edm = 14.5    real::b2 = 2.3    real::kb = 1.380649e-23    real::a_i = 0.08    real::a_j = 4.33d-5    real(kind=8)::h = 1.06e-34    real(kind=8)::w = 2.37e3    integer::b1 = 3390end module contants     real function gos(x) use contants implicit none real:: x real(kind=4):: m, m1, m2 real(kind=4):: em1 real::a_1, b_1 integer::j, v integer::tv, t0, trot, tread *,x,tv,t0,trot,t a_1 = a1 - (b1/trot) - (b1/t0)    a1 = -3*a_i*kb*t/(2*ed)    a2 = -3*a_j*kb*t/(2*ed) if (v >= 0 .and. v < 9) then  m1 = 0  m2 = 9  b_1 = b1  m = 1  em1 = (m1 + 0.5)*h*w  gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) else if (v >= 9 .and. v < 31) then  m1 = 9  m2 = 31  b_1 = 0.75*b1  m = 2  em1 = (m1 + 0.5)*h*w  gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) else if (v >= 31 .and. v < 55) then  m1 = 31  m2 = 55  b_1 = 0.45*b1  m = 3  em1 = (m1 + 0.5)*h*w  gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) end if! fl=q+d+c write(*,*) gos return    end function    program piok !use contants implicit none real::y,x  real,external::gos  y=gos(x) ! print*,gos  end

发自小木虫Android客户端
3楼2021-05-21 11:18:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

欢迎你来了

新虫 (小有名气)

引用回帖:
2楼: Originally posted by raymond137 at 2021-05-21 07:06:37
贴出源代码

module contants  
  implicit none  
  real::a1, a2  
  real::ed = 9.91, edm = 14.5   
real::b2 = 2.3   
real::kb = 1.380649e-23  
  real::a_i = 0.08   
real::a_j = 4.33d-5   
real(kind=8)::h = 1.06e-34   
real(kind=8)::w = 2.37e3  
  integer::b1 = 3390
end module contants   
real function gos(x)
use contants
implicit none real:: x real(kind=4):: m, m1, m2 real(kind=4):: em1
real::a_1, b_1
integer::j, v
integer::tv, t0, trot, t
read *,x,tv,t0,trot,t a_1 = a1 - (b1/trot) - (b1/t0)   
a1 = -3*a_i*kb*t/(2*ed)
   a2 = -3*a_j*kb*t/(2*ed)
if (v >= 0 .and. v < 9) then
m1 = 0  m2 = 9  b_1 = b1  m = 1  em1 = (m1 + 0.5)*h*w
gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1)))
else if (v >= 9 .and. v < 31) then
m1 = 9  m2 = 31  
b_1 = 0.75*b1
m = 2  
em1 = (m1 + 0.5)*h*w  
gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) else if (v >= 31 .and. v < 55) then  
m1 = 31
m2 = 55
  b_1 = 0.45*b1
  m = 3
em1 = (m1 + 0.5)*h*w
gos= exp(x*em1 + a_1*m1)*((1 - exp((m2 - m1)*(x*b_1*kb + a_1)))/(1 - exp(x*b_1*kb + a_1))) end if
! fl=q+d+c
write(*,*) gos
return   
end function   

program piok
!use contants
implicit none
real::y,x  
real,external::gos
y=gos(x)
! print*,gos  
end

发自小木虫Android客户端
4楼2021-05-21 11:26:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

raymond137

金虫 (正式写手)

引用回帖:
4楼: Originally posted by 欢迎你来了 at 2021-05-21 11:26:32
module contants  
  implicit none  
  real::a1, a2  
  real::ed = 9.91, edm = 14.5   
real::b2 = 2.3   
real::kb = 1.380649e-23  
  real::a_i = 0.08   
real::a_j = 4.33d-5   
real(kind=8)::h = 1 ...

我是做C++的。你这个是什么语言代码?Fortran?建议你用debug模式,加几个断点,分步跟踪。

发自小木虫Android客户端
5楼2021-05-21 11:53:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

欢迎你来了

新虫 (小有名气)

6楼2021-05-21 11:59:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 欢迎你来了 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 每到中夜,情难自抑 +37 sioc-sunj 2024-05-28 61/3050 2024-06-03 15:16 by EllieWU
[硕博家园] 讨厌酒桌文化 +9 hahamyid 2024-05-31 9/450 2024-06-03 14:16 by 北京炸酱面
[硕博家园] 实验室太吵闹,无法安静学习,怎么办? +11 utahh 2024-05-31 22/1100 2024-06-03 11:18 by utahh
[基金申请] 化学B02口青基 代表作都是什么水平的?向大佬求助 +14 arthas_007 2024-06-01 19/950 2024-06-03 11:14 by prisonrabbit
[基金申请] 离职后国自然项目不能变更单位,在新单位还能申请新的国自然项目吗 +3 宜兰 2024-05-31 11/550 2024-06-03 10:26 by lilyxjw
[论文投稿] Scientific reports 投稿 +4 chencome12 2024-06-02 6/300 2024-06-03 09:13 by chencome12
[论文投稿] 求助大神,Fe和Al离子对MOF都有淬灭,当两种离子共存时,怎么区分两种离子? 10+3 maoxiao 2024-06-02 3/150 2024-06-03 01:43 by nono2009
[基金申请] +5 河马の史诗 2024-06-02 5/250 2024-06-02 19:05 by 328838485
[教师之家] 选择 +8 459582015 2024-05-31 11/550 2024-06-02 18:30 by auvauv
[硕博家园] 各位同学能否分享一下实验室的学生劳务发放标准呀? +19 ma3252788 2024-05-30 19/950 2024-06-02 18:15 by Lcy6666
[基金申请] 讨论贴:好本子一定会被打招呼的淘汰吗 +40 ioe3937 2024-05-27 51/2550 2024-06-02 13:52 by Howard28
[论文投稿] 编辑对文章提的意见 40+3 qa8386136 2024-06-01 3/150 2024-06-02 13:28 by bobvan
[论文投稿] 求Sci期刊推荐 10+4 甄小鱼 2024-05-30 7/350 2024-06-01 10:41 by bobvan
[基金申请] 博后特别资助状态变化 +24 随梦而飞2017 2024-05-30 35/1750 2024-06-01 10:10 by 青岛阳仔
[教师之家] 中年 (金币+3) +18 459582015 2024-05-28 19/950 2024-06-01 00:41 by 沈婉婷.Girl
[文学芳草园] 物是人非 +4 myrtle 2024-05-30 4/200 2024-05-30 15:05 by mapenggao
[论文投稿] 审稿专家比较坚定的让补充实验,但实在没法补充实验,修回还有希望吗? (EPI+1) 3+3 qweasd12345 2024-05-29 6/300 2024-05-30 08:11 by qweasd12345
[有机交流] 液相纯度高,但产品析不太出来 10+4 cui19236 2024-05-27 9/450 2024-05-30 07:45 by yuanjijoy
[论文投稿] 核心初审被拒,理由是“选题的意义不明确,文章写得不像是科技论文”,怎么改 5+3 工藤雷花樱 2024-05-27 8/400 2024-05-29 10:09 by topedit
[基金申请] 如果您是国自然评审专家 +3 丁香园账户 2024-05-28 3/150 2024-05-29 06:44 by gaohui8888
信息提示
请填处理意见