/[MITgcm]/MITgcm/pkg/autodiff/autodiff_ini_model_io.F
ViewVC logotype

Diff of /MITgcm/pkg/autodiff/autodiff_ini_model_io.F

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

revision 1.6 by jmc, Wed Apr 27 14:10:06 2005 UTC revision 1.10 by heimbach, Fri Sep 2 23:21:22 2005 UTC
# Line 131  CEOP Line 131  CEOP
131       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
132  #endif  #endif
133    
134    #ifdef ALLOW_EXF
135    c
136          CALL MNC_CW_ADD_VNAME('adustress', 'Cen_xy_Hn__-__t', 3,4, myThid)
137          CALL MNC_CW_ADD_VATTR_TEXT('adustress',
138         &     'units','[cost]/[N/m^2]', myThid)
139          CALL MNC_CW_ADD_VATTR_TEXT('adustress','long_name',
140         &     'adjoint zonal wind stress', myThid)
141          CALL MNC_CW_ADD_VATTR_TEXT('adustress',
142         &     'coordinates','XC YC RC iter', myThid)
143    c
144          CALL MNC_CW_ADD_VNAME('advstress', 'Cen_xy_Hn__-__t', 3,4, myThid)
145          CALL MNC_CW_ADD_VATTR_TEXT('advstress',
146         &     'units','[cost]/[N/m^2]', myThid)
147          CALL MNC_CW_ADD_VATTR_TEXT('advstress','long_name',
148         &     'adjoint meridional wind stress', myThid)
149          CALL MNC_CW_ADD_VATTR_TEXT('advstress',
150         &     'coordinates','XC YC RC iter', myThid)
151    c
152          CALL MNC_CW_ADD_VNAME('adhflux', 'Cen_xy_Hn__-__t', 3,4, myThid)
153          CALL MNC_CW_ADD_VATTR_TEXT('adhflux',
154         &     'units','[cost]/[W/m^2]', myThid)
155          CALL MNC_CW_ADD_VATTR_TEXT('adhflux','long_name',
156         &     'adjoint net upward heat flux', myThid)
157          CALL MNC_CW_ADD_VATTR_TEXT('adhflux',
158         &     'coordinates','XC YC RC iter', myThid)
159    c
160          CALL MNC_CW_ADD_VNAME('adsflux', 'Cen_xy_Hn__-__t', 3,4, myThid)
161          CALL MNC_CW_ADD_VATTR_TEXT('adsflux',
162         &     'units','[cost]/[m/s]', myThid)
163          CALL MNC_CW_ADD_VATTR_TEXT('adsflux','long_name',
164         &     'adjoint net upward freshwater flux', myThid)
165          CALL MNC_CW_ADD_VATTR_TEXT('adsflux',
166         &     'coordinates','XC YC RC iter', myThid)
167    c
168    # ifdef ALLOW_ATM_TEMP
169          CALL MNC_CW_ADD_VNAME('adatemp', 'Cen_xy_Hn__-__t', 3,4, myThid)
170          CALL MNC_CW_ADD_VATTR_TEXT('adatemp',
171         &     'units','[cost]/[degK]', myThid)
172          CALL MNC_CW_ADD_VATTR_TEXT('adatemp','long_name',
173         &     'adjoint surface air temperature', myThid)
174          CALL MNC_CW_ADD_VATTR_TEXT('adatemp',
175         &     'coordinates','XC YC RC iter', myThid)
176    c
177          CALL MNC_CW_ADD_VNAME('adaqh', 'Cen_xy_Hn__-__t', 3,4, myThid)
178          CALL MNC_CW_ADD_VATTR_TEXT('adaqh',
179         &     'units','[cost]/[kg/kg]', myThid)
180          CALL MNC_CW_ADD_VATTR_TEXT('adaqh','long_name',
181         &     'adjoint specific humidity', myThid)
182          CALL MNC_CW_ADD_VATTR_TEXT('adaqh',
183         &     'coordinates','XC YC RC iter', myThid)
184    c
185          CALL MNC_CW_ADD_VNAME('adprecip', 'Cen_xy_Hn__-__t', 3,4, myThid)
186          CALL MNC_CW_ADD_VATTR_TEXT('adprecip',
187         &     'units','[cost]/[m/s]', myThid)
188          CALL MNC_CW_ADD_VATTR_TEXT('adprecip','long_name',
189         &     'adjoint precipitation', myThid)
190          CALL MNC_CW_ADD_VATTR_TEXT('adprecip',
191         &     'coordinates','XC YC RC iter', myThid)
192    # endif
193    
194    # ifdef ALLOW_ATM_WIND
195          CALL MNC_CW_ADD_VNAME('aduwind', 'Cen_xy_Hn__-__t', 3,4, myThid)
196          CALL MNC_CW_ADD_VATTR_TEXT('aduwind',
197         &     'units','[cost]/[m/s]', myThid)
198          CALL MNC_CW_ADD_VATTR_TEXT('aduwind','long_name',
199         &     'adjoint zonal wind speed', myThid)
200          CALL MNC_CW_ADD_VATTR_TEXT('aduwind',
201         &     'coordinates','XC YC RC iter', myThid)
202    c
203          CALL MNC_CW_ADD_VNAME('advwind', 'Cen_xy_Hn__-__t', 3,4, myThid)
204          CALL MNC_CW_ADD_VATTR_TEXT('advwind',
205         &     'units','[cost]/[m/s]', myThid)
206          CALL MNC_CW_ADD_VATTR_TEXT('advwind','long_name',
207         &     'adjoint meridional wind speed', myThid)
208          CALL MNC_CW_ADD_VATTR_TEXT('advwind',
209         &     'coordinates','XC YC RC iter', myThid)
210    # endif
211    # ifdef ALLOW_DOWNWARD_RADIATION
212          CALL MNC_CW_ADD_VNAME('adswdown', 'Cen_xy_Hn__-__t', 3,4, myThid)
213          CALL MNC_CW_ADD_VATTR_TEXT('adswdown',
214         &     'units','[cost]/[W/m^2]', myThid)
215          CALL MNC_CW_ADD_VATTR_TEXT('adswdown','long_name',
216         &     'adjoint downward shortwave radiation', myThid)
217          CALL MNC_CW_ADD_VATTR_TEXT('adswdown',
218         &     'coordinates','XC YC RC iter', myThid)
219    # endif
220    # ifdef ALLOW_CLIMSST_RELAXATION
221          CALL MNC_CW_ADD_VNAME('adclimsst', 'Cen_xy_Hn__-__t', 3,4, myThid)
222          CALL MNC_CW_ADD_VATTR_TEXT('adclimsst',
223         &     'units','[cost]/[degC]', myThid)
224          CALL MNC_CW_ADD_VATTR_TEXT('adclimsst','long_name',
225         &     'adjoint sea surface temperature', myThid)
226          CALL MNC_CW_ADD_VATTR_TEXT('adclimsst',
227         &     'coordinates','XC YC RC iter', myThid)
228    # endif
229    # ifdef ALLOW_CLIMSSS_RELAXATION
230          CALL MNC_CW_ADD_VNAME('adclimsss', 'Cen_xy_Hn__-__t', 3,4, myThid)
231          CALL MNC_CW_ADD_VATTR_TEXT('adclimsss',
232         &     'units','[cost]/[ppt]', myThid)
233          CALL MNC_CW_ADD_VATTR_TEXT('adclimsss','long_name',
234         &     'adjoint sea surface salinity', myThid)
235          CALL MNC_CW_ADD_VATTR_TEXT('adclimsss',
236         &     'coordinates','XC YC RC iter', myThid)
237    # endif
238    c
239    #endif
240    
241    #ifdef ALLOW_SEAICE
242    c
243          CALL MNC_CW_ADD_VNAME('adarea', 'Cen_xy_Hn__-__t', 3,4, myThid)
244          CALL MNC_CW_ADD_VATTR_TEXT('adarea',
245         &     'units','[cost]/[m^2/m^2]', myThid)
246          CALL MNC_CW_ADD_VATTR_TEXT('adarea','long_name',
247         &     'adjoint fractional ice-covered area', myThid)
248          CALL MNC_CW_ADD_VATTR_TEXT('adarea',
249         &     'coordinates','XC YC RC iter', myThid)
250    c
251          CALL MNC_CW_ADD_VNAME('adheff', 'Cen_xy_Hn__-__t', 3,4, myThid)
252          CALL MNC_CW_ADD_VATTR_TEXT('adheff',
253         &     'units','[cost]/[m]', myThid)
254          CALL MNC_CW_ADD_VATTR_TEXT('adheff','long_name',
255         &     'adjoint effective ice thickness', myThid)
256          CALL MNC_CW_ADD_VATTR_TEXT('adheff',
257         &     'coordinates','XC YC RC iter', myThid)
258    c
259          CALL MNC_CW_ADD_VNAME('adhsnow', 'Cen_xy_Hn__-__t', 3,4, myThid)
260          CALL MNC_CW_ADD_VATTR_TEXT('adhsnow',
261         &     'units','[cost]/[m]', myThid)
262          CALL MNC_CW_ADD_VATTR_TEXT('adhsnow','long_name',
263         &     'adjoint snow thickness', myThid)
264          CALL MNC_CW_ADD_VATTR_TEXT('adhsnow',
265         &     'coordinates','XC YC RC iter', myThid)
266    c
267    # ifdef ALLOW_DYNAMICS
268          CALL MNC_CW_ADD_VNAME('aduice', 'U_xy_Hn__-__t', 3,4, myThid)
269          CALL MNC_CW_ADD_VATTR_TEXT('aduice',
270         &     'units','[cost]/[m/s]', myThid)
271          CALL MNC_CW_ADD_VATTR_TEXT('aduice','long_name',
272         &     'adjoint zonal ice velocity', myThid)
273          CALL MNC_CW_ADD_VATTR_TEXT('aduice',
274         &     'coordinates','XG YG RC iter', myThid)
275    c
276          CALL MNC_CW_ADD_VNAME('advice', 'V_xy_Hn__-__t', 3,4, myThid)
277          CALL MNC_CW_ADD_VATTR_TEXT('advice',
278         &     'units','[cost]/[m/s]', myThid)
279          CALL MNC_CW_ADD_VATTR_TEXT('advice','long_name',
280         &     'adjoint meridional ice velocity', myThid)
281          CALL MNC_CW_ADD_VATTR_TEXT('advice',
282         &     'coordinates','XG YG RC iter', myThid)
283    # endif
284    c
285    #endif
286    
287  C     Write coordinates to "adstate" file  C     Write coordinates to "adstate" file
288        CALL MNC_CW_SET_UDIM('adstate', 0, myThid)        CALL MNC_CW_SET_UDIM('adstate', 0, myThid)
289        CALL MNC_CW_RS_W('R','adstate',0,0,'XC',xC, myThid)        CALL MNC_CW_RS_W('R','adstate',0,0,'XC',xC, myThid)
# Line 144  C     Write coordinates to "adstate" fil Line 297  C     Write coordinates to "adstate" fil
297        CALL MNC_CW_RS_W('R','adstate',0,0,'RC',rC, myThid)        CALL MNC_CW_RS_W('R','adstate',0,0,'RC',rC, myThid)
298        CALL MNC_CW_RS_W('R','adstate',0,0,'RF',rF, myThid)        CALL MNC_CW_RS_W('R','adstate',0,0,'RF',rF, myThid)
299    
300    #ifdef ALLOW_EXF
301          CALL MNC_CW_SET_UDIM('adexf', 0, myThid)
302          CALL MNC_CW_RS_W('R','adexf',0,0,'XC',xC, myThid)
303          CALL MNC_CW_RS_W('R','adexf',0,0,'YC',yC, myThid)
304          CALL MNC_CW_RS_W('R','adexf',0,0,'XU',xG, myThid)
305          CALL MNC_CW_RS_W('R','adexf',0,0,'YU',yC, myThid)
306          CALL MNC_CW_RS_W('R','adexf',0,0,'XV',xC, myThid)
307          CALL MNC_CW_RS_W('R','adexf',0,0,'YV',yG, myThid)
308          CALL MNC_CW_RS_W('R','adexf',0,0,'XG',xG, myThid)
309          CALL MNC_CW_RS_W('R','adexf',0,0,'YG',yG, myThid)
310          CALL MNC_CW_RS_W('R','adexf',0,0,'RC',rC, myThid)
311          CALL MNC_CW_RS_W('R','adexf',0,0,'RF',rF, myThid)
312    #endif
313    
314    #ifdef ALLOW_SEAICE
315          CALL MNC_CW_SET_UDIM('adseaice', 0, myThid)
316          CALL MNC_CW_RS_W('R','adseaice',0,0,'XC',xC, myThid)
317          CALL MNC_CW_RS_W('R','adseaice',0,0,'YC',yC, myThid)
318          CALL MNC_CW_RS_W('R','adseaice',0,0,'XU',xG, myThid)
319          CALL MNC_CW_RS_W('R','adseaice',0,0,'YU',yC, myThid)
320          CALL MNC_CW_RS_W('R','adseaice',0,0,'XV',xC, myThid)
321          CALL MNC_CW_RS_W('R','adseaice',0,0,'YV',yG, myThid)
322          CALL MNC_CW_RS_W('R','adseaice',0,0,'XG',xG, myThid)
323          CALL MNC_CW_RS_W('R','adseaice',0,0,'YG',yG, myThid)
324          CALL MNC_CW_RS_W('R','adseaice',0,0,'RC',rC, myThid)
325          CALL MNC_CW_RS_W('R','adseaice',0,0,'RF',rF, myThid)
326    #endif
327    
328        ENDIF        ENDIF
329  #endif  #endif
330    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22