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

Contents 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 - (show 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 C $Header: /u/gcmpack/MITgcm/pkg/autodiff/autodiff_ini_model_io.F,v 1.2 2004/10/10 06:08:48 edhill Exp $
2 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 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
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