/[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.30 by edhill, Thu Feb 5 00:13:47 2004 UTC revision 1.34 by adcroft, Mon Mar 8 16:47:25 2004 UTC
# Line 81  C     == Local variables == Line 81  C     == Local variables ==
81        CHARACTER*(MAX_LEN_MBUF) suff        CHARACTER*(MAX_LEN_MBUF) suff
82        INTEGER iRec        INTEGER iRec
83  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
84        _RL mnc_d_iter        _RL mnc_iter
85  #endif  #endif
86  CEOP  CEOP
87    
# Line 168  c       CALL WRITE_FLD_XYZ_RS('hFacS.',s Line 168  c       CALL WRITE_FLD_XYZ_RS('hFacS.',s
168           _BARRIER           _BARRIER
169    
170  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
171  C     Test the MNC package  
172        mnc_d_iter = myIter  C     Write dynvars using the MNC package
173        write(*,*) 'mnc_d_iter, myIter = ', mnc_d_iter, myIter        mnc_iter = myIter
174        CALL MNC_CW_RL_W_R(myThid,'state',0,0,'iter',-1,mnc_d_iter)        CALL MNC_CW_RL_W_D(myThid,'state',0,0,'iter',-1,mnc_iter)
175        CALL MNC_CW_RL_W_D(myThid,'state',0,0,'U',0,    uVel)        CALL MNC_CW_RL_W_D(myThid,'state',0,0,'model_time',0,myTime)
176  #endif        CALL MNC_CW_RL_W_D(myThid,'state',0,0,'U', 0, uVel)
177          CALL MNC_CW_RL_W_D(myThid,'state',0,0,'V', 0, vVel)
178          CALL MNC_CW_RL_W_D(myThid,'state',0,0,'T', 0, theta)
179          CALL MNC_CW_RL_W_D(myThid,'state',0,0,'S', 0, salt)
180          CALL MNC_CW_RL_W_D(myThid,'state',0,0,'Eta', 0, etaN)
181          CALL MNC_CW_RL_W_D(myThid,'state',0,0,'W', 0, wVel)
182          IF ( useDynP_inEos_Zc .OR. myIter.NE.nIter0 ) THEN
183            CALL MNC_CW_RL_W_R(myThid,'totPhiHyd',0,0,'iter',-1,mnc_iter)
184            CALL MNC_CW_RL_W_D(myThid,'totPhiHyd',0,0,'totPhiHyd',
185         &       0,totPhiHyd)
186          ENDIF
187          IF ((buoyancyRelation .NE. 'ATMOSPHERIC')
188         &     .AND. (myIter .NE. nIter0)) THEN
189            CALL MNC_CW_RL_W_D(myThid,'phiHydLow',0,0,'iter',-1,mnc_iter)
190            CALL MNC_CW_RL_W_D(myThid,'phiHydLow',0,0,'phiHydLow',
191         &       0,phiHydLow)
192          ENDIF
193    
194    #ifdef ALLOW_PASSIVE_TRACER
195          CALL MNC_CW_RL_W_D(myThid,'state',0,0,'tr1', 0, tr1)
196    #endif /*  ALLOW_PASSIVE_TRACER  */
197    
198    #ifdef ALLOW_NONHYDROSTATIC
199          IF (nonHydroStatic) THEN
200            CALL MNC_CW_RL_W_D(myThid,'state',0,0,'phi_nh', 0, phi_nh)
201          ENDIF
202    #endif /* ALLOW_NONHYDROSTATIC */
203    
204    #endif /*  ALLOW_MNC  */
205    
206        ENDIF        ENDIF
207    
Legend:
Removed from v.1.30  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.22