setup_domain.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	setup_domain

Function/Subroutine Documentation

setup_domain()

subroutine setup_domain

(

)

Definition at line 42 of file setup_domain.f90.

!==============================================================
! IF(PAR): 
!          Uses METIS to determine local domains
!          Create global to local and local to global maps
!==============================================================
  USE mod_par
  USE all_vars
  USE mod_ncdio
  USE mod_setup
  
  !==============================================================================|
  !  GENERATE LOCAL NODE CONNECTIVITY (NV) FROM GLOBAL NODE CONNECTIVITY (NVG)   |
  !  USING LOCAL TO GLOBAL MAPPING FOR INTERIOR ELEMENTS (EGID)                  |
  !  AND LOCAL TO GLOBAL MAPPING FOR HALO ELEMENTS (HE_LST)                      |
  !  TRANSFER GLOBAL COORDINATES TO LOCAL DOMAINS
  !==============================================================================|

  IMPLICIT NONE
  INTEGER :: SENDER,STATUS
  INTEGER :: nvals, IERR, I, EGL
  

     if (dbg_set(dbg_log)) &
          & write(ipt,*) "! SETTING UP FVCOM PROCESSOR ENVIRONMENT"
  
  ! SETUP OUTPUT MODE OR TURN IT OFF
  CALL setup_mpi_io_mode(use_mpi_io_mode,mpi_fvcom_group)
  
  IF(serial) THEN
     
     ! DO NOT ALLOCATE EL_PID FOR NON-MULTIPROCESSOR CASES
          
     ! SHIFT GRID DIMENSIONS TO LOCAL FOR ANY SERIAL CASE
     mt=mgl
     nt=ngl
     
     m = mgl
     n = ngl
     
     ! TRANSFER CONNECTIVTY TO LOCAL DOMAIN NAME VARIABLES
     ALLOCATE(nv(0:nt,4))
     nv   = nvg
     nv(0,:) = 0
     !DEALLOCATE(NVG)
     

     ! MUST ALLOCATE EL_PID FOR OUTPUT OF DATA
     ALLOCATE(el_pid(ngl))  ;  el_pid = 1


     ! SETUP NGID
     ALLOCATE(ngid(0:mgl))  ; 
     DO i=0,mgl
        ngid(i)=i
     END DO

     ngid_x => ngid

     ! SETUP NLID
     ALLOCATE(nlid(0:mgl))
     DO i=0,mgl
        nlid(i) = i
     END DO

     nlid_x => nlid

     ! SETUP EGID
     ALLOCATE(egid(0:ngl))  ; 
     DO i=0,ngl
        egid(i)=i
     END DO

     egid_x => egid

     ! SETUP ELID
     ALLOCATE(elid(0:ngl))
     DO i=0,ngl
        elid(i) = i
     END DO

     elid_x => elid

     ALLOCATE(nde_id(mgl)); nde_id=0

  ELSE ! IF PAR

  END IF
  

  ! PASS OBC TO LOCAL VARIABLES AND GENERATE MAP FROM GLOBAL

  IF(obc_longshore_flow_on) CALL genmap_lsf

  IF(obc_on) CALL genmap_obc


  IF (use_mpi_io_mode) then
     CALL fatal_error("'USE_MPI_IO_MODE' SHOULD NEVER BE TRUE WHEN C&
          &OMPILED WITHOUT MULTIPROCESSOR")
  END IF
  
  if(dbg_set(dbg_log))  write(ipt,*) "! Finished SETUP_LOCAL_DOMAINS"
  

Here is the call graph for this function:

Here is the caller graph for this function: