/[MITgcm]/MITgcm/pkg/mnc/mnc_cw_write_grid_info.F
ViewVC logotype

Diff of /MITgcm/pkg/mnc/mnc_cw_write_grid_info.F

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

revision 1.12 by jmc, Fri Aug 28 19:48:07 2009 UTC revision 1.13 by mlosch, Tue Nov 24 14:53:50 2009 UTC
# Line 159  C     Write the GRID.h variables Line 159  C     Write the GRID.h variables
159        END        END
160    
161  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
162    
163    CBOP
164    C !ROUTINE: MNC_CW_WRITE_GRID_COORD
165    
166    C !INTERFACE:
167          SUBROUTINE MNC_CW_WRITE_GRID_COORD(
168         I     fname,
169         I     myThid )
170    
171    C     !DESCRIPTION:
172    C     This subroutine writes coordinates to one netcdf file
173    C     (\texttt{fname.xxxxx.nc}) per tile.
174    
175    C     !USES:
176          implicit none
177    #include "SIZE.h"
178    #include "GRID.h"
179    
180    C     !INPUT PARAMETERS:
181          integer myThid
182          character*(*) fname
183    CEOP
184    C     Write coordinates to file "fname"
185          _BEGIN_MASTER( myThid )
186    
187          CALL MNC_CW_SET_UDIM(fname, 0, myThid)
188          CALL MNC_CW_RS_W('R',fname,0,0,'XC',xC, myThid)
189          CALL MNC_CW_RS_W('R',fname,0,0,'YC',yC, myThid)
190          CALL MNC_CW_RS_W('R',fname,0,0,'XG',xG, myThid)
191          CALL MNC_CW_RS_W('R',fname,0,0,'YG',yG, myThid)
192    C--   Here we can add as much grid information as we like, but this
193    C--   will increase the size of each of the affected files. It is
194    C--   a matter of taste, but here are some candidates:
195    C--   Probably it is enough to write the coordinates of the
196    C--   center points and the corners.
197    C     CALL MNC_CW_RS_W('R',fname,0,0,'XU',xG, myThid)
198    C     CALL MNC_CW_RS_W('R',fname,0,0,'YU',yC, myThid)
199    C     CALL MNC_CW_RS_W('R',fname,0,0,'XV',xC, myThid)
200    C     CALL MNC_CW_RS_W('R',fname,0,0,'YV',yG, myThid)
201    C--   Only required for diagnostics files, but unclear what the vertical
202    C--   coordinate means in this case, so do not write anything for now
203    C     CALL MNC_CW_RS_W('R',fname,0,0,'RC',rC, myThid)
204    C     CALL MNC_CW_RS_W('R',fname,0,0,'RF',rF, myThid)
205    C--   If required write mask information.
206    C     CALL MNC_CW_RS_W('R',fname,0,0,'HFacC',HFacC,myThid)
207    
208          _END_MASTER(myThid)
209    
210          RETURN
211          END
212    
213    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13

  ViewVC Help
Powered by ViewVC 1.1.22