2 |
|
|
3 |
#include "COST_CPPOPTIONS.h" |
#include "COST_CPPOPTIONS.h" |
4 |
|
|
|
|
|
5 |
subroutine cost_averagesinit( mythid ) |
subroutine cost_averagesinit( mythid ) |
6 |
|
|
7 |
c ================================================================== |
c ================================================================== |
32 |
|
|
33 |
#include "ecco_cost.h" |
#include "ecco_cost.h" |
34 |
#include "ctrl_dummy.h" |
#include "ctrl_dummy.h" |
35 |
|
#ifdef ALLOW_SEAICE |
36 |
|
# include "SEAICE_COST.h" |
37 |
|
#endif |
38 |
|
|
39 |
c == routine arguments == |
c == routine arguments == |
40 |
|
|
76 |
#endif |
#endif |
77 |
|
|
78 |
#if (defined (ALLOW_THETA_COST_CONTRIBUTION) || \ |
#if (defined (ALLOW_THETA_COST_CONTRIBUTION) || \ |
|
defined (ALLOW_SST_COST_CONTRIBUTION) || \ |
|
79 |
defined (ALLOW_CTDT_COST_CONTRIBUTION) || \ |
defined (ALLOW_CTDT_COST_CONTRIBUTION) || \ |
|
defined (ALLOW_CTDTCLIM_COST_CONTRIBUTION) || \ |
|
80 |
defined (ALLOW_XBT_COST_CONTRIBUTION) || \ |
defined (ALLOW_XBT_COST_CONTRIBUTION) || \ |
81 |
defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \ |
defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \ |
82 |
defined (ALLOW_OBCS_COST_CONTRIBUTION)) |
defined (ALLOW_OBCS_COST_CONTRIBUTION)) |
89 |
enddo |
enddo |
90 |
enddo |
enddo |
91 |
xx_tbar_mean_dummy = 0. _d 0 |
xx_tbar_mean_dummy = 0. _d 0 |
92 |
|
#else |
93 |
|
#ifdef ALLOW_SST_COST_CONTRIBUTION |
94 |
|
c-- Initialise temperature average (2d). |
95 |
|
k = 1 |
96 |
|
do j = jmin,jmax |
97 |
|
do i = imin,imax |
98 |
|
tbar(i,j,bi,bj) = 0. _d 0 |
99 |
|
enddo |
100 |
|
enddo |
101 |
|
xx_tbar_mean_dummy = 0. _d 0 |
102 |
|
#endif |
103 |
#endif |
#endif |
104 |
|
|
105 |
#if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \ |
#if (defined (ALLOW_SALT_COST_CONTRIBUTION) || \ |
|
defined (ALLOW_SSS_COST_CONTRIBUTION) || \ |
|
106 |
defined (ALLOW_CTDS_COST_CONTRIBUTION) || \ |
defined (ALLOW_CTDS_COST_CONTRIBUTION) || \ |
|
defined (ALLOW_CTDSCLIM_COST_CONTRIBUTION) || \ |
|
107 |
defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \ |
defined (ALLOW_DRIFT_COST_CONTRIBUTION) || \ |
108 |
defined (ALLOW_OBCS_COST_CONTRIBUTION)) |
defined (ALLOW_OBCS_COST_CONTRIBUTION)) |
109 |
c-- Initialise salt average. |
c-- Initialise salt average. |
115 |
enddo |
enddo |
116 |
enddo |
enddo |
117 |
xx_sbar_mean_dummy = 0. _d 0 |
xx_sbar_mean_dummy = 0. _d 0 |
118 |
|
#else |
119 |
|
#ifdef ALLOW_SSS_COST_CONTRIBUTION |
120 |
|
c-- Initialise salinity average (2d). |
121 |
|
k = 1 |
122 |
|
do j = jmin,jmax |
123 |
|
do i = imin,imax |
124 |
|
sbar(i,j,bi,bj) = 0. _d 0 |
125 |
|
enddo |
126 |
|
enddo |
127 |
|
xx_sbar_mean_dummy = 0. _d 0 |
128 |
|
#endif |
129 |
#endif |
#endif |
130 |
|
|
131 |
#if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \ |
#if (defined (ALLOW_DRIFTER_COST_CONTRIBUTION) || \ |
132 |
defined (ALLOW_CURMTR_COST_CONTRIBUTION) || \ |
defined (ALLOW_OBCS_COST_CONTRIBUTION)) |
|
defined (OBCS_AGEOS_COST_CONTRIBUTION)) |
|
133 |
c-- Initialise uvel, vvel average. |
c-- Initialise uvel, vvel average. |
134 |
do k = 1,nr |
do k = 1,nr |
135 |
do j = jmin,jmax |
do j = jmin,jmax |
185 |
enddo |
enddo |
186 |
enddo |
enddo |
187 |
xx_sflux_mean_dummy = 0. _d 0 |
xx_sflux_mean_dummy = 0. _d 0 |
188 |
#endif |
#endif |
189 |
|
|
190 |
|
#ifdef ALLOW_SEAICE |
191 |
|
|
192 |
|
# ifdef ALLOW_SEAICE_COST_SMR_AREA |
193 |
|
c-- Initialise salinity average (2d). |
194 |
|
do j = jmin,jmax |
195 |
|
do i = imin,imax |
196 |
|
smrareabar(i,j,bi,bj) = 0. _d 0 |
197 |
|
enddo |
198 |
|
enddo |
199 |
|
xx_smrareabar_mean_dummy = 0. _d 0 |
200 |
|
# endif |
201 |
|
|
202 |
|
#endif |
203 |
|
|
204 |
enddo |
enddo |
205 |
enddo |
enddo |
206 |
|
|
|
#ifdef ALLOW_COST_INI_FIN |
|
|
xx_theta_ini_fin_dummy = 0. _d 0 |
|
|
xx_salt_ini_fin_dummy = 0. _d 0 |
|
|
#endif |
|
|
|
|
207 |
return |
return |
208 |
end |
end |
209 |
|
|