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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

xiuxialy

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 0.5
帖子: 6
在線: 6.6小時
蟲號: 1997605
注冊: 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ù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 已經(jīng)有4人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有6人回復(fù)
081700化工學(xué)碩調(diào)劑已經(jīng)有3人回復(fù)

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

FORTRAN新手求助主程序循環(huán)問題已經(jīng)有10人回復(fù)
【求助】請問如何解決insufficient virtual memory 已經(jīng)有10人回復(fù)
【求助】Python 中求解微分方程 odeint 的用法已經(jīng)有11人回復(fù)
【求助】ANSYS 中命令流編程（fortran）已經(jīng)有11人回復(fù)
【求助】fortran程序未出錯，但運(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ù)值計算咨詢

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

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

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

因為 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小時
蟲號: 1997605
注冊: 2012-09-13
專業(yè): 光學(xué)

引用回帖:

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

謝謝，昨天發(fā)帖后我編了個小程序發(fā)現(xiàn)這個問題了，確實(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小時
蟲號: 1997605
注冊: 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)行屏幕會停止，而加在最后一行后運(yùn)行屏幕會一閃而過，帖子可能有點(diǎn)長，但都是定義部分，請各位大牛指點(diǎn)一下

贊一下

回復(fù)此樓

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

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

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

返回列表

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

[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 生物學(xué)調(diào)劑招人！�。� +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +8	zlingli 2026-03-13	8/400	2026-03-19 16:32 by 輕松不少隨
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 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
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7	手機(jī)用戶 2026-03-17	7/350	2026-03-18 14:34 by vgtyfty
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky

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

24小時熱門版塊排行榜

xiuxialy

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

» 猜你喜歡

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

pippi6

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

xiuxialy

xiuxialy

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

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