檢視 Subroutine compositionmapgen 的原始碼

[[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 <br><br>
    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   <br><br>
    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) <br>
  end subroutine