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

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

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

revision 1.24 by jmc, Wed Oct 8 18:55:29 2003 UTC revision 1.25 by edhill, Thu Oct 9 04:19:18 2003 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  C--  File read_write.F: Routines to handle mid-level I/O interface.  C--  File read_write.F: Routines to handle mid-level I/O interface.
# Line 97  C       Raw fields Line 98  C       Raw fields
98          CALL READ_REC_XYZ_RL( 'gsNm1',  gSNm1, 1,myIter, myThid)          CALL READ_REC_XYZ_RL( 'gsNm1',  gSNm1, 1,myIter, myThid)
99          CALL READ_REC_XY_RL ('etaN', etaN, 1,myIter, myThid)          CALL READ_REC_XY_RL ('etaN', etaN, 1,myIter, myThid)
100  #ifdef INCLUDE_CD_CODE  #ifdef INCLUDE_CD_CODE
101            IF (useCDscheme) THEN
102          CALL READ_REC_XY_RL ('etaNm1', etaNm1, 1,myIter, myThid)          CALL READ_REC_XY_RL ('etaNm1', etaNm1, 1,myIter, myThid)
103          CALL READ_REC_XYZ_RL(  'uVelD', uVelD, 1,myIter, myThid)          CALL READ_REC_XYZ_RL(  'uVelD', uVelD, 1,myIter, myThid)
104          CALL READ_REC_XYZ_RL(  'vVelD', vVelD, 1,myIter, myThid)          CALL READ_REC_XYZ_RL(  'vVelD', vVelD, 1,myIter, myThid)
# Line 104  C       Raw fields Line 106  C       Raw fields
106          CALL READ_REC_XYZ_RL(  'vNm1',  vNM1, 1,myIter, myThid)          CALL READ_REC_XYZ_RL(  'vNm1',  vNM1, 1,myIter, myThid)
107  c       CALL READ_REC_XYZ_RL(  'guCD',  guCD, 1,myIter, myThid)  c       CALL READ_REC_XYZ_RL(  'guCD',  guCD, 1,myIter, myThid)
108  c       CALL READ_REC_XYZ_RL(  'gvCD',  gvCD, 1,myIter, myThid)  c       CALL READ_REC_XYZ_RL(  'gvCD',  gvCD, 1,myIter, myThid)
109            ENDIF
110  #endif  #endif
111    
112  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
# Line 140  C--     Read model fields Line 143  C--     Read model fields
143          CALL MDSREADFIELD(fn,prec,'RL',Nr,gSnm1, 12,myThid)          CALL MDSREADFIELD(fn,prec,'RL',Nr,gSnm1, 12,myThid)
144          CALL MDSREADFIELD(fn,prec,'RL', 1,etaN,12*Nr+1,myThid)          CALL MDSREADFIELD(fn,prec,'RL', 1,etaN,12*Nr+1,myThid)
145  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
146          IF ( nonlinFreeSurf.GE.0)          IF ( nonlinFreeSurf.GT.0)
147       &  CALL MDSREADFIELD(fn,prec,'RL',1,etaH,12*Nr+2,myThid)       &  CALL MDSREADFIELD(fn,prec,'RL',1,etaH,12*Nr+2,myThid)
148  #endif  #endif
149    
# Line 150  C--     Read model fields Line 153  C--     Read model fields
153          ENDIF          ENDIF
154    
155  #ifdef INCLUDE_CD_CODE  #ifdef INCLUDE_CD_CODE
156            IF (useCDscheme) THEN
157          WRITE(fn,'(A,A10)') 'pickup_cd.',suff          WRITE(fn,'(A,A10)') 'pickup_cd.',suff
158          CALL MDSREADFIELD(fn,prec,'RL',Nr,uVelD,    1,myThid)          CALL MDSREADFIELD(fn,prec,'RL',Nr,uVelD,    1,myThid)
159          CALL MDSREADFIELD(fn,prec,'RL',Nr,vVelD,    2,myThid)          CALL MDSREADFIELD(fn,prec,'RL',Nr,vVelD,    2,myThid)
# Line 158  C--     Read model fields Line 162  C--     Read model fields
162  c       CALL MDSREADFIELD(fn,prec,'RL',Nr,guCD,     5,myThid)  c       CALL MDSREADFIELD(fn,prec,'RL',Nr,guCD,     5,myThid)
163  c       CALL MDSREADFIELD(fn,prec,'RL',Nr,gvCD,     6,myThid)  c       CALL MDSREADFIELD(fn,prec,'RL',Nr,gvCD,     6,myThid)
164          CALL MDSREADFIELD(fn,prec,'RL', 1,etaNm1,6*Nr+1,myThid)          CALL MDSREADFIELD(fn,prec,'RL', 1,etaNm1,6*Nr+1,myThid)
165            ENDIF
166  #endif /* INCLUDE_CD_CODE */  #endif /* INCLUDE_CD_CODE */
167    
168  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
# Line 216  c     _EXCH_XYZ_R8(gvNM1 , myThid ) Line 221  c     _EXCH_XYZ_R8(gvNM1 , myThid )
221       & _EXCH_XYZ_RL( totPhiHyd, myThid )       & _EXCH_XYZ_RL( totPhiHyd, myThid )
222    
223  #ifdef INCLUDE_CD_CODE  #ifdef INCLUDE_CD_CODE
224           IF (useCDscheme) THEN
225  c**** CALL EXCH_DUV_XYZ_RL(uVelD,vVelD,.TRUE.,myThid)  c**** CALL EXCH_DUV_XYZ_RL(uVelD,vVelD,.TRUE.,myThid)
226  c**** CALL EXCH_DUV_XYZ_RL(guCD,gvCD,.TRUE.,myThid)  c**** CALL EXCH_DUV_XYZ_RL(guCD,gvCD,.TRUE.,myThid)
227        _EXCH_XYZ_R8( uVelD,    myThid )        _EXCH_XYZ_R8( uVelD,    myThid )
# Line 226  c     _EXCH_XYZ_R8( vNM1,     myThid ) Line 232  c     _EXCH_XYZ_R8( vNM1,     myThid )
232  c     _EXCH_XYZ_R8( guCD,     myThid )  c     _EXCH_XYZ_R8( guCD,     myThid )
233  c     _EXCH_XYZ_R8( gvCD,     myThid )  c     _EXCH_XYZ_R8( gvCD,     myThid )
234        _EXCH_XY_R8( etaNm1, myThid )        _EXCH_XY_R8( etaNm1, myThid )
235  #endif        ENDIF
236    #endif /* INCLUDE_CD_CODE */
237  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
238          IF ( nonHydrostatic ) THEN          IF ( nonHydrostatic ) THEN
239           _EXCH_XYZ_R8(phi_nh, myThid )           _EXCH_XYZ_R8(phi_nh, myThid )
# Line 358  C       Raw fields Line 365  C       Raw fields
365          CALL WRITE_REC_XYZ_RL( 'gSNm1',  gSNm1,  1,myIter, myThid)          CALL WRITE_REC_XYZ_RL( 'gSNm1',  gSNm1,  1,myIter, myThid)
366          CALL WRITE_REC_XY_RL ('etaN', etaN,  1,myIter, myThid)          CALL WRITE_REC_XY_RL ('etaN', etaN,  1,myIter, myThid)
367  #ifdef INCLUDE_CD_CODE  #ifdef INCLUDE_CD_CODE
368            IF (useCDscheme) THEN
369          CALL WRITE_REC_XY_RL          CALL WRITE_REC_XY_RL
370       &   ( 'etaNm1', etaNm1,  1,myIter, myThid)       &   ( 'etaNm1', etaNm1,  1,myIter, myThid)
371          CALL WRITE_REC_XYZ_RL(  'uVelD', uVelD,  1,myIter, myThid)          CALL WRITE_REC_XYZ_RL(  'uVelD', uVelD,  1,myIter, myThid)
# Line 366  C       Raw fields Line 374  C       Raw fields
374          CALL WRITE_REC_XYZ_RL(  'vNM1',  vNM1,  1,myIter, myThid)          CALL WRITE_REC_XYZ_RL(  'vNM1',  vNM1,  1,myIter, myThid)
375  c       CALL WRITE_REC_XYZ_RL(  'guCD',  guCD,  1,myIter, myThid)  c       CALL WRITE_REC_XYZ_RL(  'guCD',  guCD,  1,myIter, myThid)
376  c       CALL WRITE_REC_XYZ_RL(  'gvCD',  gvCD,  1,myIter, myThid)  c       CALL WRITE_REC_XYZ_RL(  'gvCD',  gvCD,  1,myIter, myThid)
377  #endif          ENDIF
378    #endif /* INCLUDE_CD_CODE */
379    
380    
381  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
# Line 418  C--     Write model fields Line 426  C--     Write model fields
426       &                      1,myIter,myThid)       &                      1,myIter,myThid)
427          ENDIF          ENDIF
428  #ifdef INCLUDE_CD_CODE  #ifdef INCLUDE_CD_CODE
429            IF (useCDscheme) THEN
430          IF ( permCheckPoint ) THEN          IF ( permCheckPoint ) THEN
431           WRITE(fn,'(A,I10.10)') 'pickup_cd.',myIter           WRITE(fn,'(A,I10.10)') 'pickup_cd.',myIter
432          ELSE          ELSE
# Line 435  c       CALL MDSWRITEFIELD(fn,prec,lgf,' Line 444  c       CALL MDSWRITEFIELD(fn,prec,lgf,'
444          CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,vNM1, 6,myIter,myThid)          CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,vNM1, 6,myIter,myThid)
445          CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,etaNm1,6*Nr+1,          CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,etaNm1,6*Nr+1,
446       &                     myIter,myThid)       &                     myIter,myThid)
447            ENDIF
448  #endif /* INCLUDE_CD_CODE */  #endif /* INCLUDE_CD_CODE */
449  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
450          IF ( nonHydrostatic ) THEN          IF ( nonHydrostatic ) THEN
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22