var/da/da_tools/da_llxy_kma_global.inc

   1 subroutine da_llxy_kma_global(lat,lon,x,y)
   2
   3    !----------------------------------------------------------------------------
   4    ! Purpose:  calculates the(x,y) location(dot) in the global grids
   5    !           from latitudes and longitudes
   6    !----------------------------------------------------------------------------
   7
   8    implicit none
   9
  10    real, intent(in)  :: lat, lon
  11    real, intent(out) :: x, y
  12
  13    real              :: xlat, xlon
  14
  15    if (trace_use_frequent) call da_trace_entry("da_llxy_kma_global")
  16
  17    xlat = lat - start_lat
  18    xlon = lon - start_lon
  19
  20    if (xlat < 0.0) xlat = xlat + 180.0
  21    if (xlon < 0.0) xlon = xlon + 360.0
  22
  23    x = start_x + xlon/delt_lon
  24    y = start_y + xlat/delt_lat
  25
  26    if (trace_use_frequent) call da_trace_exit("da_llxy_kma_global")
  27
  28 end subroutine da_llxy_kma_global
  29
  30