/[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.3 by jmc, Tue Jul 6 00:54:07 2004 UTC revision 1.6 by edhill, Wed Sep 22 21:19:44 2004 UTC
# Line 29  C       |-- MNC_INIT Line 29  C       |-- MNC_INIT
29  #include "GRID.h"  #include "GRID.h"
30  #include "EEPARAMS.h"  #include "EEPARAMS.h"
31  #include "PARAMS.h"  #include "PARAMS.h"
32    #ifdef ALLOW_MNC
33    #include "MNC_PARAMS.h"
34    #endif
35    
36  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
37  C     myThid -  Number of this instances  C     myThid -  Number of this instances
# Line 50  C     Set globalFiles flag for READ_WRIT Line 53  C     Set globalFiles flag for READ_WRIT
53  C     Define coordinates for all MNC files ---------------------------------  C     Define coordinates for all MNC files ---------------------------------
54  C     XC  C     XC
55        CALL MNC_CW_ADD_VNAME('XC', 'Cen_xy_Hn__-__-', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('XC', 'Cen_xy_Hn__-__-', 3,4, myThid)
56        CALL MNC_CW_ADD_VATTR_TEXT('XC',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('XC','description',
57       &     'X coordinate of cell center (T-P point)',myThid)       &     'X coordinate of cell center (T-P point)',myThid)
58        CALL MNC_CW_ADD_VATTR_TEXT('XC',1,        CALL MNC_CW_ADD_VATTR_TEXT('XC',
59       &     'units', 'degree_east', myThid)       &     'units', 'degree_east', myThid)
60    
61  C     YC  C     YC
62        CALL MNC_CW_ADD_VNAME('YC', 'Cen_xy_Hn__-__-', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('YC', 'Cen_xy_Hn__-__-', 3,4, myThid)
63        CALL MNC_CW_ADD_VATTR_TEXT('YC',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('YC','description',
64       &     'Y coordinate of cell center (T-P point)',myThid)       &     'Y coordinate of cell center (T-P point)',myThid)
65        CALL MNC_CW_ADD_VATTR_TEXT('YC',1,        CALL MNC_CW_ADD_VATTR_TEXT('YC',
66       &     'units', 'degree_north', myThid)       &     'units', 'degree_north', myThid)
67    
68  C     XU  C     XU
69        CALL MNC_CW_ADD_VNAME('XU', 'U_xy_Hn__-__-', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('XU', 'U_xy_Hn__-__-', 3,4, myThid)
70        CALL MNC_CW_ADD_VATTR_TEXT('XU',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('XU','description',
71       &     'X coordinate of U point',myThid)       &     'X coordinate of U point',myThid)
72        CALL MNC_CW_ADD_VATTR_TEXT('XU',1,        CALL MNC_CW_ADD_VATTR_TEXT('XU',
73       &     'units', 'degree_east', myThid)       &     'units', 'degree_east', myThid)
74    
75  C     YU  C     YU
76        CALL MNC_CW_ADD_VNAME('YU', 'U_xy_Hn__-__-', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('YU', 'U_xy_Hn__-__-', 3,4, myThid)
77        CALL MNC_CW_ADD_VATTR_TEXT('YU',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('YU','description',
78       &     'Y coordinate of U point',myThid)       &     'Y coordinate of U point',myThid)
79        CALL MNC_CW_ADD_VATTR_TEXT('YU',1,        CALL MNC_CW_ADD_VATTR_TEXT('YU',
80       &     'units', 'degree_north', myThid)       &     'units', 'degree_north', myThid)
81    
82  C     XV  C     XV
83        CALL MNC_CW_ADD_VNAME('XV', 'V_xy_Hn__-__-', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('XV', 'V_xy_Hn__-__-', 3,4, myThid)
84        CALL MNC_CW_ADD_VATTR_TEXT('XV',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('XV','description',
85       &     'X coordinate of V point',myThid)       &     'X coordinate of V point',myThid)
86        CALL MNC_CW_ADD_VATTR_TEXT('XV',1,        CALL MNC_CW_ADD_VATTR_TEXT('XV',
87       &     'units', 'degree_east', myThid)       &     'units', 'degree_east', myThid)
88    
89  C     YV  C     YV
90        CALL MNC_CW_ADD_VNAME('YV', 'V_xy_Hn__-__-', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('YV', 'V_xy_Hn__-__-', 3,4, myThid)
91        CALL MNC_CW_ADD_VATTR_TEXT('YV',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('YV','description',
92       &     'Y coordinate of V point',myThid)       &     'Y coordinate of V point',myThid)
93        CALL MNC_CW_ADD_VATTR_TEXT('YV',1,        CALL MNC_CW_ADD_VATTR_TEXT('YV',
94       &     'units', 'degree_north', myThid)       &     'units', 'degree_north', myThid)
95    
96  C     XG  C     XG
97        CALL MNC_CW_ADD_VNAME('XG', 'Cor_xy_Hn__-__-', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('XG', 'Cor_xy_Hn__-__-', 3,4, myThid)
98        CALL MNC_CW_ADD_VATTR_TEXT('XG',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('XG','description',
99       &     'X coordinate of cell corner (Vorticity point)',myThid)       &     'X coordinate of cell corner (Vorticity point)',myThid)
100        CALL MNC_CW_ADD_VATTR_TEXT('XG',1,        CALL MNC_CW_ADD_VATTR_TEXT('XG',
101       &     'units', 'degree_east', myThid)       &     'units', 'degree_east', myThid)
102    
103  C     YG  C     YG
104        CALL MNC_CW_ADD_VNAME('YG', 'Cor_xy_Hn__-__-', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('YG', 'Cor_xy_Hn__-__-', 3,4, myThid)
105        CALL MNC_CW_ADD_VATTR_TEXT('YG',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('YG','description',
106       &     'Y coordinate of cell corner (Vorticity point)',myThid)       &     'Y coordinate of cell corner (Vorticity point)',myThid)
107        CALL MNC_CW_ADD_VATTR_TEXT('YG',1,        CALL MNC_CW_ADD_VATTR_TEXT('YG',
108       &     'units', 'degree_north', myThid)       &     'units', 'degree_north', myThid)
109    
110  C     RC  C     RC
111        CALL MNC_CW_ADD_VNAME('RC', '-_-_--__C__-', 0,0, myThid)        CALL MNC_CW_ADD_VNAME('RC', '-_-_--__C__-', 0,0, myThid)
112        CALL MNC_CW_ADD_VATTR_TEXT('RC',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('RC','description',
113       &     'R coordinate of cell center',myThid)       &     'R coordinate of cell center',myThid)
114        CALL MNC_CW_ADD_VATTR_TEXT('RC',1,        CALL MNC_CW_ADD_VATTR_TEXT('RC',
115       &     'units', 'm', myThid)       &     'units', 'm', myThid)
116    
117  C     RF  C     RF
118        CALL MNC_CW_ADD_VNAME('RF', '-_-_--__I__-', 0,0, myThid)        CALL MNC_CW_ADD_VNAME('RF', '-_-_--__I__-', 0,0, myThid)
119        CALL MNC_CW_ADD_VATTR_TEXT('RF',1,'description',        CALL MNC_CW_ADD_VATTR_TEXT('RF','description',
120       &     'R coordinate of cell interface',myThid)       &     'R coordinate of cell interface',myThid)
121        CALL MNC_CW_ADD_VATTR_TEXT('RF',1,        CALL MNC_CW_ADD_VATTR_TEXT('RF',
122       &     'units', 'm', myThid)       &     'units', 'm', myThid)
123  C     End define coordinates for all MNC files ----------------------------  C     End define coordinates for all MNC files ----------------------------
124    
125    
126  C     Create MNC definitions for DYNVARS.h variables  C     Create MNC definitions for DYNVARS.h variables
127        CALL MNC_CW_ADD_VNAME('iter', '-_-_--__-__t', 0,0, myThid)        CALL MNC_CW_ADD_VNAME('iter', '-_-_--__-__t', 0,0, myThid)
128        CALL MNC_CW_ADD_VATTR_TEXT('iter',1,        CALL MNC_CW_ADD_VATTR_TEXT('iter',
129       &     'long_name','iteration_count', myThid)       &     'long_name','iteration_count', myThid)
130    
131        CALL MNC_CW_ADD_VNAME('model_time', '-_-_--__-__t', 0,0, myThid)        CALL MNC_CW_ADD_VNAME('model_time', '-_-_--__-__t', 0,0, myThid)
132        CALL MNC_CW_ADD_VATTR_TEXT('model_time',1,        CALL MNC_CW_ADD_VATTR_TEXT('model_time',
133       &     'long_name','Model Time', myThid)       &     'long_name','Model Time', myThid)
134        CALL MNC_CW_ADD_VATTR_TEXT('model_time',1,'units','s', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('model_time','units','s', myThid)
135    
136        CALL MNC_CW_ADD_VNAME('U', 'U_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('U', 'U_xy_Hn__C__t', 4,5, myThid)
137        CALL MNC_CW_ADD_VATTR_TEXT('U',1,'units','m/s', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('U','units','m/s', myThid)
138        CALL MNC_CW_ADD_VATTR_TEXT('U',1,        CALL MNC_CW_ADD_VATTR_TEXT('U',
139       &     'coordinates','XU YU RC iter', myThid)       &     'coordinates','XU YU RC iter', myThid)
140    
141        CALL MNC_CW_ADD_VNAME('V', 'V_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('V', 'V_xy_Hn__C__t', 4,5, myThid)
142        CALL MNC_CW_ADD_VATTR_TEXT('V',1,'units','m/s', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('V','units','m/s', myThid)
143        CALL MNC_CW_ADD_VATTR_TEXT('V',1,        CALL MNC_CW_ADD_VATTR_TEXT('V',
144       &     'coordinates','XV YV RC iter', myThid)       &     'coordinates','XV YV RC iter', myThid)
145    
146        CALL MNC_CW_ADD_VNAME('T', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('T', 'Cen_xy_Hn__C__t', 4,5, myThid)
147        CALL MNC_CW_ADD_VATTR_TEXT('T',1,'units','degC', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('T','units','degC', myThid)
148        CALL MNC_CW_ADD_VATTR_TEXT('T',1,'long_name',        CALL MNC_CW_ADD_VATTR_TEXT('T','long_name',
149       &     'potential_temperature', myThid)       &     'potential_temperature', myThid)
150        CALL MNC_CW_ADD_VATTR_TEXT('T',1,        CALL MNC_CW_ADD_VATTR_TEXT('T',
151       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
152    
153        CALL MNC_CW_ADD_VNAME('S', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('S', 'Cen_xy_Hn__C__t', 4,5, myThid)
154        CALL MNC_CW_ADD_VATTR_TEXT('S',1,'long_name',        CALL MNC_CW_ADD_VATTR_TEXT('S','long_name',
155       &     'salinity', myThid)       &     'salinity', myThid)
156        CALL MNC_CW_ADD_VATTR_TEXT('S',1,        CALL MNC_CW_ADD_VATTR_TEXT('S',
157       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
158    
159        CALL MNC_CW_ADD_VNAME('gUnm1', 'U_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('gUnm1', 'U_xy_Hn__C__t', 4,5, myThid)
# Line 159  C     Create MNC definitions for DYNVARS Line 162  C     Create MNC definitions for DYNVARS
162        CALL MNC_CW_ADD_VNAME('gSnm1', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('gSnm1', 'Cen_xy_Hn__C__t', 4,5, myThid)
163    
164        CALL MNC_CW_ADD_VNAME('Eta', 'Cen_xy_Hn__-__t', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('Eta', 'Cen_xy_Hn__-__t', 3,4, myThid)
165        CALL MNC_CW_ADD_VATTR_TEXT('Eta',1,'long_name',        CALL MNC_CW_ADD_VATTR_TEXT('Eta','long_name',
166       &     'free-surface_r-anomaly', myThid)       &     'free-surface_r-anomaly', myThid)
167        CALL MNC_CW_ADD_VATTR_TEXT('Eta',1,'units','m', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('Eta','units','m', myThid)
168        CALL MNC_CW_ADD_VATTR_TEXT('Eta',1,        CALL MNC_CW_ADD_VATTR_TEXT('Eta',
169       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
170    
171        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)
172        CALL MNC_CW_ADD_VATTR_TEXT('EtaH',1,'long_name',        CALL MNC_CW_ADD_VATTR_TEXT('EtaH','long_name',
173       &     'column-thickness_r-anomaly', myThid)       &     'column-thickness_r-anomaly', myThid)
174        CALL MNC_CW_ADD_VATTR_TEXT('EtaH',1,'units','m', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('EtaH','units','m', myThid)
175        CALL MNC_CW_ADD_VATTR_TEXT('EtaH',1,        CALL MNC_CW_ADD_VATTR_TEXT('EtaH',
176       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
177        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)
178    
179        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)
180        CALL MNC_CW_ADD_VATTR_TEXT('W',1,'units','m/s', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('W','units','m/s', myThid)
181        CALL MNC_CW_ADD_VATTR_TEXT('W',1,        CALL MNC_CW_ADD_VATTR_TEXT('W',
182       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
183    
184        CALL MNC_CW_ADD_VNAME('phiHyd', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('phiHyd', 'Cen_xy_Hn__C__t', 4,5, myThid)
# Line 183  C     Create MNC definitions for DYNVARS Line 186  C     Create MNC definitions for DYNVARS
186        CALL MNC_CW_ADD_VNAME('phi_nh', 'Cen_xy_Hn__C__t', 4,5, myThid)        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)        CALL MNC_CW_ADD_VNAME('gW', 'Cen_xy_Hn__C__t', 4,5, myThid)
188    
189        CALL MNC_CW_ADD_VNAME('tr1', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VATTR_TEXT('T',
       CALL MNC_CW_ADD_VATTR_TEXT('tr1',1,  
      &     'long_name','passive_tracer_1', myThid)  
       CALL MNC_CW_ADD_VATTR_TEXT('T',1,  
190       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
191    
192  C     Write coordinates to "state" file  C     Write coordinates to "state" file
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22