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

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

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

revision 1.16 by jmc, Sun Jan 13 22:46:38 2013 UTC revision 1.17 by torge, Fri May 3 19:29:24 2013 UTC
# Line 79  C     record number < 0 : a hack not to Line 79  C     record number < 0 : a hack not to
79    
80  C--   write Sea-Ice Thermodynamics State variables, starting with 3-D fields:  C--   write Sea-Ice Thermodynamics State variables, starting with 3-D fields:
81         IF ( .NOT.useThSIce ) THEN         IF ( .NOT.useThSIce ) THEN
82            j = j + 1
83            nj = nj-1
84          IF (SEAICE_multDim.GT.1) THEN          IF (SEAICE_multDim.GT.1) THEN
          j = j + 1  
          nj = nj-1  
85           CALL WRITE_REC_3D_RL(fn,fp,MULTDIM,TICES, nj, myIter, myThid )           CALL WRITE_REC_3D_RL(fn,fp,MULTDIM,TICES, nj, myIter, myThid )
86           IF (j.LE.listDim) wrFldList(j) = 'siTICES '           IF (j.LE.listDim) wrFldList(j) = 'siTICES '
87  C-    switch to 2-D fields:  C-    switch to 2-D fields:
88           nj = nj*MULTDIM  c         nj = nj*MULTDIM
89             nj = nj-MULTDIM+1
90          ELSE          ELSE
          j = j + 1  
          nj = nj-1  
91           CALL WRITE_REC_3D_RL( fn, fp,  1, TICE  , nj, myIter, myThid )           CALL WRITE_REC_3D_RL( fn, fp,  1, TICE  , nj, myIter, myThid )
92           IF (j.LE.listDim) wrFldList(j) = 'siTICE  '           IF (j.LE.listDim) wrFldList(j) = 'siTICE  '
93          ENDIF          ENDIF
94    
95  C---  continue to write 2-D fields:  C---  continue to write 2-D fields:
96          j = j + 1          j = j + 1
97            nj = nj-1
98  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
         nj = nj-nITD  
99          CALL WRITE_REC_3D_RL(fn, fp, nITD, AREAITD , nj, myIter, myThid)          CALL WRITE_REC_3D_RL(fn, fp, nITD, AREAITD , nj, myIter, myThid)
100          IF (j.LE.listDim) wrFldList(j) = 'siAREAn '          IF (j.LE.listDim) wrFldList(j) = 'siAREAn '
101            nj = nj-nITD+1
102  #else  #else
         nj = nj-1  
103          CALL WRITE_REC_3D_RL( fn, fp,  1, AREA , nj, myIter, myThid )          CALL WRITE_REC_3D_RL( fn, fp,  1, AREA , nj, myIter, myThid )
104          IF (j.LE.listDim) wrFldList(j) = 'siAREA  '          IF (j.LE.listDim) wrFldList(j) = 'siAREA  '
105  #endif  #endif
106    
107          j = j + 1          j = j + 1
108            nj = nj-1
109  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
         nj = nj-nITD  
110          CALL WRITE_REC_3D_RL(fn, fp, nITD, HEFFITD , nj, myIter, myThid)          CALL WRITE_REC_3D_RL(fn, fp, nITD, HEFFITD , nj, myIter, myThid)
111          IF (j.LE.listDim) wrFldList(j) = 'siHEFFn '          IF (j.LE.listDim) wrFldList(j) = 'siHEFFn '
112            nj = nj-nITD+1
113  #else  #else
         nj = nj-1  
114          CALL WRITE_REC_3D_RL( fn, fp,  1, HEFF , nj, myIter, myThid )          CALL WRITE_REC_3D_RL( fn, fp,  1, HEFF , nj, myIter, myThid )
115          IF (j.LE.listDim) wrFldList(j) = 'siHEFF  '          IF (j.LE.listDim) wrFldList(j) = 'siHEFF  '
116  #endif  #endif
117    
118          j = j + 1          j = j + 1
119            nj = nj-1
120  #ifdef SEAICE_ITD  #ifdef SEAICE_ITD
         nj = nj-nITD  
121          CALL WRITE_REC_3D_RL(fn, fp, nITD, HSNOWITD, nj, myIter, myThid)          CALL WRITE_REC_3D_RL(fn, fp, nITD, HSNOWITD, nj, myIter, myThid)
122          IF (j.LE.listDim) wrFldList(j) = 'siHSNOWn'          IF (j.LE.listDim) wrFldList(j) = 'siHSNOWn'
123            nj = nj-nITD+1
124  #else  #else
         nj = nj-1  
125          CALL WRITE_REC_3D_RL( fn, fp,  1, HSNOW , nj, myIter, myThid )          CALL WRITE_REC_3D_RL( fn, fp,  1, HSNOW , nj, myIter, myThid )
126          IF (j.LE.listDim) wrFldList(j) = 'siHSNOW '          IF (j.LE.listDim) wrFldList(j) = 'siHSNOW '
127  #endif  #endif
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22