プログラム中で同じ処理を何度も行いたい場合がよくある．そのたびごとに同 じプログラムを何度も書くのは馬鹿げている．同じ処理はサブルーチンとして 外部手続きとしておくと都合がよい．
サブルーチンを呼び出す側のプログラムは「主プログラム(main program)」と呼ばれる．サブルーチンは主プログラムのend文以降に，主プログラムと独立に次のように書く. subroutine文以下endまでがサブルーチンである.
c main program
implicit none
real*8 d(10),f,g
......
call S(f,g,d)
......
end
subroutine S(a,b,c)
implicit none
real*8 a,b,c(10)
......
......
return
end
| サブルーチン文 | S(a,b,c) | サブルーチンの定義をおこなう |
| S | サブルーチン名：英文字で始まる6文字以内の英数字 |
|---|---|
| a,b,c | 仮引数：使う変数は型宣言すること |
| call 文 | call S(f,g,d) | サブルーチンを呼び出す |
| f,g,d | 実引数：仮引数と順序, 型を一致させる必要あり |
|---|
サブルーチンはend文で終わらなければならない．サブルーチンから主プログラムに戻る場所では
return
を書く．return文は何個あってもよい(endは一個：end文のみでreturn文がない場合はendで呼ばれた場所に戻る)．それでは次のプログラム
を入力して実行してみよう．
プログラム１
c main program
implicit none
real*8 a,b,c
a=1.D0
b=2.D0
call wa(a,b,c)
write(6,*) c
end
c subroutine
subroutine wa(a,b,c)
implicit none
real*8 a,b,c
c=a+b
return
end
呼ぶ方(call wa(a,b))と呼ばれる方(subroutine wa(a,b))の引数は順序，個数，型を一致させなければならない．逆に言うと，引数の個数，型が一致していれば名前は違っていてもよい．サブルーチンを使う場合は, サブルーチンとの変数, 配列のやりとりに注意すること．どれが入力でどれが出力かを頭にいれてプログラムを作成すること.
次のプログラムは何をするのか考えてみよう．
c main program
implicit none
real*8 a,b,c,d,e
a=1.D0
b=2.D0
d=3.D0
e=4.D0
call wa(a,b,c)
write(6,*) c
call wa(d,e,c)
write(6,*) c
end
subroutine wa(a,b,c)
implicit none
real*8 a,b,c
c=a+b
return
end
プログラム1
c 1-1
implicit none
real*8 a,b
a=1.D0
b=2.D0
call wa(a,b)
write(6,*) b
end
subroutine wa(a,b,c)
implicit none
real*8 a,b
b=a+b
return
end
プログラム2
c 1-2
implicit none
real*8 a,b
a=1.D0
b=2.D0
call wa(a,b)
write(6,*) b
end
subroutine wa(a,b)
implicit none
integer a,b
b=a+b
return
end
プログラム3
c 1-3
implicit none
real*8 a(2),b
a(1)=1.D0
a(2)=2.D0
b=0.D0
call wa(a,b)
write(6,*) b
end
subroutine wa(a,b)
implicit none
real*8 a,b
b=a(1)+a(2)
return
end
プログラム1-2を修正し, 二つの数の積を計算するサブルーチンを書き加えよ. これを用いて変数a,bの積（答えは2.D0）を算出せよ.
変数a,bに代入された値を入れ替えるサブルーチンswap(a,b)を作成せよ．サブルーチンとは変数a,bがやりとりされる．
実数が10個格納された配列a(10)を，a(1)から大きい順に並べ替えるサブルーチンsort(n,a)を作成せよ．サブルーチンとは整数nおよび配列aがやりとりされる．