sumplatf_mod.F90

Go to the documentation of this file.

MODULE sumplatf_mod
CONTAINS
SUBROUTINE sumplatf(KDGL,KPROCA,KMYSETA,&
                   &kultpp,kprocl,kptrls)

!**** *SUMPLATF * - Initialize fourier space distibution in N-S direction

!     Purpose.
!     --------


!**   Interface.
!     ----------
!        *CALL* *SUMPLATF *

!     Explicit arguments - input :
!     -------------------- 
!                          KDGL       -last  latitude
!                          KPROCA     -number of processors in A direction
!                          KMYSETA    -process number in A direction

!     Explicit arguments - output:
!     -------------------- 

!                          KULTPP     -number of latitudes in process
!                                      (in Fourier space)
!                          KPROCL     -process responsible for latitude
!                                      (in Fourier space)
!                          KPTRLS     -pointer to first global latitude
!                                      of process (in Fourier space)

!        Implicit arguments :
!        --------------------


!     Method.
!     -------
!        See documentation

!     Externals.   SUMPLATB and SUEMPLATB.
!     ----------

!     Reference.
!     ----------
!        ECMWF Research Department documentation of the IFS

!     Author.
!     -------
!        MPP Group *ECMWF*

!     Modifications.
!     --------------
!        Original : 95-10-01
!        David Dent:97-06-02 parameters KFRSTLAT etc added
!        JF. Estrade:97-11-13 Adaptation to ALADIN case
!        J.Boutahar: 98-07-06  phasing with CY19
!        Modified 98-08-10 by K. YESSAD: removal of LRPOLE option + cleanings
!         (correct computation of extrapolar latitudes for KPROCL).
!        Modified 98-12-07 by K. YESSAD and C. FISCHER: cleaning.
!         - merge old sumplat.F and suemplat.F
!         - gather 'lelam' code and 'not lelam' code.
!         - clean (useless duplication of variables, non doctor features).
!         - remodularise according to lelam/not lelam
!           -> lelam features in new routine suemplatb.F,
!              not lelam features in new routine sumplatb.F
!     ------------------------------------------------------------------

USE parkind1  ,ONLY : jpim     ,jprb

use tpm_geometry

use sumplatb_mod


IMPLICIT NONE


!     * DUMMY:
INTEGER(KIND=JPIM),INTENT(IN)  :: KDGL
INTEGER(KIND=JPIM),INTENT(IN)  :: KPROCA
INTEGER(KIND=JPIM),INTENT(IN)  :: KMYSETA
INTEGER(KIND=JPIM),INTENT(OUT) :: KULTPP(:)
INTEGER(KIND=JPIM),INTENT(OUT) :: KPROCL(:)
INTEGER(KIND=JPIM),INTENT(OUT) :: KPTRLS(:)

!     * LOCAL:
INTEGER(KIND=JPIM) :: INDIC(kproca),ILAST(kproca)

!     LOCAL INTEGER SCALARS
INTEGER(KIND=JPIM) :: IA, ILAT, ISTART, IMEDIAP,IRESTM, JA,  JLTLOC
LOGICAL :: LLSPLIT

!      -----------------------------------------------------------------

!*       1.    CODE DEPENDING ON 'LELAM': COMPUTATION OF
!              KMEDIAP, KRESTM, INDIC, ILAST.
!              -----------------------------------------

llsplit = .false.

CALL sumplatb(1,kdgl,kproca,g%NLOEN,llsplit,&
 &imediap,irestm,indic,ilast)

!      -----------------------------------------------------------------

!*       2.    CODE NOT DEPENDING ON 'LELAM': 
!              ------------------------------



!     * Definitions related to distribution of latitudes along sets
!       ------------ in fourier-space -----------------------------
istart = 0
kultpp(1) = ilast(1)
DO ja=1,kproca
  IF(ja > 1) THEN
    IF(ilast(ja) /= 0) THEN
      kultpp(ja) = ilast(ja)-ilast(ja-1)
    ELSE
      kultpp(ja) = 0
    ENDIF
  ENDIF
  DO jltloc=1,kultpp(ja)
    ilat = istart + jltloc
    kprocl(ilat) = ja
  ENDDO
  istart = istart + kultpp(ja)
ENDDO

!     * Computes KPTRLS.

ia = kprocl(1)
kptrls(ia) = 1
DO ja=ia+1,kproca
  kptrls(ja) = kptrls(ja-1) + kultpp(ja-1)
ENDDO

END SUBROUTINE sumplatf
END MODULE sumplatf_mod