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

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

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

Revision 1.3 - (hide annotations) (download)
Tue Nov 16 05:42:11 2004 UTC (19 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint56
Changes since 1.2: +16 -1 lines 
More on dsvd vs. MITgcm interfacing
o handling of g_, ad, via admtlm_vector (mds...vector)
o use ctrl_pack/unpack for admtlm_vector I/O
o use optimcycle for dsvd iteration
o make sure norm is w.r.t. derived quantities
1 heimbach 1.3 C $Header: /u/gcmpack/MITgcm/pkg/autodiff/autodiff_ini_model_io.F,v 1.2 2004/10/10 06:08:48 edhill Exp $
2 heimbach 1.1 C $Name: $
3    
4     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
5    
6     #include "PACKAGES_CONFIG.h"
7     #include "CPP_OPTIONS.h"
8    
9     CBOP
10     C !ROUTINE: AUTODIFF_INI_MODEL_IO
11     C !INTERFACE:
12     SUBROUTINE AUTODIFF_INI_MODEL_IO( myThid )
13    
14     C !DESCRIPTION: \bv
15     C autodiff_ini_model_io() is where run-time/experiment specific data are
16     C passed to any I/O packages ready that will be used for I/O of model state
17     C variables.
18     C This is specifically for setting up once only information such as
19     C shape/size of variables, units, etc... and is only for state variables.
20     C \ev
21    
22     IMPLICIT NONE
23     #include "SIZE.h"
24     #include "GRID.h"
25     #include "EEPARAMS.h"
26     #include "PARAMS.h"
27     #ifdef ALLOW_AUTODIFF_MONITOR
28     #include "adcommon.h"
29     #endif
30    
31     C !INPUT/OUTPUT PARAMETERS:
32     C myThid - Number of this instances
33     INTEGER myThid
34     CEOP
35    
36    
37     #if (defined (ALLOW_MNC) && defined (ALLOW_AUTODIFF_MONITOR))
38    
39     IF (useMNC) THEN
40    
41     C Create MNC definitions for DYNVARS.h variables
42     c CALL MNC_CW_ADD_VNAME('iter', '-_-_--__-__t', 0,0, myThid)
43     c CALL MNC_CW_ADD_VATTR_TEXT('iter',
44     c & 'long_name','iteration_count', myThid)
45    
46     c CALL MNC_CW_ADD_VNAME('model_time', '-_-_--__-__t', 0,0, myThid)
47     c CALL MNC_CW_ADD_VATTR_TEXT('model_time',
48     c & 'long_name','Model Time', myThid)
49     c CALL MNC_CW_ADD_VATTR_TEXT('model_time','units','s', myThid)
50    
51     CALL MNC_CW_ADD_VNAME('adU', 'U_xy_Hn__C__t', 4,5, myThid)
52     CALL MNC_CW_ADD_VATTR_TEXT('adU','units','[cost]/[m/s]', myThid)
53     CALL MNC_CW_ADD_VATTR_TEXT('adU','long_name',
54     & 'adj. zonal velocity', myThid)
55     CALL MNC_CW_ADD_VATTR_TEXT('adU',
56     & 'coordinates','XU YU RC iter', myThid)
57    
58     CALL MNC_CW_ADD_VNAME('adV', 'V_xy_Hn__C__t', 4,5, myThid)
59     CALL MNC_CW_ADD_VATTR_TEXT('adV','units','[cost]/[m/s]', myThid)
60     CALL MNC_CW_ADD_VATTR_TEXT('adV','long_name',
61     & 'adj. merid. velocity', myThid)
62     CALL MNC_CW_ADD_VATTR_TEXT('adV',
63     & 'coordinates','XV YV RC iter', myThid)
64    
65     CALL MNC_CW_ADD_VNAME('adT', 'Cen_xy_Hn__C__t', 4,5, myThid)
66     CALL MNC_CW_ADD_VATTR_TEXT('adT','units','[cost]/[degC]', myThid)
67     CALL MNC_CW_ADD_VATTR_TEXT('adT','long_name',
68     & 'adj. potential_temperature', myThid)
69     CALL MNC_CW_ADD_VATTR_TEXT('adT',
70     & 'coordinates','XC YC RC iter', myThid)
71    
72     CALL MNC_CW_ADD_VNAME('adS', 'Cen_xy_Hn__C__t', 4,5, myThid)
73     CALL MNC_CW_ADD_VATTR_TEXT('adS','units','[cost]/[psu]', myThid)
74     CALL MNC_CW_ADD_VATTR_TEXT('adS','long_name',
75     & 'adj. salinity', myThid)
76     CALL MNC_CW_ADD_VATTR_TEXT('adS',
77     & 'coordinates','XC YC RC iter', myThid)
78    
79     CALL MNC_CW_ADD_VNAME('adEta', 'Cen_xy_Hn__-__t', 3,4, myThid)
80     CALL MNC_CW_ADD_VATTR_TEXT('adEta','units','[cost]/[m]', myThid)
81     CALL MNC_CW_ADD_VATTR_TEXT('adEta','long_name',
82     & 'adj. free-surface_r-anomaly', myThid)
83     CALL MNC_CW_ADD_VATTR_TEXT('adEta',
84     & 'coordinates','XC YC RC iter', myThid)
85    
86     CALL MNC_CW_ADD_VNAME('adW', 'Cen_xy_Hn__C__t', 4,5, myThid)
87     CALL MNC_CW_ADD_VATTR_TEXT('adW','units','[cost]/[m/s]', myThid)
88     CALL MNC_CW_ADD_VATTR_TEXT('adW','long_name',
89     & 'adj. vertical velocity', myThid)
90     CALL MNC_CW_ADD_VATTR_TEXT('adW',
91     & 'coordinates','XC YC RC iter', myThid)
92    
93 heimbach 1.3 CALL MNC_CW_ADD_VNAME('adSST', 'Cen_xy_Hn__-__t', 3,4, myThid)
94     CALL MNC_CW_ADD_VATTR_TEXT('adSST',
95     & 'units','[cost]/[degC]', myThid)
96     CALL MNC_CW_ADD_VATTR_TEXT('adSST','long_name',
97     & 'adj. sea_surface_temperature', myThid)
98     CALL MNC_CW_ADD_VATTR_TEXT('adSST',
99     & 'coordinates','XC YC RC iter', myThid)
100    
101     CALL MNC_CW_ADD_VNAME('adSSS', 'Cen_xy_Hn__-__t', 3,4, myThid)
102     CALL MNC_CW_ADD_VATTR_TEXT('adSSS',
103     & 'units','[cost]/[degC]', myThid)
104     CALL MNC_CW_ADD_VATTR_TEXT('adSSS','long_name',
105     & 'adj. sea_surface_temperature', myThid)
106     CALL MNC_CW_ADD_VATTR_TEXT('adSSS',
107     & 'coordinates','XC YC RC iter', myThid)
108 heimbach 1.1
109     C Write coordinates to "adstate" file
110     CALL MNC_CW_SET_UDIM('adstate', 0, myThid)
111     CALL MNC_CW_RS_W('R','adstate',0,0,'XC',xC, myThid)
112     CALL MNC_CW_RS_W('R','adstate',0,0,'YC',yC, myThid)
113     CALL MNC_CW_RS_W('R','adstate',0,0,'XU',xG, myThid)
114     CALL MNC_CW_RS_W('R','adstate',0,0,'YU',yC, myThid)
115     CALL MNC_CW_RS_W('R','adstate',0,0,'XV',xC, myThid)
116     CALL MNC_CW_RS_W('R','adstate',0,0,'YV',yG, myThid)
117     CALL MNC_CW_RS_W('R','adstate',0,0,'XG',xG, myThid)
118     CALL MNC_CW_RS_W('R','adstate',0,0,'YG',yG, myThid)
119     CALL MNC_CW_RS_W('R','adstate',0,0,'RC',rC, myThid)
120     CALL MNC_CW_RS_W('R','adstate',0,0,'RF',rF, myThid)
121    
122     ENDIF
123     #endif
124    
125     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
126    
127     RETURN
128     END
  ViewVC Help
Powered by ViewVC 1.1.22