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