File: C:\NOAA\NEMS_11731\src\atmos\gfs\phys\lon_lat_para.f

1           SUBROUTINE LONLAT_PARA(global_lats_r,XLON,XLAT,lonsperlar)
2     !
3     c***********************************************************************
4     !
5           USE module_gfs_machine, ONLY : kind_grid
6     
7           use resol_def
8           use layout1
9           use gg_def
10           use physcons, pi => con_pi
11           implicit none
12           integer i,j,lat
13           integer                 lonsperlar(latr)
14           real (kind=kind_grid) tpi,hpi,bphi
15           PARAMETER (TPI=2.E0*PI,HPI=0.5E0*PI)
16           integer              global_lats_r(latr)
17           real (kind=kind_grid) XLON(lonr,lats_node_r)
18           real (kind=kind_grid) XLAT(lonr,lats_node_r)
19     !
20           xlon=0.
21           xlat=0.
22      
23           DO j=1,lats_node_r
24             lat = global_lats_r(ipt_lats_node_r-1+j)
25             BPHI = TPI/lonsperlar(lat)
26             if (lat.le.latr2) then
27               DO i=1,lonsperlar(lat)
28                 XLON(I,J) = (i-1) * BPHI
29                 XLAT(I,J) = HPI - colrad_r(lat)
30               ENDDO
31             else
32               DO i=1,lonsperlar(lat)
33                 XLON(I,J) =  (i-1) * BPHI
34                 XLAT(I,J) = colrad_r(lat)-HPI
35               ENDDO
36             endif
37           ENDDO
38      
39           RETURN
40           END
41      
42