edge_len.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	edge_len

Function/Subroutine Documentation

edge_len()

subroutine edge_len

(

)

Definition at line 41 of file edge_len.f90.

  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

Here is the call graph for this function:

Here is the caller graph for this function: