/[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.3 by heimbach, Tue Nov 16 05:42:11 2004 UTC revision 1.5 by edhill, Sat Dec 18 02:18:55 2004 UTC
# Line 38  CEOP Line 38  CEOP
38    
39        IF (useMNC) THEN        IF (useMNC) THEN
40    
 C     Create MNC definitions for DYNVARS.h variables  
 c      CALL MNC_CW_ADD_VNAME('iter', '-_-_--__-__t', 0,0, myThid)  
 c      CALL MNC_CW_ADD_VATTR_TEXT('iter',  
 c     &     'long_name','iteration_count', myThid)  
   
 c      CALL MNC_CW_ADD_VNAME('model_time', '-_-_--__-__t', 0,0, myThid)  
 c      CALL MNC_CW_ADD_VATTR_TEXT('model_time',  
 c     &     'long_name','Model Time', myThid)  
 c      CALL MNC_CW_ADD_VATTR_TEXT('model_time','units','s', myThid)  
   
41        CALL MNC_CW_ADD_VNAME('adU', 'U_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('adU', 'U_xy_Hn__C__t', 4,5, myThid)
42        CALL MNC_CW_ADD_VATTR_TEXT('adU','units','[cost]/[m/s]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adU','units','[cost]/[m/s]', myThid)
43        CALL MNC_CW_ADD_VATTR_TEXT('adU','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adU','long_name',
44       &     'adj. zonal velocity', myThid)       &     'adjoint zonal velocity', myThid)
45        CALL MNC_CW_ADD_VATTR_TEXT('adU',        CALL MNC_CW_ADD_VATTR_TEXT('adU',
46       &     'coordinates','XU YU RC iter', myThid)       &     'coordinates','XU YU RC iter', myThid)
47    
48        CALL MNC_CW_ADD_VNAME('adV', 'V_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('adV', 'V_xy_Hn__C__t', 4,5, myThid)
49        CALL MNC_CW_ADD_VATTR_TEXT('adV','units','[cost]/[m/s]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adV','units','[cost]/[m/s]', myThid)
50        CALL MNC_CW_ADD_VATTR_TEXT('adV','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adV','long_name',
51       &     'adj. merid. velocity', myThid)       &     'adjoint merid. velocity', myThid)
52        CALL MNC_CW_ADD_VATTR_TEXT('adV',        CALL MNC_CW_ADD_VATTR_TEXT('adV',
53       &     'coordinates','XV YV RC iter', myThid)       &     'coordinates','XV YV RC iter', myThid)
54    
55        CALL MNC_CW_ADD_VNAME('adT', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('adT', 'Cen_xy_Hn__C__t', 4,5, myThid)
56        CALL MNC_CW_ADD_VATTR_TEXT('adT','units','[cost]/[degC]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adT','units','[cost]/[degC]', myThid)
57        CALL MNC_CW_ADD_VATTR_TEXT('adT','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adT','long_name',
58       &     'adj. potential_temperature', myThid)       &     'adjoint potential_temperature', myThid)
59        CALL MNC_CW_ADD_VATTR_TEXT('adT',        CALL MNC_CW_ADD_VATTR_TEXT('adT',
60       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
61    
62        CALL MNC_CW_ADD_VNAME('adS', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('adS', 'Cen_xy_Hn__C__t', 4,5, myThid)
63        CALL MNC_CW_ADD_VATTR_TEXT('adS','units','[cost]/[psu]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adS','units','[cost]/[psu]', myThid)
64        CALL MNC_CW_ADD_VATTR_TEXT('adS','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adS','long_name',
65       &     'adj. salinity', myThid)       &     'adjoint salinity', myThid)
66        CALL MNC_CW_ADD_VATTR_TEXT('adS',        CALL MNC_CW_ADD_VATTR_TEXT('adS',
67       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
68    
69        CALL MNC_CW_ADD_VNAME('adEta', 'Cen_xy_Hn__-__t', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('adEta', 'Cen_xy_Hn__-__t', 3,4, myThid)
70        CALL MNC_CW_ADD_VATTR_TEXT('adEta','units','[cost]/[m]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adEta','units','[cost]/[m]', myThid)
71        CALL MNC_CW_ADD_VATTR_TEXT('adEta','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adEta','long_name',
72       &     'adj. free-surface_r-anomaly', myThid)       &     'adjoint free-surface_r-anomaly', myThid)
73        CALL MNC_CW_ADD_VATTR_TEXT('adEta',        CALL MNC_CW_ADD_VATTR_TEXT('adEta',
74       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
75    
76        CALL MNC_CW_ADD_VNAME('adW', 'Cen_xy_Hn__C__t', 4,5, myThid)        CALL MNC_CW_ADD_VNAME('adW', 'Cen_xy_Hn__C__t', 4,5, myThid)
77        CALL MNC_CW_ADD_VATTR_TEXT('adW','units','[cost]/[m/s]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adW','units','[cost]/[m/s]', myThid)
78        CALL MNC_CW_ADD_VATTR_TEXT('adW','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adW','long_name',
79       &     'adj. vertical velocity', myThid)       &     'adjoint vertical velocity', myThid)
80        CALL MNC_CW_ADD_VATTR_TEXT('adW',        CALL MNC_CW_ADD_VATTR_TEXT('adW',
81       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
82    
83        CALL MNC_CW_ADD_VNAME('adSST', 'Cen_xy_Hn__-__t', 3,4, myThid)        CALL MNC_CW_ADD_VNAME('adQnet', 'Cen_xy_Hn__-__t', 3,4, myThid)
84          CALL MNC_CW_ADD_VATTR_TEXT('adQnet',
85         &     'units','[cost]/[W/m^2]', myThid)
86          CALL MNC_CW_ADD_VATTR_TEXT('adQnet','long_name',
87         &     'adjoint net upward heat flux', myThid)
88          CALL MNC_CW_ADD_VATTR_TEXT('adQnet',
89         &     'coordinates','XC YC RC iter', myThid)
90    
91          CALL MNC_CW_ADD_VNAME('adEmpmr', 'Cen_xy_Hn__-__t', 3,4, myThid)
92          CALL MNC_CW_ADD_VATTR_TEXT('adEmpmr',
93         &     'units','[cost]/[m/s]', myThid)
94          CALL MNC_CW_ADD_VATTR_TEXT('adEmpmr','long_name',
95         &     'adjoint net upward freshwater flux', myThid)
96          CALL MNC_CW_ADD_VATTR_TEXT('adEmpmr',
97         &     'coordinates','XC YC RC iter', myThid)
98    
99          CALL MNC_CW_ADD_VNAME(    'adFu', 'Cen_xy_Hn__-__t', 3,4, myThid)
100          CALL MNC_CW_ADD_VATTR_TEXT('adFu',
101         &     'units','[cost]/[N/m^2]', myThid)
102          CALL MNC_CW_ADD_VATTR_TEXT('adFu','long_name',
103         &     'adjoint zonal wind stress', myThid)
104          CALL MNC_CW_ADD_VATTR_TEXT('adFu',
105         &     'coordinates','XC YC RC iter', myThid)
106    
107          CALL MNC_CW_ADD_VNAME(    'adFv', 'Cen_xy_Hn__-__t', 3,4, myThid)
108          CALL MNC_CW_ADD_VATTR_TEXT('adFv',
109         &     'units','[cost]/[N/m^2]', myThid)
110          CALL MNC_CW_ADD_VATTR_TEXT('adFv','long_name',
111         &     'adjoint zonal meridional stress', myThid)
112          CALL MNC_CW_ADD_VATTR_TEXT('adFv',
113         &     'coordinates','XC YC RC iter', myThid)
114    
115    #ifdef ALLOW_SST0_CONTROL
116          CALL MNC_CW_ADD_VNAME(    'adSST', 'Cen_xy_Hn__-__t', 3,4, myThid)
117        CALL MNC_CW_ADD_VATTR_TEXT('adSST',        CALL MNC_CW_ADD_VATTR_TEXT('adSST',
118       &     'units','[cost]/[degC]', myThid)       &     'units','[cost]/[degC]', myThid)
119        CALL MNC_CW_ADD_VATTR_TEXT('adSST','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adSST','long_name',
120       &     'adj. sea_surface_temperature', myThid)       &     'adjoint sea_surface_temperature', myThid)
121        CALL MNC_CW_ADD_VATTR_TEXT('adSST',        CALL MNC_CW_ADD_VATTR_TEXT('adSST',
122       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
123    #endif
124    
125        CALL MNC_CW_ADD_VNAME('adSSS', 'Cen_xy_Hn__-__t', 3,4, myThid)  #ifdef ALLOW_SSS0_CONTROL
126          CALL MNC_CW_ADD_VNAME(    'adSSS', 'Cen_xy_Hn__-__t', 3,4, myThid)
127        CALL MNC_CW_ADD_VATTR_TEXT('adSSS',        CALL MNC_CW_ADD_VATTR_TEXT('adSSS',
128       &     'units','[cost]/[degC]', myThid)       &     'units','[cost]/[psu]', myThid)
129        CALL MNC_CW_ADD_VATTR_TEXT('adSSS','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adSSS','long_name',
130       &     'adj. sea_surface_temperature', myThid)       &     'adjoint sea_surface_salinity', myThid)
131        CALL MNC_CW_ADD_VATTR_TEXT('adSSS',        CALL MNC_CW_ADD_VATTR_TEXT('adSSS',
132       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
133    #endif
134    
135  C     Write coordinates to "adstate" file  C     Write coordinates to "adstate" file
136        CALL MNC_CW_SET_UDIM('adstate', 0, myThid)        CALL MNC_CW_SET_UDIM('adstate', 0, myThid)
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22