[問題] 跑迴圈

[nestea911399 (卡進無底洞)]

PARAMETER (RC=3.0D-4, RN=2.0D-4)

N=0  ! Shower counter.
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C  ------------------------  Shower simulation starts here.
C
  101 CONTINUE
      N=N+1.0D0
C
C  **********  Set the initial state of the primary particle.
C
      KPAR=KPARP
        RY=RC*RAND(1.0D0)**(1.0/3.0)
        IF(RY.LT.2.0D-4) GO TO 101
        PTHETA=ACOS(2*RAND(1.0D0)-1)
        PPHI=2*PI*RAND(1.0D0)
      X=RY*SIN(PTHETA)*COS(PPHI)
      Y=RY*SIN(PTHETA)*SIN(PPHI)
      Z=RY*COS(PTHETA)

請問我現在假設有一個球的模型,而每一個亂數得到的座標會是空間上的一點

我現在是想點數均勻分布在大球跟小球之間

所以我假設IF(RY.LT.2.0D-4) GO TO 101,若是RY小於小球的半徑就重新跑

但是有個問題是跑完了就會由下一個接下去

例如

1符合但2不符合就從3開始

但是2就沒辦法記錄到就變成只有兩個結果,

有辦法解決這個問題嗎

謝謝大家


--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.114.106.110
→ latinboy: 把N=N+1移到IF下一行...  問問題之前你有沒有想過@@
→ nestea911399:N=N+1不能移到下面,因為下面還有整個迴圈像上ㄧ篇那
→ nestea911399:要是移到下面KPAR就沒辦法紀錄粒子數
推 meteorology:所以你是想要RY小於小球的半徑也要輸出嗎？
→ nestea911399:可以想像只要把點分布在甜甜圈就中間空心的不要
→ meteorology:可是就你的例如感覺在空心部份你是要粒子數資料
→ nestea911399:就是當RY小於R時我不要執行,RY大於RN才執行
→ nestea911399:像我的寫法好像可以達到我的目的但是有些粒就紀錄不
→ nestea911399:到,因為若RY<R時,N就會再+1就變成第三顆繼續執行
推 meteorology:恩那你就在if下先把你要的粒子write到檔案 再GOTO
→ meteorology:到101
→ nestea911399:可是那是亂數產生我沒辦法知道哪一次要紀錄,我只能把
→ nestea911399:不要的丟掉
推 meteorology:那必須看你要的變數是什麼了～
推 meteorology:因為每個人專業不同  所以這些輸出可能跟理論有關
→ meteorology:建議你先去翻翻理論～或許就可以知道要的是什麼～
推 latinboy: 你不PO出完整的程式碼 或是完整的do block 很難回答
→ latinboy: 感覺就是if處裡一下就能解決的問題
→ latinboy: 你還是問實驗室學長姐比較快...
→ nestea911399:完整CODE在上ㄧ篇已經有了,不好意思,我只是想問說,如
→ nestea911399:何讓每一顆粒子都是在小球以外出現,謝謝
→ meteorology:你不是因為沒有辦法看到半徑小於小球的數據？
→ meteorology:如果是這樣你原本的程式不是就可以了？
→ meteorology:因為你已經做了大於小球的資料輸出
→ nestea911399:可是因為是亂數產生有可能會有小於RN的情況發生,所以
→ nestea911399:若是用我的方法當小於RN就重新跑那會有沒記錄到的情
→ nestea911399:況,謝謝
推 meteorology:你亂數產生是產生RY嗎？
→ nestea911399:對
推 meteorology:你的問題有點亂～因為這樣看下來的確小於RN的半徑就
→ meteorology:會不見 如果你要小於RN的也要列出來  有兩個辦法
→ meteorology:一個是把IF去掉  讓迴圈跑完  並且輸出RY的值 從中挑
→ meteorology:選大於RN的RY的結果
→ meteorology:另外就是在程式中小於RN的值做另外的檔案輸出
→ irlmpo:code若不是自己寫的,建議你不要看程式碼,再重寫一次
→ irlmpo:ps.僅是個人建議,無惡意
推 meteorology:幫原po說一下～～這支程式是上萬行～XDD
→ nestea911399:謝謝M大幫我說明謝謝大家
推 noonee:那我的推荐是  把每個獨立的地方都儘量拆成subroutine
→ noonee:或許很蠢  不過釐清問題很方便
推 meteorology:樓上  這不蠢阿  這是好方法！y