<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_XYLL_MERC'><A href='../../html_code/tools/da_xyll_merc.inc.html#DA_XYLL_MERC' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
subroutine da_xyll_merc(x, y, proj, lat, lon) 2,2
!-----------------------------------------------------------------------
! Compute the lat/lon from i/j for mercator projection
!-----------------------------------------------------------------------
implicit none
real,intent(in) :: x
real,intent(in) :: y
type(proj_info),intent(in) :: proj
real, intent(out) :: lat
real, intent(out) :: lon
if (trace_use_frequent) call da_trace_entry
("da_xyll_merc")
lat = 2.0*ATAN(EXP(proj%dlon*(proj%rsw + y-1.0)))*deg_per_rad - 90.0
lon = (x-1.0)*proj%dlon*deg_per_rad + proj%lon1
if (lon > 180.0) lon = lon - 360.0
if (lon < -180.0) lon = lon + 360.0
if (trace_use_frequent) call da_trace_exit("da_xyll_merc")
end subroutine da_xyll_merc