29 |
#include "ctrl.h" |
#include "ctrl.h" |
30 |
#include "ctrl_dummy.h" |
#include "ctrl_dummy.h" |
31 |
#include "optim.h" |
#include "optim.h" |
32 |
|
#ifdef ALLOW_ECCO |
33 |
|
# include "ecco_cost.h" |
34 |
|
#endif |
35 |
|
|
36 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
37 |
c == routine arguments == |
c == routine arguments == |
62 |
character*( 80) fnamebottomdrag |
character*( 80) fnamebottomdrag |
63 |
|
|
64 |
_RL fac |
_RL fac |
65 |
|
_RL tmptest |
66 |
|
|
67 |
c == external == |
c == external == |
68 |
integer ilnblnk |
integer ilnblnk |
105 |
do k = 1,nr |
do k = 1,nr |
106 |
do j = jmin,jmax |
do j = jmin,jmax |
107 |
do i = imin,imax |
do i = imin,imax |
108 |
|
#ifdef ALLOW_ECCO |
109 |
|
IF (abs(tmpfld3d(i,j,k,bi,bj)).gt. |
110 |
|
$ 2.0/sqrt(wtheta(k,bi,bj))) |
111 |
|
$ tmpfld3d(i,j,k,bi,bj)= |
112 |
|
$ sign(2.0/sqrt(wtheta(k,bi,bj)),tmpfld3d(i,j,k,bi,bj)) |
113 |
|
#endif |
114 |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
115 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
116 |
if(theta(i,j,k,bi,bj).lt.-2.0) |
if(theta(i,j,k,bi,bj).lt.-2.0) |
117 |
& theta(i,j,k,bi,bj)= -2.0 |
& theta(i,j,k,bi,bj)= -2.0 |
118 |
enddo |
enddo |
119 |
enddo |
enddo |
120 |
enddo |
enddo |
121 |
enddo |
enddo |
122 |
enddo |
enddo |
123 |
|
|
124 |
#endif |
#endif |
125 |
|
|
126 |
#ifdef ALLOW_SALT0_CONTROL |
#ifdef ALLOW_SALT0_CONTROL |
137 |
do k = 1,nr |
do k = 1,nr |
138 |
do j = jmin,jmax |
do j = jmin,jmax |
139 |
do i = imin,imax |
do i = imin,imax |
140 |
|
#ifdef ALLOW_ECCO |
141 |
|
IF (abs(tmpfld3d(i,j,k,bi,bj)).gt. |
142 |
|
$ 2.0/sqrt(wsalt(k,bi,bj))) |
143 |
|
$ tmpfld3d(i,j,k,bi,bj)= |
144 |
|
$ sign(2.0/sqrt(wsalt(k,bi,bj)),tmpfld3d(i,j,k,bi,bj)) |
145 |
|
#endif |
146 |
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) + |
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) + |
147 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
148 |
|
|
149 |
enddo |
enddo |
150 |
enddo |
enddo |
151 |
enddo |
enddo |
152 |
enddo |
enddo |
153 |
enddo |
enddO |
154 |
#endif |
#endif |
155 |
|
|
156 |
#ifdef ALLOW_TR10_CONTROL |
#ifdef ALLOW_TR10_CONTROL |