SUBROUTINE LONLAT_PARA(global_lats_r,XLON,XLAT,lonsperlar) ! c*********************************************************************** ! USE MACHINE , ONLY : kind_rad,kind_phys use resol_def use layout1 use gg_def use physcons, pi => con_pi implicit none integer i,j,lat integer lonsperlar(latr) real (kind=kind_rad) tpi,hpi,bphi PARAMETER (TPI=2.E0*PI,HPI=0.5E0*PI) integer global_lats_r(latr) real (kind=kind_rad) XLON(lonr,lats_node_r) real (kind=kind_rad) XLAT(lonr,lats_node_r) ! xlon=0. xlat=0. DO j=1,lats_node_r lat = global_lats_r(ipt_lats_node_r-1+j) BPHI = TPI/lonsperlar(lat) if (lat.le.latr2) then DO i=1,lonsperlar(lat) XLON(I,J) = (i-1) * BPHI XLAT(I,J) = HPI - colrad_r(lat) ENDDO else DO i=1,lonsperlar(lat) XLON(I,J) = (i-1) * BPHI XLAT(I,J) = colrad_r(lat)-HPI ENDDO endif ENDDO RETURN END