module finalist_internal_11
  use finalist_header
contains
  subroutine Formura_internal_258 ()
    double precision  :: a
    double precision  :: a_0
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(0), iz+(0))=Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_m1_Om_87_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(0), iz+(36))=Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_87_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(0), iz+(2))=Facet_m1_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_87_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(36), iz+(0))=Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_87_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(2), iz+(0))=Facet_m1_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_87_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(36), iz+(36))=Facet_m1_0_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_87_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(2), iz+(36))=Facet_m1_0_0_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_87_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(36), iz+(2))=Facet_m1_0_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_87_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(2), iz+(2))=Facet_m1_0_0_src__1_0_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_87_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(0), iz+(0))=Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_87_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc0(ix+(2), iy+(0), iz+(0))=Facet_0_m1_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_87_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(0), iz+(36))=Facet_0_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_87_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc0(ix+(2), iy+(0), iz+(36))=Facet_0_m1_0_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_87_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(0), iz+(2))=Facet_0_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_87_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc0(ix+(2), iy+(0), iz+(2))=Facet_0_m1_0_src__0_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_87_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(36), iz+(0))=Facet_0_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_87_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc0(ix+(2), iy+(36), iz+(0))=Facet_0_0_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_87_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(2), iz+(0))=Facet_0_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_87_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  258
             Rsc0(ix+(2), iy+(2), iz+(0))=Facet_0_0_m1_src__0_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_87_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(36), iz+(36))=Ridge_0_0_0_Om_87_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc0(ix+(2), iy+(36), iz+(36))=Ridge_0_0_0_Om_87_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(2), iz+(36))=Ridge_0_0_0_Om_87_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  258
             Rsc0(ix+(2), iy+(2), iz+(36))=Ridge_0_0_0_Om_87_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(36), iz+(2))=Ridge_0_0_0_Om_87_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc0(ix+(2), iy+(36), iz+(2))=Ridge_0_0_0_Om_87_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(2), iz+(2))=Ridge_0_0_0_Om_87_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(0), iz+(0))=Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_m1_Om_88_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(0), iz+(36))=Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_88_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(0), iz+(2))=Facet_m1_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_88_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(36), iz+(0))=Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_88_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(2), iz+(0))=Facet_m1_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_88_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(36), iz+(36))=Facet_m1_0_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_88_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(2), iz+(36))=Facet_m1_0_0_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_88_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(36), iz+(2))=Facet_m1_0_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_88_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(2), iz+(2))=Facet_m1_0_0_src__1_0_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_88_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(0), iz+(0))=Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_88_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc1(ix+(2), iy+(0), iz+(0))=Facet_0_m1_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_88_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(0), iz+(36))=Facet_0_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_88_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc1(ix+(2), iy+(0), iz+(36))=Facet_0_m1_0_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_88_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(0), iz+(2))=Facet_0_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_88_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc1(ix+(2), iy+(0), iz+(2))=Facet_0_m1_0_src__0_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_88_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(36), iz+(0))=Facet_0_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_88_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc1(ix+(2), iy+(36), iz+(0))=Facet_0_0_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_88_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(2), iz+(0))=Facet_0_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_88_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  258
             Rsc1(ix+(2), iy+(2), iz+(0))=Facet_0_0_m1_src__0_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_88_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(36), iz+(36))=Ridge_0_0_0_Om_88_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc1(ix+(2), iy+(36), iz+(36))=Ridge_0_0_0_Om_88_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(2), iz+(36))=Ridge_0_0_0_Om_88_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  258
             Rsc1(ix+(2), iy+(2), iz+(36))=Ridge_0_0_0_Om_88_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(36), iz+(2))=Ridge_0_0_0_Om_88_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  258
             Rsc1(ix+(2), iy+(36), iz+(2))=Ridge_0_0_0_Om_88_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(2), iz+(2))=Ridge_0_0_0_Om_88_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz, a, a_0)
    do  iz = 3 ,  38 ,  1

       do  iy = 3 ,  38 ,  1

          do  ix = 3 ,  262 ,  1

             a=Rsc0(ix+(-1), iy+(-1), iz+(-1))
             a_0=Rsc1(ix+(-1), iy+(-1), iz+(-1))

 Rsc3(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( (-( ( ( ( 1.0 / 900.0 ) * a ) * a_0 ) * a_0 )) + ( ( 1.0 / 86400.0 ) * ( 1.0 - a ) ) ) + ( ( ( 0.1 * 2.3e-9 ) * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc0(ix+(0), iy+(-1), iz+(-1)) - a ) - ( a - &
Rsc0(ix+(-2), iy+(-1), iz+(-1)) ) ) + ( ( Rsc0(ix+(-1), iy+(0), iz+(-1)) - a ) - ( a - Rsc0(ix+(-1), iy+(-2), iz+(-1)) ) ) ) + ( ( Rsc0(ix+(-1), iy+(-1), iz+(0)) - a ) - ( a - Rsc0(ix+(-1), iy+(-1), iz+(-2)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_258

  subroutine Formura_internal_259 ()
    double precision  :: a
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz, a)
    do  iz = 3 ,  38 ,  1

       do  iy = 3 ,  38 ,  1

          do  ix = 3 ,  262 ,  1

             a=Rsc1(ix+(-1), iy+(-1), iz+(-1))

 Rsc2(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( ( ( ( ( 1.0 / 900.0 ) * Rsc0(ix+(-1), iy+(-1), iz+(-1)) ) * a ) * a ) - ( ( 6.0 / 86400.0 ) * a ) ) + ( ( 6.1e-11 * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc1(ix+(0), iy+(-1), iz+(-1)) - a ) - ( a - &
Rsc1(ix+(-2), iy+(-1), iz+(-1)) ) ) + ( ( Rsc1(ix+(-1), iy+(0), iz+(-1)) - a ) - ( a - Rsc1(ix+(-1), iy+(-2), iz+(-1)) ) ) ) + ( ( Rsc1(ix+(-1), iy+(-1), iz+(0)) - a ) - ( a - Rsc1(ix+(-1), iy+(-1), iz+(-2)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_259

  subroutine Formura_internal_26 ()
    double precision  :: a
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz, a)
    do  iz = 1 ,  36 ,  1

       do  iy = 1 ,  36 ,  1

          do  ix = 3 ,  254 ,  1

             a=Rsc1(ix+(-1), iy+(1), iz+(1))

 Rsc2(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( ( ( ( ( 1.0 / 900.0 ) * Rsc0(ix+(-1), iy+(1), iz+(1)) ) * a ) * a ) - ( ( 6.0 / 86400.0 ) * a ) ) + ( ( 6.1e-11 * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc1(ix+(0), iy+(1), iz+(1)) - a ) - ( a - &
Rsc1(ix+(-2), iy+(1), iz+(1)) ) ) + ( ( Rsc1(ix+(-1), iy+(2), iz+(1)) - a ) - ( a - Rsc1(ix+(-1), iy+(0), iz+(1)) ) ) ) + ( ( Rsc1(ix+(-1), iy+(1), iz+(2)) - a ) - ( a - Rsc1(ix+(-1), iy+(1), iz+(0)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  252
             Ridge_0_0_0_Om_82_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc2(ix+(2), iy+(34), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  252
             Ridge_0_0_0_Om_82_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(34))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  252
             Ridge_0_0_0_Om_82_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(2), iy+(34), iz+(34))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  252
             Facet_0_0_1_src__1_0_0_dest__1_0_1_Send%Ridge_0_0_1_Om_82_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  252
             Facet_0_0_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_0_1_Om_82_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(2), iy+(34), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  252
             Facet_0_1_0_src__1_0_0_dest__1_1_0_Send%Ridge_0_1_0_Om_82_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  252
             Facet_0_1_0_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_0_Om_82_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(34))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  252
             Facet_0_1_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_1_Om_82_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_26

  subroutine Formura_internal_260 ()
    double precision  :: a
    double precision  :: a_0
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(0), iz+(0))=Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_m1_Om_89_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(0), iz+(38))=Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_89_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(0), iz+(2))=Facet_m1_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_89_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(38), iz+(0))=Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_89_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(2), iz+(0))=Facet_m1_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_89_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(38), iz+(38))=Facet_m1_0_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_89_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(2), iz+(38))=Facet_m1_0_0_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_89_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(38), iz+(2))=Facet_m1_0_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_89_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(2), iz+(2))=Facet_m1_0_0_src__1_0_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_89_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(262), iy+(0), iz+(0))=Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_89_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc3(ix+(2), iy+(0), iz+(0))=Facet_0_m1_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_89_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(262), iy+(0), iz+(38))=Facet_0_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_89_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc3(ix+(2), iy+(0), iz+(38))=Facet_0_m1_0_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_89_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(262), iy+(0), iz+(2))=Facet_0_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_89_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc3(ix+(2), iy+(0), iz+(2))=Facet_0_m1_0_src__0_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_89_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(262), iy+(38), iz+(0))=Facet_0_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_89_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc3(ix+(2), iy+(38), iz+(0))=Facet_0_0_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_89_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc3(ix+(262), iy+(2), iz+(0))=Facet_0_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_89_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  260
             Rsc3(ix+(2), iy+(2), iz+(0))=Facet_0_0_m1_src__0_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_89_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(262), iy+(38), iz+(38))=Ridge_0_0_0_Om_89_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc3(ix+(2), iy+(38), iz+(38))=Ridge_0_0_0_Om_89_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc3(ix+(262), iy+(2), iz+(38))=Ridge_0_0_0_Om_89_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  260
             Rsc3(ix+(2), iy+(2), iz+(38))=Ridge_0_0_0_Om_89_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc3(ix+(262), iy+(38), iz+(2))=Ridge_0_0_0_Om_89_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc3(ix+(2), iy+(38), iz+(2))=Ridge_0_0_0_Om_89_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc3(ix+(262), iy+(2), iz+(2))=Ridge_0_0_0_Om_89_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(0), iz+(0))=Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_m1_Om_90_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(0), iz+(38))=Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_90_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(0), iz+(2))=Facet_m1_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_90_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(38), iz+(0))=Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_90_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(2), iz+(0))=Facet_m1_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_90_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(38), iz+(38))=Facet_m1_0_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_90_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(2), iz+(38))=Facet_m1_0_0_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_90_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(38), iz+(2))=Facet_m1_0_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_90_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(2), iz+(2))=Facet_m1_0_0_src__1_0_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_90_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(262), iy+(0), iz+(0))=Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_90_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc2(ix+(2), iy+(0), iz+(0))=Facet_0_m1_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_90_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(262), iy+(0), iz+(38))=Facet_0_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_90_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc2(ix+(2), iy+(0), iz+(38))=Facet_0_m1_0_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_90_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(262), iy+(0), iz+(2))=Facet_0_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_90_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc2(ix+(2), iy+(0), iz+(2))=Facet_0_m1_0_src__0_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_90_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(262), iy+(38), iz+(0))=Facet_0_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_90_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc2(ix+(2), iy+(38), iz+(0))=Facet_0_0_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_90_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc2(ix+(262), iy+(2), iz+(0))=Facet_0_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_90_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  260
             Rsc2(ix+(2), iy+(2), iz+(0))=Facet_0_0_m1_src__0_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_90_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(262), iy+(38), iz+(38))=Ridge_0_0_0_Om_90_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc2(ix+(2), iy+(38), iz+(38))=Ridge_0_0_0_Om_90_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc2(ix+(262), iy+(2), iz+(38))=Ridge_0_0_0_Om_90_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  36
          do ix = 1 ,  260
             Rsc2(ix+(2), iy+(2), iz+(38))=Ridge_0_0_0_Om_90_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc2(ix+(262), iy+(38), iz+(2))=Ridge_0_0_0_Om_90_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  2
          do ix = 1 ,  260
             Rsc2(ix+(2), iy+(38), iz+(2))=Ridge_0_0_0_Om_90_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc2(ix+(262), iy+(2), iz+(2))=Ridge_0_0_0_Om_90_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz, a, a_0)
    do  iz = 3 ,  40 ,  1

       do  iy = 3 ,  40 ,  1

          do  ix = 3 ,  264 ,  1

             a=Rsc3(ix+(-1), iy+(-1), iz+(-1))
             a_0=Rsc2(ix+(-1), iy+(-1), iz+(-1))

 Rsc1(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( (-( ( ( ( 1.0 / 900.0 ) * a ) * a_0 ) * a_0 )) + ( ( 1.0 / 86400.0 ) * ( 1.0 - a ) ) ) + ( ( ( 0.1 * 2.3e-9 ) * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc3(ix+(0), iy+(-1), iz+(-1)) - a ) - ( a - &
Rsc3(ix+(-2), iy+(-1), iz+(-1)) ) ) + ( ( Rsc3(ix+(-1), iy+(0), iz+(-1)) - a ) - ( a - Rsc3(ix+(-1), iy+(-2), iz+(-1)) ) ) ) + ( ( Rsc3(ix+(-1), iy+(-1), iz+(0)) - a ) - ( a - Rsc3(ix+(-1), iy+(-1), iz+(-2)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_260

  subroutine Formura_internal_261 ()
    double precision  :: a
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz, a)
    do  iz = 3 ,  40 ,  1

       do  iy = 3 ,  40 ,  1

          do  ix = 3 ,  264 ,  1

             a=Rsc2(ix+(-1), iy+(-1), iz+(-1))

 Rsc0(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( ( ( ( ( 1.0 / 900.0 ) * Rsc3(ix+(-1), iy+(-1), iz+(-1)) ) * a ) * a ) - ( ( 6.0 / 86400.0 ) * a ) ) + ( ( 6.1e-11 * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc2(ix+(0), iy+(-1), iz+(-1)) - a ) - ( a - &
Rsc2(ix+(-2), iy+(-1), iz+(-1)) ) ) + ( ( Rsc2(ix+(-1), iy+(0), iz+(-1)) - a ) - ( a - Rsc2(ix+(-1), iy+(-2), iz+(-1)) ) ) ) + ( ( Rsc2(ix+(-1), iy+(-1), iz+(0)) - a ) - ( a - Rsc2(ix+(-1), iy+(-1), iz+(-2)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_261

  subroutine Formura_internal_262 ()
    double precision  :: a
    double precision  :: a_0
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(0), iz+(0))=Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_m1_Om_91_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(0), iz+(40))=Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_91_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(0), iz+(2))=Facet_m1_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_91_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(40), iz+(0))=Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_91_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(2), iz+(0))=Facet_m1_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_91_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(40), iz+(40))=Facet_m1_0_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_91_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(2), iz+(40))=Facet_m1_0_0_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_91_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(40), iz+(2))=Facet_m1_0_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_91_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(2), iz+(2))=Facet_m1_0_0_src__1_0_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_91_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(264), iy+(0), iz+(0))=Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_91_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc1(ix+(2), iy+(0), iz+(0))=Facet_0_m1_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_91_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(264), iy+(0), iz+(40))=Facet_0_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_91_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc1(ix+(2), iy+(0), iz+(40))=Facet_0_m1_0_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_91_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(264), iy+(0), iz+(2))=Facet_0_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_91_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc1(ix+(2), iy+(0), iz+(2))=Facet_0_m1_0_src__0_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_91_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(264), iy+(40), iz+(0))=Facet_0_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_91_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc1(ix+(2), iy+(40), iz+(0))=Facet_0_0_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_91_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc1(ix+(264), iy+(2), iz+(0))=Facet_0_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_91_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  262
             Rsc1(ix+(2), iy+(2), iz+(0))=Facet_0_0_m1_src__0_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_91_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(264), iy+(40), iz+(40))=Ridge_0_0_0_Om_91_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc1(ix+(2), iy+(40), iz+(40))=Ridge_0_0_0_Om_91_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc1(ix+(264), iy+(2), iz+(40))=Ridge_0_0_0_Om_91_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  262
             Rsc1(ix+(2), iy+(2), iz+(40))=Ridge_0_0_0_Om_91_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc1(ix+(264), iy+(40), iz+(2))=Ridge_0_0_0_Om_91_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc1(ix+(2), iy+(40), iz+(2))=Ridge_0_0_0_Om_91_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc1(ix+(264), iy+(2), iz+(2))=Ridge_0_0_0_Om_91_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(0), iz+(0))=Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_m1_Om_92_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(0), iz+(40))=Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_92_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(0), iz+(2))=Facet_m1_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_m1_0_Om_92_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(40), iz+(0))=Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_92_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(2), iz+(0))=Facet_m1_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_m1_Om_92_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(40), iz+(40))=Facet_m1_0_0_src__1_1_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_92_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(2), iz+(40))=Facet_m1_0_0_src__1_0_1_dest__0_0_0_Recv%Ridge_m1_0_0_Om_92_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(40), iz+(2))=Facet_m1_0_0_src__1_1_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_92_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(2), iz+(2))=Facet_m1_0_0_src__1_0_0_dest__0_0_0_Recv%Ridge_m1_0_0_Om_92_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(264), iy+(0), iz+(0))=Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_92_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc0(ix+(2), iy+(0), iz+(0))=Facet_0_m1_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_m1_Om_92_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(264), iy+(0), iz+(40))=Facet_0_m1_0_src__1_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_92_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc0(ix+(2), iy+(0), iz+(40))=Facet_0_m1_0_src__0_1_1_dest__0_0_0_Recv%Ridge_0_m1_0_Om_92_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(264), iy+(0), iz+(2))=Facet_0_m1_0_src__1_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_92_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc0(ix+(2), iy+(0), iz+(2))=Facet_0_m1_0_src__0_1_0_dest__0_0_0_Recv%Ridge_0_m1_0_Om_92_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(264), iy+(40), iz+(0))=Facet_0_0_m1_src__1_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_92_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc0(ix+(2), iy+(40), iz+(0))=Facet_0_0_m1_src__0_1_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_92_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc0(ix+(264), iy+(2), iz+(0))=Facet_0_0_m1_src__1_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_92_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  262
             Rsc0(ix+(2), iy+(2), iz+(0))=Facet_0_0_m1_src__0_0_1_dest__0_0_0_Recv%Ridge_0_0_m1_Om_92_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(264), iy+(40), iz+(40))=Ridge_0_0_0_Om_92_r__1_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc0(ix+(2), iy+(40), iz+(40))=Ridge_0_0_0_Om_92_r__0_1_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc0(ix+(264), iy+(2), iz+(40))=Ridge_0_0_0_Om_92_r__1_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  38
          do ix = 1 ,  262
             Rsc0(ix+(2), iy+(2), iz+(40))=Ridge_0_0_0_Om_92_r__0_0_1_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  2
             Rsc0(ix+(264), iy+(40), iz+(2))=Ridge_0_0_0_Om_92_r__1_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  2
          do ix = 1 ,  262
             Rsc0(ix+(2), iy+(40), iz+(2))=Ridge_0_0_0_Om_92_r__0_1_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  38
       do iy = 1 ,  38
          do ix = 1 ,  2
             Rsc0(ix+(264), iy+(2), iz+(2))=Ridge_0_0_0_Om_92_r__1_0_0_r__0_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz, a, a_0)
    do  iz = -5 ,  34 ,  1

       do  iy = -5 ,  34 ,  1

          do  ix = -5 ,  258 ,  1

             a=Rsc1(ix+(7), iy+(7), iz+(7))
             a_0=Rsc0(ix+(7), iy+(7), iz+(7))

 U(ix+(8)+(0), iy+(8)+(0), iz+(8)+(0))=( a + ( 200.0 * ( ( (-( ( ( ( 1.0 / 900.0 ) * a ) * a_0 ) * a_0 )) + ( ( 1.0 / 86400.0 ) * ( 1.0 - a ) ) ) + ( ( ( 0.1 * 2.3e-9 ) * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc1(ix+(8), iy+(7), iz+(7)) - a ) - ( a - &
Rsc1(ix+(6), iy+(7), iz+(7)) ) ) + ( ( Rsc1(ix+(7), iy+(8), iz+(7)) - a ) - ( a - Rsc1(ix+(7), iy+(6), iz+(7)) ) ) ) + ( ( Rsc1(ix+(7), iy+(7), iz+(8)) - a ) - ( a - Rsc1(ix+(7), iy+(7), iz+(6)) ) ) ) ) ) ) )

          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_262

  subroutine Formura_internal_263 ()
    double precision  :: a
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz, a)
    do  iz = -5 ,  34 ,  1

       do  iy = -5 ,  34 ,  1

          do  ix = -5 ,  258 ,  1

             a=Rsc0(ix+(7), iy+(7), iz+(7))

 V(ix+(8)+(0), iy+(8)+(0), iz+(8)+(0))=( a + ( 200.0 * ( ( ( ( ( ( 1.0 / 900.0 ) * Rsc1(ix+(7), iy+(7), iz+(7)) ) * a ) * a ) - ( ( 6.0 / 86400.0 ) * a ) ) + ( ( 6.1e-11 * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc0(ix+(8), iy+(7), iz+(7)) - a ) - ( a - &
Rsc0(ix+(6), iy+(7), iz+(7)) ) ) + ( ( Rsc0(ix+(7), iy+(8), iz+(7)) - a ) - ( a - Rsc0(ix+(7), iy+(6), iz+(7)) ) ) ) + ( ( Rsc0(ix+(7), iy+(7), iz+(8)) - a ) - ( a - Rsc0(ix+(7), iy+(7), iz+(6)) ) ) ) ) ) ) )

          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_263

  subroutine Formura_internal_27 ()
    double precision  :: a
    double precision  :: a_0
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(0), iz+(0))=Ridge_0_0_0_Om_81_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc3(ix+(254), iy+(0), iz+(0))=Facet_1_0_0_src__0_0_0_dest__1_0_0_Recv%Ridge_1_0_0_Om_81_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(0), iz+(0))=Ridge_0_0_0_Om_82_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  36
       do iy = 1 ,  36
          do ix = 1 ,  2
             Rsc2(ix+(254), iy+(0), iz+(0))=Facet_1_0_0_src__0_0_0_dest__1_0_0_Recv%Ridge_1_0_0_Om_82_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz, a, a_0)
    do  iz = 1 ,  34 ,  1

       do  iy = 1 ,  34 ,  1

          do  ix = 3 ,  256 ,  1

             a=Rsc3(ix+(-1), iy+(1), iz+(1))
             a_0=Rsc2(ix+(-1), iy+(1), iz+(1))

 Rsc1(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( (-( ( ( ( 1.0 / 900.0 ) * a ) * a_0 ) * a_0 )) + ( ( 1.0 / 86400.0 ) * ( 1.0 - a ) ) ) + ( ( ( 0.1 * 2.3e-9 ) * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc3(ix+(0), iy+(1), iz+(1)) - a ) - ( a - &
Rsc3(ix+(-2), iy+(1), iz+(1)) ) ) + ( ( Rsc3(ix+(-1), iy+(2), iz+(1)) - a ) - ( a - Rsc3(ix+(-1), iy+(0), iz+(1)) ) ) ) + ( ( Rsc3(ix+(-1), iy+(1), iz+(2)) - a ) - ( a - Rsc3(ix+(-1), iy+(1), iz+(0)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  254
             Ridge_0_0_0_Om_83_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc1(ix+(2), iy+(32), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  254
             Ridge_0_0_0_Om_83_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(32))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  254
             Ridge_0_0_0_Om_83_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc1(ix+(2), iy+(32), iz+(32))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  254
             Facet_0_0_1_src__1_0_0_dest__1_0_1_Send%Ridge_0_0_1_Om_83_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  254
             Facet_0_0_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_0_1_Om_83_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc1(ix+(2), iy+(32), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  254
             Facet_0_1_0_src__1_0_0_dest__1_1_0_Send%Ridge_0_1_0_Om_83_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  254
             Facet_0_1_0_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_0_Om_83_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(32))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  254
             Facet_0_1_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_1_Om_83_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_27

  subroutine Formura_internal_28 ()
    double precision  :: a
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz, a)
    do  iz = 1 ,  34 ,  1

       do  iy = 1 ,  34 ,  1

          do  ix = 3 ,  256 ,  1

             a=Rsc2(ix+(-1), iy+(1), iz+(1))

 Rsc0(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( ( ( ( ( 1.0 / 900.0 ) * Rsc3(ix+(-1), iy+(1), iz+(1)) ) * a ) * a ) - ( ( 6.0 / 86400.0 ) * a ) ) + ( ( 6.1e-11 * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc2(ix+(0), iy+(1), iz+(1)) - a ) - ( a - &
Rsc2(ix+(-2), iy+(1), iz+(1)) ) ) + ( ( Rsc2(ix+(-1), iy+(2), iz+(1)) - a ) - ( a - Rsc2(ix+(-1), iy+(0), iz+(1)) ) ) ) + ( ( Rsc2(ix+(-1), iy+(1), iz+(2)) - a ) - ( a - Rsc2(ix+(-1), iy+(1), iz+(0)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  254
             Ridge_0_0_0_Om_84_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc0(ix+(2), iy+(32), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  254
             Ridge_0_0_0_Om_84_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(32))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  254
             Ridge_0_0_0_Om_84_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc0(ix+(2), iy+(32), iz+(32))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  34
          do ix = 1 ,  254
             Facet_0_0_1_src__1_0_0_dest__1_0_1_Send%Ridge_0_0_1_Om_84_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  254
             Facet_0_0_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_0_1_Om_84_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc0(ix+(2), iy+(32), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  2
          do ix = 1 ,  254
             Facet_0_1_0_src__1_0_0_dest__1_1_0_Send%Ridge_0_1_0_Om_84_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  254
             Facet_0_1_0_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_0_Om_84_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(32))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  254
             Facet_0_1_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_1_Om_84_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_28

  subroutine Formura_internal_29 ()
    double precision  :: a
    double precision  :: a_0
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(0), iz+(0))=Ridge_0_0_0_Om_83_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc1(ix+(256), iy+(0), iz+(0))=Facet_1_0_0_src__0_0_0_dest__1_0_0_Recv%Ridge_1_0_0_Om_83_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(0), iz+(0))=Ridge_0_0_0_Om_84_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  34
       do iy = 1 ,  34
          do ix = 1 ,  2
             Rsc0(ix+(256), iy+(0), iz+(0))=Facet_1_0_0_src__0_0_0_dest__1_0_0_Recv%Ridge_1_0_0_Om_84_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz, a, a_0)
    do  iz = 1 ,  32 ,  1

       do  iy = 1 ,  32 ,  1

          do  ix = 3 ,  258 ,  1

             a=Rsc1(ix+(-1), iy+(1), iz+(1))
             a_0=Rsc0(ix+(-1), iy+(1), iz+(1))

 Rsc2(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( (-( ( ( ( 1.0 / 900.0 ) * a ) * a_0 ) * a_0 )) + ( ( 1.0 / 86400.0 ) * ( 1.0 - a ) ) ) + ( ( ( 0.1 * 2.3e-9 ) * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc1(ix+(0), iy+(1), iz+(1)) - a ) - ( a - &
Rsc1(ix+(-2), iy+(1), iz+(1)) ) ) + ( ( Rsc1(ix+(-1), iy+(2), iz+(1)) - a ) - ( a - Rsc1(ix+(-1), iy+(0), iz+(1)) ) ) ) + ( ( Rsc1(ix+(-1), iy+(1), iz+(2)) - a ) - ( a - Rsc1(ix+(-1), iy+(1), iz+(0)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  32
       do iy = 1 ,  2
          do ix = 1 ,  256
             Ridge_0_0_0_Om_85_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc2(ix+(2), iy+(30), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  32
          do ix = 1 ,  256
             Ridge_0_0_0_Om_85_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(30))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  256
             Ridge_0_0_0_Om_85_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(2), iy+(30), iz+(30))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  32
          do ix = 1 ,  256
             Facet_0_0_1_src__1_0_0_dest__1_0_1_Send%Ridge_0_0_1_Om_85_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  256
             Facet_0_0_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_0_1_Om_85_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(2), iy+(30), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  32
       do iy = 1 ,  2
          do ix = 1 ,  256
             Facet_0_1_0_src__1_0_0_dest__1_1_0_Send%Ridge_0_1_0_Om_85_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  256
             Facet_0_1_0_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_0_Om_85_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(30))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  256
             Facet_0_1_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_1_Om_85_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_29

  subroutine Formura_internal_30 ()
    double precision  :: a
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz, a)
    do  iz = 1 ,  32 ,  1

       do  iy = 1 ,  32 ,  1

          do  ix = 3 ,  258 ,  1

             a=Rsc0(ix+(-1), iy+(1), iz+(1))

 Rsc3(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( ( ( ( ( 1.0 / 900.0 ) * Rsc1(ix+(-1), iy+(1), iz+(1)) ) * a ) * a ) - ( ( 6.0 / 86400.0 ) * a ) ) + ( ( 6.1e-11 * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc0(ix+(0), iy+(1), iz+(1)) - a ) - ( a - &
Rsc0(ix+(-2), iy+(1), iz+(1)) ) ) + ( ( Rsc0(ix+(-1), iy+(2), iz+(1)) - a ) - ( a - Rsc0(ix+(-1), iy+(0), iz+(1)) ) ) ) + ( ( Rsc0(ix+(-1), iy+(1), iz+(2)) - a ) - ( a - Rsc0(ix+(-1), iy+(1), iz+(0)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  32
       do iy = 1 ,  2
          do ix = 1 ,  256
             Ridge_0_0_0_Om_86_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc3(ix+(2), iy+(30), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  32
          do ix = 1 ,  256
             Ridge_0_0_0_Om_86_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(30))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  256
             Ridge_0_0_0_Om_86_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc3(ix+(2), iy+(30), iz+(30))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  32
          do ix = 1 ,  256
             Facet_0_0_1_src__1_0_0_dest__1_0_1_Send%Ridge_0_0_1_Om_86_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  256
             Facet_0_0_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_0_1_Om_86_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc3(ix+(2), iy+(30), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  32
       do iy = 1 ,  2
          do ix = 1 ,  256
             Facet_0_1_0_src__1_0_0_dest__1_1_0_Send%Ridge_0_1_0_Om_86_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  256
             Facet_0_1_0_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_0_Om_86_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(30))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  256
             Facet_0_1_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_1_Om_86_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_30

  subroutine Formura_internal_31 ()
    double precision  :: a
    double precision  :: a_0
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  32
       do iy = 1 ,  32
          do ix = 1 ,  2
             Rsc2(ix+(0), iy+(0), iz+(0))=Ridge_0_0_0_Om_85_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  32
       do iy = 1 ,  32
          do ix = 1 ,  2
             Rsc2(ix+(258), iy+(0), iz+(0))=Facet_1_0_0_src__0_0_0_dest__1_0_0_Recv%Ridge_1_0_0_Om_85_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  32
       do iy = 1 ,  32
          do ix = 1 ,  2
             Rsc3(ix+(0), iy+(0), iz+(0))=Ridge_0_0_0_Om_86_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  32
       do iy = 1 ,  32
          do ix = 1 ,  2
             Rsc3(ix+(258), iy+(0), iz+(0))=Facet_1_0_0_src__0_0_0_dest__1_0_0_Recv%Ridge_1_0_0_Om_86_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz, a, a_0)
    do  iz = 1 ,  30 ,  1

       do  iy = 1 ,  30 ,  1

          do  ix = 3 ,  260 ,  1

             a=Rsc2(ix+(-1), iy+(1), iz+(1))
             a_0=Rsc3(ix+(-1), iy+(1), iz+(1))

 Rsc0(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( (-( ( ( ( 1.0 / 900.0 ) * a ) * a_0 ) * a_0 )) + ( ( 1.0 / 86400.0 ) * ( 1.0 - a ) ) ) + ( ( ( 0.1 * 2.3e-9 ) * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc2(ix+(0), iy+(1), iz+(1)) - a ) - ( a - &
Rsc2(ix+(-2), iy+(1), iz+(1)) ) ) + ( ( Rsc2(ix+(-1), iy+(2), iz+(1)) - a ) - ( a - Rsc2(ix+(-1), iy+(0), iz+(1)) ) ) ) + ( ( Rsc2(ix+(-1), iy+(1), iz+(2)) - a ) - ( a - Rsc2(ix+(-1), iy+(1), iz+(0)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  30
       do iy = 1 ,  2
          do ix = 1 ,  258
             Ridge_0_0_0_Om_87_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc0(ix+(2), iy+(28), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  30
          do ix = 1 ,  258
             Ridge_0_0_0_Om_87_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(28))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Ridge_0_0_0_Om_87_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc0(ix+(2), iy+(28), iz+(28))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  30
          do ix = 1 ,  258
             Facet_0_0_1_src__1_0_0_dest__1_0_1_Send%Ridge_0_0_1_Om_87_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Facet_0_0_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_0_1_Om_87_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc0(ix+(2), iy+(28), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  30
       do iy = 1 ,  2
          do ix = 1 ,  258
             Facet_0_1_0_src__1_0_0_dest__1_1_0_Send%Ridge_0_1_0_Om_87_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Facet_0_1_0_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_0_Om_87_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(28))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Facet_0_1_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_1_Om_87_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc0(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_31

  subroutine Formura_internal_32 ()
    double precision  :: a
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz, a)
    do  iz = 1 ,  30 ,  1

       do  iy = 1 ,  30 ,  1

          do  ix = 3 ,  260 ,  1

             a=Rsc3(ix+(-1), iy+(1), iz+(1))

 Rsc1(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( ( ( ( ( 1.0 / 900.0 ) * Rsc2(ix+(-1), iy+(1), iz+(1)) ) * a ) * a ) - ( ( 6.0 / 86400.0 ) * a ) ) + ( ( 6.1e-11 * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc3(ix+(0), iy+(1), iz+(1)) - a ) - ( a - &
Rsc3(ix+(-2), iy+(1), iz+(1)) ) ) + ( ( Rsc3(ix+(-1), iy+(2), iz+(1)) - a ) - ( a - Rsc3(ix+(-1), iy+(0), iz+(1)) ) ) ) + ( ( Rsc3(ix+(-1), iy+(1), iz+(2)) - a ) - ( a - Rsc3(ix+(-1), iy+(1), iz+(0)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  30
       do iy = 1 ,  2
          do ix = 1 ,  258
             Ridge_0_0_0_Om_88_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc1(ix+(2), iy+(28), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  30
          do ix = 1 ,  258
             Ridge_0_0_0_Om_88_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(28))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Ridge_0_0_0_Om_88_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc1(ix+(2), iy+(28), iz+(28))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  30
          do ix = 1 ,  258
             Facet_0_0_1_src__1_0_0_dest__1_0_1_Send%Ridge_0_0_1_Om_88_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Facet_0_0_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_0_1_Om_88_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc1(ix+(2), iy+(28), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  30
       do iy = 1 ,  2
          do ix = 1 ,  258
             Facet_0_1_0_src__1_0_0_dest__1_1_0_Send%Ridge_0_1_0_Om_88_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Facet_0_1_0_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_0_Om_88_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(28))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  258
             Facet_0_1_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_1_Om_88_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc1(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_32

  subroutine Formura_internal_33 ()
    double precision  :: a
    double precision  :: a_0
    integer :: ix
    integer :: iy
    integer :: iz
    !$omp parallel

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  30
       do iy = 1 ,  30
          do ix = 1 ,  2
             Rsc0(ix+(0), iy+(0), iz+(0))=Ridge_0_0_0_Om_87_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  30
       do iy = 1 ,  30
          do ix = 1 ,  2
             Rsc0(ix+(260), iy+(0), iz+(0))=Facet_1_0_0_src__0_0_0_dest__1_0_0_Recv%Ridge_1_0_0_Om_87_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  30
       do iy = 1 ,  30
          do ix = 1 ,  2
             Rsc1(ix+(0), iy+(0), iz+(0))=Ridge_0_0_0_Om_88_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  30
       do iy = 1 ,  30
          do ix = 1 ,  2
             Rsc1(ix+(260), iy+(0), iz+(0))=Facet_1_0_0_src__0_0_0_dest__1_0_0_Recv%Ridge_1_0_0_Om_88_r__0_0_0_r__1_0_0(ix, iy, iz)
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz, a, a_0)
    do  iz = 1 ,  28 ,  1

       do  iy = 1 ,  28 ,  1

          do  ix = 3 ,  262 ,  1

             a=Rsc0(ix+(-1), iy+(1), iz+(1))
             a_0=Rsc1(ix+(-1), iy+(1), iz+(1))

 Rsc3(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( (-( ( ( ( 1.0 / 900.0 ) * a ) * a_0 ) * a_0 )) + ( ( 1.0 / 86400.0 ) * ( 1.0 - a ) ) ) + ( ( ( 0.1 * 2.3e-9 ) * ( 1000.0 * 1000.0 ) ) * ( ( ( ( Rsc0(ix+(0), iy+(1), iz+(1)) - a ) - ( a - &
Rsc0(ix+(-2), iy+(1), iz+(1)) ) ) + ( ( Rsc0(ix+(-1), iy+(2), iz+(1)) - a ) - ( a - Rsc0(ix+(-1), iy+(0), iz+(1)) ) ) ) + ( ( Rsc0(ix+(-1), iy+(1), iz+(2)) - a ) - ( a - Rsc0(ix+(-1), iy+(1), iz+(0)) ) ) ) ) ) ) ))

          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  28
       do iy = 1 ,  2
          do ix = 1 ,  260
             Ridge_0_0_0_Om_89_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc3(ix+(2), iy+(26), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  28
          do ix = 1 ,  260
             Ridge_0_0_0_Om_89_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(26))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Ridge_0_0_0_Om_89_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc3(ix+(2), iy+(26), iz+(26))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  28
          do ix = 1 ,  260
             Facet_0_0_1_src__1_0_0_dest__1_0_1_Send%Ridge_0_0_1_Om_89_r__1_0_0_r__1_0_1(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Facet_0_0_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_0_1_Om_89_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc3(ix+(2), iy+(26), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  28
       do iy = 1 ,  2
          do ix = 1 ,  260
             Facet_0_1_0_src__1_0_0_dest__1_1_0_Send%Ridge_0_1_0_Om_89_r__1_0_0_r__1_1_0(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Facet_0_1_0_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_0_Om_89_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(26))
          end do
       end do
    end do

    !$omp do collapse(2) private(ix, iy, iz)
    do iz = 1 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  260
             Facet_0_1_1_src__1_0_0_dest__1_1_1_Send%Ridge_0_1_1_Om_89_r__1_0_0_r__1_1_1(ix, iy, iz)=Rsc3(ix+(2), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_33

end module finalist_internal_11


c***FROM: http://www.netlib.org/blas/d1mach.f
      DOUBLE PRECISION FUNCTION D1MACH(I)
      INTEGER I
C
C  DOUBLE-PRECISION MACHINE CONSTANTS
C  D1MACH( 1) = B**(EMIN-1), THE SMALLEST POSITIVE MAGNITUDE.
C  D1MACH( 2) = B**EMAX*(1 - B**(-T)), THE LARGEST MAGNITUDE.
C  D1MACH( 3) = B**(-T), THE SMALLEST RELATIVE SPACING.
C  D1MACH( 4) = B**(1-T), THE LARGEST RELATIVE SPACING.
C  D1MACH( 5) = LOG10(B)
C
      INTEGER SMALL(2)
      INTEGER LARGE(2)
      INTEGER RIGHT(2)
      INTEGER DIVER(2)
      INTEGER LOG10(2)
      INTEGER SC, CRAY1(38), J
      COMMON /D9MACH/ CRAY1
      SAVE SMALL, LARGE, RIGHT, DIVER, LOG10, SC
      DOUBLE PRECISION DMACH(5)
      EQUIVALENCE (DMACH(1),SMALL(1))
      EQUIVALENCE (DMACH(2),LARGE(1))
      EQUIVALENCE (DMACH(3),RIGHT(1))
      EQUIVALENCE (DMACH(4),DIVER(1))
      EQUIVALENCE (DMACH(5),LOG10(1))
C  THIS VERSION ADAPTS AUTOMATICALLY TO MOST CURRENT MACHINES.
C  R1MACH CAN HANDLE AUTO-DOUBLE COMPILING, BUT THIS VERSION OF
C  D1MACH DOES NOT, BECAUSE WE DO NOT HAVE QUAD CONSTANTS FOR
C  MANY MACHINES YET.
C  TO COMPILE ON OLDER MACHINES, ADD A C IN COLUMN 1
C  ON THE NEXT LINE
      DATA SC/0/
C  AND REMOVE THE C FROM COLUMN 1 IN ONE OF THE SECTIONS BELOW.
C  CONSTANTS FOR EVEN OLDER MACHINES CAN BE OBTAINED BY
C          mail netlib@research.bell-labs.com
C          send old1mach from blas
C  PLEASE SEND CORRECTIONS TO dmg OR ehg@bell-labs.com.
C
C     MACHINE CONSTANTS FOR THE HONEYWELL DPS 8/70 SERIES.
C      DATA SMALL(1),SMALL(2) / O402400000000, O000000000000 /
C      DATA LARGE(1),LARGE(2) / O376777777777, O777777777777 /
C      DATA RIGHT(1),RIGHT(2) / O604400000000, O000000000000 /
C      DATA DIVER(1),DIVER(2) / O606400000000, O000000000000 /
C      DATA LOG10(1),LOG10(2) / O776464202324, O117571775714 /, SC/987/
C
C     MACHINE CONSTANTS FOR PDP-11 FORTRANS SUPPORTING
C     32-BIT INTEGERS.
C      DATA SMALL(1),SMALL(2) /    8388608,           0 /
C      DATA LARGE(1),LARGE(2) / 2147483647,          -1 /
C      DATA RIGHT(1),RIGHT(2) /  612368384,           0 /
C      DATA DIVER(1),DIVER(2) /  620756992,           0 /
C      DATA LOG10(1),LOG10(2) / 1067065498, -2063872008 /, SC/987/
C
C     MACHINE CONSTANTS FOR THE UNIVAC 1100 SERIES.
C      DATA SMALL(1),SMALL(2) / O000040000000, O000000000000 /
C      DATA LARGE(1),LARGE(2) / O377777777777, O777777777777 /
C      DATA RIGHT(1),RIGHT(2) / O170540000000, O000000000000 /
C      DATA DIVER(1),DIVER(2) / O170640000000, O000000000000 /
C      DATA LOG10(1),LOG10(2) / O177746420232, O411757177572 /, SC/987/
C
C     ON FIRST CALL, IF NO DATA UNCOMMENTED, TEST MACHINE TYPES.
      IF (SC .NE. 987) THEN
         DMACH(1) = 1.D13
         IF (      SMALL(1) .EQ. 1117925532
     *       .AND. SMALL(2) .EQ. -448790528) THEN
*           *** IEEE BIG ENDIAN ***
            SMALL(1) = 1048576
            SMALL(2) = 0
            LARGE(1) = 2146435071
            LARGE(2) = -1
            RIGHT(1) = 1017118720
            RIGHT(2) = 0
            DIVER(1) = 1018167296
            DIVER(2) = 0
            LOG10(1) = 1070810131
            LOG10(2) = 1352628735
         ELSE IF ( SMALL(2) .EQ. 1117925532
     *       .AND. SMALL(1) .EQ. -448790528) THEN
*           *** IEEE LITTLE ENDIAN ***
            SMALL(2) = 1048576
            SMALL(1) = 0
            LARGE(2) = 2146435071
            LARGE(1) = -1
            RIGHT(2) = 1017118720
            RIGHT(1) = 0
            DIVER(2) = 1018167296
            DIVER(1) = 0
            LOG10(2) = 1070810131
            LOG10(1) = 1352628735
         ELSE IF ( SMALL(1) .EQ. -2065213935
     *       .AND. SMALL(2) .EQ. 10752) THEN
*               *** VAX WITH D_FLOATING ***
            SMALL(1) = 128
            SMALL(2) = 0
            LARGE(1) = -32769
            LARGE(2) = -1
            RIGHT(1) = 9344
            RIGHT(2) = 0
            DIVER(1) = 9472
            DIVER(2) = 0
            LOG10(1) = 546979738
            LOG10(2) = -805796613
         ELSE IF ( SMALL(1) .EQ. 1267827943
     *       .AND. SMALL(2) .EQ. 704643072) THEN
*               *** IBM MAINFRAME ***
            SMALL(1) = 1048576
            SMALL(2) = 0
            LARGE(1) = 2147483647
            LARGE(2) = -1
            RIGHT(1) = 856686592
            RIGHT(2) = 0
            DIVER(1) = 873463808
            DIVER(2) = 0
            LOG10(1) = 1091781651
            LOG10(2) = 1352628735
         ELSE IF ( SMALL(1) .EQ. 1120022684
     *       .AND. SMALL(2) .EQ. -448790528) THEN
*           *** CONVEX C-1 ***
            SMALL(1) = 1048576
            SMALL(2) = 0
            LARGE(1) = 2147483647
            LARGE(2) = -1
            RIGHT(1) = 1019215872
            RIGHT(2) = 0
            DIVER(1) = 1020264448
            DIVER(2) = 0
            LOG10(1) = 1072907283
            LOG10(2) = 1352628735
         ELSE IF ( SMALL(1) .EQ. 815547074
     *       .AND. SMALL(2) .EQ. 58688) THEN
*           *** VAX G-FLOATING ***
            SMALL(1) = 16
            SMALL(2) = 0
            LARGE(1) = -32769
            LARGE(2) = -1
            RIGHT(1) = 15552
            RIGHT(2) = 0
            DIVER(1) = 15568
            DIVER(2) = 0
            LOG10(1) = 1142112243
            LOG10(2) = 2046775455
         ELSE
            DMACH(2) = 1.D27 + 1
            DMACH(3) = 1.D27
            LARGE(2) = LARGE(2) - RIGHT(2)
            IF (LARGE(2) .EQ. 64 .AND. SMALL(2) .EQ. 0) THEN
               CRAY1(1) = 67291416
               DO 10 J = 1, 20
                  CRAY1(J+1) = CRAY1(J) + CRAY1(J)
 10               CONTINUE
               CRAY1(22) = CRAY1(21) + 321322
               DO 20 J = 22, 37
                  CRAY1(J+1) = CRAY1(J) + CRAY1(J)
 20               CONTINUE
               IF (CRAY1(38) .EQ. SMALL(1)) THEN
*                  *** CRAY ***
                  CALL I1MCRY(SMALL(1), J, 8285, 8388608, 0)
                  SMALL(2) = 0
                  CALL I1MCRY(LARGE(1), J, 24574, 16777215, 16777215)
                  CALL I1MCRY(LARGE(2), J, 0, 16777215, 16777214)
                  CALL I1MCRY(RIGHT(1), J, 16291, 8388608, 0)
                  RIGHT(2) = 0
                  CALL I1MCRY(DIVER(1), J, 16292, 8388608, 0)
                  DIVER(2) = 0
                  CALL I1MCRY(LOG10(1), J, 16383, 10100890, 8715215)
                  CALL I1MCRY(LOG10(2), J, 0, 16226447, 9001388)
               ELSE
                  WRITE(*,9000)
                  STOP 779
                  END IF
            ELSE
               WRITE(*,9000)
               STOP 779
               END IF
            END IF
         SC = 987
         END IF
*    SANITY CHECK
      IF (DMACH(4) .GE. 1.0D0) STOP 778
      IF (I .LT. 1 .OR. I .GT. 5) THEN
         WRITE(*,*) 'D1MACH(I): I =',I,' is out of bounds.'
         STOP
         END IF
      D1MACH = DMACH(I)
      RETURN
 9000 FORMAT(/' Adjust D1MACH by uncommenting data statements'/
     *' appropriate for your machine.')
* /* Standard C source for D1MACH -- remove the * in column 1 */
*#include <stdio.h>
*#include <float.h>
*#include <math.h>
*double d1mach_(long *i)
*{
*	switch(*i){
*	  case 1: return DBL_MIN;
*	  case 2: return DBL_MAX;
*	  case 3: return DBL_EPSILON/FLT_RADIX;
*	  case 4: return DBL_EPSILON;
*	  case 5: return log10((double)FLT_RADIX);
*	  }
*	fprintf(stderr, "invalid argument: d1mach(%ld)\n", *i);
*	exit(1); return 0; /* some compilers demand return values */
*}
      END
      SUBROUTINE I1MCRY(A, A1, B, C, D)
**** SPECIAL COMPUTATION FOR OLD CRAY MACHINES ****
      INTEGER A, A1, B, C, D
      A1 = 16777216*B + C
      A = 16777216*A1 + D
      END


      SUBROUTINE G1CHAR(NCHAR)
C
C          ------------------------------------------------
C          ROUTINE NO. (1032)   VERSION (A8.7)    14:JUL:97
C          ------------------------------------------------
C
C          THIS DRAWS OUT THE CHARACTER NO. <NCHAR>.
C
C
      INTEGER IXPTS(30),IYPTS(30),IXUNDL(7),IYUNDL(7)
      LOGICAL DRAW,DRAWN
      LOGICAL TYPMOD,KHRDW1,OPEN,CURTRA
C
      COMMON /T1CATT/ IUNDL1,ITAL1
      COMMON /T1CDCC/ SINCHR,COSCHR
      COMMON /T1CDIM/ MAGN1,OBLAT1
      COMMON /T1CFIN/ INDLEN(2,256),WIDCHR(256)
      COMMON /T1CFON/ KFONT1
      COMMON /T1CMOD/ TYPMOD
      COMMON /T1CPOS/ XCHAR,YCHAR
      COMMON /T1DENA/ OPEN
      COMMON /T1DLIM/ DLIMX,DLIMY
      COMMON /T1DRES/ DRESX,DRESY
      COMMON /T1HRDC/ KHRDW1
      COMMON /T1HRDM/ MAPCHR(256)
      COMMON /T1LIND/ PATREM
      COMMON /T1PPAT/ XPLOTA,YPLOTA
      COMMON /T1PPBT/ XPLOTB,YPLOTB
      COMMON /T1SPEC/ XSPEC(4600),YSPEC(4600),ICHPTR(128),NEXTPT,NSTORE
      COMMON /T1TCOB/ XCBX,XCBY,YCBX,YCBY,RPICXB,RPICYB,RDEVXB,RDEVYB
      COMMON /T1TCOC/ XCCX,XCCY,YCCX,YCCY,RPICXC,RPICYC,RDEVXC,RDEVYC,
     &                CURTRA
      COMMON /T1TSPA/ X1TYPW,X2TYPW,Y1TYPW,Y2TYPW
C
      DATA IXUNDL /7,7,6,0,7,0,6/, IYUNDL /0,2,1,1,0,0,0/,
     &     CHRSIZ /1.428571E-4/, DOTSIZ /0.0005/
C
C
C          IF NO DEVICE CHANNEL IS ON, NOTHING CAN BE DONE.
C
C          IF TYPING MODE HAS BEEN SET, THE CHAR. IS NOT
C          DRAWN IF IT LIES OUTSIDE THE TYPEWRITER WINDOW.
C
      IF (.NOT.OPEN)                            RETURN
      IF (TYPMOD) THEN
        IF ((XCHAR-X1TYPW)*(XCHAR-X2TYPW).GT.0.0) RETURN
        IF ((YCHAR-Y1TYPW)*(YCHAR-Y2TYPW).GT.0.0) RETURN
      ENDIF
C
C          IF HARDWARE CHARS. ARE BEING DONE, A CHECK IS MADE
C          TO ENSURE THE POSITION LIES WITHIN THE DEVICE WINDOW.
C
      IF (KHRDW1) THEN
C
        XCHARB= XCBX*(XCHAR-RPICXB)+XCBY*(YCHAR-RPICYB)+RDEVXB
        YCHARB= YCBX*(XCHAR-RPICXB)+YCBY*(YCHAR-RPICYB)+RDEVYB
        IF (.NOT.CURTRA) THEN
          XCHARA= XCHARB
          YCHARA= YCHARB
        ELSE
          XCHARA= XCCX*(XCHARB-RPICXC)+XCCY*(YCHARB-RPICYC)+RDEVXC
          YCHARA= YCCX*(XCHARB-RPICXC)+YCCY*(YCHARB-RPICYC)+RDEVYC
        ENDIF
C
        IF (XCHARA.LT.0.0.OR.XCHARA.GT.DLIMX) RETURN
        IF (YCHARA.LT.0.0.OR.YCHARA.GT.DLIMY) RETURN
C
        IXCHAR= DRESX*XCHARA+0.5
        IYCHAR= DRESY*YCHARA+0.5
        CALL G1FILB(IXCHAR,IYCHAR,-2,MAPCHR(NCHAR+1))
        RETURN
C
      ELSE
C
        XSAVEB= XPLOTB
        YSAVEB= YPLOTB
        XSAVEA= XPLOTA
        YSAVEA= YPLOTA
        REMSAV= PATREM
        IF (KFONT1.LT.11) THEN
C
C          FOR SOFTWARE CHARS., FIRST THE SPECIFICATION IS
C          FETCHED. IF UNDERLINE IS REQUIRED, THE NECESSARY
C          EXTRA SPECIFICATION IS APPENDED USING LOOP-100.
C
          CALL G1SPEC(NCHAR,KFONT1,IXPTS,IYPTS,LENGTH)
          IF (IUNDL1.GT.0) THEN
C
            NADD= IUNDL1*3+1
C
            DO 100 IADD= 1,NADD
              IAPPND= LENGTH+IADD
              IXPTS(IAPPND)= IXUNDL(IADD)
              IYPTS(IAPPND)= IYUNDL(IADD)
  100       CONTINUE
C
            LENGTH= LENGTH+NADD
          ENDIF
C
C          IF THE FINAL SPEC. HAS ZERO LENGTH, NOTHING IS DONE.
C
C          LOOP-200 SCANS THROUGH THE COORDINATE PAIRS,
C          DRAWING THE ELEMENTAL LINES OF THE CHARACTER.
C          SPECIAL FUNCTIONS ARE GIVEN BY <IXV>= 7, THE
C          ACTION BEING DETERMINED BY <IYV> AS FOLLOWS:
C
C          = 0, THE NEXT ELEMENTAL LINE IS NOT VISIBLE,
C          = 1, THE Y-OFFSET IS SET TO -3,
C          = 2, THE Y-OFFSET IS SET TO -2,
C          = 3, THE Y-OFFSET IS SET TO -1,
C          = 4, THE Y-OFFSET IS SET TO  0,
C          = 5, THE Y-OFFSET IS SET TO +1,
C          = 6, THE Y-OFFSET IS SET TO +2,
C          = 7, (UNDEFINED).
C
C          A LOZENGE IS DRAWN WHEN TWO SUCCESSIVE VECTORS ARE
C          NOT VISIBLE, OR IF THE LAST VECTOR IS NOT VISIBLE.
C
          IF (LENGTH.LE.0) RETURN
C
          IYOFF= 0
          DRAW= .FALSE.
          DRAWN= .TRUE.
C
          DO 200 IPT= 1,LENGTH
            IXV= IXPTS(IPT)
            IYV= IYPTS(IPT)
            IF (IXV.LE.6) GO TO 1
            IF (IYV.NE.0) THEN
              IYOFF= IYV-4
              GO TO 200
            ENDIF
C
            DRAW= .FALSE.
            IF (.NOT.DRAWN) GO TO 2
            GO TO 200
C
    1       IYV= IYV+IYOFF
            X1= (IXV-3)*MAGN1*CHRSIZ
            Y1= (IYV-2)*MAGN1*CHRSIZ
            X1= (X1+Y1*ITAL1*0.25)*OBLAT1
            XNOW= XCHAR+X1*COSCHR-Y1*SINCHR
            YNOW= YCHAR+X1*SINCHR+Y1*COSCHR
            IF (.NOT.DRAW.AND.IPT.EQ.LENGTH) GO TO 2
C
            CALL G1LINE(XNOW,YNOW,DRAW)
            DRAWN= DRAW
            DRAW= .TRUE.
            GO TO 200
C
    2       CALL G1LINE(XNOW,       YNOW+DOTSIZ,.FALSE.)
            CALL G1LINE(XNOW-DOTSIZ,YNOW,       .TRUE.)
            CALL G1LINE(XNOW,       YNOW-DOTSIZ,.TRUE.)
            CALL G1LINE(XNOW+DOTSIZ,YNOW,       .TRUE.)
            CALL G1LINE(XNOW,       YNOW+DOTSIZ,.TRUE.)
  200     CONTINUE
C
        ELSE
          NCHAR1= NCHAR+1
          IF (KFONT1.LT.101) THEN
C
C          FONTS 11-24 ARE HIGH-QUALITY FONTS.
C
            IF (NCHAR1.GT.128) THEN
              IPOS= 1
              IF (NCHAR1.NE.NSTORE) CALL G1CHIO(2,NCHAR1)
C
              NSTORE= NCHAR1
            ELSE
              IF (ICHPTR(NCHAR1).EQ.0) CALL G1CHIO(2,NCHAR1)
C
              IPOS= ICHPTR(NCHAR1)
            ENDIF
C
            DRAW= .FALSE.
            IP= IPOS-1
C
            DO 400 IPT= 1,INDLEN(2,NCHAR1)
              IP= IP+1
              XC= XSPEC(IP)
              IF(XC.GE.-100.0) THEN
                X1= (XC-3.0)*MAGN1*CHRSIZ*OBLAT1
                Y1= (YSPEC(IP)-2.0)*MAGN1*CHRSIZ
                XNOW= XCHAR+X1*COSCHR-Y1*SINCHR
                YNOW= YCHAR+X1*SINCHR+Y1*COSCHR
                CALL G1LINE(XNOW,YNOW,DRAW)
                DRAW= .TRUE.
              ELSE
                DRAW= .FALSE.
              ENDIF
C
  400       CONTINUE
C
            IF (IUNDL1.GT.0) THEN
C
C          UNDERLINE HIGH QUALITY CHARACTERS.
C
              XPOS= 3.0*WIDCHR(NCHAR1)
              YPOS= -3.0
              DRAW= .FALSE.
C
              DO 500 LIN= 1,IUNDL1*2
                XNOW= XCHAR+(XPOS*COSCHR*OBLAT1-YPOS*SINCHR)*
     &                MAGN1*CHRSIZ
                YNOW= YCHAR+(XPOS*SINCHR*OBLAT1+YPOS*COSCHR)*
     &                MAGN1*CHRSIZ
                CALL G1LINE(XNOW,YNOW,DRAW)
                IF (MOD(LIN,2).EQ.1) XPOS= -XPOS
                IF (LIN.EQ.2) YPOS= YPOS-1.0
C
                DRAW= .NOT.DRAW
  500         CONTINUE
C
            ENDIF
C
          ELSE
C
C          FONTS 101 ONWARDS ARE POSTSCRIPT FONTS.
C
C          A CHECK IS MADE TO ENSURE THE POSITION LIES WITHIN
C          THE DEVICE WINDOW.
C
            X1= -3.0*MAGN1*CHRSIZ*OBLAT1*WIDCHR(NCHAR1)
            Y1= -2.0*MAGN1*CHRSIZ
            XNOW= XCHAR+X1*COSCHR-Y1*SINCHR
            YNOW= YCHAR+X1*SINCHR+Y1*COSCHR
            XCHARB= XCBX*(XNOW-RPICXB)+XCBY*(YNOW-RPICYB)+RDEVXB
            YCHARB= YCBX*(XNOW-RPICXB)+YCBY*(YNOW-RPICYB)+RDEVYB
            IF (.NOT.CURTRA) THEN
              XCHARA= XCHARB
              YCHARA= YCHARB
            ELSE
              XCHARA= XCCX*(XCHARB-RPICXC)+XCCY*(YCHARB-RPICYC)+RDEVXC
              YCHARA= YCCX*(XCHARB-RPICXC)+YCCY*(YCHARB-RPICYC)+RDEVYC
            ENDIF
C
            IF (XCHARA.LT.0.0.OR.XCHARA.GT.DLIMX) RETURN
            IF (YCHARA.LT.0.0.OR.YCHARA.GT.DLIMY) RETURN
C
            IXCHAR= DRESX*XCHARA+0.5
            IYCHAR= DRESY*YCHARA+0.5
            CALL G1FILB(IXCHAR,IYCHAR,-2,NCHAR)
          ENDIF
        ENDIF
C
        XPLOTB= XSAVEB
        YPLOTB= YSAVEB
        XPLOTA= XSAVEA
        YPLOTA= YSAVEA
        PATREM= REMSAV
      ENDIF
C
      RETURN
      END