local_slope_param.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
SUBROUTINE local_slope_param(JPINDX1,JPINDX2)

!!                03/14  : M Lafaysse, modifs for optimization and parallelization

USE modd_slope_effect, ONLY: xzsl,xzs_xy,xslopang,xslopazi,xsurf_triangle,&
                  xzs_thread,xzs_xy_thread,xslopang_thread,xslopazi_thread,xsurf_triangle_thread,&
                  xxhat,xyhat,xxhat_thread,xyhat_thread,nnxloc,nnyloc,nixloc,niyloc,&
                  lrevertgrid,nix,niy,nindx1_x,nindx2_x,nindx1_y,nindx2_y

USE modd_surfex_mpi, ONLY : nrank, npio, nproc, ncomm, nsize, nindex, nsize_task

IMPLICIT NONE

INTEGER,INTENT(IN)::jpindx1,jpindx2
INTEGER,SAVE::nindx1_mpi_omp,nindx2_mpi_omp,iminmpi
INTEGER::jj


! Il faut déterminer la position de début  du thread MPI par rapport au ZZS1D_FULL de init_slope_param :
!$OMP SINGLE
IF (nproc>1) THEN
  DO jj=1,SIZE(nindex)
    IF(nindex(jj)==nrank) THEN
      iminmpi=jj
      EXIT
    ENDIF
  END DO
ELSE
  iminmpi=1
END IF
!$OMP END SINGLE

nindx1_mpi_omp=iminmpi+jpindx1-1
nindx2_mpi_omp=iminmpi+jpindx2-1


IF (lrevertgrid) THEN
  nindx1_x=mod(nindx1_mpi_omp-1,nix)+1
  nindx2_x=mod(nindx2_mpi_omp-1,nix)+3
  nindx1_y=niy-(nindx2_mpi_omp-1)/nix
  nindx2_y=niy-(nindx1_mpi_omp-1)/nix+2
ELSE
  nindx1_x=mod(nindx1_mpi_omp-1,nix)+1
  nindx2_x=mod(nindx2_mpi_omp-1,nix)+3
  nindx1_y=(nindx1_mpi_omp-1)/nix+1
  nindx2_y=(nindx2_mpi_omp-1)/nix+3
ENDIF

nnxloc=nindx2_x-nindx1_x+1
nnyloc=nindx2_y-nindx1_y+1
nixloc=nnxloc-2
niyloc=nnyloc-2

ALLOCATE(xzs_thread(nnxloc,nnyloc))
ALLOCATE(xzs_xy_thread(nnxloc,nnyloc))
ALLOCATE(xslopang_thread(nnxloc,nnyloc,4))
ALLOCATE(xslopazi_thread(nnxloc,nnyloc,4))
ALLOCATE(xsurf_triangle_thread(nnxloc,nnyloc,4))

ALLOCATE(xxhat_thread(nnxloc))
ALLOCATE(xyhat_thread(nnyloc))


xxhat_thread(:)=xxhat(nindx1_x:nindx2_x)
xyhat_thread(:)=xyhat(nindx1_y:nindx2_y)

xzs_thread=xzsl(nindx1_x:nindx2_x,nindx1_y:nindx2_y)
xzs_xy_thread=xzs_xy(nindx1_x:nindx2_x,nindx1_y:nindx2_y)
xslopang_thread=xslopang(nindx1_x:nindx2_x,nindx1_y:nindx2_y,:)
xslopazi_thread=xslopazi(nindx1_x:nindx2_x,nindx1_y:nindx2_y,:)
xsurf_triangle_thread=xsurf_triangle(nindx1_x:nindx2_x,nindx1_y:nindx2_y,:)


END SUBROUTINE local_slope_param