/[MITgcm]/MITgcm/eesupp/src/print.F
ViewVC logotype

Diff of /MITgcm/eesupp/src/print.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.10 by adcroft, Wed Jul 15 22:16:10 1998 UTC revision 1.11 by cnh, Sat Sep 5 17:52:13 1998 UTC
# Line 564  C--   Calculate field range Line 564  C--   Calculate field range
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 )
# Line 601  C--   Write field title and statistics Line 601  C--   Write field title and statistics
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)')
# Line 730  C      X across, Z down slice Line 735  C      X across, Z down slice
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
# Line 792  C      Data Line 797  C      Data
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
# Line 810  C      Data Line 819  C      Data
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,
# Line 958  C--   Calculate field range Line 967  C--   Calculate field range
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
# Line 995  C--   Write field title and statistics Line 1004  C--   Write field title and statistics
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)')
# Line 1124  C      X across, Z down slice Line 1138  C      X across, Z down slice
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
# Line 1184  C      Data Line 1198  C      Data
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
# Line 1201  C      Data Line 1219  C      Data
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,
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22