/[MITgcm]/MITgcm/model/src/write_state.F
ViewVC logotype

Diff of /MITgcm/model/src/write_state.F

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

revision 1.1 by adcroft, Wed May 5 18:32:35 1999 UTC revision 1.17 by adcroft, Mon Mar 4 17:26:41 2002 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5  C--  File read_write.F: Routines to handle mid-level I/O interface.  #undef  MULTIPLE_RECORD_STATE_FILES
6  C--   Contents  
7  C--   o WRITE_STATE - Write out model state.  CBOP
8    C     !ROUTINE: WRITE_STATE
9  CStartofinterface  C     !INTERFACE:
10        SUBROUTINE WRITE_STATE ( forceOutput, myCurrentTime,        SUBROUTINE WRITE_STATE ( myTime, myIter, myThid )
11       &    myIter, myThid )  C     !DESCRIPTION: \bv
12  C     /==========================================================\  C     *==========================================================*
13  C     | SUBROUTINE WRITE_STATE                                   |  C     | SUBROUTINE WRITE_STATE                                    
14  C     | o Controlling routine for IO to dump model state.        |  C     | o Controlling routine for IO to dump model state.        
15  C     |==========================================================|  C     *==========================================================*
16  C     | Write model state files for post-processing. This file   |  C     | Write model state files for post-processing. This file  
17  C     | includes code for diagnosing W and RHO for output.       |  C     | includes code for diagnosing W and RHO for output.      
18  C     \==========================================================/  C     | File write_state.F: Routines to handle mid-level I/O
19        IMPLICIT NONE  C     | interface.
20    C     | o WRITE_STATE - Write out model state.
21    C     | The following CPP flag (MULTIPLE_RECORD_STATE_FILES) is
22    C     | #define/#undefed here since it is specific to this routine
23    C     | and very user-preference specific.
24    C     |
25    C     | If #undefed (default) the state files are written as in all versions
26    C     | prior to checkpoint32, where a file is created per variable, per time
27    C     | and per tile. This *has* to be the default because most users use this
28    C     | mode and all utilities and scripts (diagnostic) assume this form.
29    C     | It is also robust, as explained below.
30    C     |
31    C     | If #defined, subsequent snap-shots are written as records in the
32    C     | same file (no iteration number in filenames).
33    C     | Advantages: - fewer files
34    C     |       - for small problems, is easy to copy the output around
35    C     | Disadvantages:
36    C     |       - breaks a lot of diagnostic scripts
37    C     |       - for large or long problems this creates huge files
38    C     |       - is an unexpected, unsolicited change in behaviour which came
39    C     |         as a surprise (in c32) and inconvenience to several users
40    C     |       - can not accomodate changing the frequency of output
41    C     |         after a pickup (this is trivial in previous method
42    C     |         but needs new code and parameters in this new method)
43    C     | Known Bugs:
44    C     |       - if the length of integration is not exactly an integer
45    C     |         times the output frequency then the last record written
46    C     |         (at end of integration) overwrites a previously written
47    C     |         record corresponding to an earier time. *BE WARNED*
48    C     *==========================================================*
49    C     \ev
50    
51    C     !USES:
52          IMPLICIT NONE
53  C     == Global variables ===  C     == Global variables ===
54  #include "SIZE.h"  #include "SIZE.h"
55  #include "EEPARAMS.h"  #include "EEPARAMS.h"
56  #include "PARAMS.h"  #include "PARAMS.h"
57  #include "DYNVARS.h"  #include "DYNVARS.h"
58  #include "CG2D.h"  #include "GRID.h"
59  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_PASSIVE_TRACER
60  #include "CG3D.h"  #include "TR1.h"
 #include "GW.h"  
61  #endif  #endif
   
62        LOGICAL  DIFFERENT_MULTIPLE        LOGICAL  DIFFERENT_MULTIPLE
63        EXTERNAL DIFFERENT_MULTIPLE        EXTERNAL DIFFERENT_MULTIPLE
64        INTEGER  IO_ERRCOUNT        INTEGER  IO_ERRCOUNT
65        EXTERNAL IO_ERRCOUNT        EXTERNAL IO_ERRCOUNT
66    
67    C     !INPUT/OUTPUT PARAMETERS:
68  C     == Routine arguments ==  C     == Routine arguments ==
69  C     myThid - Thread number for this instance of the routine.  C     myThid - Thread number for this instance of the routine.
70  C     myIter - Iteration number  C     myIter - Iteration number
71  C     myCurrentTime - Current time of simulation ( s )  C     myTime - Current time of simulation ( s )
72        LOGICAL forceOutput        _RL     myTime
       REAL    myCurrentTime  
73        INTEGER myThid        INTEGER myThid
74        INTEGER myIter        INTEGER myIter
 CEndofinterface  
75    
76    C     !LOCAL VARIABLES:
77  C     == Local variables ==  C     == Local variables ==
 C     suff - Hold suffix part of a filename  
 C     beginIOErrCount - Begin and end IO error counts  
 C     endIOErrCount  
 C     msgBuf - Error message buffer  
       INTEGER beginIOErrCount  
       INTEGER endIOErrCount  
       CHARACTER*(MAX_LEN_MBUF) msgBuf  
