!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     
     !program main 
     !&主程序 
      program main
      REAL(8),DIMENSION(12,1000)::A
      real(8)::slide_R,nnx,mmcx,mmsx,fN,fM,shrinkstrain1,          &
               shrinkstrain,FF,F,t,t0,tt,t1,t01,Fck,RH,G,           &
               D1,D2,D3,W1,W2,W3,deflection,Qx,slide,strainc,       &
               strains,stressc,stresss,Ms,Mc,Nx,p,U,V,R,R1,R2,       &
               Ss,I1,L,L1,x,H,H1,Hcc,Hc,Ec,Es,Ac,b,b1,alpha,       &
               Hss,Hs,As,Is,Ic,Yc,Ys,ysc
      integer::n,i,ii,nn,j,jj,jjj
      integer,parameter::m=11
      A=0.0
      nn=m-1
      write(*,100)'Please Look up In_put.txt and Input Data'
100   format(//A50//)
      pause
     !Input the Data of Basic Parameter from In_put.value 
     !& 从In_put.value文件输入基本参量
      OPEN(8,FILE='in_put.txt',STATUS='old')
      read(8,'(17F20.5,I5)')slide_R,Fck,RH,t1,t01,p,R1,R2,L,L1,    &
                            Hcc,Ec,Es,b,Hss,As,Is,n
      write(*,*)slide_R
      close(8)    
     !Solve the Chatacteristic Value of Section and Stiff of 
     !Shear connectors
     !& 求解组合梁截面的特征值和剪力连接件的刚度      
      call main6(slide_R,p,U,V,R,R1,R2,Ss,I1,H,H1,Hcc,Hc,Ec,Es,    &
                 Ac,b,b1,alpha,Hss,Hs,As,Is,Ic,Yc,Ys,ysc)
     !Solve Midle varies of Elastic Solution
     !& 求解弹性解的中间变量     
           U=R*(1.0/(Ec*Ac)+1.0/(Es*As)+H1**2.0/(Ec*Ic+Es*Is))
           V=-R*(H1/(Ec*Ic+Es*Is))
           G=1.0/(Ec*Ic+Es*Is)
           D1=p*(1.0+(V*H1)/U)
           D2=(p*L*(1.0+(V*H1)/U))/2.0
           D3=(p*(L-L1)*(1.0+(V*H1)/U))/2.0
           W1=(L1*(-D1*L1+4.0*(D2-D3)))/8.0
           W2=(L*(D1*(L**2.0-3.0*L1**2.0)-6.0*D2*(L-2.0*L1)           &
              -12.0*D3*L1))/48.0
           W3=D1*(L**3.0-3.0*L*L1**2.0+2.0*L1**3.0)/48.0           &
              -(D2*(L-L1)**2.0)/8.0-(D3*L1*(2.0*L-L1))/8.0
           FF=F(RH,t01,H1,Fck,Ec,Ic,Es,Is)
     !Solve Elastic Solution 
     !& 求解弹性解
     !The Discretizations of x along the Beam Axis are Carried Out
     !& 沿梁轴线方向对梁长进行离散化   
         do 10 i=0,n
            x=-i*L/(2.0*n)
            if (x.ge.-L1/2.0)then
               Nx=-(V*p)*(L-L1)/(2.0*U)
               Mc=(Ec*Ic*p*(L-L1)*(1.0+(V*H1)/U))/(2.0*(Ec*Ic+Es*Is))
               Ms=(Es*Is*p*(L-L1)*(1.0+(V*H1)/U))/(2.0*(Ec*Ic+Es*Is))
               Qx=0.0
               deflection=G*((D3*x**2.0)/2.0+W3)
            else
               Nx=-(V*p)*(L+2.0*x)/(2.0*U)
               Mc=(Ec*Ic*p*(1.0+(V*H1)/U)*(x+L/2.0))/(Ec*Ic+Es*Is)
               Ms=(Es*Is*p*(1.0+(V*H1)/U)*(x+L/2.0))/(Ec*Ic+Es*Is)
               Qx=-(V*p)/U
               deflection=G*((D1*x**3.0)/6.0+(D2*x**2.0)/2+W1*x+W2)
            end if
            strainc=-Nx/(Ec*Ac)+(Mc*Hc)/(Ec*Ic)
            stressc=strainc*Ec
            strains=Nx/(Es*As)-(Ms*Hs)/(Es*Is)
            stresss=strains*Es
            slide=Qx/R
            A(1,1)=t01
            A(2,1)=x
            A(3,1)=nx
            A(4,1)=mc
            A(5,1)=ms
            A(6,1)=strainc
            A(7,1)=strains
            A(8,1)=stressc
            A(9,1)=stresss
            A(10,1)=qx
            A(11,1)=slide
            A(12,1)=deflection
            open(9,file='out_put.txt')
            write(9,'(12A20)')'t=','x=','nx=','mc=','ms=',           &
                              'strainc=','strains=','stressc=',       &
                              'stresss=','qx=','slide=',           &
                              'deflection='
      !The Discretizations in Time are initially Carried Out
      !对时间t进行初始离散化        
            tt=(t1-t01)/nn
            do 50 ii=1,nn
               t=t01+ii*tt
               t0=t01
      !Set initiatal Value to Unknown Increments
      !& 给未知增量赋初值           
               nx=A(3,1)
               mc=A(4,1)
               ms=A(5,1)
               strainc=A(6,1)
               strains=A(7,1)
               stressc=A(8,1)
               stresss=A(9,1)
               qx=A(10,1)
               slide=A(11,1)
               deflection=A(12,1)
               nnx=0.0
               mmcx=0.0
               mmsx=0.0
               fN=0.0
               fM=0.0
               shrinkstrain1=0.0
               shrinkstrain=0.0
      !Solve Viscoelastic Solution
      !& 求解粘弹性解
           write(*,*)
               call main2(nnx,mmcx,mmsx,fN,fM,shrinkstrain1,       &
                          shrinkstrain,x,Nx,Mc,Ms,strainc,           &
                          strains,stressc,stresss,Qx,slide,           &
                          deflection,p,U,V,R,L,L1,t,tt,t0,Ec,      &
                          Es,Is,H,Hcc,b,Hss,As,FF)
                  jjj=ii+1
               A(1,jjj)=t
               A(2,jjj)=x
               A(3,jjj)=nx
               A(4,jjj)=mc
               A(5,jjj)=ms
               A(6,jjj)=strainc
               A(7,jjj)=strains
               A(8,jjj)=stressc
               A(9,jjj)=stresss
               A(10,jjj)=qx
               A(11,jjj)=slide
               A(12,jjj)=deflection
50            continue
       !Output the Data of Elastic Solution and viscoelastic Solution 
      !to out_put.value
      !& 输出弹性和粘弹性解到out_put.value文件    
            write(9,300)((a(j,jj),j=1,12),jj=1,m)
300            format(1x,12f20.10)
10       continue

       write(9,'(A20)')'R='
       write(9,500)R
500       format(1x,f20.10)
       close(9)
      !Finish Program Main
      !& 结束主程序
       write(*,200)'Please Look up Out_put.value and Found&
                    the Data of Output'
200       format(//A60//)
       end program main