/[MITgcm]/MITgcm/pkg/seaice/seaice_pickup.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_pickup.F

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

revision 1.1 by heimbach, Mon Nov 11 22:01:21 2002 UTC revision 1.2 by heimbach, Tue Nov 12 20:47:27 2002 UTC
# Line 0  Line 1 
1    C
2    
3    #include "SEAICE_OPTIONS.h"
4    
5    C--   File seaice_pickup.F: Routines to handle reading and writing
6    C--   sea ice pickup files.
7    C--   Contents
8    C--   o SEAICE_READ_PICKUP  - Read sea ice pickup file for restarting.
9    C--   o SEAICE_WRITE_PICKUP - Write out sea ice pickup file.
10    
11    CStartofinterface
12          SUBROUTINE SEAICE_READ_PICKUP ( myThid )
13    C     /==========================================================\
14    C     | SUBROUTINE SEAICE_READ_PICKUP                            |
15    C     | o Read in sea ice pickup file for restarting.            |
16    C     |==========================================================|
17    C     \==========================================================/
18          IMPLICIT NONE
19    
20    C     == Global variables ===
21    #include "SIZE.h"
22    #include "EEPARAMS.h"
23    #include "PARAMS.h"
24    #include "SEAICE.h"
25    
26    C     == Routine arguments ==
27    C     myThid - Thread number for this instance of the routine.
28          INTEGER myThid
29    CEndofinterface
30    
31    #ifdef ALLOW_SEAICE
32    
33    C     == Local variables ==
34          CHARACTER*(MAX_LEN_FNAM) fn
35    
36    C--    Going to do some IO. Make everyone except master thread wait.
37           _BARRIER
38           _BEGIN_MASTER( myThid )
39    
40    C--    Read ice model fields
41           WRITE(fn,'(A,I10.10)') 'pickup_seaice.',nIter0
42           CALL MDSREADFIELD( fn, precFloat64, 'RL', 1, TICE , 1, myThid )
43           CALL MDSREADFIELD( fn, precFloat64, 'RL', 1, YNEG , 2, myThid )
44           CALL MDSREADFIELD( fn, precFloat64, 'RL', 1, HSNOW, 3, myThid )
45           CALL MDSREADFIELD( fn, precFloat64, 'RL', 3, UICE , 2, myThid )
46           CALL MDSREADFIELD( fn, precFloat64, 'RL', 3, VICE , 3, myThid )
47           CALL MDSREADFIELD( fn, precFloat64, 'RL', 3, HEFF , 4, myThid )
48           CALL MDSREADFIELD( fn, precFloat64, 'RL', 3, AREA , 5, myThid )
49    
50           _END_MASTER( myThid )
51           _BARRIER
52    
53    C--    Update overlap regions
54           CALL EXCH_RL( UICE, OLx, OLx, OLy, OLy, 3, OLx, OLy,
55         I      FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
56           CALL EXCH_RL( VICE, OLx, OLx, OLy, OLy, 3, OLx, OLy,
57         I      FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
58           CALL EXCH_RL( HEFF, OLx, OLx, OLy, OLy, 3, OLx, OLy,
59         I      FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
60           CALL EXCH_RL( AREA, OLx, OLx, OLy, OLy, 3, OLx, OLy,
61         I      FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid )
62           _EXCH_XY_R4(TICE , myThid )
63           _EXCH_XY_R4(YNEG , myThid )
64           _EXCH_XY_R4(HSNOW, myThid )
65    
66    #endif ALLOW_SEAICE
67    
68          RETURN
69          END
70    
71    CStartofinterface
72          SUBROUTINE SEAICE_WRITE_PICKUP ( modelEnd, myCurrentTime,
73         &                              myIt, myThid )
74    C     /==========================================================\
75    C     | SUBROUTINE SEAICE_WRITE_PICKUP                           |
76    C     | o Write sea ice pickup file for restarting.              |
77    C     |==========================================================|
78    C     \==========================================================/
79          IMPLICIT NONE
80    
81    C     == Global variables ===
82    #include "SIZE.h"
83    #include "EEPARAMS.h"
84    #include "PARAMS.h"
85    #include "SEAICE.h"
86    
87          LOGICAL  DIFFERENT_MULTIPLE
88          EXTERNAL DIFFERENT_MULTIPLE
89          INTEGER  IO_ERRCOUNT
90          EXTERNAL IO_ERRCOUNT
91    
92    C     == Routine arguments ==
93    C     modelEnd - Checkpoint call at end of model run.
94    C     myCurrentTime - Current time of simulation ( s )
95    C     myIt     - Iteration number
96    C     myThid   - Thread number for this instance of the routine.
97          LOGICAL modelEnd    
98          INTEGER myThid
99          INTEGER myIt
100          _RL     myCurrentTime
101    CEndofinterface
102    
103    #ifdef ALLOW_SEAICE
104    
105    C     == Common blocks ==
106          COMMON /PCKP_GBLFLS/ globalFile
107          LOGICAL globalFile
108    
109    C     == Local variables ==
110    C     permCheckPoint - Flag indicating whether a permanent checkpoint
111    C                      will be written.
112          LOGICAL permCheckPoint  
113          CHARACTER*(MAX_LEN_FNAM) fn
114          INTEGER prec
115          LOGICAL lgf
116    
117          permCheckPoint = .FALSE.
118          permCheckPoint=
119         & DIFFERENT_MULTIPLE(pChkptFreq,myCurrentTime,
120         &                    myCurrentTime-deltaTClock)
121    
122          IF (
123         &    (.NOT. modelEnd .AND. (
124         &     permCheckPoint
125         &     .OR.
126         &     DIFFERENT_MULTIPLE(chkptFreq,
127         &    myCurrentTime,myCurrentTime-deltaTClock)
128         &     )
129         &    )
130         &     .OR.
131         &    (
132         &     modelEnd
133         &     .AND. .NOT.
134         &     permCheckPoint
135         &     .AND. .NOT.
136         &     DIFFERENT_MULTIPLE(chkptFreq,
137         &    myCurrentTime,myCurrentTime-deltaTClock)
138         &    )
139         & ) THEN
140    
141    C--    Going to do some IO. Make everyone except master thread wait.
142           _BARRIER
143           _BEGIN_MASTER( myThid )
144    
145           prec = precFloat64
146           lgf = globalFile
147    
148    C--    Write model fields
149           IF ( permCheckPoint ) THEN
150            WRITE(fn,'(A,I10.10)') 'pickup_seaice.',myIt
151           ELSE
152            WRITE(fn,'(A,A)') 'pickup_seaice.',checkPtSuff(nCheckLev)
153           ENDIF
154           CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,  TICE, 1, myIt,myThid)
155           CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,  YNEG, 2, myIt,myThid)
156           CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1, HSNOW, 3, myIt,myThid)
157           CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 3,  UICE, 2, myIt,myThid)
158           CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 3,  VICE, 3, myIt,myThid)
159           CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 3,  HEFF, 4, myIt,myThid)
160           CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 3,  AREA, 5, myIt,myThid)
161    
162           _END_MASTER( myThid )
163           _BARRIER
164    
165          ENDIF
166    
167    #endif ALLOW_SEAICE
168    
169          RETURN
170          END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22