!______________________________________________________________________________ 
! SUBROUTINE PlotNetwork2D(heads_tails, node_coords, 
! arc_values, arc_labels, node_values, node_labels, 
! number_arcs, number_nodes, 
! color_nodes, resize_nodes, color_arcs, resize_arcs, 
! arc_size_range, arc_color_range, node_size_range, node_color_range, 
! arc_color_table, node_color_table, vector_type, axis, my_func)
!______________________________________________________________________________ 
PROGRAM Network_Plot_test
  USE Precision, ONLY: wp=>r_wp 
  USE Network_Graphics    
  IMPLICIT NONE 
   
  INTEGER, PARAMETER :: n_nodes=5, n_arcs=2*n_nodes, n_dim=2 
  REAL(KIND=wp), DIMENSION(n_nodes) :: node_values
  REAL(KIND=wp), DIMENSION(n_arcs) :: arc_values
  REAL, DIMENSION(n_dim,n_nodes) :: node_coords  
  INTEGER, DIMENSION(2,n_arcs) :: heads_tails 
  INTEGER :: i,j,k, node, arc  
  REAL, DIMENSION(n_dim,n_arcs) :: random_array
   
  !__________________________________ 
  CALL RANDOM_SEED()
  CALL RANDOM_NUMBER(random_array)
  heads_tails=INT(REAL(n_nodes)*random_array(1:2,:))+1
  CALL RANDOM_NUMBER(random_array)
  node_coords=random_array(:,1:n_nodes)  
  node_values=random_array(1,1:n_nodes)
  arc_values=random_array(2,:)*10.0
  ! WRITE(*,"(A,25I5)") "Heads: ",heads_tails(1,:)
  ! WRITE(*,"(A,25I5)") "Tails: ",heads_tails(2,:)  
  !__________________________________ 

  DO i=1,1
	  ! This is an example of using PlotNetwork2D plot: 
	  CALL InitNetworkGraphics(file="TestNetwork2D.ps",file_type="CONS",  &
	                    plot_title=(/"Example of a network plot                   ", &
	                                 "                                            ", &  
	                                 "$v(i)=i R \[1 + \tanh(\frac{i-i_c}{\xi})\]  "/), &
	                    x_label="Current flow \Huge{$\rightarrow$}", &	                       
                       page_size=(/5000,3000/) ) !, label_format="9F1")
	  CALL PlotNetwork2D(heads_tails=heads_tails,node_coords=node_coords, &
	                     node_values=node_values, arc_values=arc_values, &
	                     color_nodes=.TRUE.,resize_nodes=.TRUE., &
	                     color_arcs=.TRUE.,resize_arcs=.TRUE., &	                                          
	                     vector_type=1311, axis=(/0.0,0.0,1.0,1.0/), &
	                     color_table="RAIN", colorbar_position="VERTICAL", &
	                     node_colorbar_format="5D1",arc_colorbar_format="5E1", &
	                     node_labels=(/(node**2,node=1,n_nodes)/), &
	                     number_arcs=.TRUE., &
	                     axis_labels_format=(/"NAME","LINE","NONE","LINE"/), &	                  
                        arc_size_range=(/0.0,25.0/),arc_color_range=(/0.0,1.0/), &
                        node_size_range=(/-1.0,2.0/),node_color_range=(/-1.0,2.0/))	  	                                         
	  CALL EndNetworkGraphics()  
  END DO
   
END PROGRAM Network_Plot_test