yomprad.F90

Go to the documentation of this file.

MODULE yomprad

USE parkind1  ,ONLY : jpim     ,jprb

IMPLICIT NONE

SAVE

!----------------------------------------------------------------
!*    new radiation grid data structures
!----------------------------------------------------------------

! NRESOL_ID : transformation package resolution identifier
! NSMAX     : spectral truncation
! NASM0     : address in a spectral array of (m, n=m)
! MYMS      : actual wave numbers handled by this processor
! NDGLG     : number of latitude rows
! NDLON     : length of a row of latitude near equator
! NDGSAL    : Local version of NDGSA , always 1
! NDGENL    : Number of latitude rows for which this process has grid points
! NDGSAH    : 1-NROWIDE
! NDGENH    : NDGENL+NROWIDE
! NDGSAG    : Global version of NDGSA
! NDGENG    : Global version of NDGEN
! NDSUR1    : over dimensioning of NDLON for technical reasons (at least 2)
! NDLSUR    : NDLON+NDSUR1
! NDGSUR    : number of additional rows at each pole for semi-lagrangian
!           : interpolation
! MYFRSTACTLAT : first actual lat on this PE in grid-point space,
!              : it is nfrstlat(my_region_ns)
! MYLSTACTLAT  : last actual lat on this PE in grid-point space,
!              : it is nlstlat(my_region_ns)
! NGPTOT    : number of local grid points
! NGPTOTG   : total number of grid points
! NRGRI     : number of grid points on a latitude row
! NLOENG    : as above but extended at poles for SL interpolation
! GELAM     : radiation grid geographic longitude "lambda"
! GELAT     : radiation grid geographic latitude "theta"
! GESLO     : sine of geographic longitude "sin(lambda)"
! GECLO     : cosine of geographic longitude "cos(lambda)
! GEMU      : sine of geographic latitude "sin(theta)"
! RMU       : mu              sin(theta)
! RSQM2     : SQRT(R1MU2)     cos(theta)
! RLATIG    : arcsin(mu)      theta  GLOBAL VIEW
! RLATI     : arcsin(mu)      theta
! RIPI0     : bi-cubic interpolation coefficients
! RIPI1     : bi-cubic interpolation coefficients
! RIPI2     : bi-cubic interpolation coefficients

TYPE radiation_grid_struct
INTEGER(KIND=JPIM)                            :: nresol_id, ngptot, ngptotg, &
 & NGPTOTMX, NSPEC2, NSMAX, &
 & NPTRFLOFF, NUMP, NDLON, &
 & NDGSAL, NDGENL, NDGSAH, NDGENH, &
 & NDGLG, NDGSAG, NDGENG, NDLSUR, &
 & NFRSTLOFF, NDSUR1,NDGSUR, &
 & MYFRSTACTLAT, MYLSTACTLAT  
INTEGER(KIND=JPIM), POINTER, DIMENSION(:)     :: nrgri, nloeng, nptrfrstlat, nfrstlat, &
 & NLSTLAT, MYMS, NASM0  
INTEGER(KIND=JPIM), POINTER, DIMENSION(:,:)   :: nsta, nonl
REAL(KIND=JPRB),    POINTER, DIMENSION(:)     :: gelam, gelat, geclo, geslo, gemu, &
 & RMU, RSQM2, RLATIG, RLATI, &
 & RIPI0, RIPI1, RIPI2  
end TYPE radiation_grid_struct

TYPE(radiation_grid_struct) :: radgrid

!----------------------------------------------------------------
!*    radiation on demand comms data structures
!----------------------------------------------------------------

LOGICAL :: lradondem
LOGICAL :: lradondem_active
INTEGER(KIND=JPIM) :: nfixradfld(2)
INTEGER(KIND=JPIM),ALLOCATABLE :: mask_ri1(:)
INTEGER(KIND=JPIM),ALLOCATABLE :: mask_ri2(:)
INTEGER(KIND=JPIM),ALLOCATABLE :: mask_ro1(:)
INTEGER(KIND=JPIM),ALLOCATABLE :: mask_ro2(:)

!----------------------------------------------------------------
!*    radiation interpolation/load balancing data structures
!----------------------------------------------------------------

