var/da/da_tools/da_llxy_latlon.inc

   1 subroutine da_llxy_latlon(lat, lon, proj, x, y)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: Compute the x/y location of a lat/lon on a LATLON
   5    !          (cylindrical equidistant) grid.
   6    !-----------------------------------------------------------------------
   7
   8    implicit none
   9
  10    real, intent(in)             :: lat
  11    real, intent(in)             :: lon
  12    type(proj_info), intent(in)  :: proj
  13    real, intent(out)            :: x
  14    real, intent(out)            :: y
  15
  16    real                         :: deltalat
  17    real                         :: deltalon
  18    real                         :: lon360
  19    real                         :: latinc
  20    real                         :: loninc
  21
  22    if (trace_use_frequent) call da_trace_entry("da_llxy_latlon")
  23
  24    ! To account for issues around the dateline, convert the incoming
  25    ! longitudes to be 0->360.0
  26    if (lon < 0) then
  27       lon360 = lon + 360.0
  28    else
  29       lon360 = lon
  30    end if
  31
  32    deltalat = lat - proj%lat1
  33    deltalon = lon360 - proj%lon1
  34
  35    !For cylindrical equidistant, dx == dy
  36    loninc = proj%dx*360.0/(2.0*EARTH_RADIUS_M*PI)
  37    latinc = proj%dx*360.0/(2.0*EARTH_RADIUS_M*PI)
  38
  39    ! Compute x/y
  40    x = deltalon/loninc
  41    y = deltalat/latinc
  42
  43    x = x + proj%knowni
  44    y = y + proj%knownj
  45
  46
  47    if (trace_use_frequent) call da_trace_exit("da_llxy_latlon")
  48
  49 end subroutine da_llxy_latlon
  50
  51