/[MITgcm]/MITgcm/pkg/autodiff/adcommon.h
ViewVC logotype

Diff of /MITgcm/pkg/autodiff/adcommon.h

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

revision 1.3 by heimbach, Fri Jan 11 18:52:13 2002 UTC revision 1.31 by jmc, Thu Aug 2 01:13:29 2012 UTC
# Line 1  Line 1 
1     These lines are uncommented on purpose to provoke  C $Header$
2     a model crash. You need to customize this header file  C $Name$
    to your problem.  
    These common blocks are extracted from the  
    automatically created tangent linear code.  
    You need to make sure that they are up-to-date  
    (i.e. in right order), and customize them  
    accordingly.  
3    
4     heimbach@mit.edu 11-Jan-2001  C--   These common blocks are extracted from the
5    C--   automatically created tangent linear code.
6    C--   You need to make sure that they are up-to-date
7    C--   (i.e. in right order), and customize them accordingly.
8    C--
9    C--   heimbach@mit.edu 11-Jan-2001
10    
11        common /addynvars_r/  #ifdef ALLOW_AUTODIFF_MONITOR
12    
13          common /addynvars_r/
14       &                     adetan,       &                     adetan,
15       &                     aduvel, advvel, adwvel,       &                     aduvel, advvel, adwvel,
16       &                     adtheta, adsalt,       &                     adtheta, adsalt,
17       &                     adgu, adgv, adgt, adgs,       &                     adgu, adgv, adgt, adgs,
18    #ifdef ALLOW_ADAMSBASHFORTH_3
19         &                     adgunm, adgvnm, adgtnm, adgsnm
20    #else
21       &                     adgunm1, adgvnm1, adgtnm1, adgsnm1       &                     adgunm1, adgvnm1, adgtnm1, adgsnm1
22    #endif
23        _RL adetan(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)        _RL adetan(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
24        _RL adgs(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL adgs(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
       _RL adgsnm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
25        _RL adgt(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL adgt(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
       _RL adgtnm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
26        _RL adgu(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL adgu(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
       _RL adgunm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
27        _RL adgv(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL adgv(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
       _RL adgvnm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
28        _RL adsalt(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL adsalt(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
29        _RL adtheta(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL adtheta(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
30        _RL aduvel(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL aduvel(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
31        _RL advvel(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL advvel(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
32        _RL adwvel(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL adwvel(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
33    #ifdef ALLOW_ADAMSBASHFORTH_3
34          _RL adgtnm(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy,2)
35          _RL adgsnm(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy,2)
36          _RL adgunm(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy,2)
37          _RL adgvnm(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy,2)
38    #else
39          _RL adgtnm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
40          _RL adgsnm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
41          _RL adgunm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
42          _RL adgvnm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
43    #endif
44    
45        common /addynvars_cd/        common /addynvars_r_2/
46         &                     adetah
47          _RL adetah(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
48    
49    #ifdef ALLOW_AUTODIFF_MONITOR_DIAG
50          common /addynvars_diag/
51         &                     adtotphihyd, adrhoinsitu
52          _RL adrhoinsitu(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
53          _RL adtotphihyd(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
54    #endif
55    
56    #ifdef ALLOW_CD_CODE
57          common /addynvars_cd/
58       &                      aduveld, advveld,       &                      aduveld, advveld,
59       &                       adetanm1,       &                       adetanm1,
60       &                      adunm1, advnm1,       &                      adunm1, advnm1
      &                      adgucd, adgvcd  
61        _RL aduveld(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL aduveld(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
62        _RL advveld(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL advveld(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
63        _RL adetanm1(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)        _RL adetanm1(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
64        _RL adunm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL adunm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
65        _RL advnm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)        _RL advnm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
66        _RL adgucd(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  #endif
67        _RL adgvcd(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
68          COMMON /adffields_fu/ adfu
69          COMMON /adffields_fv/ adfv
70          COMMON /adffields_Qnet/ adQnet
71          COMMON /adffields_Qsw/ adQsw
72          COMMON /adffields_dQdT/ addQdT
73          COMMON /adffields_EmPmR/ adEmPmR
74          COMMON /adffields_saltFlux/ adsaltFlux
75          COMMON /adffields_SST/ adSST
76          COMMON /adffields_SSS/ adSSS
77          COMMON /adffields_lambdaThetaClimRelax/ adlambdaThetaClimRelax
78          COMMON /adffields_lambdaSaltClimRelax/ adlambdaSaltClimRelax
79    #ifdef ATMOSPHERIC_LOADING
80          COMMON /adffields_pload/ adpload
81          COMMON /adffields_sIceLoad/ adsIceLoad
82    #endif
83    
84          _RS  adfu       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
85          _RS  adfv       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86          _RS  adQnet     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87          _RS  adQsw      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88          _RS  addQdT     (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
89          _RS  adEmPmR    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90          _RS  adsaltFlux (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
91          _RS  adSST      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92          _RS  adSSS      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
93          _RS  adlambdaThetaClimRelax
94         &    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
95          _RS  adlambdaSaltClimRelax
96         &    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
97    #ifdef ATMOSPHERIC_LOADING
98          _RS  adpload    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
99          _RS  adsIceLoad (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
100    #endif
101    
102    #ifdef ALLOW_DIFFKR_CONTROL
103          COMMON /ADDYNVARS_DIFFKR/
104         &                       addiffKr
105          _RL  addiffKr (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
106    #endif
107    #ifdef ALLOW_KAPGM_CONTROL
108          COMMON /ADDYNVARS_KAPGM/
109         &                       adkapgm
110          _RL  adkapgm (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
111    #endif
112    #ifdef ALLOW_KAPREDI_CONTROL
113          COMMON /ADDYNVARS_KAPREDI/
114         &                       adkapredi
115          _RL  adkapredi (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
116    #endif
117    #ifdef ALLOW_BOTTOMDRAG_CONTROL
118          COMMON /ADDYNVARS_BOTTOMDRAG/
119         &                adbottomdragfld
120          _RL  adbottomdragfld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
121    #endif
122    #ifdef ALLOW_EDDYPSI_CONTROL
123          _RS adEddyPsiX(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
124          _RS adEddyPsiY(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
125          common /adeddypsiffields/ adeddypsix, adeddypsiy
126    #endif
127    
128    #ifdef ALLOW_EXF
129    
130          _RL adhflux(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
131          _RL adsflux(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
132          common /adexf_hsflux_r/ adhflux, adsflux
133    
134          _RL adustress(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
135          _RL advstress(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
136          common /adexf_stress_r/ adustress, advstress
137    
138          _RL adwspeed(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
139          common /adexf_wspeed_r/ adwspeed
140    
141    # ifdef ALLOW_ATM_TEMP
142          _RL adatemp     (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
143          _RL adaqh       (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
144          _RL adlwflux    (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
145          _RL adprecip    (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
146          common /adexf_atm_temp_r/ adatemp, adaqh, adlwflux,
147         & adprecip
148    #  ifdef SHORTWAVE_HEATING
149          _RL adswflux    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
150          common /adexf_swflux_r/ adswflux
151    #  endif
152    # endif
153    # ifdef ALLOW_ATM_WIND
154          _RL aduwind     (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
155          _RL advwind     (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
156          common /adexf_atm_wind_r/ aduwind, advwind
157    # endif
158    # ifdef ALLOW_DOWNWARD_RADIATION
159          _RL adswdown    (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
160          _RL adlwdown    (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
161          common /adexf_rad_down_r/
162         &     adswdown, adlwdown
163    # endif
164    # ifdef ALLOW_CLIMSST_RELAXATION
165          _RL adclimsst(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
166          common /adexf_clim_sst_r/ adclimsst
167    # endif
168    # ifdef ALLOW_CLIMSSS_RELAXATION
169          _RL adclimsss(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
170          common /adexf_clim_sss_r/ adclimsss
171    # endif
172    
173    #endif /* ALLOW_EXF */
174    
175    #ifdef ALLOW_SEAICE
176          _RL adarea  (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
177          _RL adheff  (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
178          _RL adhsnow (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
179          _RL aduice  (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
180          _RL advice  (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
181          common /adseaice_dynvars_1/
182         &     adarea, adheff, adhsnow, aduice, advice
183    # ifdef SEAICE_VARIABLE_SALINITY
184          _RL adhsalt (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
185          common /seaice_salinity_r/
186         &     adhsalt
187    # endif
188    #endif /* ALLOW_SEAICE */
189    
190    #ifdef ALLOW_GGL90
191          _RL adggl90tke     (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
192          _RL adggl90diffkr  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
193          common /adggl90_fields/
194         &     adggl90tke, adggl90diffkr
195    #endif
196    
197    #ifdef ALLOW_DEPTH_CONTROL
198          _RL adr_low_control(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
199          _RL adhfacc(1-olx:snx+olx,1-oly:sny+oly,1:nr,nsx,nsy)
200          _RL adhfacs(1-olx:snx+olx,1-oly:sny+oly,1:nr,nsx,nsy)
201          _RL adhfacw(1-olx:snx+olx,1-oly:sny+oly,1:nr,nsx,nsy)
202          _RL adrecip_hfacc(1-olx:snx+olx,1-oly:sny+oly,1:nr,nsx,nsy)
203          _RL adrecip_hfacs(1-olx:snx+olx,1-oly:sny+oly,1:nr,nsx,nsy)
204          _RL adrecip_hfacw(1-olx:snx+olx,1-oly:sny+oly,1:nr,nsx,nsy)
205          common /adgrid_r/ adr_low_control,
206         &                  adrecip_hfacc, adrecip_hfacw, adrecip_hfacs
207          common /adgrid_r_c/ adhfacc
208          common /adgrid_r_s/ adhfacs
209          common /adgrid_r_w/ adhfacw
210    #endif /* ALLOW_DEPTH_CONTROL */
211    
212  cph      common /adtr1_r/  #endif /* ALLOW_AUTODIFF_MONITOR */
 cph     &                 adtr1, adgtr1, adgtr1nm1  
 cph      _RL adgtr1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
 cph      _RL adgtr1nm1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
 cph      _RL adtr1(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
   
       common /adffields/  
      &                   adfu, adfv  
      &                 , adqnet, adempmr  
 cph     &                 , adsst, adsss  
       _RL adfu(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)  
       _RL adfv(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)  
       _RL adqnet(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)  
       _RL adempmr(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)  
 cph      _RL adsst(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)  
 cph      _RL adsss(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)  
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.31
  ViewVC Help
Powered by ViewVC 1.1.22