564 |
DO K=kMin, kMax, kStr |
DO K=kMin, kMax, kStr |
565 |
DO J=jMin, jMax, jStr |
DO J=jMin, jMax, jStr |
566 |
DO I=iMin, iMax, iStr |
DO I=iMin, iMax, iStr |
567 |
IF ( fld(I,J,K,bi,bj) .NE. 0. ) THEN |
IF (printMapIncludesZeros .OR. fld(I,J,K,bi,bj) .NE. 0.) THEN |
568 |
IF ( fld(I,J,K,bi,bj) .LT. fMin ) |
IF ( fld(I,J,K,bi,bj) .LT. fMin ) |
569 |
& fMin = fld(I,J,K,bi,bj) |
& fMin = fld(I,J,K,bi,bj) |
570 |
IF ( fld(I,J,K,bi,bj) .GT. fMax ) |
IF ( fld(I,J,K,bi,bj) .GT. fMax ) |
601 |
& '// CMAX = ', fMax |
& '// CMAX = ', fMax |
602 |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
603 |
& SQUEEZE_RIGHT, 1) |
& SQUEEZE_RIGHT, 1) |
604 |
WRITE(msgBuf,'(A,1PE30.15)') |
IF ( validRange ) THEN |
605 |
& '// CINT = ', fRange/FLOAT(lChlist-1) |
WRITE(msgBuf,'(A,1PE30.15)') |
606 |
|
& '// CINT = ', fRange/FLOAT(lChlist-1) |
607 |
|
ELSE |
608 |
|
WRITE(msgBuf,'(A,1PE30.15)') |
609 |
|
& '// CINT = ', 0. |
610 |
|
ENDIF |
611 |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
612 |
& SQUEEZE_RIGHT, 1) |
& SQUEEZE_RIGHT, 1) |
613 |
WRITE(msgBuf,'(A,1024A1)') |
WRITE(msgBuf,'(A,1024A1)') |
735 |
pltStep = sNy |
pltStep = sNy |
736 |
pltLab = 'J =' |
pltLab = 'J =' |
737 |
ENDIF |
ENDIF |
738 |
IF ( validRange ) THEN |
C IF ( validRange ) THEN |
739 |
C Header |
C Header |
740 |
C Data |
C Data |
741 |
DO bk=pltBlo, pltBhi, pltBstr |
DO bk=pltBlo, pltBhi, pltBstr |
797 |
ELSEIF ( plotMode .EQ. PRINT_MAP_YZ ) THEN |
ELSEIF ( plotMode .EQ. PRINT_MAP_YZ ) THEN |
798 |
val = fld(K,I,J,bk,bi) |
val = fld(K,I,J,bk,bi) |
799 |
ENDIF |
ENDIF |
800 |
IDX = NINT( |
IF ( validRange ) THEN |
801 |
|
IDX = NINT( |
802 |
& FLOAT( lChList-1 )*( val-fMin ) / (fRange) |
& FLOAT( lChList-1 )*( val-fMin ) / (fRange) |
803 |
& )+1 |
& )+1 |
804 |
|
ELSE |
805 |
|
IDX = 1 |
806 |
|
ENDIF |
807 |
IF ( iBuf .LE. MAX_LEN_PLOTBUF ) |
IF ( iBuf .LE. MAX_LEN_PLOTBUF ) |
808 |
& plotBuf(iBuf:iBuf) = chList(IDX:IDX) |
& plotBuf(iBuf:iBuf) = chList(IDX:IDX) |
809 |
IF ( val .EQ. 0. ) THEN |
IF ( val .EQ. 0. ) THEN |
819 |
ENDDO |
ENDDO |
820 |
ENDDO |
ENDDO |
821 |
ENDDO |
ENDDO |
822 |
ENDIF |
C ENDIF |
823 |
C-- Write delimiter |
C-- Write delimiter |
824 |
msgBuf = '// =======================================================' |
msgBuf = '// =======================================================' |
825 |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
967 |
DO K=kMin, kMax, kStr |
DO K=kMin, kMax, kStr |
968 |
DO J=jMin, jMax, jStr |
DO J=jMin, jMax, jStr |
969 |
DO I=iMin, iMax, iStr |
DO I=iMin, iMax, iStr |
970 |
C IF ( fld(I,J,K,bi,bj) .NE. 0. ) THEN |
IF ( printMapIncludesZeros .OR. fld(I,J,K,bi,bj) .NE. 0. ) THEN |
971 |
IF ( fld(I,J,K,bi,bj) .LT. fMin ) |
IF ( fld(I,J,K,bi,bj) .LT. fMin ) |
972 |
& fMin = fld(I,J,K,bi,bj) |
& fMin = fld(I,J,K,bi,bj) |
973 |
IF ( fld(I,J,K,bi,bj) .GT. fMax ) |
IF ( fld(I,J,K,bi,bj) .GT. fMax ) |
974 |
& fMax = fld(I,J,K,bi,bj) |
& fMax = fld(I,J,K,bi,bj) |
975 |
C ENDIF |
ENDIF |
976 |
ENDDO |
ENDDO |
977 |
ENDDO |
ENDDO |
978 |
ENDDO |
ENDDO |
1004 |
& '// CMAX = ', fMax |
& '// CMAX = ', fMax |
1005 |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
1006 |
& SQUEEZE_RIGHT, 1) |
& SQUEEZE_RIGHT, 1) |
1007 |
WRITE(msgBuf,'(A,1PE30.15)') |
IF ( validRange ) THEN |
1008 |
|
WRITE(msgBuf,'(A,1PE30.15)') |
1009 |
& '// CINT = ', fRange/FLOAT(lChlist-1) |
& '// CINT = ', fRange/FLOAT(lChlist-1) |
1010 |
|
ELSE |
1011 |
|
WRITE(msgBuf,'(A,1PE30.15)') |
1012 |
|
& '// CINT = ', 0. |
1013 |
|
ENDIF |
1014 |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
1015 |
& SQUEEZE_RIGHT, 1) |
& SQUEEZE_RIGHT, 1) |
1016 |
WRITE(msgBuf,'(A,1024A1)') |
WRITE(msgBuf,'(A,1024A1)') |
1138 |
pltStep = sNy |
pltStep = sNy |
1139 |
pltLab = 'J =' |
pltLab = 'J =' |
1140 |
ENDIF |
ENDIF |
1141 |
IF ( validRange ) THEN |
C IF ( validRange ) THEN |
1142 |
C Header |
C Header |
1143 |
C Data |
C Data |
1144 |
DO bk=pltBlo, pltBhi, pltBstr |
DO bk=pltBlo, pltBhi, pltBstr |
1198 |
ELSEIF ( plotMode .EQ. PRINT_MAP_YZ ) THEN |
ELSEIF ( plotMode .EQ. PRINT_MAP_YZ ) THEN |
1199 |
val = fld(K,I,J,bk,bi) |
val = fld(K,I,J,bk,bi) |
1200 |
ENDIF |
ENDIF |
1201 |
IDX = NINT( |
IF ( validRange ) THEN |
1202 |
& FLOAT( lChList-1 )*( val-fMin ) / (fRange) |
IDX = NINT( |
1203 |
& )+1 |
& FLOAT( lChList-1 )*( val-fMin ) / (fRange) |
1204 |
|
& )+1 |
1205 |
|
ELSE |
1206 |
|
IDX = 1 |
1207 |
|
ENDIF |
1208 |
IF ( iBuf .LE. MAX_LEN_PLOTBUF ) |
IF ( iBuf .LE. MAX_LEN_PLOTBUF ) |
1209 |
& plotBuf(iBuf:iBuf) = chList(IDX:IDX) |
& plotBuf(iBuf:iBuf) = chList(IDX:IDX) |
1210 |
IF ( val .EQ. 0. ) THEN |
IF ( val .EQ. 0. ) THEN |
1219 |
ENDDO |
ENDDO |
1220 |
ENDDO |
ENDDO |
1221 |
ENDDO |
ENDDO |
1222 |
ENDIF |
C ENDIF |
1223 |
C-- Write delimiter |
C-- Write delimiter |
1224 |
msgBuf = '// =======================================================' |
msgBuf = '// =======================================================' |
1225 |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |