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

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

       do  iy = 1 ,  38 ,  1

          do  ix = 1 ,  262 ,  1

             a=U(ix+(251), iy+(27), iz+(27))
             a_0=V(ix+(251), iy+(27), iz+(27))

 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 ) ) * ( ( ( ( U(ix+(252), iy+(27), iz+(27)) - a ) - ( a - &
U(ix+(250), iy+(27), iz+(27)) ) ) + ( ( U(ix+(251), iy+(28), iz+(27)) - a ) - ( a - U(ix+(251), iy+(26), iz+(27)) ) ) ) + ( ( U(ix+(251), iy+(27), iz+(28)) - a ) - ( a - U(ix+(251), iy+(27), 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 ,  2
             Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_m1_Om_79_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(260), iy+(36), iz+(36))
          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
             Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_0_Om_79_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(260), iy+(36), iz+(0))
          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
             Facet_m1_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_m1_0_Om_79_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc0(ix+(260), iy+(36), 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 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_m1_Om_79_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(260), iy+(0), iz+(36))
          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
             Facet_m1_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_m1_Om_79_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc0(ix+(260), iy+(0), iz+(36))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_0_Om_79_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(260), 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 ,  38
          do ix = 1 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_0_Om_79_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc0(ix+(260), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_0_0_Om_79_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc0(ix+(260), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_1_1_Send%Ridge_m1_0_0_Om_79_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc0(ix+(260), 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 ,  2
             Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_m1_Om_79_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(0), iy+(36), iz+(36))
          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
             Facet_0_m1_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_m1_Om_79_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc0(ix+(0), iy+(36), iz+(36))
          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
             Facet_0_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_0_Om_79_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(0), iy+(36), 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 ,  262
             Facet_0_m1_0_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_0_Om_79_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc0(ix+(0), iy+(36), iz+(0))
          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
             Facet_0_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_0_m1_0_Om_79_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc0(ix+(0), iy+(36), iz+(0))
          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
             Facet_0_m1_0_src__1_1_1_dest__1_0_1_Send%Ridge_0_m1_0_Om_79_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc0(ix+(0), iy+(36), 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 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_0_m1_Om_79_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(36))
          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
             Facet_0_0_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_0_m1_Om_79_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(36))
          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
             Facet_0_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_0_0_m1_Om_79_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(36))
          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
             Facet_0_0_m1_src__1_1_1_dest__1_1_0_Send%Ridge_0_0_m1_Om_79_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(36))
          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
             Ridge_0_0_0_Om_79_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  262
             Ridge_0_0_0_Om_79_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  38
          do ix = 1 ,  2
             Ridge_0_0_0_Om_79_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  38
          do ix = 1 ,  262
             Ridge_0_0_0_Om_79_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_79_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_79_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_79_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_136

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

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

       do  iy = 1 ,  38 ,  1

          do  ix = 1 ,  262 ,  1

             a=V(ix+(251), iy+(27), iz+(27))

 Rsc1(ix+(0), iy+(0), iz+(0))=(+( a + ( 200.0 * ( ( ( ( ( ( 1.0 / 900.0 ) * U(ix+(251), iy+(27), iz+(27)) ) * a ) * a ) - ( ( 6.0 / 86400.0 ) * a ) ) + ( ( 6.1e-11 * ( 1000.0 * 1000.0 ) ) * ( ( ( ( V(ix+(252), iy+(27), iz+(27)) - a ) - ( a - &
V(ix+(250), iy+(27), iz+(27)) ) ) + ( ( V(ix+(251), iy+(28), iz+(27)) - a ) - ( a - V(ix+(251), iy+(26), iz+(27)) ) ) ) + ( ( V(ix+(251), iy+(27), iz+(28)) - a ) - ( a - V(ix+(251), iy+(27), 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 ,  2
             Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_m1_Om_80_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(260), iy+(36), iz+(36))
          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
             Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_0_Om_80_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(260), iy+(36), iz+(0))
          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
             Facet_m1_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_m1_0_Om_80_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc1(ix+(260), iy+(36), 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 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_m1_Om_80_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(260), iy+(0), iz+(36))
          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
             Facet_m1_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_m1_Om_80_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc1(ix+(260), iy+(0), iz+(36))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_0_Om_80_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(260), 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 ,  38
          do ix = 1 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_0_Om_80_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc1(ix+(260), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_0_0_Om_80_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc1(ix+(260), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_1_1_Send%Ridge_m1_0_0_Om_80_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc1(ix+(260), 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 ,  2
             Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_m1_Om_80_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(0), iy+(36), iz+(36))
          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
             Facet_0_m1_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_m1_Om_80_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc1(ix+(0), iy+(36), iz+(36))
          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
             Facet_0_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_0_Om_80_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(0), iy+(36), 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 ,  262
             Facet_0_m1_0_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_0_Om_80_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc1(ix+(0), iy+(36), iz+(0))
          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
             Facet_0_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_0_m1_0_Om_80_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc1(ix+(0), iy+(36), iz+(0))
          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
             Facet_0_m1_0_src__1_1_1_dest__1_0_1_Send%Ridge_0_m1_0_Om_80_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc1(ix+(0), iy+(36), 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 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_0_m1_Om_80_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(36))
          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
             Facet_0_0_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_0_m1_Om_80_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(36))
          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
             Facet_0_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_0_0_m1_Om_80_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(36))
          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
             Facet_0_0_m1_src__1_1_1_dest__1_1_0_Send%Ridge_0_0_m1_Om_80_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(36))
          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
             Ridge_0_0_0_Om_80_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  262
             Ridge_0_0_0_Om_80_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  38
          do ix = 1 ,  2
             Ridge_0_0_0_Om_80_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  38
          do ix = 1 ,  262
             Ridge_0_0_0_Om_80_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_80_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_80_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_80_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_137

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

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

       do  iy = 1 ,  36 ,  1

          do  ix = 1 ,  260 ,  1

             a=Rsc0(ix+(1), iy+(1), iz+(1))
             a_0=Rsc1(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 ) ) * ( ( ( ( Rsc0(ix+(2), iy+(1), iz+(1)) - a ) - ( a - &
Rsc0(ix+(0), 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 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_m1_Om_81_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(258), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_0_Om_81_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(258), 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 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_m1_0_Om_81_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc2(ix+(258), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_m1_Om_81_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(258), 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 ,  36
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_m1_Om_81_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc2(ix+(258), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_0_Om_81_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(258), 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 ,  36
          do ix = 1 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_0_Om_81_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc2(ix+(258), iy+(0), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_0_0_Om_81_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc2(ix+(258), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_1_1_Send%Ridge_m1_0_0_Om_81_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc2(ix+(258), 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 ,  2
             Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_m1_Om_81_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  260
             Facet_0_m1_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_m1_Om_81_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_0_Om_81_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  260
             Facet_0_m1_0_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_0_Om_81_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_0_m1_0_Om_81_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  260
             Facet_0_m1_0_src__1_1_1_dest__1_0_1_Send%Ridge_0_m1_0_Om_81_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_0_m1_Om_81_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  260
             Facet_0_0_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_0_m1_Om_81_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  36
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_0_0_m1_Om_81_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  36
          do ix = 1 ,  260
             Facet_0_0_m1_src__1_1_1_dest__1_1_0_Send%Ridge_0_0_m1_Om_81_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Ridge_0_0_0_Om_81_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(0), 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
             Ridge_0_0_0_Om_81_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  36
          do ix = 1 ,  2
             Ridge_0_0_0_Om_81_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  36
          do ix = 1 ,  260
             Ridge_0_0_0_Om_81_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc2(ix+(0), iy+(0), 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 ,  2
             Ridge_0_0_0_Om_81_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc2(ix+(0), iy+(0), 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 ,  260
             Ridge_0_0_0_Om_81_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc2(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_81_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc2(ix+(0), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_138

  subroutine Formura_internal_139 ()
    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 = 1 ,  260 ,  1

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

 Rsc3(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+(2), iy+(1), iz+(1)) - a ) - ( a - &
Rsc1(ix+(0), 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 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_m1_Om_82_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(258), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_0_Om_82_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(258), 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 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_m1_0_Om_82_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc3(ix+(258), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_m1_Om_82_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(258), 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 ,  36
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_m1_Om_82_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc3(ix+(258), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_0_Om_82_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(258), 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 ,  36
          do ix = 1 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_0_Om_82_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc3(ix+(258), iy+(0), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_0_0_Om_82_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc3(ix+(258), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_1_1_Send%Ridge_m1_0_0_Om_82_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc3(ix+(258), 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 ,  2
             Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_m1_Om_82_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
          do ix = 1 ,  260
             Facet_0_m1_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_m1_Om_82_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_0_Om_82_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
          do ix = 1 ,  260
             Facet_0_m1_0_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_0_Om_82_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_0_m1_0_Om_82_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc3(ix+(0), 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 ,  260
             Facet_0_m1_0_src__1_1_1_dest__1_0_1_Send%Ridge_0_m1_0_Om_82_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_0_m1_Om_82_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  260
             Facet_0_0_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_0_m1_Om_82_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  36
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_0_0_m1_Om_82_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  36
          do ix = 1 ,  260
             Facet_0_0_m1_src__1_1_1_dest__1_1_0_Send%Ridge_0_0_m1_Om_82_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
             Ridge_0_0_0_Om_82_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(0), 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
             Ridge_0_0_0_Om_82_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  36
          do ix = 1 ,  2
             Ridge_0_0_0_Om_82_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  36
          do ix = 1 ,  260
             Ridge_0_0_0_Om_82_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc3(ix+(0), iy+(0), 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 ,  2
             Ridge_0_0_0_Om_82_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc3(ix+(0), iy+(0), 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 ,  260
             Ridge_0_0_0_Om_82_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc3(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_82_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc3(ix+(0), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_139

  subroutine Formura_internal_14 ()
    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 = 1 ,  256 ,  1

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

 Rsc2(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+(2), iy+(1), iz+(1)) - a ) - ( a - &
Rsc0(ix+(0), 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 ,  32
          do ix = 1 ,  2
             Ridge_0_0_0_Om_86_r__0_0_0_r__1_0_0(ix, iy, iz)=Rsc2(ix+(254), iy+(0), 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__0_0_0_r__0_1_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Ridge_0_0_0_Om_86_r__0_0_0_r__1_1_0(ix, iy, iz)=Rsc2(ix+(254), 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__0_0_0_r__0_0_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  32
          do ix = 1 ,  2
             Ridge_0_0_0_Om_86_r__0_0_0_r__1_0_1(ix, iy, iz)=Rsc2(ix+(254), 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__0_0_0_r__0_1_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Ridge_0_0_0_Om_86_r__0_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(254), 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__0_0_0_dest__0_0_1_Send%Ridge_0_0_1_Om_86_r__0_0_0_r__0_0_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  32
          do ix = 1 ,  2
             Facet_0_0_1_src__0_0_0_dest__1_0_1_Send%Ridge_0_0_1_Om_86_r__0_0_0_r__1_0_1(ix, iy, iz)=Rsc2(ix+(254), 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__0_0_0_dest__0_1_1_Send%Ridge_0_0_1_Om_86_r__0_0_0_r__0_1_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_0_1_src__0_0_0_dest__1_1_1_Send%Ridge_0_0_1_Om_86_r__0_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(254), 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__0_0_0_dest__0_1_0_Send%Ridge_0_1_0_Om_86_r__0_0_0_r__0_1_0(ix, iy, iz)=Rsc2(ix+(0), iy+(0), 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 ,  2
             Facet_0_1_0_src__0_0_0_dest__1_1_0_Send%Ridge_0_1_0_Om_86_r__0_0_0_r__1_1_0(ix, iy, iz)=Rsc2(ix+(254), 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__0_0_0_dest__0_1_1_Send%Ridge_0_1_0_Om_86_r__0_0_0_r__0_1_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_0_1_0_src__0_0_0_dest__1_1_1_Send%Ridge_0_1_0_Om_86_r__0_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(254), 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__0_0_0_dest__0_1_1_Send%Ridge_0_1_1_Om_86_r__0_0_0_r__0_1_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_0_1_1_src__0_0_0_dest__1_1_1_Send%Ridge_0_1_1_Om_86_r__0_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(254), iy+(0), iz+(0))
          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
             Facet_1_0_0_src__0_0_0_dest__1_0_0_Send%Ridge_1_0_0_Om_86_r__0_0_0_r__1_0_0(ix, iy, iz)=Rsc2(ix+(0), iy+(0), 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 ,  2
             Facet_1_0_0_src__0_0_0_dest__1_1_0_Send%Ridge_1_0_0_Om_86_r__0_0_0_r__1_1_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_1_0_0_src__0_0_0_dest__1_0_1_Send%Ridge_1_0_0_Om_86_r__0_0_0_r__1_0_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_1_0_0_src__0_0_0_dest__1_1_1_Send%Ridge_1_0_0_Om_86_r__0_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_1_0_1_src__0_0_0_dest__1_0_1_Send%Ridge_1_0_1_Om_86_r__0_0_0_r__1_0_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_1_0_1_src__0_0_0_dest__1_1_1_Send%Ridge_1_0_1_Om_86_r__0_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_1_1_0_src__0_0_0_dest__1_1_0_Send%Ridge_1_1_0_Om_86_r__0_0_0_r__1_1_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_1_1_0_src__0_0_0_dest__1_1_1_Send%Ridge_1_1_0_Om_86_r__0_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_1_1_1_src__0_0_0_dest__1_1_1_Send%Ridge_1_1_1_Om_86_r__0_0_0_r__1_1_1(ix, iy, iz)=Rsc2(ix+(0), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_14

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

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

       do  iy = 1 ,  34 ,  1

          do  ix = 1 ,  258 ,  1

             a=Rsc2(ix+(1), iy+(1), iz+(1))
             a_0=Rsc3(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 ) ) * ( ( ( ( Rsc2(ix+(2), iy+(1), iz+(1)) - a ) - ( a - &
Rsc2(ix+(0), 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 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_m1_Om_83_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(256), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_0_Om_83_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(256), 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 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_m1_0_Om_83_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc1(ix+(256), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_m1_Om_83_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(256), 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 ,  34
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_m1_Om_83_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc1(ix+(256), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_0_Om_83_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(256), 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 ,  34
          do ix = 1 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_0_Om_83_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc1(ix+(256), iy+(0), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_0_0_Om_83_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc1(ix+(256), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_1_1_Send%Ridge_m1_0_0_Om_83_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc1(ix+(256), 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 ,  2
             Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_m1_Om_83_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  2
          do ix = 1 ,  258
             Facet_0_m1_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_m1_Om_83_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_0_Om_83_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  2
          do ix = 1 ,  258
             Facet_0_m1_0_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_0_Om_83_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_0_m1_0_Om_83_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc1(ix+(0), 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 ,  258
             Facet_0_m1_0_src__1_1_1_dest__1_0_1_Send%Ridge_0_m1_0_Om_83_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc1(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_0_m1_Om_83_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  258
             Facet_0_0_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_0_m1_Om_83_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  34
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_0_0_m1_Om_83_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  34
          do ix = 1 ,  258
             Facet_0_0_m1_src__1_1_1_dest__1_1_0_Send%Ridge_0_0_m1_Om_83_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  2
             Ridge_0_0_0_Om_83_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc1(ix+(0), 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
             Ridge_0_0_0_Om_83_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  34
          do ix = 1 ,  2
             Ridge_0_0_0_Om_83_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc1(ix+(0), 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 ,  34
          do ix = 1 ,  258
             Ridge_0_0_0_Om_83_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc1(ix+(0), iy+(0), 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 ,  2
             Ridge_0_0_0_Om_83_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc1(ix+(0), iy+(0), 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 ,  258
             Ridge_0_0_0_Om_83_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_83_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc1(ix+(0), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_140

  subroutine Formura_internal_141 ()
    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 = 1 ,  258 ,  1

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

 Rsc0(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+(2), iy+(1), iz+(1)) - a ) - ( a - &
Rsc3(ix+(0), 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 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_m1_Om_84_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(256), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_0_Om_84_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(256), 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 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_m1_0_Om_84_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc0(ix+(256), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_m1_Om_84_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(256), 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 ,  34
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_m1_Om_84_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc0(ix+(256), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_0_Om_84_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(256), 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 ,  34
          do ix = 1 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_0_Om_84_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc0(ix+(256), iy+(0), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_0_0_Om_84_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc0(ix+(256), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_1_1_Send%Ridge_m1_0_0_Om_84_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc0(ix+(256), 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 ,  2
             Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_m1_Om_84_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  2
          do ix = 1 ,  258
             Facet_0_m1_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_m1_Om_84_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_0_Om_84_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  2
          do ix = 1 ,  258
             Facet_0_m1_0_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_0_Om_84_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_0_m1_0_Om_84_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc0(ix+(0), 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 ,  258
             Facet_0_m1_0_src__1_1_1_dest__1_0_1_Send%Ridge_0_m1_0_Om_84_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc0(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_0_m1_Om_84_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  258
             Facet_0_0_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_0_m1_Om_84_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  34
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_0_0_m1_Om_84_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  34
          do ix = 1 ,  258
             Facet_0_0_m1_src__1_1_1_dest__1_1_0_Send%Ridge_0_0_m1_Om_84_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  2
             Ridge_0_0_0_Om_84_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc0(ix+(0), 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
             Ridge_0_0_0_Om_84_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  34
          do ix = 1 ,  2
             Ridge_0_0_0_Om_84_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc0(ix+(0), 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 ,  34
          do ix = 1 ,  258
             Ridge_0_0_0_Om_84_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc0(ix+(0), iy+(0), 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 ,  2
             Ridge_0_0_0_Om_84_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc0(ix+(0), iy+(0), 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 ,  258
             Ridge_0_0_0_Om_84_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_84_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc0(ix+(0), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_141

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

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

       do  iy = 1 ,  32 ,  1

          do  ix = 1 ,  256 ,  1

             a=Rsc1(ix+(1), iy+(1), iz+(1))
             a_0=Rsc0(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 ) ) * ( ( ( ( Rsc1(ix+(2), iy+(1), iz+(1)) - a ) - ( a - &
Rsc1(ix+(0), 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 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_m1_Om_85_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(254), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_0_Om_85_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(254), 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 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_m1_0_Om_85_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc3(ix+(254), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_m1_Om_85_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(254), 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 ,  32
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_m1_Om_85_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc3(ix+(254), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_0_Om_85_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(254), 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 ,  32
          do ix = 1 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_0_Om_85_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc3(ix+(254), iy+(0), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_0_0_Om_85_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc3(ix+(254), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_1_1_Send%Ridge_m1_0_0_Om_85_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc3(ix+(254), 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 ,  2
             Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_m1_Om_85_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
          do ix = 1 ,  256
             Facet_0_m1_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_m1_Om_85_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_0_Om_85_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
          do ix = 1 ,  256
             Facet_0_m1_0_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_0_Om_85_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_0_m1_0_Om_85_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc3(ix+(0), 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_m1_0_src__1_1_1_dest__1_0_1_Send%Ridge_0_m1_0_Om_85_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_0_m1_Om_85_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(0), 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_0_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_0_m1_Om_85_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  32
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_0_0_m1_Om_85_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  32
          do ix = 1 ,  256
             Facet_0_0_m1_src__1_1_1_dest__1_1_0_Send%Ridge_0_0_m1_Om_85_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  2
             Ridge_0_0_0_Om_85_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc3(ix+(0), 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
             Ridge_0_0_0_Om_85_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  32
          do ix = 1 ,  2
             Ridge_0_0_0_Om_85_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc3(ix+(0), 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 ,  32
          do ix = 1 ,  256
             Ridge_0_0_0_Om_85_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc3(ix+(0), iy+(0), 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 ,  2
             Ridge_0_0_0_Om_85_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc3(ix+(0), iy+(0), 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_1_1_r__1_0_1(ix, iy, iz)=Rsc3(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_85_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc3(ix+(0), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_142

  subroutine Formura_internal_143 ()
    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 = 1 ,  256 ,  1

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

 Rsc2(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+(2), iy+(1), iz+(1)) - a ) - ( a - &
Rsc0(ix+(0), 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 ,  2
       do iy = 1 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_m1_Om_86_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(254), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_m1_0_Om_86_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(254), 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 ,  2
             Facet_m1_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_m1_0_Om_86_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc2(ix+(254), 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 ,  2
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_m1_Om_86_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(254), 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 ,  32
          do ix = 1 ,  2
             Facet_m1_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_m1_Om_86_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc2(ix+(254), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_0_Send%Ridge_m1_0_0_Om_86_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(254), 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 ,  32
          do ix = 1 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_1_0_Send%Ridge_m1_0_0_Om_86_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc2(ix+(254), iy+(0), 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 ,  2
             Facet_m1_0_0_src__1_1_1_dest__0_0_1_Send%Ridge_m1_0_0_Om_86_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc2(ix+(254), iy+(0), iz+(0))
          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
             Facet_m1_0_0_src__1_1_1_dest__0_1_1_Send%Ridge_m1_0_0_Om_86_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc2(ix+(254), 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 ,  2
             Facet_0_m1_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_m1_Om_86_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  256
             Facet_0_m1_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_m1_Om_86_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_0_Send%Ridge_0_m1_0_Om_86_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  256
             Facet_0_m1_0_src__1_1_1_dest__1_0_0_Send%Ridge_0_m1_0_Om_86_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Facet_0_m1_0_src__1_1_1_dest__0_0_1_Send%Ridge_0_m1_0_Om_86_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc2(ix+(0), 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_m1_0_src__1_1_1_dest__1_0_1_Send%Ridge_0_m1_0_Om_86_r__1_1_1_r__1_0_1(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_0_0_Send%Ridge_0_0_m1_Om_86_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(0), 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_0_m1_src__1_1_1_dest__1_0_0_Send%Ridge_0_0_m1_Om_86_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  32
          do ix = 1 ,  2
             Facet_0_0_m1_src__1_1_1_dest__0_1_0_Send%Ridge_0_0_m1_Om_86_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  32
          do ix = 1 ,  256
             Facet_0_0_m1_src__1_1_1_dest__1_1_0_Send%Ridge_0_0_m1_Om_86_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  2
             Ridge_0_0_0_Om_86_r__1_1_1_r__0_0_0(ix, iy, iz)=Rsc2(ix+(0), 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
             Ridge_0_0_0_Om_86_r__1_1_1_r__1_0_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  32
          do ix = 1 ,  2
             Ridge_0_0_0_Om_86_r__1_1_1_r__0_1_0(ix, iy, iz)=Rsc2(ix+(0), 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 ,  32
          do ix = 1 ,  256
             Ridge_0_0_0_Om_86_r__1_1_1_r__1_1_0(ix, iy, iz)=Rsc2(ix+(0), iy+(0), 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 ,  2
             Ridge_0_0_0_Om_86_r__1_1_1_r__0_0_1(ix, iy, iz)=Rsc2(ix+(0), iy+(0), 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_1_1_r__1_0_1(ix, iy, iz)=Rsc2(ix+(0), iy+(0), iz+(0))
          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
             Ridge_0_0_0_Om_86_r__1_1_1_r__0_1_1(ix, iy, iz)=Rsc2(ix+(0), iy+(0), iz+(0))
          end do
       end do
    end do

    !$omp end parallel

  end subroutine Formura_internal_143

end module finalist_internal_3

max_issues_repo_name: formura/formura


text:
      dimension p1(2),t(6,4),b1(2),b2(2),al1(2),al2(2),g1(2),g2(2)
      save
      if(nlin.eq.0) then
        do 20 l=1,2
          ll=2*l
          b2(l)=t(6-ll,ll-1)*t(6-ll,ll-1)+t(7-ll,ll-1)*t(7-ll,ll-1)
          write(*,*) b2(l)
   20   continue
      endif
      end

! { dg-final { scan-tree-dump-times "number of SCoPs: 1" 1 "graphite" { xfail *-*-* } } }
! { dg-final { cleanup-tree-dump "graphite" } }

ext: f

max_stars_repo_name: maurizioabba/rose

max_issues_repo_name: sujankh/rose-matlab

max_forks_repo_name: sujankh/rose-matlab


text:
C-------------------------------------------------------------C
C  An Extended Split-Radix DIF FFT                            C
C  Complex input and output in data arrays X and Y            C
C  Length is N = 2**M                                         C
C                                                             C
C  D. Takahashi, University of Tsukuba, Sep. 2002             C
C                                                             C
C  Reference: D. Takahashi, "An extended split-radix FFT      C
C             algorithm," IEEE Signal Processing Lett.,       C
C             vol. 8, pp. 145-147, May 2001.                  C
C-------------------------------------------------------------C
      SUBROUTINE ESRFFT(X,Y,IBETA,N,M,IV)
      IMPLICIT REAL*8 (A-H,O-Z)
      REAL*8 X(N),Y(N)
      INTEGER*4 IBETA(N/2)
      PARAMETER (C21=0.70710678118654752D0,
     +           TWOPI=6.28318530717958647D0)
C-------Beta vector table-------
      IBETA(1)=1
      IBETA(2)=0
      IBETA(3)=0
      IBETA(4)=0
      ID=1
      DO 30 K=3,M-1
        DO 20 J=1,4
          IS=(J+3)*ID+1
          DO 10 I=IS,IS+ID-1
            IBETA(I)=IBETA(I-IS+1)
   10     CONTINUE
   20   CONTINUE
        ID=2*ID
   30 CONTINUE
C
      IF (IV .EQ. -1) THEN
        DO 40 I=1,N
          Y(I)=-Y(I)
   40   CONTINUE
      END IF
C-------L shaped butterflies-------
      L=1
      N2=N
      DO 70 K=1,M-2
        N8=N2/8
        E=-TWOPI/DBLE(N2)
        A=0.0D0
        DO 60 J=1,N8
          A3=3.0D0*A
          A5=5.0D0*A
          A7=7.0D0*A
          CC1=DCOS(A)
          SS1=DSIN(A)
          CC3=DCOS(A3)
          SS3=DSIN(A3)
          CC5=DCOS(A5)
          SS5=DSIN(A5)
          CC7=DCOS(A7)
          SS7=DSIN(A7)
          A=DBLE(J)*E
          DO 50 I=1,L
            IF (IBETA(I) .EQ. 1) THEN
              I0=(I-1)*N2+J
              I1=I0+N8
              I2=I1+N8
              I3=I2+N8
              I4=I3+N8
              I5=I4+N8
              I6=I5+N8
              I7=I6+N8
              X0=X(I0)-X(I4)
              Y0=Y(I0)-Y(I4)
              X1=X(I1)-X(I5)
              Y1=Y(I1)-Y(I5)
              X2=Y(I2)-Y(I6)
              Y2=X(I6)-X(I2)
              X3=X(I3)-X(I7)
              Y3=Y(I3)-Y(I7)
              X(I0)=X(I0)+X(I4)
              Y(I0)=Y(I0)+Y(I4)
              X(I1)=X(I1)+X(I5)
              Y(I1)=Y(I1)+Y(I5)
              X(I2)=X(I2)+X(I6)
              Y(I2)=Y(I2)+Y(I6)
              X(I3)=X(I3)+X(I7)
              Y(I3)=Y(I3)+Y(I7)
              U0=X0+C21*(X1-X3)
              V0=Y0+C21*(Y1-Y3)
              U1=X0-C21*(X1-X3)
              V1=Y0-C21*(Y1-Y3)
              U2=X2+C21*(Y1+Y3)
              V2=Y2-C21*(X1+X3)
              U3=X2-C21*(Y1+Y3)
              V3=Y2+C21*(X1+X3)
              X(I4)=CC1*(U0+U2)-SS1*(V0+V2)
              Y(I4)=CC1*(V0+V2)+SS1*(U0+U2)
              X(I5)=CC5*(U1+U3)-SS5*(V1+V3)
              Y(I5)=CC5*(V1+V3)+SS5*(U1+U3)
              X(I6)=CC3*(U1-U3)-SS3*(V1-V3)
              Y(I6)=CC3*(V1-V3)+SS3*(U1-U3)
              X(I7)=CC7*(U0-U2)-SS7*(V0-V2)
              Y(I7)=CC7*(V0-V2)+SS7*(U0-U2)
            END IF
   50     CONTINUE
   60   CONTINUE
        L=2*L
        N2=N2/2
   70 CONTINUE
C-------Length four butterflies-------
      DO 80 I=1,N/4
        IF (IBETA(I) .EQ. 1) THEN
          I0=4*I-3
          I1=I0+1
          I2=I1+1
          I3=I2+1
          X0=X(I0)-X(I2)
          Y0=Y(I0)-Y(I2)
          X1=Y(I1)-Y(I3)
          Y1=X(I3)-X(I1)
          X(I0)=X(I0)+X(I2)
          Y(I0)=Y(I0)+Y(I2)
          X(I1)=X(I1)+X(I3)
          Y(I1)=Y(I1)+Y(I3)
          X(I2)=X0+X1
          Y(I2)=Y0+Y1
          X(I3)=X0-X1
          Y(I3)=Y0-Y1
        END IF
   80 CONTINUE
C-------Length two butterflies-------
      DO 90 I=1,N/2
        IF (IBETA(I) .EQ. 1) THEN
          I0=2*I-1
          I1=I0+1
          X0=X(I0)-X(I1)
          Y0=Y(I0)-Y(I1)
          X(I0)=X(I0)+X(I1)
          Y(I0)=Y(I0)+Y(I1)
          X(I1)=X0
          Y(I1)=Y0
        END IF
   90 CONTINUE
C-------Digit reverse counter-------
      J=1
      DO 110 I=1,N-1
        IF (I .LT. J) THEN
          X0=X(I)
          Y0=Y(I)
          X(I)=X(J)
          Y(I)=Y(J)
          X(J)=X0
          Y(J)=Y0
        END IF
        K=N/2
  100   IF (K .LT. J) THEN
          J=J-K
          K=K/2
          GO TO 100
        END IF
        J=J+K
  110 CONTINUE
C
      IF (IV .EQ. -1) THEN
        DO 120 I=1,N
          X(I)=X(I)/DBLE(N)
          Y(I)=-Y(I)/DBLE(N)
  120   CONTINUE
      END IF
      RETURN
      END

max_stars_repo_name: mreineck/benchfft

max_issues_repo_name: mreineck/benchfft

max_forks_repo_name: mreineck/benchfft


text:
      SUBROUTINE CALERF(ARG,RESULT,JINT)
C------------------------------------------------------------------
C
C This packet evaluates  erf(x),  erfc(x),  and  exp(x*x)*erfc(x)
C   for a real argument  x.  It contains three FUNCTION type
C   subprograms: ERF, ERFC, and ERFCX (or DERF, DERFC, and DERFCX),
C   and one SUBROUTINE type subprogram, CALERF.  The calling
C   statements for the primary entries are:
C
C                   Y=ERF(X)     (or   Y=DERF(X)),
C
C                   Y=ERFC(X)    (or   Y=DERFC(X)),
C   and
C                   Y=ERFCX(X)   (or   Y=DERFCX(X)).
C
C   The routine  CALERF  is intended for internal packet use only,
C   all computations within the packet being concentrated in this
C   routine.  The function subprograms invoke  CALERF  with the
C   statement
C
C          CALL CALERF(ARG,RESULT,JINT)
C
C   where the parameter usage is as follows
C
C      Function                     Parameters for CALERF
C       call              ARG                  Result          JINT
C
C     ERF(ARG)      ANY REAL ARGUMENT         ERF(ARG)          0
C     ERFC(ARG)     ABS(ARG) .LT. XBIG        ERFC(ARG)         1
C     ERFCX(ARG)    XNEG .LT. ARG .LT. XMAX   ERFCX(ARG)        2
C
C   The main computation evaluates near-minimax approximations
C   from "Rational Chebyshev approximations for the error function"
C   by W. J. Cody, Math. Comp., 1969, PP. 631-638.  This
C   transportable program uses rational functions that theoretically
C   approximate  erf(x)  and  erfc(x)  to at least 18 significant
C   decimal digits.  The accuracy achieved depends on the arithmetic
C   system, the compiler, the intrinsic functions, and proper
C   selection of the machine-dependent constants.
C
C*******************************************************************
C*******************************************************************
C
C Explanation of machine-dependent constants.  Let
C
C   XMIN   = the smallest positive floating-point number.
C
C Then the following machine-dependent constants must be declared
C   in DATA statements.  IEEE values are provided as a default.
C
C   XINF   = the largest positive finite floating-point number.
C   XNEG   = the largest negative argument acceptable to ERFCX;
C            the negative of the solution to the equation
C            2*exp(x*x) = XINF.
C   XSMALL = argument below which erf(x) may be represented by
C            2*x/sqrt(pi)  and above which  x*x  will not underflow.
C            A conservative value is the largest machine number X
C            such that   1.0 + X = 1.0   to machine precision.
C   XBIG   = largest argument acceptable to ERFC;  solution to
C            the equation:  W(x) * (1-0.5/x**2) = XMIN,  where
C            W(x) = exp(-x*x)/[x*sqrt(pi)].
C   XHUGE  = argument above which  1.0 - 1/(2*x*x) = 1.0  to
C            machine precision.  A conservative value is
C            1/[2*sqrt(XSMALL)]
C   XMAX   = largest acceptable argument to ERFCX; the minimum
C            of XINF and 1/[sqrt(pi)*XMIN].
C
C   Approximate values for some important machines are:
C
C                          XMIN       XINF        XNEG     XSMALL
C
C  CDC 7600      (S.P.)  3.13E-294   1.26E+322   -27.220  7.11E-15
C  CRAY-1        (S.P.)  4.58E-2467  5.45E+2465  -75.345  7.11E-15
C  IEEE (IBM/XT,
C    SUN, etc.)  (S.P.)  1.18E-38    3.40E+38     -9.382  5.96E-8
C  IEEE (IBM/XT,
C    SUN, etc.)  (D.P.)  2.23D-308   1.79D+308   -26.628  1.11D-16
C  IBM 195       (D.P.)  5.40D-79    7.23E+75    -13.190  1.39D-17
C  UNIVAC 1108   (D.P.)  2.78D-309   8.98D+307   -26.615  1.73D-18
C  VAX D-Format  (D.P.)  2.94D-39    1.70D+38     -9.345  1.39D-17
C  VAX G-Format  (D.P.)  5.56D-309   8.98D+307   -26.615  1.11D-16
C
C
C                          XBIG       XHUGE       XMAX
C
C  CDC 7600      (S.P.)  25.922      8.39E+6     1.80X+293
C  CRAY-1        (S.P.)  75.326      8.39E+6     5.45E+2465
C  IEEE (IBM/XT,
C    SUN, etc.)  (S.P.)   9.194      2.90E+3     4.79E+37
C  IEEE (IBM/XT,
C    SUN, etc.)  (D.P.)  26.543      6.71D+7     2.53D+307
C  IBM 195       (D.P.)  13.306      1.90D+8     7.23E+75
C  UNIVAC 1108   (D.P.)  26.582      5.37D+8     8.98D+307
C  VAX D-Format  (D.P.)   9.269      1.90D+8     1.70D+38
C  VAX G-Format  (D.P.)  26.569      6.71D+7     8.98D+307
C
C*******************************************************************
C*******************************************************************
C
C Error returns
C
C  The program returns  ERFC = 0      for  ARG .GE. XBIG;
C
C                       ERFCX = XINF  for  ARG .LT. XNEG;
C      and
C                       ERFCX = 0     for  ARG .GE. XMAX.
C
C
C Intrinsic functions required are:
C
C     ABS, AINT, EXP
C
C
C  Author: W. J. Cody
C          Mathematics and Computer Science Division
C          Argonne National Laboratory
C          Argonne, IL 60439
C
C  Latest modification: March 12, 1992
C
C------------------------------------------------------------------
      INTEGER I,JINT
CS    REAL
      DOUBLE PRECISION
     1     A,ARG,B,C,D,DEL,FOUR,HALF,P,ONE,Q,RESULT,SIXTEN,SQRPI,
     2     TWO,THRESH,X,XBIG,XDEN,XHUGE,XINF,XMAX,XNEG,XNUM,XSMALL,
     3     Y,YSQ,ZERO
      DIMENSION A(5),B(4),C(9),D(8),P(6),Q(5)
C------------------------------------------------------------------
C  Mathematical constants
C------------------------------------------------------------------
CS    DATA FOUR,ONE,HALF,TWO,ZERO/4.0E0,1.0E0,0.5E0,2.0E0,0.0E0/,
CS   1     SQRPI/5.6418958354775628695E-1/,THRESH/0.46875E0/,
CS   2     SIXTEN/16.0E0/
      DATA FOUR,ONE,HALF,TWO,ZERO/4.0D0,1.0D0,0.5D0,2.0D0,0.0D0/,
     1     SQRPI/5.6418958354775628695D-1/,THRESH/0.46875D0/,
     2     SIXTEN/16.0D0/
C------------------------------------------------------------------
C  Machine-dependent constants
C------------------------------------------------------------------
CS    DATA XINF,XNEG,XSMALL/3.40E+38,-9.382E0,5.96E-8/,
CS   1     XBIG,XHUGE,XMAX/9.194E0,2.90E3,4.79E37/
      DATA XINF,XNEG,XSMALL/1.79D308,-26.628D0,1.11D-16/,
     1     XBIG,XHUGE,XMAX/26.543D0,6.71D7,2.53D307/
C------------------------------------------------------------------
C  Coefficients for approximation to  erf  in first interval
C------------------------------------------------------------------
CS    DATA A/3.16112374387056560E00,1.13864154151050156E02,
CS   1       3.77485237685302021E02,3.20937758913846947E03,
CS   2       1.85777706184603153E-1/
CS    DATA B/2.36012909523441209E01,2.44024637934444173E02,
CS   1       1.28261652607737228E03,2.84423683343917062E03/
      DATA A/3.16112374387056560D00,1.13864154151050156D02,
     1       3.77485237685302021D02,3.20937758913846947D03,
     2       1.85777706184603153D-1/
      DATA B/2.36012909523441209D01,2.44024637934444173D02,
     1       1.28261652607737228D03,2.84423683343917062D03/
C------------------------------------------------------------------
C  Coefficients for approximation to  erfc  in second interval
C------------------------------------------------------------------
CS    DATA C/5.64188496988670089E-1,8.88314979438837594E0,
CS   1       6.61191906371416295E01,2.98635138197400131E02,
CS   2       8.81952221241769090E02,1.71204761263407058E03,
CS   3       2.05107837782607147E03,1.23033935479799725E03,
CS   4       2.15311535474403846E-8/
CS    DATA D/1.57449261107098347E01,1.17693950891312499E02,
CS   1       5.37181101862009858E02,1.62138957456669019E03,
CS   2       3.29079923573345963E03,4.36261909014324716E03,
CS   3       3.43936767414372164E03,1.23033935480374942E03/
      DATA C/5.64188496988670089D-1,8.88314979438837594D0,
     1       6.61191906371416295D01,2.98635138197400131D02,
     2       8.81952221241769090D02,1.71204761263407058D03,
     3       2.05107837782607147D03,1.23033935479799725D03,
     4       2.15311535474403846D-8/
      DATA D/1.57449261107098347D01,1.17693950891312499D02,
     1       5.37181101862009858D02,1.62138957456669019D03,
     2       3.29079923573345963D03,4.36261909014324716D03,
     3       3.43936767414372164D03,1.23033935480374942D03/
C------------------------------------------------------------------
C  Coefficients for approximation to  erfc  in third interval
C------------------------------------------------------------------
CS    DATA P/3.05326634961232344E-1,3.60344899949804439E-1,
CS   1       1.25781726111229246E-1,1.60837851487422766E-2,
CS   2       6.58749161529837803E-4,1.63153871373020978E-2/
CS    DATA Q/2.56852019228982242E00,1.87295284992346047E00,
CS   1       5.27905102951428412E-1,6.05183413124413191E-2,
CS   2       2.33520497626869185E-3/
      DATA P/3.05326634961232344D-1,3.60344899949804439D-1,
     1       1.25781726111229246D-1,1.60837851487422766D-2,
     2       6.58749161529837803D-4,1.63153871373020978D-2/
      DATA Q/2.56852019228982242D00,1.87295284992346047D00,
     1       5.27905102951428412D-1,6.05183413124413191D-2,
     2       2.33520497626869185D-3/
C------------------------------------------------------------------
      X = ARG
      Y = ABS(X)
      IF (Y .LE. THRESH) THEN
C------------------------------------------------------------------
C  Evaluate  erf  for  |X| <= 0.46875
C------------------------------------------------------------------
            YSQ = ZERO
            IF (Y .GT. XSMALL) YSQ = Y * Y
            XNUM = A(5)*YSQ
            XDEN = YSQ
            DO 20 I = 1, 3
               XNUM = (XNUM + A(I)) * YSQ
               XDEN = (XDEN + B(I)) * YSQ
   20       CONTINUE
            RESULT = X * (XNUM + A(4)) / (XDEN + B(4))
            IF (JINT .NE. 0) RESULT = ONE - RESULT
            IF (JINT .EQ. 2) RESULT = EXP(YSQ) * RESULT
            GO TO 800
C------------------------------------------------------------------
C  Evaluate  erfc  for 0.46875 <= |X| <= 4.0
C------------------------------------------------------------------
         ELSE IF (Y .LE. FOUR) THEN
            XNUM = C(9)*Y
            XDEN = Y
            DO 120 I = 1, 7
               XNUM = (XNUM + C(I)) * Y
               XDEN = (XDEN + D(I)) * Y
  120       CONTINUE
            RESULT = (XNUM + C(8)) / (XDEN + D(8))
            IF (JINT .NE. 2) THEN
               YSQ = AINT(Y*SIXTEN)/SIXTEN
               DEL = (Y-YSQ)*(Y+YSQ)
               RESULT = EXP(-YSQ*YSQ) * EXP(-DEL) * RESULT
            END IF
C------------------------------------------------------------------
C  Evaluate  erfc  for |X| > 4.0
C------------------------------------------------------------------
         ELSE
            RESULT = ZERO
            IF (Y .GE. XBIG) THEN
               IF ((JINT .NE. 2) .OR. (Y .GE. XMAX)) GO TO 300
               IF (Y .GE. XHUGE) THEN
                  RESULT = SQRPI / Y
                  GO TO 300
               END IF
            END IF
            YSQ = ONE / (Y * Y)
            XNUM = P(6)*YSQ
            XDEN = YSQ
            DO 240 I = 1, 4
               XNUM = (XNUM + P(I)) * YSQ
               XDEN = (XDEN + Q(I)) * YSQ
  240       CONTINUE
            RESULT = YSQ *(XNUM + P(5)) / (XDEN + Q(5))
            RESULT = (SQRPI -  RESULT) / Y
            IF (JINT .NE. 2) THEN
               YSQ = AINT(Y*SIXTEN)/SIXTEN
               DEL = (Y-YSQ)*(Y+YSQ)
               RESULT = EXP(-YSQ*YSQ) * EXP(-DEL) * RESULT
            END IF
      END IF
C------------------------------------------------------------------
C  Fix up for negative argument, erf, etc.
C------------------------------------------------------------------
  300 IF (JINT .EQ. 0) THEN
            RESULT = (HALF - RESULT) + HALF
            IF (X .LT. ZERO) RESULT = -RESULT
         ELSE IF (JINT .EQ. 1) THEN
            IF (X .LT. ZERO) RESULT = TWO - RESULT
         ELSE
            IF (X .LT. ZERO) THEN
               IF (X .LT. XNEG) THEN
                     RESULT = XINF
                  ELSE
                     YSQ = AINT(X*SIXTEN)/SIXTEN
                     DEL = (X-YSQ)*(X+YSQ)
                     Y = EXP(YSQ*YSQ) * EXP(DEL)
                     RESULT = (Y+Y) - RESULT
               END IF
            END IF
      END IF
  800 RETURN
C---------- Last card of CALERF ----------
      END

max_stars_repo_name: sbwhitecap/clocc-hg

max_issues_repo_name: sbwhitecap/clocc-hg

max_forks_repo_name: sbwhitecap/clocc-hg