78        CHARACTER*(MAX_LEN_MBUF) suff        CHARACTER*(MAX_LEN_MBUF) suff
79        CHARACTER*(80) fName        INTEGER iRec
80    CEOP
81    
82        IF (        IF (
83       &  ( DIFFERENT_MULTIPLE(dumpFreq,myCurrentTime,       &     DIFFERENT_MULTIPLE(dumpFreq,myTime,
84       &    myCurrentTime-deltaTClock) .AND. myCurrentTime.NE.startTime )       &     myTime-deltaTClock) .OR.
85       &  .OR. forceOutput       &     myTime.EQ.endTime .OR.
86       & ) THEN       &     myTime.EQ.startTime
87         &     ) THEN
88  C--    Going to really do some IO. Make everyone except master thread wait.           iRec = int ( (myTime-startTime) / dumpFreq +1.5 )
89         _BARRIER          
90         _BEGIN_MASTER( myThid )  C--   Going to really do some IO. Make everyone except master thread wait.
91             _BARRIER
92          WRITE(suff,'(I10.10)') myIter           _BEGIN_MASTER( myThid )
93    
94  C--     Write model fields  C--     Write model fields
95    
96    #ifdef MULTIPLE_RECORD_STATE_FILES
97    
98    C       Write each snap-shot as a new record in one file per variable
99    C       - creates relatively few files but these files can become huge
100            CALL WRITE_REC_XYZ_RL( 'U',  uVel,iRec,myIter,myThid)
101            CALL WRITE_REC_XYZ_RL( 'V',  vVel,iRec,myIter,myThid)
102            CALL WRITE_REC_XYZ_RL( 'T', theta,iRec,myIter,myThid)
103            CALL WRITE_REC_XYZ_RL( 'S',  salt,iRec,myIter,myThid)
104            CALL WRITE_REC_XY_RL('Eta',etaN,iRec,myIter,myThid)
105            CALL WRITE_REC_XYZ_RL( 'W',wVel,iRec,myIter,myThid)
106    #ifdef ALLOW_PASSIVE_TRACER
107            CALL WRITE_REC_XYZ_RL( 'TR1',tr1,iRec,myIter,myThid)
108    #endif
109    #ifdef ALLOW_NONHYDROSTATIC
110            IF (nonHydroStatic) THEN
111             CALL WRITE_REC_XYZ_RL( 'PNH',phi_nh,iRec,myIter,myThid)
112            ENDIF
113    #endif /* ALLOW_NONHYDROSTATIC */
114    #ifdef NONLIN_FRSURF
115            CALL WRITE_REC_XYZ_RS('hFacC.',hFacC,iRec,myIter,myThid)
116            CALL WRITE_FLD_XYZ_RS('hFacW.',hFacW,iRec,myIter,myThid)
117            CALL WRITE_FLD_XYZ_RS('hFacS.',hFacS,iRec,myIter,myThid)
118    #endif /* NONLIN_FRSURF */
119    
120    #else /* MULTIPLE_RECORD_STATE_FILES */
121    
122    C       Write each snap-shot as a new file (original and default method)
123    C       - creates many files but for large configurations is easier to
124    C         transfer analyse a particular snap-shots
125            WRITE(suff,'(I10.10)') myIter
126          CALL WRITE_FLD_XYZ_RL( 'U.',suff,uVel,myIter,myThid)          CALL WRITE_FLD_XYZ_RL( 'U.',suff,uVel,myIter,myThid)
127          CALL WRITE_FLD_XYZ_RL( 'V.',suff,vVel,myIter,myThid)          CALL WRITE_FLD_XYZ_RL( 'V.',suff,vVel,myIter,myThid)
128          CALL WRITE_FLD_XYZ_RL( 'T.',suff,theta,myIter,myThid)          CALL WRITE_FLD_XYZ_RL( 'T.',suff,theta,myIter,myThid)
129          CALL WRITE_FLD_XYZ_RL( 'S.',suff,salt,myIter,myThid)          CALL WRITE_FLD_XYZ_RL( 'S.',suff,salt,myIter,myThid)
130          CALL WRITE_FLD_XY_RL( 'PS.',suff,cg2d_x,myIter,myThid)          CALL WRITE_FLD_XY_RL('Eta.',suff,etaN,myIter,myThid)
131            CALL WRITE_FLD_XYZ_RL( 'W.',suff,wVel,myIter,myThid)
132    #ifdef ALLOW_PASSIVE_TRACER
133            CALL WRITE_FLD_XYZ_RL( 'TR1.',suff,tr1,myIter,myThid)
134    #endif
135  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
136          IF (nonHydroStatic) THEN          IF (nonHydroStatic) THEN
137           CALL WRITE_FLD_XYZ_RL( 'PNH.',suff,cg3d_x,myIter,myThid)           CALL WRITE_FLD_XYZ_RL( 'PNH.',suff,phi_nh,myIter,myThid)
138          ENDIF          ENDIF
139          CALL WRITE_FLD_XYZ_RL( 'W.',suff,wVel,myIter,myThid)  #endif /* ALLOW_NONHYDROSTATIC */
140  #endif  #ifdef NONLIN_FRSURF
141            CALL WRITE_FLD_XYZ_RS('hFacC.',suff,hFacC,myIter,myThid)
142         _END_MASTER( myThid )          CALL WRITE_FLD_XYZ_RS('hFacW.',suff,hFacW,myIter,myThid)
143         _BARRIER          CALL WRITE_FLD_XYZ_RS('hFacS.',suff,hFacS,myIter,myThid)
144    #endif /* NONLIN_FRSURF */
145    
146    #endif /* MULTIPLE_RECORD_STATE_FILES */
147    
148    #ifdef ALLOW_PTRACERS
149           IF (usePTRACERS) CALL PTRACERS_WRITE_STATE(myIter,myTime,myThid)
150    #endif /* ALLOW_PTRACERS */
151             _END_MASTER( myThid )
152             _BARRIER
153    
154        ENDIF        ENDIF
155    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22