"Subroutine compositionmapgen" 修訂間的差異

於 2017年12月5日 (二) 09:54 的最新修訂

Subroutine compositionmapgen is a subroutine used to called after the program read-in the gmsh file. This routine is designed for user to generate composition map for the random alloy fluctuation. However, it is not limited on this. This one usually goes with subsroutine readEg. And this subroutine is activated by the command $gencompositionmap. The example is as follows

 Subroutine compositionmapgen
   use meshelement
   use global3d
   use compositionvar
   implicit none
   real(dp) :: x,y,z,qw_begin,qw_end
   integer(i8b) :: nodetype ,i ,j, k 


   open(155,FILE="In.out")
   write(155,"(a11)")"$MeshFormat"
   write(155,"(a5)")"2 0 8"
   write(155,"(a9)")"$NodeData"
   write(155,"(a1)")"1"
   write(155,"(a8)")"Al_check"
   write(155,"(a1)")"1"
   write(155,"(a3)")"1.0"
   write(155,"(a1)")"3"
   write(155,"(a1)")"0"
   write(155,"(a1)")"1"
   write(155,*) mshnd%n   


   do i=1,mshnd%n ! Eg_check.out
       nodetype = mshnd%element(i) ! Get the node number
       x=mshnd%x(i)
       y=mshnd%y(i)
       z=mshnd%z(i)
       Incom_inter =0.0
       if(z>=qw_begin .and. z<= qw_end) then
           Incom_inter=randm()        
       end if
       write(155,*)  i ,Incom_inter 
    end do
    write(155,*) "$EndNodeData"
    close(155) 

 end subroutine