Update fn_green3d.f90