FEATURE - Add sol_num parameter to match routine. Adding sol_num parameter and set the wanted number of solution component then this will make Solution{sol_num}.bin file containing the {sol_num}th solution on Wt. If there aren't any sol_num then this will work just as same as before.

Author: satelite2517

match/ideal.f

@@ -305,7 +305,9 @@ END SUBROUTINE ideal_transform
 c-----------------------------------------------------------------------
 c     declarations.
 c-----------------------------------------------------------------------
-      SUBROUTINE ideal_build
+      SUBROUTINE ideal_build(sol_num)
+
+      INTEGER, INTENT(IN) :: sol_num
 
       INTEGER :: istep,ifix,jfix,kfix,ieq,info
       INTEGER, DIMENSION(mpert) :: ipiv
@@ -319,7 +321,8 @@ SUBROUTINE ideal_build
          uedge=0
          uedge(mripple-mlow+1)=1
       ELSE
-         uedge=wt(:,1)
+         WRITE(*,*)"sol_num = ",sol_num
+         uedge=wt(:,sol_num)
       ENDIF
       temp2=soltype(mstep)%u(:,1:mpert,1)
       CALL zgetrf(mpert,mpert,temp2,mpert,ipiv,info)
@@ -357,12 +360,15 @@ END SUBROUTINE ideal_build
 c-----------------------------------------------------------------------
 c     declarations.
 c-----------------------------------------------------------------------
-      SUBROUTINE ideal_write
+      SUBROUTINE ideal_write(sol_num)
+
+      INTEGER, INTENT(IN) :: sol_num
 
       LOGICAL, PARAMETER :: diagnose=.FALSE.
       INTEGER :: ipert,istep,m
       REAL(r8), DIMENSION(mpert) :: singfac
       COMPLEX(r8), DIMENSION(mpert,0:mstep) :: b
+      CHARACTER(128) :: outputname
 c-----------------------------------------------------------------------
 c     compute normal perturbed magnetic field.
 c-----------------------------------------------------------------------
@@ -375,7 +381,14 @@ SUBROUTINE ideal_write
 c     write binary output for graphs.
 c-----------------------------------------------------------------------
       WRITE(*,*)"Write binary output for graphs."
-      CALL bin_open(bin_unit,"solutions.bin","UNKNOWN","REWIND","none")
+
+      IF (sol_num == 1) THEN
+         outputname = "solutions.bin"
+      ELSE
+         WRITE(outputname, '(A,I0,A)') "solutions", sol_num, ".bin"
+      END IF
+
+      CALL bin_open(bin_unit, outputname, "UNKNOWN", "REWIND", "none")
       DO ipert=1,mpert
          DO istep=0,mstep
             WRITE(bin_unit)REAL(psifac(istep),4),REAL(rho(istep),4),

match/match.f

@@ -29,7 +29,7 @@ PROGRAM match_main
      $     big,small,ns,root,scan_flag,root_flag,contour_flag,full,
      $     z_level_type,phi_level,transform_flag,poly_form_type,
      $     poly_test_type,scan0,scan1,nscan,delta0,delta1,condense,
-     $     ripple_flag,mripple
+     $     ripple_flag,mripple,sol_num
 c-----------------------------------------------------------------------
 c     format statements.
 c-----------------------------------------------------------------------
@@ -55,8 +55,8 @@ PROGRAM match_main
 c-----------------------------------------------------------------------
       IF(ideal_flag)THEN
          CALL ideal_transform
-         CALL ideal_build
-         CALL ideal_write
+         CALL ideal_build(sol_num)
+         CALL ideal_write(sol_num)
          CALL ideal_chord
          IF(contour_flag)CALL ideal_contour
       ENDIF

match/match_mod.f

@@ -42,6 +42,7 @@ MODULE match_mod
      $     kin_flag=.FALSE.,con_flag=.FALSE.
       LOGICAL, DIMENSION(:), POINTER :: sing_flag
       INTEGER :: mfix,mhigh,mlow,mpert,mstep,nn,msing,mmatch,mbit,mterm
+      INTEGER :: sol_num = 1
       INTEGER, DIMENSION(:), POINTER :: fixstep
       REAL(r8) :: sfac0=1
       REAL(r8), DIMENSION(:), POINTER :: psifac,rho,q