pd2 = abs(rtod(asin(sinh)))
do case
case msphe > 0 .and. msphe < 90 mspc = msph -pd2 case msphe > 90 .and. msphe < 180 mspc = msph +pd2 case msphe > 180 .and. msphe < 270 mspc = msph -pd2 case msphe > 270 .and. msphe < 360
mspc = msph + pd2
endcase
close databases

