subroutine da_llxy_merc(lat, lon, proj, x, y) !----------------------------------------------------------------------- ! 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