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

Annotation of /MITgcm/model/src/write_grid.F

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

Revision 1.2 - (hide annotations) (download)
Fri Sep 10 12:19:30 2004 UTC (19 years, 8 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint55b_post, checkpoint55, checkpoint54f_post, checkpoint55a_post
Changes since 1.1: +152 -185 lines 
 o overhaul of IO so that we now have flags for MDSIO and/or MNC
   - all verification tests compile and run with linux_ia32_g77
   - defaults are compatible with current input files--nothing
     should change if you were not previously using MNC
   - MNC output has been added in numerous places (eg. timeave)
     but there are still a few writes not yet do-able with MNC
     (this is in progress)
   - flags now allow for either/or/both use of MDSIO and MNC and
     documentation will soon follow
   - numerous small formatting cleanups for ProTeX
1 edhill 1.2 C $Header: /u/gcmpack/MITgcm/model/src/write_grid.F,v 1.1 2004/05/13 15:40:53 adcroft Exp $
2     C $Name: $
3 adcroft 1.1
4     #include "PACKAGES_CONFIG.h"
5     #include "CPP_OPTIONS.h"
6    
7     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
8     CBOP
9     C !ROUTINE: WRITE_GRID
10     C !INTERFACE:
11     SUBROUTINE WRITE_GRID(
12     I myThid )
13    
14 edhill 1.2 C !DESCRIPTION:
15     C Writes the model geometry and grid arrays to file(s) using which
16     C ever I/O package is available/active. This includes horizontal
17     C and vertical grid and finite volume modulating factors (hFacs).
18     C Essentailly, it dumps the contents of GRID.h to disk.
19 adcroft 1.1
20     C !CALLING SEQUENCE:
21     C INITIALIZE_FIXED
22     C |
23     C |-- WRITE_GRID
24    
25     IMPLICIT NONE
26     #include "SIZE.h"
27     #include "EEPARAMS.h"
28     #include "PARAMS.h"
29 edhill 1.2 #ifdef ALLOW_MNC
30     #include "MNC_PARAMS.h"
31     #endif
32 adcroft 1.1 #include "GRID.h"
33    
34     C !INPUT/OUTPUT PARAMETERS:
35     C myThid - Number of this instances
36     INTEGER myThid
37    
38     C !LOCAL VARIABLES:
39     character*(4) bfname
40     C tmpfld - Temporary array used to compute & write Total Depth
41     C has to be in common for multi threading
42     C ?aja: not sure why the COMMON block is necessary - should ask JMC?
43     COMMON / LOCAL_INI_MASKS_ETC / tmpfld
44     _RS tmpfld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
45     INTEGER i,j,bi,bj
46     CEOP
47    
48 edhill 1.2 C Calculate the fluid thickness in R coordinates as seen by the
49     C model. This may differ from what the user specified due to
50     C partial cells, etc.
51 adcroft 1.1 DO bj = myByLo(myThid), myByHi(myThid)
52 edhill 1.2 DO bi = myBxLo(myThid), myBxHi(myThid)
53     DO j=1-Oly,sNy+Oly
54     DO i=1-Olx,sNx+Olx
55     C Total fluid column thickness (r_unit) :
56     tmpfld(i,j,bi,bj) = Ro_surf(i,j,bi,bj) - R_low(i,j,bi,bj)
57     ENDDO
58     ENDDO
59 adcroft 1.1 ENDDO
60     ENDDO
61    
62     _BEGIN_MASTER( myThid )
63 edhill 1.2
64     IF (snapshot_mdsio) THEN
65     C Write horizontal grid arrays
66     CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid)
67     CALL WRITE_FLD_XY_RS( 'YC',' ',YC,0,myThid)
68     CALL WRITE_FLD_XY_RS( 'XG',' ',XG,0,myThid)
69     CALL WRITE_FLD_XY_RS( 'YG',' ',YG,0,myThid)
70     CALL WRITE_FLD_XY_RS( 'RAC',' ',rA,0,myThid)
71     CALL WRITE_FLD_XY_RS( 'RAW',' ',rAw,0,myThid)
72     CALL WRITE_FLD_XY_RS( 'RAS',' ',rAs,0,myThid)
73     CALL WRITE_FLD_XY_RS( 'DXG',' ',DXG,0,myThid)
74     CALL WRITE_FLD_XY_RS( 'DYG',' ',DYG,0,myThid)
75     CALL WRITE_FLD_XY_RS( 'DXC',' ',DXC,0,myThid)
76     CALL WRITE_FLD_XY_RS( 'DYC',' ',DYC,0,myThid)
77     C Write 3D geometry arrays
78     CALL WRITE_FLD_XY_RS( 'Depth',' ',tmpfld,0,myThid)
79     CALL WRITE_FLD_XYZ_RS( 'hFacC',' ',hFacC,0,myThid)
80     CALL WRITE_FLD_XYZ_RS( 'hFacW',' ',hFacW,0,myThid)
81     CALL WRITE_FLD_XYZ_RS( 'hFacS',' ',hFacS,0,myThid)
82     IF (buoyancyRelation .EQ. 'ATMOSPHERIC')
83     & CALL WRITE_FLD_XY_RS( 'topo_P',' ',Ro_surf,0,myThid)
84    
85     ENDIF
86 adcroft 1.1 _END_MASTER(myThid)
87    
88     #ifdef ALLOW_MNC
89 edhill 1.2 IF (useMNC .AND. snapshot_mnc) THEN
90     bfname='grid';
91 adcroft 1.1
92 edhill 1.2 C Grid spacing, areas and hFacs
93     CALL MNC_CW_ADD_VNAME('drC', '-_-_--__C__-', 0,0, myThid)
94     CALL MNC_CW_ADD_VNAME('drF', '-_-_--__C__-', 0,0, myThid)
95     CALL MNC_CW_ADD_VNAME('saFac', '-_-_--__C__-', 0,0, myThid)
96    
97     CALL MNC_CW_ADD_VNAME('dxC', 'Cen_xy_Hn__-__-', 3,4, myThid)
98     CALL MNC_CW_ADD_VNAME('dyC', 'Cen_xy_Hn__-__-', 3,4, myThid)
99     CALL MNC_CW_ADD_VNAME('dxF', 'Cen_xy_Hn__-__-', 3,4, myThid)
100     CALL MNC_CW_ADD_VNAME('dyF', 'Cen_xy_Hn__-__-', 3,4, myThid)
101     CALL MNC_CW_ADD_VNAME('dxG', 'Cen_xy_Hn__-__-', 3,4, myThid)
102     CALL MNC_CW_ADD_VNAME('dyG', 'Cen_xy_Hn__-__-', 3,4, myThid)
103     CALL MNC_CW_ADD_VNAME('dxV', 'U_xy_Hn__-__-', 3,4, myThid)
104     CALL MNC_CW_ADD_VNAME('dyU', 'V_xy_Hn__-__-', 3,4, myThid)
105     CALL MNC_CW_ADD_VNAME('rA', 'Cen_xy_Hn__-__-', 3,4, myThid)
106     CALL MNC_CW_ADD_VNAME('rAw', 'Cen_xy_Hn__-__-', 3,4, myThid)
107     CALL MNC_CW_ADD_VNAME('rAs', 'Cen_xy_Hn__-__-', 3,4, myThid)
108     CALL MNC_CW_ADD_VNAME('rAz', 'Cen_xy_Hn__-__-', 3,4, myThid)
109     CALL MNC_CW_ADD_VNAME('fCori', 'Cen_xy_Hn__-__-', 3,4,myThid)
110     CALL MNC_CW_ADD_VNAME('fCoriG','Cen_xy_Hn__-__-', 3,4,myThid)
111     CALL MNC_CW_ADD_VNAME('R_low','Cen_xy_Hn__-__-', 3,4, myThid)
112     CALL MNC_CW_ADD_VNAME('Ro_surf','Cen_xy_Hn__-__-',3,4,myThid)
113     CALL MNC_CW_ADD_VNAME('Depth', 'Cen_xy_Hn__-__-', 3,4,myThid)
114     CALL MNC_CW_ADD_VNAME('HFacC', 'Cen_xy_Hn__C__-', 4,5,myThid)
115     CALL MNC_CW_ADD_VNAME('HFacW', 'Cen_xy_Hn__C__-', 4,5,myThid)
116     CALL MNC_CW_ADD_VNAME('HFacS', 'Cen_xy_Hn__C__-', 4,5,myThid)
117    
118     CALL MNC_CW_ADD_VATTR_TEXT('drC',1,'description',
119     & 'r cell center separation',myThid)
120     CALL MNC_CW_ADD_VATTR_TEXT('drF',1,'description',
121     & 'r cell face separation',myThid)
122     CALL MNC_CW_ADD_VATTR_TEXT('saFac',1,'description',
123     & 'shallow atmosphere factor',myThid)
124     CALL MNC_CW_ADD_VATTR_TEXT('dxC',1,'description',
125     & 'x cell center separation',myThid)
126     CALL MNC_CW_ADD_VATTR_TEXT('dyC',1,'description',
127     & 'y cell center separation',myThid)
128     CALL MNC_CW_ADD_VATTR_TEXT('dxF',1,'description',
129     & 'x cell face separation',myThid)
130     CALL MNC_CW_ADD_VATTR_TEXT('dyF',1,'description',
131 adcroft 1.1 & 'y cell face separation',myThid)
132 edhill 1.2 CALL MNC_CW_ADD_VATTR_TEXT('dxG',1,'description',
133 adcroft 1.1 & 'x cell face separation',myThid)
134 edhill 1.2 CALL MNC_CW_ADD_VATTR_TEXT('dyG',1,'description',
135     & 'y cell face separation',myThid)
136     CALL MNC_CW_ADD_VATTR_TEXT('dxV',1,'description',
137     & 'x v-velocity separation',myThid)
138     CALL MNC_CW_ADD_VATTR_TEXT('dyU',1,'description',
139     & 'y u-velocity separation',myThid)
140     CALL MNC_CW_ADD_VATTR_TEXT('rA',1,'description',
141     & 'r-face area',myThid)
142     CALL MNC_CW_ADD_VATTR_TEXT('rAw',1,'description',
143     & '',myThid)
144     CALL MNC_CW_ADD_VATTR_TEXT('rAs',1,'description',
145     & '',myThid)
146     CALL MNC_CW_ADD_VATTR_TEXT('rAz',1,'description',
147     & '',myThid)
148     CALL MNC_CW_ADD_VATTR_TEXT('fCori',1,'description',
149     & 'Coriolis f',myThid)
150     CALL MNC_CW_ADD_VATTR_TEXT('fCoriG',1,'description',
151     & 'Coriolis f',myThid)
152     CALL MNC_CW_ADD_VATTR_TEXT('R_low',1,'description',
153     & 'base of fluid in r-units',myThid)
154     CALL MNC_CW_ADD_VATTR_TEXT('Ro_surf',1,'description',
155     & 'surface reference (at rest) position',myThid)
156     CALL MNC_CW_ADD_VATTR_TEXT('Depth',1,'description',
157     & 'fluid thickness in r coordinates (at rest)',myThid)
158     CALL MNC_CW_ADD_VATTR_TEXT('HFacC',1,'description',
159     & 'C-facing vertical fraction of open cell',myThid)
160     CALL MNC_CW_ADD_VATTR_TEXT('HFacW',1,'description',
161     & 'W-facing vertical fraction of open cell',myThid)
162     CALL MNC_CW_ADD_VATTR_TEXT('HFacS',1,'description',
163     & 'S-facing vertical fraction of open cell',myThid)
164    
165     C Write the GRID.h variables to a file
166     CALL MNC_CW_SET_UDIM(bfname, 0, myThid)
167     CALL MNC_CW_RL_W('D',bfname,0,0,'RC',rC,myThid)
168     CALL MNC_CW_RL_W('D',bfname,0,0,'RF',rF,myThid)
169     CALL MNC_CW_RL_W('D',bfname,0,0,'drC',drC,myThid)
170     CALL MNC_CW_RL_W('D',bfname,0,0,'drF',drF,myThid)
171     CALL MNC_CW_RL_W('D',bfname,0,0,'saFac',saFac,myThid)
172     CALL MNC_CW_RL_W('D',bfname,0,0,'XC',xC,myThid)
173     CALL MNC_CW_RL_W('D',bfname,0,0,'YC',yC,myThid)
174     CALL MNC_CW_RL_W('D',bfname,0,0,'XU',xG,myThid)
175     CALL MNC_CW_RL_W('D',bfname,0,0,'YU',yC,myThid)
176     CALL MNC_CW_RL_W('D',bfname,0,0,'XV',xC,myThid)
177     CALL MNC_CW_RL_W('D',bfname,0,0,'YV',yG,myThid)
178     CALL MNC_CW_RL_W('D',bfname,0,0,'XG',xG,myThid)
179     CALL MNC_CW_RL_W('D',bfname,0,0,'YG',yG,myThid)
180     CALL MNC_CW_RL_W('D',bfname,0,0,'dxC',dxC,myThid)
181     CALL MNC_CW_RL_W('D',bfname,0,0,'dyC',dyC,myThid)
182     CALL MNC_CW_RL_W('D',bfname,0,0,'dxF',dxF,myThid)
183     CALL MNC_CW_RL_W('D',bfname,0,0,'dyF',dyF,myThid)
184     CALL MNC_CW_RL_W('D',bfname,0,0,'dxG',dxG,myThid)
185     CALL MNC_CW_RL_W('D',bfname,0,0,'dyG',dyG,myThid)
186     CALL MNC_CW_RL_W('D',bfname,0,0,'dxV',dxV,myThid)
187     CALL MNC_CW_RL_W('D',bfname,0,0,'dyU',dyU,myThid)
188     CALL MNC_CW_RL_W('D',bfname,0,0,'rA',rA,myThid)
189     CALL MNC_CW_RL_W('D',bfname,0,0,'rAw',rAw,myThid)
190     CALL MNC_CW_RL_W('D',bfname,0,0,'rAs',rAs,myThid)
191     CALL MNC_CW_RL_W('D',bfname,0,0,'rAz',rAz,myThid)
192     CALL MNC_CW_RL_W('D',bfname,0,0,'fCori',fCori,myThid)
193     CALL MNC_CW_RL_W('D',bfname,0,0,'fCoriG',fCoriG,myThid)
194     CALL MNC_CW_RL_W('D',bfname,0,0,'R_low',R_low,myThid)
195     CALL MNC_CW_RL_W('D',bfname,0,0,'Ro_surf',Ro_surf,myThid)
196     CALL MNC_CW_RL_W('D',bfname,0,0,'Depth',tmpfld,myThid)
197     CALL MNC_CW_RL_W('D',bfname,0,0,'HFacC',HFacC,myThid)
198     CALL MNC_CW_RL_W('D',bfname,0,0,'HFacW',HFacW,myThid)
199     CALL MNC_CW_RL_W('D',bfname,0,0,'HFacS',HFacS,myThid)
200    
201 adcroft 1.1 ENDIF
202     #endif /* ALLOW_MNC */
203    
204     RETURN
205     END
206    
207     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
208    
  ViewVC Help
Powered by ViewVC 1.1.22