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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

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ù)此樓

» 猜你喜歡

一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑已經(jīng)有6人回復(fù)
一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑已經(jīng)有12人回復(fù)
289求調(diào)劑已經(jīng)有11人回復(fù)
食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有7人回復(fù)
生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 已經(jīng)有6人回復(fù)
321求調(diào)劑已經(jīng)有7人回復(fù)
343求調(diào)劑已經(jīng)有4人回復(fù)
305求調(diào)劑已經(jīng)有3人回復(fù)
本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑已經(jīng)有6人回復(fù)
332求調(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ù)
【求助】求指點，在Fortran里面，怎樣實現(xiàn)“數(shù)組維度可調(diào)”的數(shù)組? 已經(jīng)有5人回復(fù)
【求助】三維矩陣運(yùn)算已經(jīng)有8人回復(fù)
【求助】Fortran程序，如何讓exe文件可以單獨運(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的回帖

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)行屏幕會一閃而過，帖子可能有點長，但都是定義部分，請各位大牛指點一下

贊一下

回復(fù)此樓

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

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

查看全部 4 個回答

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)這個問題了，確實是這樣，不過仍然很感謝！

贊一下

回復(fù)此樓

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

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

查看全部 4 個回答

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

[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +5	YYYYYNNNNN 2026-03-26	6/300	2026-03-28 17:50 by 狂炫麥當(dāng)當(dāng)
[考研] 材料與化工272求調(diào)劑 +9	阿斯蒂芬2004 2026-03-28	9/450	2026-03-28 15:21 by 1018329917
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 化學(xué)308分求調(diào)劑 +8	你好明天你好 2026-03-23	9/450	2026-03-27 14:01 by 楊光于青云
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6	kotoko_ik 2026-03-23	7/350	2026-03-27 12:29 by 惠州彭于晏
[考研] 359求調(diào)劑 +4	王了個楠 2026-03-25	4/200	2026-03-27 08:43 by 不吃魚的貓
[考研] 341求調(diào)劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid
[考研] 081700 調(diào)劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 各位老師您好：本人初試372分 +5	jj涌77 2026-03-25	6/300	2026-03-25 14:15 by mapenggao
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5	超贊的 2026-03-24	5/250	2026-03-24 21:07 by 星空星月
[考研] 化工專碩求調(diào)劑 +3	question挽風(fēng) 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky
[考研] 305分求調(diào)劑（食品工程） +5	Sxy112 2026-03-21	7/350	2026-03-24 12:27 by 544594351
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106
[考研] 070300，一志愿北航320求調(diào)劑 +3	Jerry0216 2026-03-22	5/250	2026-03-23 09:16 by 。。堂堂

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

24小時熱門版塊排行榜

xiuxialy

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

» 猜你喜歡

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

xiuxialy

pippi6

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

xiuxialy

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

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