setup_domain.f90

Go to the documentation of this file.

!/===========================================================================/
! Copyright (c) 2007, The University of Massachusetts Dartmouth 
! Produced at the School of Marine Science & Technology 
! Marine Ecosystem Dynamics Modeling group
! All rights reserved.
!
! FVCOM has been developed by the joint UMASSD-WHOI research team. For 
! details of authorship and attribution of credit please see the FVCOM
! technical manual or contact the MEDM group.
!
! 
! This file is part of FVCOM. For details, see http://fvcom.smast.umassd.edu 
! The full copyright notice is contained in the file COPYRIGHT located in the 
! root directory of the FVCOM code. This original header must be maintained
! in all distributed versions.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" 
! AND ANY EXPRESS OR  IMPLIED WARRANTIES, INCLUDING,  BUT NOT  LIMITED TO,
! THE IMPLIED WARRANTIES OF MERCHANTABILITY AND  FITNESS FOR A PARTICULAR
! PURPOSE ARE DISCLAIMED.  
!
!/---------------------------------------------------------------------------/
! CVS VERSION INFORMATION
! $Id$
! $Name$
! $Revision$
!/===========================================================================/

!==============================================================
SUBROUTINE setup_domain
!==============================================================
! 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"
  

END SUBROUTINE setup_domain