/[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.13 by cnh, Fri Nov 6 22:44:42 1998 UTC revision 1.19 by adcroft, Fri Sep 21 14:31:12 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_EEOPTIONS.h"  #include "CPP_EEOPTIONS.h"
5    
# Line 18  C--    o print_mapr8    Formats ABCD... Line 19  C--    o print_mapr8    Formats ABCD...
19  C--                     Uses print_message for writing  C--                     Uses print_message for writing
20  C--    o print_message  Does IO with unhighlighted header  C--    o print_message  Does IO with unhighlighted header
21    
22  CStartOfInterface  CBOP              
23    
24    C     !ROUTINE: PRINT_ERROR
25    
26    C     !INTERFACE:
27        SUBROUTINE PRINT_ERROR( message , myThid )        SUBROUTINE PRINT_ERROR( message , myThid )
28  C     /============================================================\        IMPLICIT NONE
29  C     | SUBROUTINE PRINT_ERROR                                     |  
30  C     | o Write out error message using "standard" format.         |  C     !DESCRIPTION:
31  C     | Notes                                                      |  C     *============================================================*
32  C     | =====                                                      |  C     | SUBROUTINE PRINT_ERROR                                      
33  C     | o Some system   I/O is not "thread-safe". For this reason  |  C     | o Write out error message using "standard" format.          
34  C     |   without the FMTFTN_IO_THREAD_SAFE directive set a        |  C     *============================================================*
35  C     |   critical region is defined around the write here. In some|  C     | Notes                                                      
36  C     |   cases  BEGIN_CRIT() is approximated by only doing writes |  C     | =====                                                      
37  C     |   for thread number 1 - writes for other threads are       |  C     | o Some system   I/O is not "thread-safe". For this reason  
38  C     |   ignored!                                                 |  C     |   without the FMTFTN_IO_THREAD_SAFE directive set a        
39  C     | o In a non-parallel form these routines can still be used. |  C     |   critical region is defined around the write here. In some
40  C     |   to produce pretty printed output!                        |  C     |   cases  BEGIN_CRIT() is approximated by only doing writes  
41  C     \============================================================/  C     |   for thread number 1 - writes for other threads are        
42    C     |   ignored!                                                  
43    C     | o In a non-parallel form these routines are still used  
44    C     |   to produce pretty printed output. The process and thread
45    C     |   id prefix is omitted in this case.
46    C     *============================================================*
47    
48    C     !USES:
49  C     == Global data ==  C     == Global data ==
50  #include "SIZE.h"  #include "SIZE.h"
51  #include "EEPARAMS.h"  #include "EEPARAMS.h"
52  #include "EESUPPORT.h"  #include "EESUPPORT.h"
 C     == Routine arguments ==  
       CHARACTER*(*) message  
       INTEGER       myThid  
 CEndOfInterface  
53        INTEGER  IFNBLNK        INTEGER  IFNBLNK
54        EXTERNAL IFNBLNK        EXTERNAL IFNBLNK
55        INTEGER  ILNBLNK        INTEGER  ILNBLNK
56        EXTERNAL ILNBLNK        EXTERNAL ILNBLNK
57    
58    C     !INPUT/OUTPUT PARAMETERS:
59    C     == Routine arguments ==
60    C     message :: Text string to print
61    C     myThid  :: Thread number of this instance
62          CHARACTER*(*) message
63          INTEGER       myThid
64    
65    C     !LOCAL VARIABLES:
66  C     == Local variables ==  C     == Local variables ==
67    C     iStart, iEnd :: Temps. for string indexing
68    C     idString     :: Temp. for building message prefix
69        INTEGER iStart        INTEGER iStart
70        INTEGER iEnd        INTEGER iEnd
71        CHARACTER*9 idString        CHARACTER*9 idString
72    CEOP
73    
74  C--   Find beginning and end of message  C--   Find beginning and end of message
75        iStart = IFNBLNK( message )        iStart = IFNBLNK( message )
76        iEnd   = ILNBLNK( message )        iEnd   = ILNBLNK( message )
# Line 105  C Line 126  C
126        GOTO 1000        GOTO 1000
127        END        END
128    
129  CStartofinterface  CBOP
130    C     !ROUTINE: PRINT_LIST_I
131    
132    C     !INTERFACE:
133        SUBROUTINE PRINT_LIST_I( fld, lFld, index_type,        SUBROUTINE PRINT_LIST_I( fld, lFld, index_type,
134       &                         markEnd, compact, ioUnit )       &                         markEnd, compact, ioUnit )
135  C     /==========================================================\        IMPLICIT NONE
136  C     | o SUBROUTINE PRINT_LIST_I                                |  C     !DESCRIPTION:
137  C     |==========================================================|  C     *==========================================================*
138  C     | Routine for producing list of values for a field with    |  C     | o SUBROUTINE PRINT_LIST_I                                
139  C     | duplicate values collected into                          |  C     *==========================================================*
140  C     |    n @ value                                             |  C     | Routine for producing list of values for a field with    
141  C     | record.                                                  |  C     | duplicate values collected into                          
142  C     \==========================================================/  C     |    n @ value                                              
143    C     | record.                                                  
144    C     *==========================================================*
145    
146    C     !USES:
147  C     == Global data ==    C     == Global data ==  
148  #include "SIZE.h"  #include "SIZE.h"
149  #include "EEPARAMS.h"  #include "EEPARAMS.h"
150    
151    C     !INPUT/OUTPUT PARAMETERS:
152  C     == Routine arguments ==  C     == Routine arguments ==
153  C     fld    -  Data to be printed  C     fld    ::  Data to be printed
154  C     lFld   -  Number of elements to be printed  C     lFld   ::  Number of elements to be printed
155  C     index_type - Flag indicating which type of index to print  C     index_type :: Flag indicating which type of index to print
156  C                  INDEX_K    => /* K = nnn */  C                   INDEX_K    => /* K = nnn */
157  C                  INDEX_I    => /* I = nnn */  C                   INDEX_I    => /* I = nnn */
158  C                  INDEX_J    => /* J = nnn */  C                   INDEX_J    => /* J = nnn */
159  C                  INDEX_NONE =>  C                   INDEX_NONE =>
160  C     compact -  Flag to control use of repeat symbol for same valued  C     compact ::  Flag to control use of repeat symbol for same valued
161  C                fields.  C                 fields.
162  C     markEnd -  Flag to control whether there is a separator after the  C     markEnd ::  Flag to control whether there is a separator after the
163  C                last element  C                 last element
164  C     ioUnit -  Unit number for IO.  C     ioUnit ::   Unit number for IO.
165        INTEGER lFld        INTEGER lFld
166        INTEGER index_type        INTEGER index_type
167        INTEGER fld(lFld)        INTEGER fld(lFld)
168        LOGICAL markEnd        LOGICAL markEnd
169        LOGICAL compact        LOGICAL compact
170        INTEGER ioUnit        INTEGER ioUnit
 CEndifinterface  
171    
172    C     !LOCAL VARIABLES:
173  C     == Local variables ==  C     == Local variables ==
174  C     iLo  - Range index holders for selecting elements with  C     iLo  - Range index holders for selecting elements with
175  C     iHi    with the same value  C     iHi    with the same value
# Line 160  C     K    - Loop counter Line 188  C     K    - Loop counter
188        CHARACTER*2 commOpen,commClose        CHARACTER*2 commOpen,commClose
189        CHARACTER*3 index_lab        CHARACTER*3 index_lab
190        INTEGER K        INTEGER K
191    CEOP
192    
193        IF     ( index_type .EQ. INDEX_I ) THEN        IF     ( index_type .EQ. INDEX_I ) THEN
194         index_lab = 'I ='         index_lab = 'I ='
# Line 178  C     K    - Loop counter Line 207  C     K    - Loop counter
207        xOld = fld(1)        xOld = fld(1)
208        DO K=2,lFld        DO K=2,lFld
209         xNew = fld(K  )         xNew = fld(K  )
210         IF ( .NOT. compact .OR. xNew .NE. xOld ) THEN         IF ( .NOT. compact .OR. (xNew .NE. xOld) ) THEN
211          nDup = iHi-iLo+1          nDup = iHi-iLo+1
212          IF ( nDup .EQ. 1 ) THEN          IF ( nDup .EQ. 1 ) THEN
213           WRITE(msgBuf,'(A,I9,A)') '              ',xOld,punc           WRITE(msgBuf,'(A,I9,A)') '              ',xOld,punc
# Line 186  C     K    - Loop counter Line 215  C     K    - Loop counter
215       &    WRITE(msgBuf(45:),'(A,1X,A,I3,1X,A)')       &    WRITE(msgBuf(45:),'(A,1X,A,I3,1X,A)')
216       &    commOpen,index_lab,iLo,commClose       &    commOpen,index_lab,iLo,commClose
217          ELSE          ELSE
218           WRITE(msgBuf,'(I,'' '',A,I9,A)') nDup,'@',xOld,punc           WRITE(msgBuf,'(I5,'' '',A,I9,A)') nDup,'@',xOld,punc
219           IF ( index_type .NE. INDEX_NONE )           IF ( index_type .NE. INDEX_NONE )
220       &    WRITE(msgBuf(45:),'(A,1X,A,I3,A,I3,1X,A)')       &    WRITE(msgBuf(45:),'(A,1X,A,I3,A,I3,1X,A)')
221       &    commOpen,index_lab,iLo,':',iHi,commClose       &    commOpen,index_lab,iLo,':',iHi,commClose
# Line 208  C     K    - Loop counter Line 237  C     K    - Loop counter
237       &  WRITE(msgBuf(45:),'(A,1X,A,I3,1X,A)')       &  WRITE(msgBuf(45:),'(A,1X,A,I3,1X,A)')
238       &  commOpen,index_lab,iLo,commClose       &  commOpen,index_lab,iLo,commClose
239        ELSEIF( nDup .GT. 1 ) THEN        ELSEIF( nDup .GT. 1 ) THEN
240         WRITE(msgBuf,'(I,'' '',A,I9,A)') nDup,'@',xOld,punc         WRITE(msgBuf,'(I5,'' '',A,I9,A)') nDup,'@',xOld,punc
241         IF ( index_type .NE. INDEX_NONE )         IF ( index_type .NE. INDEX_NONE )
242       &  WRITE(msgBuf(45:),'(A,1X,A,I3,A,I3,1X,A)')       &  WRITE(msgBuf(45:),'(A,1X,A,I3,A,I3,1X,A)')
243       &  commOpen,index_lab,iLo,':',iHi,commClose       &  commOpen,index_lab,iLo,':',iHi,commClose
# Line 218  C     K    - Loop counter Line 247  C     K    - Loop counter
247        RETURN        RETURN
248        END        END
249    
250  CStartofinterface  CBOP
251    C     !ROUTINE: PRINT_LIST_L
252    
253    C     !INTERFACE:
254        SUBROUTINE PRINT_LIST_L( fld, lFld, index_type, markEnd,        SUBROUTINE PRINT_LIST_L( fld, lFld, index_type, markEnd,
255       &                         compact, ioUnit )       &                         compact, ioUnit )
256  C     /==========================================================\        IMPLICIT NONE
257  C     | o SUBROUTINE PRINT_LIST_L                                |  C     !DESCRIPTION:
258  C     |==========================================================|  C     *==========================================================*
259  C     | Routine for producing list of values for a field with    |  C     | o SUBROUTINE PRINT_LIST_L                                
260  C     | duplicate values collected into                          |  C     *==========================================================*
261  C     |    n @ value                                             |  C     | Routine for producing list of values for a field with    
262  C     | record.                                                  |  C     | duplicate values collected into                          
263  C     \==========================================================/  C     |    n @ value                                              
264    C     | record.                                                  
265    C     *==========================================================*
266    
267    C     !USES:
268  C     == Global data ==    C     == Global data ==  
269  #include "SIZE.h"  #include "SIZE.h"
270  #include "EEPARAMS.h"  #include "EEPARAMS.h"
271    
272    C     !INPUT/OUTPUT PARAMETERS:
273  C     == Routine arguments ==  C     == Routine arguments ==
274  C     fld    -  Data to be printed  C     fld    -  Data to be printed
275  C     lFld   -  Number of elements to be printed  C     lFld   -  Number of elements to be printed
# Line 253  C     ioUnit -  Unit number for IO. Line 289  C     ioUnit -  Unit number for IO.
289        LOGICAL markEnd        LOGICAL markEnd
290        LOGICAL compact        LOGICAL compact
291        INTEGER ioUnit        INTEGER ioUnit
 CEndifinterface  
292    
293    C     !LOCAL VARIABLES:
294  C     == Local variables ==  C     == Local variables ==
295  C     iLo  - Range index holders for selecting elements with  C     iLo  - Range index holders for selecting elements with
296  C     iHi    with the same value  C     iHi    with the same value
# Line 273  C     K    - Loop counter Line 309  C     K    - Loop counter
309        CHARACTER*2 commOpen,commClose        CHARACTER*2 commOpen,commClose
310        CHARACTER*3 index_lab        CHARACTER*3 index_lab
311        INTEGER K        INTEGER K
312    CEOP
313    
314        IF     ( index_type .EQ. INDEX_I ) THEN        IF     ( index_type .EQ. INDEX_I ) THEN
315         index_lab = 'I ='         index_lab = 'I ='
# Line 291  C     K    - Loop counter Line 328  C     K    - Loop counter
328        xOld = fld(1)        xOld = fld(1)
329        DO K=2,lFld        DO K=2,lFld
330         xNew = fld(K  )         xNew = fld(K  )
331         IF ( .NOT. compact .OR. xNew .NEQV. xOld ) THEN         IF ( .NOT. compact .OR. (xNew .NEQV. xOld) ) THEN
332          nDup = iHi-iLo+1          nDup = iHi-iLo+1
333          IF ( nDup .EQ. 1 ) THEN          IF ( nDup .EQ. 1 ) THEN
334           WRITE(msgBuf,'(A,L5,A)') '              ',xOld,punc           WRITE(msgBuf,'(A,L5,A)') '              ',xOld,punc
# Line 299  C     K    - Loop counter Line 336  C     K    - Loop counter
336       &    WRITE(msgBuf(45:),'(A,1X,A,I3,1X,A)')       &    WRITE(msgBuf(45:),'(A,1X,A,I3,1X,A)')
337       &    commOpen,index_lab,iLo,commClose       &    commOpen,index_lab,iLo,commClose
338          ELSE          ELSE
339           WRITE(msgBuf,'(I,'' '',A,L5,A)') nDup,'@',xOld,punc           WRITE(msgBuf,'(I5,'' '',A,L5,A)') nDup,'@',xOld,punc
340           IF ( index_type .NE. INDEX_NONE )           IF ( index_type .NE. INDEX_NONE )
341       &    WRITE(msgBuf(45:),'(A,1X,A,L3,A,I3,1X,A)')       &    WRITE(msgBuf(45:),'(A,1X,A,L3,A,I3,1X,A)')
342       &    commOpen,index_lab,iLo,':',iHi,commClose       &    commOpen,index_lab,iLo,':',iHi,commClose
# Line 321  C     K    - Loop counter Line 358  C     K    - Loop counter
358       &  WRITE(msgBuf(45:),'(A,1X,A,I3,1X,A)')       &  WRITE(msgBuf(45:),'(A,1X,A,I3,1X,A)')
359       &    commOpen,index_lab,iLo,commClose       &    commOpen,index_lab,iLo,commClose
360        ELSEIF( nDup .GT. 1 ) THEN        ELSEIF( nDup .GT. 1 ) THEN
361         WRITE(msgBuf,'(I,'' '',A,L5,A)') nDup,'@',xOld,punc         WRITE(msgBuf,'(I5,'' '',A,L5,A)') nDup,'@',xOld,punc
362         IF ( index_type .NE. INDEX_NONE )         IF ( index_type .NE. INDEX_NONE )
363       &  WRITE(msgBuf(45:),'(A,1X,A,I3,A,I3,1X,A)')       &  WRITE(msgBuf(45:),'(A,1X,A,I3,A,I3,1X,A)')
364       &  commOpen,index_lab,iLo,':',iHi,commClose       &  commOpen,index_lab,iLo,':',iHi,commClose
# Line 331  C     K    - Loop counter Line 368  C     K    - Loop counter
368        RETURN        RETURN
369        END        END
370    
371  CStartofinterface  CBOP
372    C     !ROUTINE: PRINT_LIST_R8
373    C     !INTERFACE:
374        SUBROUTINE PRINT_LIST_R8( fld, lFld, index_type,        SUBROUTINE PRINT_LIST_R8( fld, lFld, index_type,
375       &    markEnd, compact, ioUnit )       &    markEnd, compact, ioUnit )
376  C     /==========================================================\        IMPLICIT NONE
377  C     | o SUBROUTINE PRINT_LIST_R8                               |  C     !DESCRIPTION:
378  C     |==========================================================|  C     *==========================================================*
379  C     | Routine for producing list of values for a field with    |  C     | o SUBROUTINE PRINT_LIST_R8                                
380  C     | duplicate values collected into                          |  C     *==========================================================*
381  C     |    n @ value                                             |  C     | Routine for producing list of values for a field with    
382  C     | record.                                                  |  C     | duplicate values collected into                          
383  C     \==========================================================/  C     |    n @ value                                              
384    C     | record.                                                  
385    C     *==========================================================*
386    
387  C     == Global data ==    C     !USES:
388    C     == Global data ==
389  #include "SIZE.h"  #include "SIZE.h"
390  #include "EEPARAMS.h"  #include "EEPARAMS.h"
391    
392    C     !INPUT/OUTPUT PARAMETERS:
393  C     == Routine arguments ==  C     == Routine arguments ==
394  C     fld    -  Data to be printed  C     fld    -  Data to be printed
395  C     lFld   -  Number of elements to be printed  C     lFld   -  Number of elements to be printed
# Line 366  C     ioUnit -  Unit number for IO. Line 409  C     ioUnit -  Unit number for IO.
409        LOGICAL markEnd        LOGICAL markEnd
410        LOGICAL compact        LOGICAL compact
411        INTEGER ioUnit        INTEGER ioUnit
 CEndifinterface  
412    
413    C     !LOCA VARIABLES:
414  C     == Local variables ==  C     == Local variables ==
415  C     iLo  - Range index holders for selecting elements with  C     iLo  - Range index holders for selecting elements with
416  C     iHi    with the same value  C     iHi    with the same value
# Line 386  C     K    - Loop counter Line 429  C     K    - Loop counter
429        CHARACTER*2 commOpen,commClose        CHARACTER*2 commOpen,commClose
430        CHARACTER*3 index_lab        CHARACTER*3 index_lab
431        INTEGER K        INTEGER K
432    CEOP
433    
434        IF     ( index_type .EQ. INDEX_I ) THEN        IF     ( index_type .EQ. INDEX_I ) THEN
435         index_lab = 'I ='         index_lab = 'I ='
# Line 404  C     K    - Loop counter Line 448  C     K    - Loop counter
448        xOld = fld(1)        xOld = fld(1)
449        DO K=2,lFld        DO K=2,lFld
450         xNew = fld(K  )         xNew = fld(K  )
451         IF ( .NOT. compact .OR. xNew .NE. xOld ) THEN         IF ( .NOT. compact .OR. (xNew .NE. xOld) ) THEN
452          nDup = iHi-iLo+1          nDup = iHi-iLo+1
453          IF ( nDup .EQ. 1 ) THEN          IF ( nDup .EQ. 1 ) THEN
454           WRITE(msgBuf,'(A,1PE23.15,A)') '              ',xOld,punc           WRITE(msgBuf,'(A,1PE23.15,A)') '              ',xOld,punc
# Line 446  C     K    - Loop counter Line 490  C     K    - Loop counter
490        RETURN        RETURN
491        END        END
492    
493  CStartOfInterface  CBOP
494    C     !ROUTINE: PRINT_MAPRS
495    C     !INTERFACE:
496        SUBROUTINE PRINT_MAPRS ( fld, fldTitle, plotMode,        SUBROUTINE PRINT_MAPRS ( fld, fldTitle, plotMode,
497       I        iLo,   iHi,   jLo,   jHi,  kLo,  kHi, nBx, nBy,       I        iLo,   iHi,   jLo,   jHi,  kLo,  kHi, nBx, nBy,
498       I       iMin,  iMax,  iStr,       I       iMin,  iMax,  iStr,
# Line 454  CStartOfInterface Line 500  CStartOfInterface
500       I       kMin, kMax,   kStr,       I       kMin, kMax,   kStr,
501       I      bxMin, bxMax,  bxStr,       I      bxMin, bxMax,  bxStr,
502       I      byMin, byMax,  byStr )       I      byMin, byMax,  byStr )
503  C     /==========================================================\        IMPLICIT NONE
504  C     | SUBROUTINE PRINT_MAPR4                                   |  C     !DESCRIPTION:
505  C     | o Does textual mapping printing of a field.              |  C     *==========================================================*
506  C     |==========================================================|  C     | SUBROUTINE PRINT_MAPR4                                    
507  C     | This routine does the actual formatting of the data      |  C     | o Does textual mapping printing of a field.              
508  C     | and printing to a file. It assumes an array using the    |  C     *==========================================================*
509  C     | MITgcm UV indexing scheme and base index variables.      |  C     | This routine does the actual formatting of the data      
510  C     | User code should call an interface routine like          |  C     | and printing to a file. It assumes an array using the    
511  C     | PLOT_FIELD_XYR4( ... ) rather than this code directly.   |  C     | MITgcm UV indexing scheme and base index variables.      
512  C     | Text plots can be oriented XY, YZ, XZ. An orientation    |  C     | User code should call an interface routine like          
513  C     | is specficied through the "plotMode" argument. All the   |  C     | PLOT_FIELD_XYR4( ... ) rather than this code directly.    
514  C     | plots made by a single call to this routine will use the |  C     | Text plots can be oriented XY, YZ, XZ. An orientation    
515  C     | same contour interval. The plot range (iMin,...,byStr)   |  C     | is specficied through the "plotMode" argument. All the    
516  C     | can be three-dimensional. A separate plot is made for    |  C     | plots made by a single call to this routine will use the  
517  C     | each point in the plot range normal to the orientation.  |  C     | same contour interval. The plot range (iMin,...,byStr)    
518  C     | e.g. if the orientation is XY (plotMode = PRINT_MAP_XY). |  C     | can be three-dimensional. A separate plot is made for    
519  C     |      kMin =1, kMax = 5 and kStr = 2 will produce three XY|  C     | each point in the plot range normal to the orientation.  
520  C     |      plots - one for K=1, one for K=3 and one for K=5.   |  C     | e.g. if the orientation is XY (plotMode = PRINT_MAP_XY).  
521  C     |      Each plot would have extents iMin:iMax step iStr    |  C     |      kMin =1, kMax = 5 and kStr = 2 will produce three XY
522  C     |      and jMin:jMax step jStr.                            |  C     |      plots - one for K=1, one for K=3 and one for K=5.    
523  C     \==========================================================/  C     |      Each plot would have extents iMin:iMax step iStr    
524    C     |      and jMin:jMax step jStr.                            
525    C     *==========================================================*
526    
527    C     !USES:
528  C     == Global variables ==  C     == Global variables ==
529  #include "SIZE.h"  #include "SIZE.h"
530  #include "EEPARAMS.h"  #include "EEPARAMS.h"
531  #include "EESUPPORT.h"  #include "EESUPPORT.h"
532          INTEGER  IFNBLNK
533          EXTERNAL IFNBLNK
534          INTEGER  ILNBLNK
535          EXTERNAL ILNBLNK
536    
537    C     !INPUT/OUTPUT PARAMETERS:
538  C     == Routine arguments ==  C     == Routine arguments ==
539  C     fld        - Real*4 array holding data to be plotted  C     fld        - Real*4 array holding data to be plotted
540  C     fldTitle   - Name of field to be plotted  C     fldTitle   - Name of field to be plotted
# Line 508  C     kStr Line 562  C     kStr
562        INTEGER kMin, kMax, kStr        INTEGER kMin, kMax, kStr
563        INTEGER bxMin, bxMax, bxStr        INTEGER bxMin, bxMax, bxStr
564        INTEGER byMin, byMax, byStr        INTEGER byMin, byMax, byStr
 CEndOfInterface  
 C     == Local variables ==  
       INTEGER  IFNBLNK  
       EXTERNAL IFNBLNK  
       INTEGER  ILNBLNK  
       EXTERNAL ILNBLNK  
565    
566    C     !LOCAL VARIABLES:
567  C     == Local variables ==  C     == Local variables ==
568  C     plotBuf - Buffer for building plot record  C     plotBuf - Buffer for building plot record
569  C     chList  - Character string used for plot  C     chList  - Character string used for plot
# Line 563  C               Str  - stride within blo Line 612  C               Str  - stride within blo
612        INTEGER I, J, K, iStrngLo, iStrngHi, iBuf, iDx        INTEGER I, J, K, iStrngLo, iStrngHi, iBuf, iDx
613        INTEGER bi, bj, bk        INTEGER bi, bj, bk
614        LOGICAL validRange        LOGICAL validRange
615    CEOP
616    
617        chList = '-abcdefghijklmnopqrstuvwxyz+'        chList = '-abcdefghijklmnopqrstuvwxyz+'
618        small  =  1. _d -15        small  =  1. _d -15
# Line 854  C--   Write delimiter Line 904  C--   Write delimiter
904        RETURN        RETURN
905        END        END
906    
907  CStartOfInterface  CBOP
908    C     !ROUTINE: PRINT_MAPRL
909    
910    C     !INTERFACE:
911        SUBROUTINE PRINT_MAPRL ( fld, fldTitle, plotMode,        SUBROUTINE PRINT_MAPRL ( fld, fldTitle, plotMode,
912       I        iLo,   iHi,   jLo,   jHi,  kLo,  kHi, nBx, nBy,       I        iLo,   iHi,   jLo,   jHi,  kLo,  kHi, nBx, nBy,
913       I       iMin,  iMax,  iStr,       I       iMin,  iMax,  iStr,
# Line 862  CStartOfInterface Line 915  CStartOfInterface
915       I       kMin, kMax,   kStr,       I       kMin, kMax,   kStr,
916       I      bxMin, bxMax,  bxStr,       I      bxMin, bxMax,  bxStr,
917       I      byMin, byMax,  byStr )       I      byMin, byMax,  byStr )
918  C     /==========================================================\        IMPLICIT NONE
 C     | SUBROUTINE PRINT_MAPRL                                   |  
 C     | o Does textual mapping printing of a field.              |  
 C     |==========================================================|  
 C     | This routine does the actual formatting of the data      |  
 C     | and printing to a file. It assumes an array using the    |  
 C     | MITgcm UV indexing scheme and base index variables.      |  
 C     | User code should call an interface routine like          |  
 C     | PLOT_FIELD_XYR8( ... ) rather than this code directly.   |  
 C     | Text plots can be oriented XY, YZ, XZ. An orientation    |  
 C     | is specficied through the "plotMode" argument. All the   |  
 C     | plots made by a single call to this routine will use the |  
 C     | same contour interval. The plot range (iMin,...,byStr)   |  
 C     | can be three-dimensional. A separate plot is made for    |  
 C     | each point in the plot range normal to the orientation.  |  
 C     | e.g. if the orientation is XY (plotMode = PRINT_MAP_XY). |  
 C     |      kMin =1, kMax = 5 and kStr = 2 will produce three XY|  
 C     |      plots - one for K=1, one for K=3 and one for K=5.   |  
 C     |      Each plot would have extents iMin:iMax step iStr    |  
 C     |      and jMin:jMax step jStr.                            |  
 C     \==========================================================/  
919    
920    C     !DESCRIPTION:
921    C     *==========================================================*
922    C     | SUBROUTINE PRINT_MAPRL                                    
923    C     | o Does textual mapping printing of a field.              
924    C     *==========================================================*
925    C     | This routine does the actual formatting of the data      
926    C     | and printing to a file. It assumes an array using the    
927    C     | MITgcm UV indexing scheme and base index variables.      
928    C     | User code should call an interface routine like          
929    C     | PLOT_FIELD_XYR8( ... ) rather than this code directly.    
930    C     | Text plots can be oriented XY, YZ, XZ. An orientation    
931    C     | is specficied through the "plotMode" argument. All the    
932    C     | plots made by a single call to this routine will use the  
933    C     | same contour interval. The plot range (iMin,...,byStr)    
934    C     | can be three-dimensional. A separate plot is made for    
935    C     | each point in the plot range normal to the orientation.  
936    C     | e.g. if the orientation is XY (plotMode = PRINT_MAP_XY).  
937    C     |      kMin =1, kMax = 5 and kStr = 2 will produce three XY
938    C     |      plots - one for K=1, one for K=3 and one for K=5.    
939    C     |      Each plot would have extents iMin:iMax step iStr    
940    C     |      and jMin:jMax step jStr.                            
941    C     *==========================================================*
942    
943    C     !USES:
944  C     == Global variables ==  C     == Global variables ==
945  #include "SIZE.h"  #include "SIZE.h"
946  #include "EEPARAMS.h"  #include "EEPARAMS.h"
947  #include "EESUPPORT.h"  #include "EESUPPORT.h"
948          INTEGER  IFNBLNK
949          EXTERNAL IFNBLNK
950          INTEGER  ILNBLNK
951          EXTERNAL ILNBLNK
952    
953    C     !INPUT/OUTPUT PARAMETERS:
954  C     == Routine arguments ==  C     == Routine arguments ==
955  C     fld        - Real*8 array holding data to be plotted  C     fld        - Real*8 array holding data to be plotted
956  C     fldTitle   - Name of field to be plotted  C     fldTitle   - Name of field to be plotted
# Line 916  C     kStr Line 978  C     kStr
978        INTEGER kMin, kMax, kStr        INTEGER kMin, kMax, kStr
979        INTEGER bxMin, bxMax, bxStr        INTEGER bxMin, bxMax, bxStr
980        INTEGER byMin, byMax, byStr        INTEGER byMin, byMax, byStr
 CEndOfInterface  
 C     == Local variables ==  
       INTEGER  IFNBLNK  
       EXTERNAL IFNBLNK  
       INTEGER  ILNBLNK  
       EXTERNAL ILNBLNK  
981    
982    C     !LOCAL VARIABLES:
983  C     == Local variables ==  C     == Local variables ==
984  C     plotBuf - Buffer for building plot record  C     plotBuf - Buffer for building plot record
985  C     chList  - Character string used for plot  C     chList  - Character string used for plot
# Line 971  C               Str  - stride within blo Line 1028  C               Str  - stride within blo
1028        INTEGER I, J, K, iStrngLo, iStrngHi, iBuf, iDx        INTEGER I, J, K, iStrngLo, iStrngHi, iBuf, iDx
1029        INTEGER bi, bj, bk        INTEGER bi, bj, bk
1030        LOGICAL validRange        LOGICAL validRange
1031    CEOP
1032    
1033        chList = '-abcdefghijklmnopqrstuvwxyz+'        chList = '-abcdefghijklmnopqrstuvwxyz+'
1034        small  = 1. _d -15        small  = 1. _d -15
# Line 1260  C--   Write delimiter Line 1318  C--   Write delimiter
1318        RETURN        RETURN
1319        END        END
1320    
1321  CStartOfInterface  CBOP
1322    C     !ROUTINE: PRINT_MESSAGE
1323    
1324    C     !INTERFACE:
1325        SUBROUTINE PRINT_MESSAGE( message, unit, sq , myThid )        SUBROUTINE PRINT_MESSAGE( message, unit, sq , myThid )
1326  C     /============================================================\        IMPLICIT NONE
1327  C     | SUBROUTINE PRINT_MESSAGE                                   |  C     !DESCRIPTION:
1328  C     | o Write out informational message using "standard" format. |  C     *============================================================*
1329  C     | Notes                                                      |  C     | SUBROUTINE PRINT_MESSAGE                                    
1330  C     | =====                                                      |  C     | o Write out informational message using "standard" format.  
1331  C     | o Some system   I/O is not "thread-safe". For this reason  |  C     *============================================================*
1332  C     |   without the FMTFTN_IO_THREAD_SAFE directive set a        |  C     | Notes                                                      
1333  C     |   critical region is defined around the write here. In some|  C     | =====                                                      
1334  C     |   cases  BEGIN_CRIT() is approximated by only doing writes |  C     | o Some system   I/O is not "thread-safe". For this reason  
1335  C     |   for thread number 1 - writes for other threads are       |  C     |   without the FMTFTN_IO_THREAD_SAFE directive set a        
1336  C     |   ignored!                                                 |  C     |   critical region is defined around the write here. In some
1337  C     | o In a non-parallel form these routines can still be used. |  C     |   cases  BEGIN_CRIT() is approximated by only doing writes  
1338  C     |   to produce pretty printed output!                        |  C     |   for thread number 1 - writes for other threads are        
1339  C     \============================================================/  C     |   ignored!                                                  
1340    C     | o In a non-parallel form these routines can still be used.  
1341    C     |   to produce pretty printed output!                        
1342    C     *============================================================*
1343    
1344    C     !USES:
1345  C     == Global data ==  C     == Global data ==
1346  #include "SIZE.h"  #include "SIZE.h"
1347  #include "EEPARAMS.h"  #include "EEPARAMS.h"
1348  #include "EESUPPORT.h"  #include "EESUPPORT.h"
1349          INTEGER  IFNBLNK
1350          EXTERNAL IFNBLNK
1351          INTEGER  ILNBLNK
1352          EXTERNAL ILNBLNK
1353    
1354    C     !INPUT/OUTPUT PARAMETERS:
1355  C     == Routine arguments ==  C     == Routine arguments ==
1356  C     message - Message to write  C     message :: Message to write
1357  C     unit    - Unit number to write to  C     unit    :: Unit number to write to
1358  C     sq      - Justification option  C     sq      :: Justification option
1359        CHARACTER*(*) message        CHARACTER*(*) message
1360        INTEGER       unit        INTEGER       unit
1361        CHARACTER*(*) sq        CHARACTER*(*) sq
1362        INTEGER  myThid        INTEGER  myThid
1363  CEndOfInterface  
1364        INTEGER  IFNBLNK  C     !LOCAL VARIABLES:
       EXTERNAL IFNBLNK  
       INTEGER  ILNBLNK  
       EXTERNAL ILNBLNK  
1365  C     == Local variables ==  C     == Local variables ==
1366    C     iStart, iEnd :: String indexing variables
1367    C     idString     :: Temp. for building prefix.
1368        INTEGER iStart        INTEGER iStart
1369        INTEGER iEnd        INTEGER iEnd
1370        CHARACTER*9 idString        CHARACTER*9 idString
1371    CEOP
1372    
1373  C--   Find beginning and end of message  C--   Find beginning and end of message
1374        IF ( sq .EQ. SQUEEZE_BOTH .OR.        IF ( sq .EQ. SQUEEZE_BOTH .OR.
1375       &     sq .EQ. SQUEEZE_LEFT ) THEN       &     sq .EQ. SQUEEZE_LEFT ) THEN
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22