/[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.2 by edhill, Sun Oct 10 06:08:48 2004 UTC revision 1.6 by jmc, Wed Apr 27 14:10:06 2005 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
5    
6  #include "PACKAGES_CONFIG.h"  #include "AUTODIFF_OPTIONS.h"
 #include "CPP_OPTIONS.h"  
7    
8  CBOP  CBOP
9  C     !ROUTINE: AUTODIFF_INI_MODEL_IO  C     !ROUTINE: AUTODIFF_INI_MODEL_IO
# Line 38  CEOP Line 37  CEOP
37    
38        IF (useMNC) THEN        IF (useMNC) THEN
39    
 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)  
   
40        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)
41        CALL MNC_CW_ADD_VATTR_TEXT('adU','units','[cost]/[m/s]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adU','units','[cost]/[m/s]', myThid)
42        CALL MNC_CW_ADD_VATTR_TEXT('adU','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adU','long_name',
43       &     'adj. zonal velocity', myThid)       &     'adjoint zonal velocity', myThid)
44        CALL MNC_CW_ADD_VATTR_TEXT('adU',        CALL MNC_CW_ADD_VATTR_TEXT('adU',
45       &     'coordinates','XU YU RC iter', myThid)       &     'coordinates','XU YU RC iter', myThid)
46    
47        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)
48        CALL MNC_CW_ADD_VATTR_TEXT('adV','units','[cost]/[m/s]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adV','units','[cost]/[m/s]', myThid)
49        CALL MNC_CW_ADD_VATTR_TEXT('adV','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adV','long_name',
50       &     'adj. merid. velocity', myThid)       &     'adjoint merid. velocity', myThid)
51        CALL MNC_CW_ADD_VATTR_TEXT('adV',        CALL MNC_CW_ADD_VATTR_TEXT('adV',
52       &     'coordinates','XV YV RC iter', myThid)       &     'coordinates','XV YV RC iter', myThid)
53    
54        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)
55        CALL MNC_CW_ADD_VATTR_TEXT('adT','units','[cost]/[degC]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adT','units','[cost]/[degC]', myThid)
56        CALL MNC_CW_ADD_VATTR_TEXT('adT','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adT','long_name',
57       &     'adj. potential_temperature', myThid)       &     'adjoint potential_temperature', myThid)
58        CALL MNC_CW_ADD_VATTR_TEXT('adT',        CALL MNC_CW_ADD_VATTR_TEXT('adT',
59       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
60    
61        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)
62        CALL MNC_CW_ADD_VATTR_TEXT('adS','units','[cost]/[psu]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adS','units','[cost]/[psu]', myThid)
63        CALL MNC_CW_ADD_VATTR_TEXT('adS','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adS','long_name',
64       &     'adj. salinity', myThid)       &     'adjoint salinity', myThid)
65        CALL MNC_CW_ADD_VATTR_TEXT('adS',        CALL MNC_CW_ADD_VATTR_TEXT('adS',
66       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
67    
68        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)
69        CALL MNC_CW_ADD_VATTR_TEXT('adEta','units','[cost]/[m]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adEta','units','[cost]/[m]', myThid)
70        CALL MNC_CW_ADD_VATTR_TEXT('adEta','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adEta','long_name',
71       &     'adj. free-surface_r-anomaly', myThid)       &     'adjoint free-surface_r-anomaly', myThid)
72        CALL MNC_CW_ADD_VATTR_TEXT('adEta',        CALL MNC_CW_ADD_VATTR_TEXT('adEta',
73       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
74    
75        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)
76        CALL MNC_CW_ADD_VATTR_TEXT('adW','units','[cost]/[m/s]', myThid)        CALL MNC_CW_ADD_VATTR_TEXT('adW','units','[cost]/[m/s]', myThid)
77        CALL MNC_CW_ADD_VATTR_TEXT('adW','long_name',        CALL MNC_CW_ADD_VATTR_TEXT('adW','long_name',
78       &     'adj. vertical velocity', myThid)       &     'adjoint vertical velocity', myThid)
79        CALL MNC_CW_ADD_VATTR_TEXT('adW',        CALL MNC_CW_ADD_VATTR_TEXT('adW',
80       &     'coordinates','XC YC RC iter', myThid)       &     'coordinates','XC YC RC iter', myThid)
81    
82          CALL MNC_CW_ADD_VNAME('adQnet', 'Cen_xy_Hn__-__t', 3,4, myThid)
83          CALL MNC_CW_ADD_VATTR_TEXT('adQnet',
84         &     'units','[cost]/[W/m^2]', myThid)
85          CALL MNC_CW_ADD_VATTR_TEXT('adQnet','long_name',
86         &     'adjoint net upward heat flux', myThid)
87          CALL MNC_CW_ADD_VATTR_TEXT('adQnet',
88         &     'coordinates','XC YC RC iter', myThid)
89    
90          CALL MNC_CW_ADD_VNAME('adEmpmr', 'Cen_xy_Hn__-__t', 3,4, myThid)
91          CALL MNC_CW_ADD_VATTR_TEXT('adEmpmr',
92         &     'units','[cost]/[m/s]', myThid)
93          CALL MNC_CW_ADD_VATTR_TEXT('adEmpmr','long_name',
94         &     'adjoint net upward freshwater flux', myThid)
95          CALL MNC_CW_ADD_VATTR_TEXT('adEmpmr',
96         &     'coordinates','XC YC RC iter', myThid)
97    
98          CALL MNC_CW_ADD_VNAME(    'adFu', 'Cen_xy_Hn__-__t', 3,4, myThid)
99          CALL MNC_CW_ADD_VATTR_TEXT('adFu',
100         &     'units','[cost]/[N/m^2]', myThid)
101          CALL MNC_CW_ADD_VATTR_TEXT('adFu','long_name',
102         &     'adjoint zonal wind stress', myThid)
103          CALL MNC_CW_ADD_VATTR_TEXT('adFu',
104         &     'coordinates','XC YC RC iter', myThid)
105    
106          CALL MNC_CW_ADD_VNAME(    'adFv', 'Cen_xy_Hn__-__t', 3,4, myThid)
107          CALL MNC_CW_ADD_VATTR_TEXT('adFv',
108         &     'units','[cost]/[N/m^2]', myThid)
109          CALL MNC_CW_ADD_VATTR_TEXT('adFv','long_name',
110         &     'adjoint zonal meridional stress', myThid)
111          CALL MNC_CW_ADD_VATTR_TEXT('adFv',
112         &     'coordinates','XC YC RC iter', myThid)
113    
114    #ifdef ALLOW_SST0_CONTROL
115          CALL MNC_CW_ADD_VNAME(    'adSST', 'Cen_xy_Hn__-__t', 3,4, myThid)
116          CALL MNC_CW_ADD_VATTR_TEXT('adSST',
117         &     'units','[cost]/[degC]', myThid)
118          CALL MNC_CW_ADD_VATTR_TEXT('adSST','long_name',
119         &     'adjoint sea_surface_temperature', myThid)
120          CALL MNC_CW_ADD_VATTR_TEXT('adSST',
121         &     'coordinates','XC YC RC iter', myThid)
122    #endif
123    
124    #ifdef ALLOW_SSS0_CONTROL
125          CALL MNC_CW_ADD_VNAME(    'adSSS', 'Cen_xy_Hn__-__t', 3,4, myThid)
126          CALL MNC_CW_ADD_VATTR_TEXT('adSSS',
127         &     'units','[cost]/[psu]', myThid)
128          CALL MNC_CW_ADD_VATTR_TEXT('adSSS','long_name',
129         &     'adjoint sea_surface_salinity', myThid)
130          CALL MNC_CW_ADD_VATTR_TEXT('adSSS',
131         &     'coordinates','XC YC RC iter', myThid)
132    #endif
133    
134  C     Write coordinates to "adstate" file  C     Write coordinates to "adstate" file
135        CALL MNC_CW_SET_UDIM('adstate', 0, myThid)        CALL MNC_CW_SET_UDIM('adstate', 0, myThid)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22