! NRIRINT  : INTERPOLATION INTERVAL FOR LOCAL LATITUDES
! NRFRSTOFF: OFFSET TO FIRST COURSE POINT FROM FIRST FINE POINT
! NRLASTOFF: OFFSET TO FIRST COURSE POINT FROM LAST  FINE POINT
! NRIMAX   : NUMBER OF COURSE POINTS FOR LOCAL LATITUDES
! NRIMAXT  : TOTAL NUMBER OF COURSE POINTS IN THIS PROCESSOR
! NRIMAXLT : TOTAL NUMBER OF COURSE POINTS TO BE PROCESSED BY THIS PROC
! NRLPRCS  : NUMBER OF PROCESSORS REQUIRED FOR LOAD BALANCING
! NRIMAXLA : LOAD BALANCING DATA - SETA OF PROCESSOR TO COMMUNICATE WITH
! NRIMAXLB : LOAD BALANCING DATA - SETB OF PROCESSOR TO COMMUNICATE WITH
! NRIMAXLN : LOAD BALANCING DATA - NUMBER OF COURSE POINTS TO SND/RCV
! NRCNEEDW : NUMBER OF COURSE POINTS NEEDED ON WESTERN END OF LOCAL LATS
! NRCNEEDE : NUMBER OF COURSE POINTS NEEDED ON EASTERN END OF LOCAL LATS
! NRCSNDW  : NUMBER OF WESTERN COURSE POINTS TO SEND TO LOCAL PROCS
! NRCSNDE  : NUMBER OF EASTERN COURSE POINTS TO SEND TO LOCAL PROCS
! NRCRCVW  : NUMBER OF WESTERN COURSE POINTS TO RCV FROM LOCAL PROCS
! NRCRCVE  : NUMBER OF WESTERN COURSE POINTS TO RCV FROM LOCAL PROCS
! NRCSNDT  : TOTAL NUMBER OF COURSE POINTS TO SEND TO LOCAL PROCESSORS
! NRCRCVT  : TOTAL NUMBER OF COURSE POINTS TO RCV FROM LOCAL PROCESSORS
! NRCRCVWO : OFFSET OF WESTERN BOUNDARY COURSE POINTS IN BOUNDARY ARRAY
! NRCRCVEO : OFFSET OF EASTERN BOUNDARY COURSE POINTS IN BOUNDARY ARRAY

! LODBGRADI: DEBUGGING FLAG, EACH PROCESSOR WILL DUMP ABOVE INFO TO A
!            PER PROCESSOR FILE CALLED 'debug_aaa_bbb' 
!            where aaa=my_region_ns and bbb=my_region_ew
! LODBGRADL: DEBUGGING FLAG, EACH PROCESSOR WILL DUMP LOAD BALANCING
!            ARRAYS as above

! NRLBCHUNKS: NUMBER OF BLOCKING LOOPS FOR LOAD-BALANCING COMMS
! NRLRCHUNKS: NUMBER OF BLOCKING LOOPS FOR LOAD-BALANCING COMMS (RESTORING)
! NRLBPOINTS: MAX NUMBER OF POINTS FOR LOAD-BALANCING COMMS
! NRLRPOINTS: MAX NUMBER OF POINTS FOR LOAD-BALANCING COMMS (RESTORING)
! NRLBDATA  : NUMBER OF DATA VALUES PER LOAD-BALANCING POINT 
! NRLRDATA  : NUMBER OF DATA VALUES PER LOAD-BALANCING POINT (RESTORING)

LOGICAL :: lodbgradi
LOGICAL :: lodbgradl
INTEGER(KIND=JPIM) :: nrlbchunks
INTEGER(KIND=JPIM) :: nrlrchunks
INTEGER(KIND=JPIM) :: nrlbpoints
INTEGER(KIND=JPIM) :: nrlrpoints
INTEGER(KIND=JPIM) :: nrlbdata
INTEGER(KIND=JPIM) :: nrlrdata
INTEGER(KIND=JPIM), PARAMETER :: jpmaxloadp=32
INTEGER(KIND=JPIM),ALLOCATABLE:: nrirint(:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrfrstoff(:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrlastoff(:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrimax(:,:)
INTEGER(KIND=JPIM) :: nrimaxt
INTEGER(KIND=JPIM) :: nrimaxlt
INTEGER(KIND=JPIM) :: nrlprcs
INTEGER(KIND=JPIM) :: nrimaxla(jpmaxloadp)
INTEGER(KIND=JPIM) :: nrimaxlb(jpmaxloadp)
INTEGER(KIND=JPIM) :: nrimaxln(jpmaxloadp)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcneedw(:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcneede(:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcsndw(:,:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcsnde(:,:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcrcvw(:,:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcrcve(:,:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcsndt(:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcrcvt(:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcrcvwo(:,:,:)
INTEGER(KIND=JPIM),ALLOCATABLE:: nrcrcveo(:,:,:)
!$OMP THREADPRIVATE(lodbgradi,lodbgradl,lradondem,lradondem_active,nfixradfld,nrimaxla,nrimaxlb)
!$OMP THREADPRIVATE(nrimaxln,nrimaxlt,nrimaxt,nrlbchunks,nrlbdata,nrlbpoints,nrlprcs,nrlrchunks)
!$OMP THREADPRIVATE(nrlrdata,nrlrpoints,radgrid)
!$OMP THREADPRIVATE(mask_ri1,mask_ri2,mask_ro1,mask_ro2,nrcneede,nrcneedw,nrcrcve,nrcrcveo,nrcrcvt)
!$OMP THREADPRIVATE(nrcrcvw,nrcrcvwo,nrcsnde,nrcsndt,nrcsndw,nrfrstoff,nrimax,nrirint,nrlastoff)
END MODULE yomprad