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