program row_tst
      
      integer ns
      parameter(ns=5)

      integer np
      parameter(np=180)

      integer specie

      real theta_i,phi_i,k_0
      common /incident/ theta_i,phi_i,k_0

      integer canopy_coho,specie_coho,plant_coho
      common /coho_flags/ canopy_coho,specie_coho,plant_coho

      complex ss_specie_inco(4,ns),ss_specie_coho(4,ns)
      complex ss_specie_cor(4,ns),ss_specie(4,ns)
      common /specie_scat_save/ss_specie_inco,ss_specie_coho,
     &     ss_specie_cor,ss_specie

      character *15 illumination
      common /antenna_pattern/ illumination

      real N_pprcm,row_spacing,row_sig
      common /row_data/ N_pprcm,row_spacing,row_sig

      real spot_size
      common /pattern_data/ spot_size

      complex ss_inco(4,ns),ss_coho(4,ns)
      complex s_plant_out(3,ns)
      
      character *15 row_plant_pdf_x,row_plant_pdf_y,row_plant_depend
      common /row_plant_dist/ row_plant_pdf_x,row_plant_pdf_y,
     &                        row_plant_depend
      
      character *15 row_space_stat
      common /row_space_model/ row_space_stat
      
      real phi_min,phi_max,phi_delta
      integer phi_sweep_flg
      common /phi_sweep_data/phi_min,phi_max,phi_delta,phi_sweep_flg

      complex ss_row_cor_save(4,ns,np),ss_row_coho_save(4,ns,np),
     &      ss_row_save(4,ns,np)
      real phi_save(np)
      common /row_phi_save/ss_row_cor_save,ss_row_coho_save,
     &      ss_row_save,phi_save

      real ss_row_cor_r,ss_row_coho_r,ss_row_r
      
      real pi,degrees

      real save(30,30,3)

      integer i,j,ii,N_phi_smple
      
      real phi_i_degrees,phi_min_degrees,phi_max_degrees,
     &     phi_delta_degrees

      pi=4.0*atan(1.0)
      degrees=pi/180.0

      specie=1

      row_plant_pdf_x='gaussian'
      row_plant_pdf_y='uniform'
      row_plant_depend='independent'
      row_space_stat='deterministic'
      

      read*, k_0
      canopy_coho=0
      read*, specie_coho
      plant_coho=0
      read*, illumination
      read*, N_pprcm
      read*, row_spacing
      read*, row_sig
      read*, spot_size
      read*, phi_min_degrees
      read*, phi_max_degrees
      read*, phi_delta_degrees
      read*, phi_sweep_flg
      read*, phi_i_degrees

      phi_min=phi_min_degrees*degrees
      phi_max=phi_max_degrees*degrees      
      phi_delta=phi_delta_degrees*degrees
      
Cc      do i=1,26
C
Cc        theta_i=pi-(5.0+i*2.5)*degrees
C
c      theta_i=pi-40.0*degrees
C
C      i=1
C
C        do j=1,25
C
C          phi_i=(-2.5+j*2.5)*degrees
      

      phi_i=phi_i_degrees*degrees
      theta_i=pi-40.0*degrees
      i=1
      j=1

          call row_power(specie,ss_inco,
     &         ss_coho,s_plant_out)

          save(i,j,1)=10.*alog10(cabs(ss_inco(1,1)))
          save(i,j,2)=10.*alog10(cabs(ss_coho(1,1)))
          save(i,j,3)=10.*alog10(cabs(ss_coho(1,1)+ss_inco(1,1)))
          
c          save(i,j,1)=ss_inco(1,1)
c          save(i,j,2)=ss_coho(1,1)
c          save(i,j,3)=ss_coho(1,1)+ss_inco(1,1)
C
C        enddo
C        
Cc      enddo
C
C          
Cc      do i=1,7
C
Cc        theta_i=(0.0+i*10.0)
C
Cc        ii=i*4.0-2.0
C
C        ii=1
C
C        do j=1,25
C
C          phi_i=(-2.5+j*2.5)

          ii=1
         print*, 'theta_i,phi_i,inco,coho,total'
         
          print*, (pi-theta_i)/degrees,phi_i/degrees,
     &         save(ii,j,1),save(ii,j,2),save(ii,j,3)
          
c        enddo
        
c      enddo

C      do i=1,7
C
C        phi_i=(-10.0+i*10.0)
C
C        ii=i*4.0-3.0
C
C        do j=1,26
C
C          theta_i=(5.0+j*2.5)
C
C          print*, theta_i,phi_i,
C     &         save(j,ii,1),save(j,ii,2),save(j,ii,3)
C          
C        enddo
C        
C      enddo

          if(phi_sweep_flg .eq. 1)then

            N_phi_smple=int((phi_max-phi_min)/phi_delta)+1

          else

            N_phi_smple=0

          endif
      
        print*,'phi_i,ss_row_cor_r,ss_row_r'
        
      do i=1,N_phi_smple+1

        ss_row_cor_r=
     &       10.*alog10(cabs(ss_row_cor_save(1,1,i)))
        ss_row_r=
     &       10.*alog10(cabs(ss_row_save(1,1,i)))
        

        
        print*,phi_save(i)/degrees,
     &       ss_row_cor_r,ss_row_r

      enddo
      
      stop
      end
      

      subroutine plant_field(s_plant)

      complex s_plant(3)

      s_plant(1)=cmplx(0.8,0.0)
      s_plant(2)=cmplx(0.8,0.0)
      s_plant(3)=cmplx(0.0,0.0)

      return
      end

      subroutine plant_power(ss_inco,ss_coho)

      complex ss_inco(4),ss_coho(4)

      ss_inco(1)=cmplx(1.0,0.0)
      ss_inco(2)=cmplx(1.0,0.0)
      ss_inco(3)=cmplx(1.0,0.0)

      ss_coho(1)=cmplx(0.0,0.0)
      ss_coho(2)=cmplx(0.0,0.0)
      ss_coho(3)=cmplx(0.0,0.0)

      return
      end