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

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

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

Revision 1.8 - (hide annotations) (download)
Wed Sep 29 03:45:14 2004 UTC (19 years, 8 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint55d_pre, checkpoint55e_post, checkpoint55d_post
Changes since 1.7: +5 -1 lines 
MNC output for adjoint state in addummy_in_stepping,
controlled via ALLOW_MNC, ALLOW_AUTODIFF_MONITOR (as is), and
useMNC, autodiff_mnc.
(TLM output t.b.d.)
1 heimbach 1.8 C $Header: /u/gcmpack/MITgcm/model/src/ini_model_io.F,v 1.7 2004/09/27 17:50:48 edhill Exp $
2 jmc 1.2 C $Name: $
3 adcroft 1.1
4     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
5    
6     #include "PACKAGES_CONFIG.h"
7     #include "CPP_OPTIONS.h"
8    
9     CBOP
10     C !ROUTINE: INI_MODEL_IO
11     C !INTERFACE:
12     SUBROUTINE INI_MODEL_IO( myThid )
13    
14     C !DESCRIPTION: \bv
15     C ini_model_io() is where run-time/experiment specific data are passed
16     C to any I/O packages ready that will be used for I/O of model state
17     C variables.
18     C This is specifically for setting up once only information such as
19     C shape/size of variables, units, etc... and is only for state variables.
20     C \ev
21    
22     C !CALLING SEQUENCE:
23     C PACKAGES_INIT_FIXED
24     C |
25     C |-- MNC_INIT
26    
27     IMPLICIT NONE
28     #include "SIZE.h"
29     #include "GRID.h"
30     #include "EEPARAMS.h"
31     #include "PARAMS.h"
32 edhill 1.4 #ifdef ALLOW_MNC
33     #include "MNC_PARAMS.h"
34     #endif
35 adcroft 1.1
36     C !INPUT/OUTPUT PARAMETERS:
37     C myThid - Number of this instances
38     INTEGER myThid
39     CEOP
40    
41    
42     C Flags specific to RW and MDSIO
43     C Set globalFiles flag for READ_WRITE_FLD package
44     CALL SET_WRITE_GLOBAL_FLD( globalFiles )
45     C Set globalFiles flag for READ_WRITE_REC package
46     CALL SET_WRITE_GLOBAL_REC( globalFiles )
47     C Set globalFiles flag for READ_WRITE_REC package
48     CALL SET_WRITE_GLOBAL_PICKUP( globalFiles )
49    
50     #ifdef ALLOW_MNC
51     IF (useMNC) THEN
52    
53     C Define coordinates for all MNC files ---------------------------------
54     C XC
55     CALL MNC_CW_ADD_VNAME('XC', 'Cen_xy_Hn__-__-', 3,4, myThid)
56 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('XC','description',
57 adcroft 1.1 & 'X coordinate of cell center (T-P point)',myThid)
58 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('XC',
59 adcroft 1.1 & 'units', 'degree_east', myThid)
60    
61     C YC
62     CALL MNC_CW_ADD_VNAME('YC', 'Cen_xy_Hn__-__-', 3,4, myThid)
63 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('YC','description',
64 adcroft 1.1 & 'Y coordinate of cell center (T-P point)',myThid)
65 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('YC',
66 adcroft 1.1 & 'units', 'degree_north', myThid)
67    
68     C XU
69     CALL MNC_CW_ADD_VNAME('XU', 'U_xy_Hn__-__-', 3,4, myThid)
70 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('XU','description',
71 adcroft 1.1 & 'X coordinate of U point',myThid)
72 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('XU',
73 adcroft 1.1 & 'units', 'degree_east', myThid)
74    
75     C YU
76     CALL MNC_CW_ADD_VNAME('YU', 'U_xy_Hn__-__-', 3,4, myThid)
77 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('YU','description',
78 adcroft 1.1 & 'Y coordinate of U point',myThid)
79 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('YU',
80 adcroft 1.1 & 'units', 'degree_north', myThid)
81    
82     C XV
83     CALL MNC_CW_ADD_VNAME('XV', 'V_xy_Hn__-__-', 3,4, myThid)
84 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('XV','description',
85 adcroft 1.1 & 'X coordinate of V point',myThid)
86 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('XV',
87 adcroft 1.1 & 'units', 'degree_east', myThid)
88    
89     C YV
90     CALL MNC_CW_ADD_VNAME('YV', 'V_xy_Hn__-__-', 3,4, myThid)
91 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('YV','description',
92 adcroft 1.1 & 'Y coordinate of V point',myThid)
93 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('YV',
94 adcroft 1.1 & 'units', 'degree_north', myThid)
95    
96     C XG
97     CALL MNC_CW_ADD_VNAME('XG', 'Cor_xy_Hn__-__-', 3,4, myThid)
98 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('XG','description',
99 adcroft 1.1 & 'X coordinate of cell corner (Vorticity point)',myThid)
100 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('XG',
101 adcroft 1.1 & 'units', 'degree_east', myThid)
102    
103     C YG
104     CALL MNC_CW_ADD_VNAME('YG', 'Cor_xy_Hn__-__-', 3,4, myThid)
105 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('YG','description',
106 adcroft 1.1 & 'Y coordinate of cell corner (Vorticity point)',myThid)
107 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('YG',
108 adcroft 1.1 & 'units', 'degree_north', myThid)
109    
110     C RC
111     CALL MNC_CW_ADD_VNAME('RC', '-_-_--__C__-', 0,0, myThid)
112 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('RC','description',
113 adcroft 1.1 & 'R coordinate of cell center',myThid)
114 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('RC',
115 adcroft 1.1 & 'units', 'm', myThid)
116    
117     C RF
118     CALL MNC_CW_ADD_VNAME('RF', '-_-_--__I__-', 0,0, myThid)
119 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('RF','description',
120 adcroft 1.1 & 'R coordinate of cell interface',myThid)
121 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('RF',
122 adcroft 1.1 & 'units', 'm', myThid)
123     C End define coordinates for all MNC files ----------------------------
124    
125    
126     C Create MNC definitions for DYNVARS.h variables
127     CALL MNC_CW_ADD_VNAME('iter', '-_-_--__-__t', 0,0, myThid)
128 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('iter',
129 adcroft 1.1 & 'long_name','iteration_count', myThid)
130    
131     CALL MNC_CW_ADD_VNAME('model_time', '-_-_--__-__t', 0,0, myThid)
132 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('model_time',
133 adcroft 1.1 & 'long_name','Model Time', myThid)
134 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('model_time','units','s', myThid)
135 adcroft 1.1
136     CALL MNC_CW_ADD_VNAME('U', 'U_xy_Hn__C__t', 4,5, myThid)
137 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('U','units','m/s', myThid)
138     CALL MNC_CW_ADD_VATTR_TEXT('U',
139 adcroft 1.1 & 'coordinates','XU YU RC iter', myThid)
140    
141     CALL MNC_CW_ADD_VNAME('V', 'V_xy_Hn__C__t', 4,5, myThid)
142 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('V','units','m/s', myThid)
143     CALL MNC_CW_ADD_VATTR_TEXT('V',
144 adcroft 1.1 & 'coordinates','XV YV RC iter', myThid)
145    
146 edhill 1.7 CALL MNC_CW_ADD_VNAME('Temp', 'Cen_xy_Hn__C__t', 4,5, myThid)
147     CALL MNC_CW_ADD_VATTR_TEXT('Temp','units','degC', myThid)
148     CALL MNC_CW_ADD_VATTR_TEXT('Temp','long_name',
149 adcroft 1.1 & 'potential_temperature', myThid)
150 edhill 1.7 CALL MNC_CW_ADD_VATTR_TEXT('Temp',
151 adcroft 1.1 & 'coordinates','XC YC RC iter', myThid)
152    
153     CALL MNC_CW_ADD_VNAME('S', 'Cen_xy_Hn__C__t', 4,5, myThid)
154 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('S','long_name',
155 adcroft 1.1 & 'salinity', myThid)
156 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('S',
157 adcroft 1.1 & 'coordinates','XC YC RC iter', myThid)
158    
159 jmc 1.2 CALL MNC_CW_ADD_VNAME('gUnm1', 'U_xy_Hn__C__t', 4,5, myThid)
160     CALL MNC_CW_ADD_VNAME('gVnm1', 'V_xy_Hn__C__t', 4,5, myThid)
161     CALL MNC_CW_ADD_VNAME('gTnm1', 'Cen_xy_Hn__C__t', 4,5, myThid)
162     CALL MNC_CW_ADD_VNAME('gSnm1', 'Cen_xy_Hn__C__t', 4,5, myThid)
163 adcroft 1.1
164     CALL MNC_CW_ADD_VNAME('Eta', 'Cen_xy_Hn__-__t', 3,4, myThid)
165 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('Eta','long_name',
166 adcroft 1.1 & 'free-surface_r-anomaly', myThid)
167 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('Eta','units','m', myThid)
168     CALL MNC_CW_ADD_VATTR_TEXT('Eta',
169 adcroft 1.1 & 'coordinates','XC YC RC iter', myThid)
170    
171 jmc 1.2 CALL MNC_CW_ADD_VNAME('EtaH', 'Cen_xy_Hn__-__t', 3,4, myThid)
172 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('EtaH','long_name',
173 jmc 1.2 & 'column-thickness_r-anomaly', myThid)
174 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('EtaH','units','m', myThid)
175     CALL MNC_CW_ADD_VATTR_TEXT('EtaH',
176 adcroft 1.1 & 'coordinates','XC YC RC iter', myThid)
177 jmc 1.3 CALL MNC_CW_ADD_VNAME('dEtaHdt', 'Cen_xy_Hn__-__t', 3,4, myThid)
178 adcroft 1.1
179     CALL MNC_CW_ADD_VNAME('W', 'Cen_xy_Hn__C__t', 4,5, myThid)
180 edhill 1.6 CALL MNC_CW_ADD_VATTR_TEXT('W','units','m/s', myThid)
181     CALL MNC_CW_ADD_VATTR_TEXT('W',
182 adcroft 1.1 & 'coordinates','XC YC RC iter', myThid)
183    
184 jmc 1.2 CALL MNC_CW_ADD_VNAME('phiHyd', 'Cen_xy_Hn__C__t', 4,5, myThid)
185 adcroft 1.1 CALL MNC_CW_ADD_VNAME('phiHydLow', 'Cen_xy_Hn__-__t', 3,4, myThid)
186     CALL MNC_CW_ADD_VNAME('phi_nh', 'Cen_xy_Hn__C__t', 4,5, myThid)
187     CALL MNC_CW_ADD_VNAME('gW', 'Cen_xy_Hn__C__t', 4,5, myThid)
188    
189     C Write coordinates to "state" file
190     CALL MNC_CW_SET_UDIM('state', 0, myThid)
191     CALL MNC_CW_RS_W('R','state',0,0,'XC',xC, myThid)
192     CALL MNC_CW_RS_W('R','state',0,0,'YC',yC, myThid)
193     CALL MNC_CW_RS_W('R','state',0,0,'XU',xG, myThid)
194     CALL MNC_CW_RS_W('R','state',0,0,'YU',yC, myThid)
195     CALL MNC_CW_RS_W('R','state',0,0,'XV',xC, myThid)
196     CALL MNC_CW_RS_W('R','state',0,0,'YV',yG, myThid)
197     CALL MNC_CW_RS_W('R','state',0,0,'XG',xG, myThid)
198     CALL MNC_CW_RS_W('R','state',0,0,'YG',yG, myThid)
199     CALL MNC_CW_RS_W('R','state',0,0,'RC',rC, myThid)
200     CALL MNC_CW_RS_W('R','state',0,0,'RF',rF, myThid)
201    
202 heimbach 1.8 #ifdef ALLOW_AUTODIFF_MONITOR
203     CALL AUTODIFF_INI_MODEL_IO( myThid )
204     #endif
205    
206 adcroft 1.1 ENDIF
207     #endif
208    
209     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
210    
211     RETURN
212     END
  ViewVC Help
Powered by ViewVC 1.1.22