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

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

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

revision 1.7 by edhill, Mon Sep 27 17:50:48 2004 UTC revision 1.17 by edhill, Sun Apr 3 05:16:43 2005 UTC
# Line 8  C---+----1----+----2----+----3----+----4 Line 8  C---+----1----+----2----+----3----+----4
8    
9  CBOP  CBOP
10  C     !ROUTINE: INI_MODEL_IO  C     !ROUTINE: INI_MODEL_IO
11    
12  C     !INTERFACE:  C     !INTERFACE:
13        SUBROUTINE INI_MODEL_IO( myThid )        SUBROUTINE INI_MODEL_IO( myThid )
14    
15  C     !DESCRIPTION: \bv  C     !DESCRIPTION:
16  C     ini_model_io() is where run-time/experiment specific data are passed  C     Pass specific setup data to any I/O packages necessary for I/O of
17  C     to any I/O packages ready that will be used for I/O of model state  C     model state variables.  This is specifically for setting up (once
18  C     variables.  C     only!)  information such as shape/size of variables, units,
19  C     This is specifically for setting up once only information such as  C     etc... and is primarily for state and snapshot variables.
 C     shape/size of variables, units, etc... and is only for state variables.  
 C     \ev  
   
 C     !CALLING SEQUENCE:  
 C     PACKAGES_INIT_FIXED  
 C       |  
 C       |-- MNC_INIT  
20    
21    C     !USES:
22        IMPLICIT NONE        IMPLICIT NONE
23  #include "SIZE.h"  #include "SIZE.h"
24  #include "GRID.h"  #include "GRID.h"
25  #include "EEPARAMS.h"  #include "EEPARAMS.h"
26  #include "PARAMS.h"  #include "PARAMS.h"
 #ifdef ALLOW_MNC  
 #include "MNC_PARAMS.h"  
 #endif  
