版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門版塊排行榜

返回列表

xiuxialy

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 0.5
帖子: 6
在線: 6.6小時(shí)
蟲號(hào): 1997605
注冊(cè): 2012-09-13
專業(yè): 光學(xué)

[求助] fortran循環(huán)求助

program ex2
  implicit none
  integer::ln,Nrank,m,Nmax,i
  complex(kind=16),allocatable::a(

,b(

,T_mat(:,

  ln=21
  Nrank=ln
  allocate(a(ln),b(ln))
  do i=1,ln
a(i)=2*i
b(i)=3*i
  end do
  do m=0,Nrank
if (m==0)then
      Nmax=Nrank
else
      Nmax=Nrank-m+1
end if
allocate(T_mat(2*Nmax,2*Nmax))
print *,m
read(*,*)
call Tmat(Nmax,ln,a,b,T_mat)
  end do
end
!*************************************************************!
subroutine Tmat(Nmax,Nrank,a,b,T)
  implicit none
  integer::Nrank,Nmax,i,j
  complex(kind=16)::a(Nrank),b(Nrank),T(2*Nmax,2*Nmax)
  do i=1,2*Nmax
do j=1,2*Nmax
   T(i,j)=0.0
   if(i==j)then
      if(j<=Nmax) then
      T(i,j)=a(i+Nrank-Nmax)
      else
      T(i,j)=b(Nrank+j-2*Nmax)
      end if
   end if
end do
  end do
end subroutine Tmat
為什么循環(huán)中只能輸出0？

回復(fù)此樓

» 猜你喜歡

329求調(diào)劑已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有11人回復(fù)
一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一）已經(jīng)有3人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有14人回復(fù)
320求調(diào)劑0856 已經(jīng)有3人回復(fù)
307求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有14人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有18人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

FORTRAN新手求助主程序循環(huán)問題已經(jīng)有10人回復(fù)
【求助】請(qǐng)問如何解決insufficient virtual memory 已經(jīng)有10人回復(fù)
【求助】Python 中求解微分方程 odeint 的用法已經(jīng)有11人回復(fù)
【求助】ANSYS 中命令流編程（fortran）已經(jīng)有11人回復(fù)
【求助】fortran程序未出錯(cuò)，但運(yùn)行不出結(jié)果，中途中斷！已經(jīng)有6人回復(fù)
【求助】求指點(diǎn)，在Fortran里面，怎樣實(shí)現(xiàn)“數(shù)組維度可調(diào)”的數(shù)組? 已經(jīng)有5人回復(fù)
【求助】三維矩陣運(yùn)算已經(jīng)有8人回復(fù)
【求助】Fortran程序，如何讓exe文件可以單獨(dú)運(yùn)行？已經(jīng)有6人回復(fù)
【求助】Fortran 如何讀取文件中特定行的內(nèi)容已經(jīng)有16人回復(fù)
【求助】復(fù)超越方程的解法已經(jīng)有19人回復(fù)
【求助】緊急求助有關(guān)MATHEMATICA的問題已經(jīng)有9人回復(fù)

1樓 2013-06-29 14:53:31

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

pippi6

鐵桿木蟲 (著名寫手)

工程和科學(xué)數(shù)值計(jì)算咨詢

應(yīng)助: 413 (碩士)
貴賓: 0.002
金幣: 7116.5
散金: 15
紅花: 63
帖子: 1639
在線: 798.9小時(shí)
蟲號(hào): 2469437
注冊(cè): 2013-05-14
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

【答案】應(yīng)助回帖

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

因?yàn)?read(*,*) 在等待數(shù)據(jù)
另外，你在do loop 里重復(fù)做 allocate 有問題，應(yīng)把
allocate(T_mat(2*Nmax,2*Nmax))
移到 do loop 之前。
后面還有數(shù)組越界之類，問題多多，需要詳細(xì)查。

贊一下(1人)

回復(fù)此樓

2樓2013-06-29 17:12:59

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

xiuxialy

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 0.5
帖子: 6
在線: 6.6小時(shí)
蟲號(hào): 1997605
注冊(cè): 2012-09-13
專業(yè): 光學(xué)

引用回帖:

2樓: Originally posted by pippi6 at 2013-06-29 17:12:59
因?yàn)?read(*,*) 在等待數(shù)據(jù)
另外，你在do loop 里重復(fù)做 allocate 有問題，應(yīng)把
allocate(T_mat(2*Nmax,2*Nmax))
移到 do loop 之前。
后面還有數(shù)組越界之類，問題多多，需要詳細(xì)查。

謝謝，昨天發(fā)帖后我編了個(gè)小程序發(fā)現(xiàn)這個(gè)問題了，確實(shí)是這樣，不過仍然很感謝！

贊一下

回復(fù)此樓

3樓2013-06-30 11:15:58

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

xiuxialy

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 0.5
帖子: 6
在線: 6.6小時(shí)
蟲號(hào): 1997605
注冊(cè): 2012-09-13
專業(yè): 光學(xué)

PROGRAM gmm01f
   implicit double precision (a-h,o-z)
   include 'gmm01f.par'
   parameter (nmp=np*(np+2),nmp0=(np+1)*(np+4)/2)
   parameter (NXMAX=3000,nangmax=181,MOR=181,ncmax=360)
   parameter (ni0=np*(np+1)*(2*np+1)/3+np*np)
   parameter (ng0=np*(2*np**3+10*np**2+19*np+5)/6)
   parameter (nrc=4*np*(np+1)*(np+2)/3+np)
   parameter (nij=nLp*(nLp-1)/2)
   integer u,v,u0,nmax(nLp),uvmax(nLp),ind(nLp),ind2(nLp)
   double precision k,lnfacd,r0(6,nLp),x(nLp),dang(nangmax),
   + r00(3,nLp),rsr0(NXMAX),rsi0(NXMAX),
   + rsx0(NXMAX),px0(NXMAX),w1(np),w2(np),w3(np),w4(np),
   + rsr(np,nLp),rsi(np,nLp),rsx(np,nLp),px(np,nLp),betar(MOR),
   + thetr(MOR),phair(MOR),smue(4,4),mue(4,4,ncmax,nangmax),
   + besj(0:2*np+1),besy(0:2*np+1),i11(nangmax),
   + i21(nangmax),i22(nangmax),i12(nangmax),inat(nangmax),
   + pol(nangmax),cscaxi(nLp),cscayi(nLp),cextxi(nLp),
   + cextyi(nLp),cabsxi(nLp),cabsyi(nLp),cexti(nLp),cabsi(nLp),
   + cscai(nLp),assymi(nLp),assymxi(nLp),assymyi(nLp),
   + cprxi(nLp),cpryi(nLp),cpri(nLp),drot(nrc,nij),
   + c0i(nLp),c1i(nLp),confg(5,nij)
   complex*16 A,B,cmz,Aj,Bj,A2,B2,Aj2,Bj2,A0,B0,ephi,ci,cin,
   + atr0(ni0,nij),btr0(ni0,nij),atr(2,np,nmp),at(nmp),bt(nmp),
   + atr1(ni0,nij),btr1(ni0,nij),ek(np,nij),ref(nLp),ref0(nLp),
   + p0(nLp,nmp),q0(nLp,nmp),an(np),bn(np),aMie(nLp,np),
   + bMie(nLp,np),B2i(nLp),
   +    as(nLp,nmp),bs(nLp,nmp),as0(nLp,nmp),bs0(nLp,nmp),
   + asc(nLp,nmp),bsc(nLp,nmp),as1(nLp,nmp),bs1(nLp,nmp),
   + ast(nLp,nmp),bst(nLp,nmp),asp(nLp,nmp),bsp(nLp,nmp),
   + asv(nLp,nmp),bsv(nLp,nmp),
   + s2x(ncmax,nangmax),s4x(ncmax,nangmax),
   +    s3y(ncmax,nangmax),s1y(ncmax,nangmax),
   +    atj(nmp),btj(nmp),py0(NXMAX),py(NXMAX),dpy(NXMAX)
   CHARACTER FLNAME*20,fileout*20,fileout1*19,fileout2*21,
   + tailn*3,cnr2*2,cnr3*3,cnr1*1,flout*22
   COMMON/MIESUB/ twopi,pih
   common/rot/bcof(0:np+2),dc(-np:np,0:nmp)
   common/fnr/fnr(0:2*(np+2))
   common/pitau/pi(nmp0),tau(nmp0)
   common/tran/atr
   common/ig0/iga0(ni0)
   common/g0/ga0(ng0)
   common/cofmnv0/cof0(ni0)
   common/crot/cofsr(nmp)
   pih = dacos(0.d0)
   twopi = 4.d0*pih
   pione  = 2.d0*pih
   ci=dcmplx(0.d0,1.d0)
   cin=dcmplx(0.d0,-1.d0)
   gcs=0.d0
   gcv=0.d0
   idpq=0
   OPEN(UNIT=1,FILE='C:\Users\lenovo\Desktop\gmm01f.in',status='old')
   READ(1,'(a20)') FLNAME
程序編譯正確，但若在最后一行前加read(*,*),運(yùn)行屏幕會(huì)停止，而加在最后一行后運(yùn)行屏幕會(huì)一閃而過，帖子可能有點(diǎn)長，但都是定義部分，請(qǐng)各位大牛指點(diǎn)一下

贊一下

回復(fù)此樓

4樓2013-09-05 21:08:03

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

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

返回列表

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

[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽陽哇塞
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[教師之家] 焦慮 +9	水冰月月野兔 2026-03-13	13/650	2026-03-19 09:50 by otani
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、

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

24小時(shí)熱門版塊排行榜

xiuxialy

[求助] fortran循環(huán)求助

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

pippi6

【答案】應(yīng)助回帖

xiuxialy

xiuxialy

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

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: