/[MITgcm]/MITgcm/pkg/diagnostics/diagstats_mnc_out.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagstats_mnc_out.F

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

revision 1.5 by edhill, Thu Jul 21 06:17:09 2005 UTC revision 1.7 by jmc, Tue Nov 1 01:53:13 2005 UTC
# Line 49  C     myThid   :: my thread Id number Line 49  C     myThid   :: my thread Id number
49  CEOP  CEOP
50    
51  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
52    #ifdef ALLOW_MNC
53        INTEGER im, ix, iv, ist        INTEGER im, ix, iv, ist
54        PARAMETER ( iv = nStats - 2 , im = nStats - 1 , ix = nStats )        PARAMETER ( iv = nStats - 2 , im = nStats - 1 , ix = nStats )
55        INTEGER i, j, k        INTEGER i, j, k
# Line 56  C     !LOCAL VARIABLES: Line 57  C     !LOCAL VARIABLES:
57        CHARACTER*(3) stat_typ(5)        CHARACTER*(3) stat_typ(5)
58        INTEGER ILNBLNK        INTEGER ILNBLNK
59        EXTERNAL ILNBLNK        EXTERNAL ILNBLNK
 #ifdef ALLOW_MNC  
60        INTEGER ii, ilen        INTEGER ii, ilen
61        CHARACTER*(MAX_LEN_FNAM) diag_mnc_bn        CHARACTER*(MAX_LEN_FNAM) diag_mnc_bn
       CHARACTER*(5) ctmp  
62        INTEGER CW_DIMS, NLEN        INTEGER CW_DIMS, NLEN
63        PARAMETER ( CW_DIMS = 10 )        PARAMETER ( CW_DIMS = 10 )
64        PARAMETER ( NLEN    = 80 )        PARAMETER ( NLEN    = 80 )
# Line 68  C     !LOCAL VARIABLES: Line 67  C     !LOCAL VARIABLES:
67        CHARACTER*(NLEN) d_cw_gname        CHARACTER*(NLEN) d_cw_gname
68        CHARACTER*(NLEN) d_cw_gname0        CHARACTER*(NLEN) d_cw_gname0
69        CHARACTER*(NLEN) dn_blnk        CHARACTER*(NLEN) dn_blnk
70    #ifdef DIAGST_MNC_NEEDSWORK
71          CHARACTER*(5) ctmp
72        _RS ztmp(Nr+Nrphys)        _RS ztmp(Nr+Nrphys)
73    #endif
74        _RL stmp(Nr+Nrphys+1,nRegions+1)        _RL stmp(Nr+Nrphys+1,nRegions+1)
75  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
76    
# Line 257  C           Copy the data into a tempora Line 259  C           Copy the data into a tempora
259                stmp(1,j+1) = statGlob(ist,0,j)                stmp(1,j+1) = statGlob(ist,0,j)
260              ENDDO              ENDDO
261                        
262              IF ((fflags(listId)(1:1) .EQ. ' ')  C-jmc: fflags is not for Statistics-Diagnostics, can be unset, and since
263       &           .OR. (fflags(listId)(1:1) .EQ. 'R')) THEN  C-     size of the output file will not be an issue here: Always write real*8
264                  c           IF ((fflags(listId)(1:1) .EQ. ' ')
265                CALL MNC_CW_RL_W('R',diag_mnc_bn,1,1,  c    &           .OR. (fflags(listId)(1:1) .EQ. 'R')) THEN
266       &             tnam, stmp, myThid)  c            
267                  c             CALL MNC_CW_RL_W('R',diag_mnc_bn,1,1,
268              ELSEIF (fflags(listId)(1:1) .EQ. 'D') THEN  c    &             tnam, stmp, myThid)
269    c            
270    c           ELSEIF (fflags(listId)(1:1) .EQ. 'D') THEN
271                                
272                CALL MNC_CW_RL_W('D',diag_mnc_bn,1,1,                CALL MNC_CW_RL_W('D',diag_mnc_bn,1,1,
273       &             tnam, stmp, myThid)       &             tnam, stmp, myThid)
274                            
275              ENDIF  c           else
276    c             write(0,*) myIter, ndId, listId
277    c             write(0,'(3A)') '>',cdiag(ndId),'<'
278    c             write(0,'(3A)') '>',fflags(listId),'<'
279    c             STOP ' in DIAGSTATS_MNC_OUT'
280    c           ENDIF
281                            
282              CALL MNC_CW_DEL_VNAME(tnam, myThid)              CALL MNC_CW_DEL_VNAME(tnam, myThid)
283                            
# Line 294  C           Copy the data into a tempora Line 303  C           Copy the data into a tempora
303                ENDDO                ENDDO
304              ENDDO              ENDDO
305                        
306              IF ((fflags(listId)(1:1) .EQ. ' ')  C-jmc: Always write real*8 (size of the output file will not be an issue here)
      &           .OR. (fflags(listId)(1:1) .EQ. 'R')) THEN  
               
               CALL MNC_CW_RL_W('R',diag_mnc_bn,1,1,  
      &             tnam, stmp, myThid)  
               
             ELSEIF (fflags(listId)(1:1) .EQ. 'D') THEN  
               
307                CALL MNC_CW_RL_W('D',diag_mnc_bn,1,1,                CALL MNC_CW_RL_W('D',diag_mnc_bn,1,1,
308       &             tnam, stmp, myThid)       &             tnam, stmp, myThid)
309                            
             ENDIF  
   
310              CALL MNC_CW_DEL_VNAME(tnam, myThid)              CALL MNC_CW_DEL_VNAME(tnam, myThid)
311                        
312            ENDIF            ENDIF
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22