32 |
#include "optim.h" |
#include "optim.h" |
33 |
#ifdef ALLOW_PTRACERS |
#ifdef ALLOW_PTRACERS |
34 |
# include "PTRACERS_SIZE.h" |
# include "PTRACERS_SIZE.h" |
35 |
# include "PTRACERS.h" |
c#include "PTRACERS_PARAMS.h" |
36 |
|
# include "PTRACERS_FIELDS.h" |
37 |
#endif |
#endif |
38 |
#ifdef ALLOW_ECCO |
#ifdef ALLOW_ECCO |
39 |
# include "ecco_cost.h" |
# include "ecco_cost.h" |
118 |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
119 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
120 |
#endif |
#endif |
121 |
|
#ifndef DISABLE_CTRL_THETA_LIMIT |
122 |
if(theta(i,j,k,bi,bj).lt.-2.0) |
if(theta(i,j,k,bi,bj).lt.-2.0) |
123 |
& theta(i,j,k,bi,bj)= -2.0 |
& theta(i,j,k,bi,bj)= -2.0 |
124 |
|
#endif |
125 |
enddo |
enddo |
126 |
enddo |
enddo |
127 |
enddo |
enddo |