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

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

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

revision 1.45 by jmc, Tue Apr 28 18:01:14 2009 UTC revision 1.46 by jmc, Sun Jun 14 21:45:12 2009 UTC
# Line 85  C------ Line 85  C------
85        ENDDO        ENDDO
86    
87  C-    Need to synchronize here before doing master-thread IO  C-    Need to synchronize here before doing master-thread IO
88        _BARRIER  C     this is done within IO routines => no longer needed
89    c     _BARRIER
90    
91  C------  C------
92  C   1) Set R_low = the Lower (in r sense) boundary of the fluid column :  C   1) Set R_low = the Lower (in r sense) boundary of the fluid column :
# Line 174  C------ read from file: Line 175  C------ read from file:
175  C- read surface topography (in m) from topoFile (case topoFile.NE.' '):  C- read surface topography (in m) from topoFile (case topoFile.NE.' '):
176    
177          CALL READ_REC_XY_RS( topoFile, topoZ, 1, 0, myThid )          CALL READ_REC_XY_RS( topoFile, topoZ, 1, 0, myThid )
         _BARRIER  
178    
179          IF (buoyancyRelation .EQ. 'ATMOSPHERIC') THEN          IF (buoyancyRelation .EQ. 'ATMOSPHERIC') THEN
180  C----  C----
# Line 187  c         CALL WRITE_FLD_XY_RS( 'topo_Z' Line 187  c         CALL WRITE_FLD_XY_RS( 'topo_Z'
187       O                       Ro_surf,       O                       Ro_surf,
188       I                       myThid )       I                       myThid )
189    
 c         _BARRIER  
190  C         This I/O is now done in write_grid.F  C         This I/O is now done in write_grid.F
191  c         CALL WRITE_FLD_XY_RS( 'topo_P',' ',Ro_surf,0,myThid)  c         CALL WRITE_FLD_XY_RS( 'topo_P',' ',Ro_surf,0,myThid)
192    
# Line 262  c          IF ( jG .EQ. Ny )  R_low(i,j, Line 261  c          IF ( jG .EQ. Ny )  R_low(i,j,
261        ENDIF        ENDIF
262    
263        IF ( debugLevel.GE.debLevB ) THEN        IF ( debugLevel.GE.debLevB ) THEN
264           _BARRIER
265         CALL PLOT_FIELD_XYRS( Ro_surf,         CALL PLOT_FIELD_XYRS( Ro_surf,
266       &           'Surface reference r-position (ini_depths)',       &           'Surface reference r-position (ini_depths)',
267       &           -1, myThid )       &           -1, myThid )
# Line 269  c          IF ( jG .EQ. Ny )  R_low(i,j, Line 269  c          IF ( jG .EQ. Ny )  R_low(i,j,
269  c     CALL WRITE_FLD_XY_RS('Ro_surf',' ',Ro_surf,0,myThid)  c     CALL WRITE_FLD_XY_RS('Ro_surf',' ',Ro_surf,0,myThid)
270    
271  C--   Everyone else must wait for the depth to be loaded  C--   Everyone else must wait for the depth to be loaded
272  C-    note: might not be necessary since all single-thread IO above  C-    note: not necessary since all single-thread IO above are followed
273  C           are followed by an EXCH (with BARRIER)  C           by an EXCH (with BARRIER) + BARRIER within IO routine
274        _BARRIER  c     _BARRIER
275    
276  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
277        IF ( useOBCS ) THEN        IF ( useOBCS ) THEN
Legend:
Removed from v.1.45  
changed lines
  Added in v.1.46
  ViewVC Help
Powered by ViewVC 1.1.22