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-|--+----| |