108 |
$ tmpfld3d(i,j,k,bi,bj)= |
$ tmpfld3d(i,j,k,bi,bj)= |
109 |
$ sign(2.0/sqrt(wtheta(k,bi,bj)),tmpfld3d(i,j,k,bi,bj)) |
$ sign(2.0/sqrt(wtheta(k,bi,bj)),tmpfld3d(i,j,k,bi,bj)) |
110 |
#endif |
#endif |
111 |
|
#ifdef ALLOW_OPENAD |
112 |
|
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
113 |
|
& fac*xx_theta(i,j,k,bi,bj) |
114 |
|
#else |
115 |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) + |
116 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
117 |
|
#endif |
118 |
if(theta(i,j,k,bi,bj).lt.-2.0) |
if(theta(i,j,k,bi,bj).lt.-2.0) |
119 |
& theta(i,j,k,bi,bj)= -2.0 |
& theta(i,j,k,bi,bj)= -2.0 |
120 |
enddo |
enddo |
145 |
$ tmpfld3d(i,j,k,bi,bj)= |
$ tmpfld3d(i,j,k,bi,bj)= |
146 |
$ sign(2.0/sqrt(wsalt(k,bi,bj)),tmpfld3d(i,j,k,bi,bj)) |
$ sign(2.0/sqrt(wsalt(k,bi,bj)),tmpfld3d(i,j,k,bi,bj)) |
147 |
#endif |
#endif |
148 |
|
#ifdef ALLOW_OPENAD |
149 |
|
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) + |
150 |
|
& fac*xx_salt(i,j,k,bi,bj) |
151 |
|
#else |
152 |
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) + |
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) + |
153 |
& fac*tmpfld3d(i,j,k,bi,bj) |
& fac*tmpfld3d(i,j,k,bi,bj) |
154 |
|
#endif |
155 |
|
|
156 |
enddo |
enddo |
157 |
enddo |
enddo |