[gradsusr] help

[furugh momenpoor]

 Hi
i read and changed GFS data with fortran program
i want to read and change geosterophic component with fortran program
 I can not write geostrophic wind and temp in fortran
 ug=-1/f   (∂h)/∂y
 vg=1/f   (∂h)/∂x
 h(z2 )-h(z1 )=R∫T(d lnp )


 ( f = coriolisparameter  , h =geopotenciel height , )



 would you please help me in this regard. what should I do?

 Thanks in advance and looking forward to your reply.

 Best Regards