fortranでサブルーチンを使って[rad]を[度]に変換するプログラムがわからない

問題
xとyを[rad]で与える。任意のxとyを与えた時、x,y,x+yを[度,分,秒]で出力するプログラムをサブルーチンを使って作成せよ。
ただし
180[°]はπ[rad]
1[°]は60[分]
1[分]は60[秒]
とする。

大学の課題で出ていて、とりあえずxだけ入力して出力するプログラムを書いているのですが、ここから先どうやって書けばいいのかいまいちわかりません。

fortran
1c xとyを[rad]で入力してx,y,x+yを[度,分,秒]で出力するプログラム
2c まずはxだけを入力して変換するプログラムを作る
3
4print *,'** Key in x[rad] **'
5read(*,*) x
6
7
8call radtodeg(s,idx,ihx,bx)
9
10stop
11end
12
13
14subroutine radtodeg(r,id,ih,b)
15data pai/3.141592/ !πを与える
16d=180*x/pai !rad->°
17id=d ![°]を整数化
18s1=d-id !小数点以下の[°]を取得
19h=s1*60 ![°]->[分]
20ih=h ![分]を整数化
21s2=h-ih !小数点以下の[分]を取得
22b=s2*60 ![分]->[秒]
23return
24end
25
26
27