/[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.7 by heimbach, Tue Aug 30 21:38:55 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    c
221    #endif
222    
223    
224  C     Write coordinates to "adstate" file  C     Write coordinates to "adstate" file
225        CALL MNC_CW_SET_UDIM('adstate', 0, myThid)        CALL MNC_CW_SET_UDIM('adstate', 0, myThid)
226        CALL MNC_CW_RS_W('R','adstate',0,0,'XC',xC, myThid)        CALL MNC_CW_RS_W('R','adstate',0,0,'XC',xC, myThid)
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22