ソースコードの動作がわからない

Question

Gfortranを使っています。
コードがどのように動作しているかがわからないです。

      do 1000 M	 =1,max
!
      resp = 0.0d0
      if  ( i <= 95 .and. i >= 33 .and. j <= 95 .and. j >=33 ) then
      do 190 i = 33,95
      do 190 j = 33,95
      pres     = ( h*s(i,j) + dx*p(i,j) ) / ( h + dx )
      P(i,j)   = P(i,j) + pres
      resp     = resp+pres*pres
  190 continue 
      else
      do 200 i = 1, 200
      do 200 j = 1, 200
      dpdxx    = (p(i+1,j) - 2.0d0*p(i,j) + p(i-1,j))/ ddx
      dpdyy    = (p(i,j+1) - 2.0d0*p(i,j) + p(i,j-1))/ ddy
      pres     = ram*( dpdxx+dpdyy - s(i,j) )
      resp     = resp+pres*pres
      P(i,j)   = P(i,j) + pres
  200 continue
      end if

33＜i,j＜95の範囲とそれ以外で分けたいのですが
上記のソースコードで正しいのでしょうか？

Answer 1

do の後の番号とcontinueの前の番号は呼応しています。CやJavaの{ }のくくりのようなものです。
したがって、do 1000 に呼応した1000 continue が末尾の方に必要です。スペースをあえて"_"を用いて段分けを強調したコードを示せば、

      do 1000 M =1,max

        resp = 0.0d0
        if  ( i <= 95 .and. i >= 33 .and. j <= 95 .and. j >=33 ) then
          do 190 i = 33,95
              do 191 j = 33,95
                pres   = ( h*s(i,j) + dx*p(i,j) ) / ( h + dx )
                P(i,j) = P(i,j) + pres
                resp   = resp+pres*pres
  191         continue 
  190     continue 

        else
          do 200 i = 1, 200
            do 201 j = 1, 200
              dpdxx   = (p(i+1,j) - 2.0d0*p(i,j) + p(i-1,j))/ ddx
              dpdyy   = (p(i,j+1) - 2.0d0*p(i,j) + p(i,j-1))/ ddy
              pres    = ram*( dpdxx+dpdyy - s(i,j) )
              resp    = resp+pres*pres
              P(i,j)  = P(i,j) + pres
  201       continue
  200     continue
        end if

 1000 continue

のように表現でき、多少は分かりやすいのではないでしょうか。

Comments

[him]

回答ありがとうございます。continueの部分を書き足すのをすっかり忘れていました。
自分としてはelse以降でi,jをともに1~200の全範囲で取っているのですが
その中にifで設定した33~95の範囲も含まれてしまっているので
上書きされてしまっていないかが心配です