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

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

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

revision 1.8 by jmc, Sun Jun 14 21:45:12 2009 UTC revision 1.9 by jmc, Fri Dec 11 13:53:07 2009 UTC
# Line 195  C--   Salinity: Line 195  C--   Salinity:
195          ENDIF          ENDIF
196         ENDIF         ENDIF
197        ENDIF        ENDIF
198    #ifdef ALLOW_NONHYDROSTATIC
199    C--   W velocity:
200          IF ( nonHydrostatic ) THEN
201           IF ( alph_AB.NE.0. .OR. beta_AB.NE.0. ) THEN
202            j = j + 1
203            CALL WRITE_REC_3D_RL( fn, fp, Nr, gwNm(1-Olx,1-Oly,1,1,1,m1),
204         &                                            -j, myIter, myThid )
205            IF (j.LE.listDim) wrFldList(j) = 'GwNm1   '
206           ENDIF
207           IF ( beta_AB.NE.0. ) THEN
208            j = j + 1
209            CALL WRITE_REC_3D_RL( fn, fp, Nr, gwNm(1-Olx,1-Oly,1,1,1,m2),
210         &                                            -j, myIter, myThid )
211            IF (j.LE.listDim) wrFldList(j) = 'GwNm2   '
212           ENDIF
213          ENDIF
214    #endif /* ALLOW_NONHYDROSTATIC */
215  #else /*  ALLOW_ADAMSBASHFORTH_3 */  #else /*  ALLOW_ADAMSBASHFORTH_3 */
216         IF ( momStepping ) THEN         IF ( momStepping ) THEN
217          j = j + 1          j = j + 1
# Line 214  C--   Salinity: Line 231  C--   Salinity:
231          CALL WRITE_REC_3D_RL( fn, fp, Nr, gsNm1,  -j, myIter, myThid )          CALL WRITE_REC_3D_RL( fn, fp, Nr, gsNm1,  -j, myIter, myThid )
232          IF (j.LE.listDim) wrFldList(j) = 'GsNm1   '          IF (j.LE.listDim) wrFldList(j) = 'GsNm1   '
233         ENDIF         ENDIF
234    #ifdef ALLOW_NONHYDROSTATIC
235           IF ( nonHydrostatic ) THEN
236            j = j + 1
237            CALL WRITE_REC_3D_RL( fn, fp, Nr, gwNm1,  -j, myIter, myThid )
238            IF (j.LE.listDim) wrFldList(j) = 'GwNm1   '
239           ENDIF
240    #endif /* ALLOW_NONHYDROSTATIC */
241  #endif /*  ALLOW_ADAMSBASHFORTH_3 */  #endif /*  ALLOW_ADAMSBASHFORTH_3 */
242    
243  C-    write Full Pressure for EOS in pressure:  C-    write Full Pressure for EOS in pressure:
# Line 228  C-    write Full Pressure for EOS in pre Line 252  C-    write Full Pressure for EOS in pre
252          CALL WRITE_REC_3D_RL( fn, fp, Nr, phi_nh, -j, myIter, myThid )          CALL WRITE_REC_3D_RL( fn, fp, Nr, phi_nh, -j, myIter, myThid )
253          IF (j.LE.listDim) wrFldList(j) = 'Phi_NHyd'          IF (j.LE.listDim) wrFldList(j) = 'Phi_NHyd'
254         ENDIF         ENDIF
        IF ( nonHydrostatic ) THEN  
         j = j + 1  
         CALL WRITE_REC_3D_RL( fn, fp, Nr, gwNm1,  -j, myIter, myThid )  
         IF (j.LE.listDim) wrFldList(j) = 'GwNm1   '  
        ENDIF  
255  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
256  #ifdef ALLOW_ADDFLUID  #ifdef ALLOW_ADDFLUID
257  C-    write mass source/sink of fluid (but not needed if selectAddFluid=-1)  C-    write mass source/sink of fluid (but not needed if selectAddFluid=-1)
# Line 249  C---  Write 2-D fields, starting with Et Line 268  C---  Write 2-D fields, starting with Et
268          nj = -( n3D*(Nr-1) + j )          nj = -( n3D*(Nr-1) + j )
269          CALL WRITE_REC_3D_RL( fn, fp, 1 , etaN,   nj, myIter, myThid )          CALL WRITE_REC_3D_RL( fn, fp, 1 , etaN,   nj, myIter, myThid )
270          IF (j.LE.listDim) wrFldList(j) = 'EtaN    '          IF (j.LE.listDim) wrFldList(j) = 'EtaN    '
271    #ifdef ALLOW_NONHYDROSTATIC
272           IF ( selectNHfreeSurf.GE.1 ) THEN
273            j = j + 1
274            nj = -( n3D*(Nr-1) + j )
275            CALL WRITE_REC_3D_RL( fn, fp, 1, dPhiNH,  nj, myIter, myThid )
276            IF (j.LE.listDim) wrFldList(j) = 'dPhiNH  '
277           ENDIF
278    #endif /* ALLOW_NONHYDROSTATIC */
279  #ifdef EXACT_CONSERV  #ifdef EXACT_CONSERV
280  c      IF ( exactConserv ) THEN  c      IF ( exactConserv ) THEN
281          j = j + 1          j = j + 1
# Line 328  C       Then set the actual unlimited di Line 355  C       Then set the actual unlimited di
355          IF ( use3Dsolver ) THEN          IF ( use3Dsolver ) THEN
356            CALL MNC_CW_RL_W('D',fn,0,0,'phi_nh', phi_nh, myThid)            CALL MNC_CW_RL_W('D',fn,0,0,'phi_nh', phi_nh, myThid)
357  c         CALL MNC_CW_RL_W('D',fn,0,0,'gW', gW, myThid)  c         CALL MNC_CW_RL_W('D',fn,0,0,'gW', gW, myThid)
358    #ifndef ALLOW_ADAMSBASHFORTH_3
359            CALL MNC_CW_RL_W('D',fn,0,0,'gWnm1', gwNm1, myThid)            CALL MNC_CW_RL_W('D',fn,0,0,'gWnm1', gwNm1, myThid)
360    #endif
361          ENDIF          ENDIF
362  #endif  #endif
363          IF ( useDynP_inEos_Zc ) THEN          IF ( useDynP_inEos_Zc ) THEN
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22