edge_len.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 edge_len
  USE all_vars
  use mod_spherical
  use mod_northpole
  use mod_utils

  implicit none
  INTEGER  :: I,IP,I1,I2,IA,IB,J,J1
  REAL(SP) :: XI, YI, X11,X33,Y11,Y33

  INTEGER :: JTMP,J2



  ! Distance between control volue edge and the node

  DO i=1,ncv_i
     ia=niec(i,1)
     ib=niec(i,2)     
     xi=0.5_sp*(xije(i,1)+xije(i,2))
     yi=0.5_sp*(yije(i,1)+yije(i,2))
     dltxncve(i,1)=xi-vx(ia)
     dltyncve(i,1)=yi-vy(ia)
     dltxncve(i,2)=xi-vx(ib)
     dltyncve(i,2)=yi-vy(ib)

  END DO


  ! Set the distance between Nodes

  IF(maxval(ntsn) > 13) CALL fatal_error &
       & ("THERE ARE MORE THAN 12 NODES AROUND ONE NODE:",&
       "PLEASE INCREASE THE SIZE OF DLTXPI AND DLTYPI IN MOD_MAIN",&
       "BUT REALLY, WHAT IS WRONG WITH YOUR MESH?")

  DO i=1,m
     DO j=1,ntsn(i)-1
        i1=nbsn(i,j)
        i2=nbsn(i,j+1)

        dltytrie(i,j) = vy(i1)-vy(i2)
        dltxtrie(i,j) = vx(i2)-vx(i1)

     END DO
  END DO

  ! Set the distance between Nodes for the North Pole region

  ! Set the distance between triangle edge centers
  IF(maxval(ntve) > 13) CALL fatal_error &
       & ("THERE ARE MORE THAN 12 CELLS AROUND ONE NODE:",&
       "PLEASE INCREASE THE SIZE OF DLVISCXPI AND DLVISCYPI IN MOD_MAIN",&
       "BUT REALLY, WHAT IS WRONG WITH YOUR MESH?")


  DO i=1,m
     DO j=1,ntve(i)
        i1=nbve(i,j)
        jtmp=nbvt(i,j)
        j1=jtmp+1-(jtmp+1)/4*3
        j2=jtmp+2-(jtmp+2)/4*3
        x11=0.5_sp*(vx(i)+vx(nv(i1,j1)))
        y11=0.5_sp*(vy(i)+vy(nv(i1,j1)))
!        X22=XC(I1)
!        Y22=YC(I1)
        x33=0.5_sp*(vx(i)+vx(nv(i1,j2)))
        y33=0.5_sp*(vy(i)+vy(nv(i1,j2)))


        dltyecec(i,j)=(y11-y33)
        dltxecec(i,j)=(x33-x11)



! Set the distance between the node and the edge Center
! NOTE: THE SIGN MATTERS!
         dltynec(i,j)=(vy(i)-y11)
         dltxnec(i,j)=(x11-vx(i))


     END DO
  END DO


   END SUBROUTINE edge_len