<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_LLXY_MERC_NEW'><A href='../../html_code/tools/da_llxy_merc_new.inc.html#DA_LLXY_MERC_NEW' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
subroutine da_llxy_merc_new(proj, info) 1,2
!-----------------------------------------------------------------------
! Purpose: Compute x/y coordinate from lat lon for mercator projection
!-----------------------------------------------------------------------
implicit none
type(proj_info), intent(in) :: proj
type(infa_type), intent(inout) :: info
real :: deltalon
integer :: n
if (trace_use) call da_trace_entry
("da_llxy_merc_new")
! FAST
! where (lon(:,:) - proj%lon1 < -180.0)
! x(:,:) = 1.0 + (lon(:,:) - proj%lon1 + 360.0)/(proj%dlon*deg_per_rad))
! elsewhere (lon(:,:) - proj%lon1 > 180.0)
! x(:,:) = 1.0 + (lon(:,:) - proj%lon1 - 360.0)/(proj%dlon*deg_per_rad))
! else
! x(:,:) = 1.0 + (lon(:,:) - proj%lon1)/(proj%dlon*deg_per_rad))
! end where
! y(:,:) = 1.0 + (ALOG(TAN(0.5*((lat(:,: + 90.0) * rad_per_deg)))) / proj%dlon - proj%rsw
! SLOW
do n=lbound(info%lat,2),ubound(info%lat,2)
deltalon = info%lon(1,n) - proj%lon1
if (deltalon < -180.0) deltalon = deltalon + 360.0
if (deltalon > 180.0) deltalon = deltalon - 360.0
info%x(:,n) = 1.0 + (deltalon/(proj%dlon*deg_per_rad))
info%y(:,n) = 1.0 + (ALOG(TAN(0.5*((info%lat(1,n) + 90.0) * rad_per_deg)))) / proj%dlon - proj%rsw
end do
if (trace_use) call da_trace_exit("da_llxy_merc_new")
end subroutine da_llxy_merc_new