mod_spherical::arcx_back Interface Reference

Public Member Functions
subroutine	arcx_back_flt (XX1, YY1, XX2, YY2, ARCX1)
subroutine	arcx_back_dbl (XX1, YY1, XX2, YY2, ARCX1)

Detailed Description

Definition at line 83 of file mod_spherical.f90.

Member Function/Subroutine Documentation

arcx_back_dbl()

subroutine mod_spherical::arcx_back::arcx_back_dbl	(	real(dp), intent(in)	XX1,
		real(dp), intent(in)	YY1,
		real(dp), intent(in)	XX2,
		real(dp), intent(in)	YY2,
		real(dp), intent(out)	ARCX1
	)

Definition at line 319 of file mod_spherical.f90.

    IMPLICIT NONE
    REAL(DP), INTENT(IN) :: XX1,YY1,XX2,YY2
    REAL(DP), INTENT(OUT) :: ARCX1

    INTEGER I
    INTEGER,PARAMETER ::NX=500
    REAL(DP) :: X1,Y1,X2,Y2,TY,A1,A2,B1,B2,C1,C2,A,B,C,X(NX+1),Y(NX+1)
    REAL(DP) :: XTMP          

    IF(xx1 == xx2)THEN
       arcx1=0.
    ELSE
       x1=xx1*deg2rad
       y1=yy1*deg2rad

       x2=xx2*deg2rad
       y2=yy2*deg2rad

       x(1)=x1
       y(1)=y1
       x(nx+1)=x2
       y(nx+1)=y2

       xtmp=x(nx+1)-x(1)
       IF(xtmp >  pi)THEN
          xtmp = real(-2*pi,dp)+xtmp
       ELSE IF(xtmp < -pi)THEN
          xtmp =  real(2*pi,dp)+xtmp
       END IF

       DO i=2,nx
          x(i)=x(i-1)+xtmp/float(nx)
          !       x(i)=x(i-1)+(x(nx+1)-x(1))/float(nx)
       END DO

       a1=dcos(y(1))*dcos(x(1))
       a2=dcos(y(nx+1))*dcos(x(nx+1))

       b1=dcos(y(1))*dsin(x(1))
       b2=dcos(y(nx+1))*dsin(x(nx+1))

       c1=dsin(y(1))
       c2=dsin(y(nx+1))

       a=a1*b2-a2*b1
       b=b1*c2-b2*c1
       c=a2*c1-a1*c2

       DO i=2,nx
          y(i)=-b*dcos(x(i))-c*dsin(x(i))
          y(i)=y(i)/a
          y(i)=datan(y(i))
       END DO

       arcx1=0.
       DO i=1,nx
          ty=0.5*(y(i)+y(i+1))
          xtmp=x(i+1)-x(i)
          IF(xtmp >  pi)THEN
             xtmp = real(-2*pi,dp)+xtmp
          ELSE IF(xtmp < -pi)THEN
             xtmp =  real(2*pi,dp)+xtmp
          END IF
          arcx1=arcx1+rearth*dcos(ty)*xtmp
          !       arcx1=arcx1+rearth*cos(ty)*(x(i+1)-x(i))
       END DO
    END IF

    RETURN

arcx_back_flt()

subroutine mod_spherical::arcx_back::arcx_back_flt	(	real(spa), intent(in)	XX1,
		real(spa), intent(in)	YY1,
		real(spa), intent(in)	XX2,
		real(spa), intent(in)	YY2,
		real(spa), intent(out)	ARCX1
	)

Definition at line 391 of file mod_spherical.f90.

    IMPLICIT NONE
    REAL(SPA), INTENT(IN) :: XX1,YY1,XX2,YY2
    REAL(SPA), INTENT(OUT)::ARCX1
    
    REAL(DP) ::ARCX_DP

    CALL arcx_back_dbl(dble(xx1),dble(yy1),dble(xx2),dble(yy2),arcx_dp)
    arcx1 = arcx_dp

---

The documentation for this interface was generated from the following file:

• FVCOMg01/mod_spherical.f90