27    
28  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
29  C     myThid -  Number of this instances  C     myThid -  Number of this instances
# Line 120  C     RF Line 112  C     RF
112       &     'R coordinate of cell interface',myThid)       &     'R coordinate of cell interface',myThid)
113        CALL MNC_CW_ADD_VATTR_TEXT('RF',        CALL MNC_CW_ADD_VATTR_TEXT('RF',
114       &     'units', 'm', myThid)       &     'units', 'm', myThid)
115    
116    C     RL
117          CALL MNC_CW_ADD_VNAME('RL', '-_-_--__L__-', 0,0, myThid)
118          CALL MNC_CW_ADD_VATTR_TEXT('RL','description',
119         &     'R coordinate of lower interface',myThid)
120          CALL MNC_CW_ADD_VATTR_TEXT('RL',
121         &     'units', 'm', myThid)
122    
123    C     RU
124          CALL MNC_CW_ADD_VNAME('RU', '-_-_--__U__-', 0,0, myThid)
125          CALL MNC_CW_ADD_VATTR_TEXT('RU','description',
126         &     'R coordinate of upper interface',myThid)
127          CALL MNC_CW_ADD_VATTR_TEXT('RU',
128         &     'units', 'm', myThid)
129  C     End define coordinates for all MNC files ----------------------------  C     End define coordinates for all MNC files ----------------------------
130    
131    
132  C     Create MNC definitions for DYNVARS.h variables  C     Create MNC definitions for DYNVARS.h variables
133        CALL MNC_CW_ADD_VNAME('iter', '-_-_--__-__t', 0,0, myThid)        CALL MNC_CW_ADD_VNAME('T', '-_-_--__-__t', 0,0, myThid)
134        CALL MNC_CW_ADD_VATTR_TEXT('iter',        CALL MNC_CW_ADD_VATTR_TEXT('T',
135       &     'long_name','iteration_count', myThid)       &     'long_name','iteration_count', myThid)
136    
137        CALL MNC_CW_ADD_VNAME('model_time', '-_-_--__-__t', 0,0, myThid)        CALL MNC_CW_ADD_VNAME('model_time', '-_-_--__-__t', 0,0, myThid)
# Line 166  C     Create MNC definitions for DYNVARS Line 172  C     Create MNC definitions for DYNVARS
172       &     'free-surface_r-anomaly', myThid)       &     'free-surface_r-anomaly', myThid)
173        CALL MNC_CW_ADD_VATTR_TEXT('Eta','units','m', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('Eta','units','m', myThid)
174        CALL MNC_CW_ADD_VATTR_TEXT('Eta',        CALL MNC_CW_ADD_VATTR_TEXT('Eta',
175       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC iter', myThid)
176    
177        CALL MNC_CW_ADD_VNAME('EtaH', 'Cen_xy_Hn__-__t', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('EtaH', 'Cen_xy_Hn__-__t', 3,4, myThid)
178        CALL MNC_CW_ADD_VATTR_TEXT('EtaH','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('EtaH','long_name',
179       &     'column-thickness_r-anomaly', myThid)       &     'column-thickness_r-anomaly', myThid)
180        CALL MNC_CW_ADD_VATTR_TEXT('EtaH','units','m', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('EtaH','units','m', myThid)
181        CALL MNC_CW_ADD_VATTR_TEXT('EtaH',        CALL MNC_CW_ADD_VATTR_TEXT('EtaH',
182       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC iter', myThid)
183        CALL MNC_CW_ADD_VNAME('dEtaHdt', 'Cen_xy_Hn__-__t', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('dEtaHdt', 'Cen_xy_Hn__-__t', 3,4, myThid)
184    
185        CALL MNC_CW_ADD_VNAME('W', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('W', 'Cen_xy_Hn__C__t', 4,5, myThid)
# Line 187  C     Create MNC definitions for DYNVARS Line 193  C     Create MNC definitions for DYNVARS
193        CALL MNC_CW_ADD_VNAME('gW', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('gW', 'Cen_xy_Hn__C__t', 4,5, myThid)
194    
195  C     Write coordinates to "state" file  C     Write coordinates to "state" file
196        CALL MNC_CW_SET_UDIM('state', 0, myThid)  C     CALL MNC_CW_SET_UDIM('state', 0, myThid)
197        CALL MNC_CW_RS_W('R','state',0,0,'XC',xC, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'XC',xC, myThid)
198        CALL MNC_CW_RS_W('R','state',0,0,'YC',yC, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'YC',yC, myThid)
199        CALL MNC_CW_RS_W('R','state',0,0,'XU',xG, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'XU',xG, myThid)
200        CALL MNC_CW_RS_W('R','state',0,0,'YU',yC, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'YU',yC, myThid)
201        CALL MNC_CW_RS_W('R','state',0,0,'XV',xC, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'XV',xC, myThid)
202        CALL MNC_CW_RS_W('R','state',0,0,'YV',yG, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'YV',yG, myThid)
203        CALL MNC_CW_RS_W('R','state',0,0,'XG',xG, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'XG',xG, myThid)
204        CALL MNC_CW_RS_W('R','state',0,0,'YG',yG, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'YG',yG, myThid)
205        CALL MNC_CW_RS_W('R','state',0,0,'RC',rC, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'RC',rC, myThid)
206        CALL MNC_CW_RS_W('R','state',0,0,'RF',rF, myThid)  C     CALL MNC_CW_RS_W('R','state',0,0,'RF',rF, myThid)
207    
208    C     Define variables used in mom_vecinv
209    C     CALL MNC_CW_ADD_VNAME('fV', 'Cen_xy_Hn__-__t', 0,0, myThid)
210    C     CALL MNC_CW_ADD_VATTR_TEXT('','units','', myThid)
211    C     CALL MNC_CW_ADD_VATTR_TEXT('','long_name',
212    C     &     'potential_temperature', myThid)
213          
214          CALL MNC_CW_ADD_VNAME('fV', 'Cen_xy_Hn__C__t', 0,0, myThid)
215          CALL MNC_CW_ADD_VNAME('fU', 'Cen_xy_Hn__C__t', 0,0, myThid)
216          CALL MNC_CW_ADD_VNAME('zV', 'Cen_xy_Hn__C__t', 0,0, myThid)
217          CALL MNC_CW_ADD_VNAME('zU', 'Cen_xy_Hn__C__t', 0,0, myThid)
218          CALL MNC_CW_ADD_VNAME('KEx','Cen_xy_Hn__C__t', 0,0, myThid)
219          CALL MNC_CW_ADD_VNAME('KEy','Cen_xy_Hn__C__t', 0,0, myThid)
220          CALL MNC_CW_ADD_VNAME('Ds', 'Cen_xy_Hn__C__t', 0,0, myThid)
221          CALL MNC_CW_ADD_VNAME('Dt', 'Cen_xy_Hn__C__t', 0,0, myThid)
222          CALL MNC_CW_ADD_VNAME('Du', 'Cen_xy_Hn__C__t', 0,0, myThid)
223          CALL MNC_CW_ADD_VNAME('Dv', 'Cen_xy_Hn__C__t', 0,0, myThid)
224          CALL MNC_CW_ADD_VNAME('Z3', 'Cen_xy_Hn__C__t', 0,0, myThid)
225          CALL MNC_CW_ADD_VNAME('W3', 'Cen_xy_Hn__C__t', 0,0, myThid)
226          CALL MNC_CW_ADD_VNAME('KE', 'Cen_xy_Hn__C__t', 0,0, myThid)
227          CALL MNC_CW_ADD_VNAME('D',  'Cen_xy_Hn__C__t', 0,0, myThid)
228    
229    C     Define variables from FFIELDS.h
230          CALL MNC_CW_ADD_VNAME('fu', 'U_xy_Hn__-__t', 3,4, myThid)
231          CALL MNC_CW_ADD_VATTR_TEXT('fu','units','N/m^2', myThid)
232          CALL MNC_CW_ADD_VATTR_TEXT('fu','description',
233         &     'Zonal surface wind stress', myThid)
234    
235          CALL MNC_CW_ADD_VNAME('fv', 'V_xy_Hn__-__t', 3,4, myThid)
236          CALL MNC_CW_ADD_VATTR_TEXT('fv','units','N/m^2', myThid)
237          CALL MNC_CW_ADD_VATTR_TEXT('fv','description',
238         &     'Meridional surface wind stress', myThid)
239    
240          CALL MNC_CW_ADD_VNAME('Qnet', 'Cen_xy_Hn__-__t', 3,4, myThid)
241          CALL MNC_CW_ADD_VATTR_TEXT('Qnet','units','W/m^2', myThid)
242          CALL MNC_CW_ADD_VATTR_TEXT('Qnet','description',
243         &     'Net upward surface heat flux (including shortwave)',
244         &     myThid)
245    
246          CALL MNC_CW_ADD_VNAME('Qsw', 'Cen_xy_Hn__-__t', 3,4, myThid)
247          CALL MNC_CW_ADD_VATTR_TEXT('Qsw','units','W/m^2', myThid)
248          CALL MNC_CW_ADD_VATTR_TEXT('Qsw','description',
249         &     'Net upward shortwave radiation', myThid)
250    
251          CALL MNC_CW_ADD_VNAME('dQdt', 'Cen_xy_Hn__-__t', 3,4, myThid)
252          CALL MNC_CW_ADD_VATTR_TEXT('dQdt','units',
253         &     'W/m^2/degrees', myThid)
254    
255          CALL MNC_CW_ADD_VNAME('EmPmR', 'Cen_xy_Hn__-__t', 3,4, myThid)
256          CALL MNC_CW_ADD_VATTR_TEXT('EmPmR','units','m/s', myThid)
257          CALL MNC_CW_ADD_VATTR_TEXT('EmPmR','description',
258         &     'Net upward freshwater flux', myThid)
259    
260          CALL MNC_CW_ADD_VNAME('saltFlux', 'Cen_xy_Hn__-__t', 3,4, myThid)
261          CALL MNC_CW_ADD_VATTR_TEXT('saltFlux','units',
262         &     'psu.kg/m^2/s', myThid)
263          CALL MNC_CW_ADD_VATTR_TEXT('saltFlux','description',
264         &     'Net upward salt flux', myThid)
265    
266          CALL MNC_CW_ADD_VNAME('SST', 'Cen_xy_Hn__-__t', 3,4, myThid)
267          CALL MNC_CW_ADD_VATTR_TEXT('SST','units','deg C', myThid)
268          CALL MNC_CW_ADD_VATTR_TEXT('SST','description',
269         &     ' Sea surface temperature for relaxation', myThid)
270    
271          CALL MNC_CW_ADD_VNAME('SSS', 'Cen_xy_Hn__-__t', 3,4, myThid)
272          CALL MNC_CW_ADD_VATTR_TEXT('SSS','units','psu', myThid)
273          CALL MNC_CW_ADD_VATTR_TEXT('SSS','description',
274         &     'Sea surface salinity for relaxation', myThid)
275    
276    
277    #ifdef ALLOW_AUTODIFF_MONITOR
278          CALL AUTODIFF_INI_MODEL_IO( myThid )
279    #endif
280    
281        ENDIF        ENDIF
282  #endif  #endif
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22