ParallelRoutineDispatchTransformation

loki/transformations/projParallelRoutineDispatch/cpphinp.intfb.h

@@ -0,0 +1,35 @@
+INTERFACE
+SUBROUTINE CPPHINP(YDGEOMETRY,YDMODEL,KIDIA,KFDIA,&
+ & PGEMU,PGELAM,&
+ & PUT0,PVT0,PQT0,PQT0L,PQT0M,PQSLT0L,PQSLT0M,&
+ & PRDELP0,PEVEL0,PCVGQSL,&
+ & PMU0,PMU0LU,PMU0M,PMU0N,PCVGQ,PCVGT)
+USE PARKIND1 , ONLY : JPIM, JPRB
+USE YOMHOOK  , ONLY : LHOOK, DR_HOOK, JPHOOK
+USE TYPE_MODEL   , ONLY : MODEL
+USE GEOMETRY_MOD , ONLY : GEOMETRY
+IMPLICIT NONE
+TYPE(GEOMETRY)    ,INTENT(IN)    :: YDGEOMETRY
+TYPE(MODEL)       ,INTENT(IN)    :: YDMODEL
+INTEGER(KIND=JPIM),INTENT(IN)    :: KIDIA 
+INTEGER(KIND=JPIM),INTENT(IN)    :: KFDIA 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PGEMU(YDGEOMETRY%YRDIM%NPROMA) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PGELAM(YDGEOMETRY%YRDIM%NPROMA) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PUT0(YDGEOMETRY%YRDIM%NPROMA,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PVT0(YDGEOMETRY%YRDIM%NPROMA,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PQT0(YDGEOMETRY%YRDIM%NPROMA,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PQT0L(YDGEOMETRY%YRDIM%NPROMA,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PQT0M(YDGEOMETRY%YRDIM%NPROMA,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PQSLT0L(YDGEOMETRY%YRDIM%NPROMA,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PQSLT0M(YDGEOMETRY%YRDIM%NPROMA,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PRDELP0(YDGEOMETRY%YRDIM%NPROMA,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PEVEL0(YDGEOMETRY%YRDIM%NPROMA,0:YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(IN)    :: PCVGQSL(YDGEOMETRY%YRDIM%NPROMA,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PMU0(YDGEOMETRY%YRDIM%NPROMA) 
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PMU0LU(YDGEOMETRY%YRDIM%NPROMA) 
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PMU0M(YDGEOMETRY%YRDIM%NPROMA) 
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PMU0N(YDGEOMETRY%YRDIM%NPROMA)
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PCVGQ(YDGEOMETRY%YRDIM%NPROMM,YDGEOMETRY%YRDIMV%NFLEVG) 
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PCVGT(YDGEOMETRY%YRDIM%NPROMM, YDGEOMETRY%YRDIMV%NFLEVG)
+END SUBROUTINE CPPHINP
+END INTERFACE

loki/transformations/projParallelRoutineDispatch/cpphinp_openacc.intfb.h

@@ -0,0 +1,40 @@
+INTERFACE
+SUBROUTINE CPPHINP (YDGEOMETRY, YDMODEL, KIDIA, KFDIA, PGEMU, PGELAM,&
+& PUT0, PVT0, PTT0L, PTT0M, PQT0, PQT0L, PQT0M, PQSLT0L, PQSLT0M, PRDELP0, PEVEL0&
+&, PCVGQSL, PMU0, PSOLO, PMU0LU, PMU0M, PMU0N, PCVGQ, PCVGT, YDSTACK)
+!$acc routine (CPPHINP_OPENACC) seq
+USE PARKIND1,ONLY:JPIM, JPRB
+USE TYPE_MODEL,ONLY:MODEL
+USE GEOMETRY_MOD,ONLY:GEOMETRY
+USE ABOR1_ACC_MOD
+USE STACK_MOD
+IMPLICIT NONE
+TYPE (GEOMETRY), INTENT (IN)::YDGEOMETRY
+TYPE (MODEL), INTENT (IN)::YDMODEL
+INTEGER (KIND=JPIM), INTENT (IN)::KIDIA
+INTEGER (KIND=JPIM), INTENT (IN)::KFDIA
+REAL (KIND=JPRB), INTENT (IN)::PGEMU (YDGEOMETRY%YRDIM%NPROMA)
+REAL (KIND=JPRB), INTENT (IN)::PGELAM (YDGEOMETRY%YRDIM%NPROMA)
+REAL (KIND=JPRB), INTENT (IN)::PUT0 (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PVT0 (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PQT0 (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PQT0L (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PQT0M (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PTT0L (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PTT0M (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PQSLT0L (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PQSLT0M (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PRDELP0 (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PEVEL0 (YDGEOMETRY%YRDIM%NPROMA, 0:YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (IN)::PCVGQSL (YDGEOMETRY%YRDIM%NPROMA, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (OUT)::PMU0 (YDGEOMETRY%YRDIM%NPROMA)
+REAL (KIND=JPRB), INTENT (OUT)::PSOLO (YDGEOMETRY%YRDIM%NPROMA)
+REAL (KIND=JPRB), INTENT (OUT)::PMU0LU (YDGEOMETRY%YRDIM%NPROMA)
+REAL (KIND=JPRB), INTENT (OUT)::PMU0M (YDGEOMETRY%YRDIM%NPROMA)
+REAL (KIND=JPRB), INTENT (OUT)::PMU0N (YDGEOMETRY%YRDIM%NPROMA)
+REAL (KIND=JPRB), INTENT (OUT)::PCVGQ (YDGEOMETRY%YRDIM%NPROMM, YDGEOMETRY%YRDIMV%NFLEVG)
+REAL (KIND=JPRB), INTENT (OUT)::PCVGT (YDGEOMETRY%YRDIM%NPROMM, YDGEOMETRY%YRDIMV%NFLEVG)
+TYPE(STACK) :: YDSTACK
+ENDSUBROUTINE CPPHINP
+
+END INTERFACE

loki/transformations/projParallelRoutineDispatch/dispatch_routine.F90

@@ -0,0 +1,71 @@
+SUBROUTINE DISPATCH_ROUTINE(YDGEOMETRY, YDCPG_BNDS, YDCPG_OPTS, &
+    & YDMF_PHYS, YDCPG_PHY0, YDCPG_DYN0, YDMF_PHYS_SURF, YDVARS, &
+    & YDMODEL)
+
+USE GEOMETRY_MOD       , ONLY : GEOMETRY
+USE MF_PHYS_TYPE_MOD   , ONLY : MF_PHYS_TYPE
+USE CPG_TYPE_MOD       , ONLY : CPG_DYN_TYPE, &
+                                & CPG_PHY_TYPE
+USE CPG_OPTS_TYPE_MOD  , ONLY : CPG_BNDS_TYPE, CPG_OPTS_TYPE
+USE MF_PHYS_SURFACE_TYPE_MOD  &
+                        & , ONLY : MF_PHYS_SURF_TYPE
+USE FIELD_VARIABLES_MOD, ONLY : FIELD_VARIABLES
+USE TYPE_MODEL         , ONLY : MODEL
+
+USE PARKIND1           , ONLY : JPIM     ,JPRB
+USE YOMHOOK            , ONLY : LHOOK    ,DR_HOOK, JPHOOK
+
+IMPLICIT NONE
+
+TYPE(GEOMETRY),                 INTENT(IN)    :: YDGEOMETRY
+TYPE(CPG_BNDS_TYPE),            INTENT(IN)    :: YDCPG_BNDS
+TYPE(CPG_OPTS_TYPE),            INTENT(IN)    :: YDCPG_OPTS
+TYPE(MF_PHYS_TYPE),             INTENT(INOUT) :: YDMF_PHYS
+TYPE(CPG_PHY_TYPE),             INTENT(IN)    :: YDCPG_PHY0
+TYPE(CPG_DYN_TYPE),             INTENT(IN)    :: YDCPG_DYN0
+TYPE(MF_PHYS_SURF_TYPE),        INTENT(INOUT) :: YDMF_PHYS_SURF
+TYPE(FIELD_VARIABLES),          INTENT(INOUT) :: YDVARS
+TYPE(MODEL),                    INTENT(IN)    :: YDMODEL
+
+#include "cpphinp.intfb.h"
+#include "mf_phys_fpl_part1.intfb.h"
+
+REAL(KIND=JPRB)     :: ZRDG_CVGQ (YDCPG_OPTS%KLON, 1:YDCPG_OPTS%KFLEVG)
+REAL(KIND=JPRB)     :: ZRDG_MU0LU (YDCPG_OPTS%KLON)
+REAL(KIND=JPRB)     :: ZRDG_MU0M (YDCPG_OPTS%KLON)
+REAL(KIND=JPRB)     :: ZRDG_MU0N (YDCPG_OPTS%KLON)
+REAL(KIND=JPRB)     :: ZRDG_MU0 (YDCPG_OPTS%KLON)
+REAL(KIND=JPRB)     :: ZPFL_FPLCH (YDCPG_OPTS%KLON, 0:YDCPG_OPTS%KFLEVG)
+REAL(KIND=JPRB)     :: ZPFL_FPLSH (YDCPG_OPTS%KLON, 0:YDCPG_OPTS%KFLEVG)
+REAL(KIND=JPRB)     :: TOTO (YDCPG_OPTS%KLON)
+
+
+INTEGER(KIND=JPIM) :: INSTEP_DEB, INSTEP_FIN
+INTEGER(KIND=JPIM) :: JLEV, JLON, JSPP
+
+REAL(KIND=JPHOOK) :: ZHOOK_HANDLE
+
+IF (LHOOK) CALL DR_HOOK('DISPATCH_ROUTINE',0,ZHOOK_HANDLE)
+
+INSTEP_DEB=1
+INSTEP_FIN=1
+!!!$ACDC PARALLEL,TARGET=OpenMP/OpenMPSingleColumn/OpenACCSingleColumn,NAME=CPPHINP {
+!$loki parallel PARALLEL,TARGET=OpenMP/OpenMPSingleColumn/OpenACCSingleColumn,NAME=CPPHINP
+
+CALL CPPHINP(YDGEOMETRY, YDMODEL, YDCPG_BNDS%KIDIA, YDCPG_BNDS%KFDIA, YDVARS%GEOMETRY%GEMU%T0,                     &
+& YDVARS%GEOMETRY%GELAM%T0, YDVARS%U%T0, YDVARS%V%T0, YDVARS%Q%T0, YDVARS%Q%DL, YDVARS%Q%DM, YDVARS%CVGQ%DL,       &
+& YDVARS%CVGQ%DM, YDCPG_PHY0%XYB%RDELP, YDCPG_DYN0%CTY%EVEL, YDVARS%CVGQ%T0, ZRDG_MU0, ZRDG_MU0LU, ZRDG_MU0M,      &
+& ZRDG_MU0N, ZRDG_CVGQ, YDMF_PHYS_SURF%GSD_VF%PZ0F)
+
+DO JLON = YDCPG_BNDS%KIDIA,YDCPG_BNDS%KFDIA
+YDMF_PHYS%OUT%CT(JLON) = TOTO(JLON)
+ENDDO
+
+!$loki end parallel
+!!!$ACDC }
+
+CALL MF_PHYS_FPL_PART1 (YDCPG_BNDS, YDCPG_OPTS, ZPFL_FPLCH, ZPFL_FPLSH, YDVARS%CPF%T0, YDVARS%SPF%T0, YDMODEL)
+
+IF (LHOOK) CALL DR_HOOK('DISPATCH_ROUTINE',1,ZHOOK_HANDLE)
+
+END SUBROUTINE DISPATCH_ROUTINE