| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in), | optional | :: | fileID_in |
SUBROUTINE Analyze(fileID_in) ! MODULES USE MODgvec_Analyze_Vars USE MODgvec_MHD3D_Vars, ONLY:which_init USE MODgvec_mhd3d_visu IMPLICIT NONE !----------------------------------------------------------------------------------------------------------------------------------- ! INPUT VARIABLES INTEGER, INTENT(IN), OPTIONAL :: fileID_in !----------------------------------------------------------------------------------------------------------------------------------- ! OUTPUT VARIABLES !----------------------------------------------------------------------------------------------------------------------------------- ! LOCAL VARIABLES LOGICAL :: vcase(4) CHARACTER(LEN=4) :: vstr INTEGER :: FileID !=================================================================================================================================== IF(.NOT.MPIroot) RETURN iAnalyze=iAnalyze+1 IF(PRESENT(FileID_in))THEN FileID=FileID_in ELSE FileID=iAnalyze END IF IF(iAnalyze.EQ.0) THEN IF(which_init.EQ.1) THEN IF(visu1D.GT.0) CALL VMEC1D_visu() IF(visu2D.GT.0) CALL VMEC3D_visu(np_visu_planes,visu_planes_minmax,.TRUE. ) IF(visu3D.GT.0) CALL VMEC3D_visu(np_visu_3D ,visu_3D_minmax ,.FALSE.) END IF END IF !iAnalyze==0 IF(visu1D.GT.0)THEN CALL visu_1d_modes(np_1d,FileID) ! END IF !visu1D IF(visu2D.GT.0)THEN vcase=.FALSE. WRITE(vstr,'(I4)')visu2D IF(INDEX(vstr,'1').NE.0) vcase(1)=.TRUE. IF(INDEX(vstr,'2').NE.0) vcase(2)=.TRUE. IF(INDEX(vstr,'3').NE.0) vcase(3)=.TRUE. IF(INDEX(vstr,'4').NE.0) vcase(4)=.TRUE. IF(vcase(1))THEN IF(iAnalyze.EQ.0) CALL visu_BC_face(np_visu_BC(1:2),visu_BC_minmax(:,:),FileID) END IF IF(vcase(2))THEN CALL visu_3D(np_visu_planes,visu_planes_minmax,.TRUE.,FileID) !only planes END IF END IF !visu2d IF(visu3D.GT.0)THEN vcase=.FALSE. WRITE(vstr,'(I4)')visu3D IF(INDEX(vstr,'1').NE.0) vcase(1)=.TRUE. IF(INDEX(vstr,'2').NE.0) vcase(2)=.TRUE. IF(INDEX(vstr,'3').NE.0) vcase(3)=.TRUE. IF(INDEX(vstr,'4').NE.0) vcase(4)=.TRUE. IF(vcase(1))THEN CALL visu_3D(np_visu_3D,visu_3D_minmax,.FALSE.,FileID) !full 3D END IF END IF !visu2d END SUBROUTINE Analyze