1 |
C $Header$ |
C $Header$ |
2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
#include "CTRL_CPPOPTIONS.h" |
#include "CTRL_OPTIONS.h" |
5 |
|
|
6 |
CBOP |
CBOP |
7 |
C !ROUTINE: ctrl_depth_ini |
C !ROUTINE: ctrl_depth_ini |
10 |
|
|
11 |
C !DESCRIPTION: \bv |
C !DESCRIPTION: \bv |
12 |
c *================================================================= |
c *================================================================= |
13 |
c | SUBROUTINE ctrl_depth_ini |
c | SUBROUTINE ctrl_depth_ini |
14 |
c | Add the depth part of the control vector to the model state |
c | Add the depth part of the control vector to the model state |
15 |
c | and update the tile halos. |
c | and update the tile halos. |
16 |
c | The control vector is defined in the header file "ctrl.h". |
c | The control vector is defined in the header file "ctrl.h". |
17 |
c *================================================================= |
c *================================================================= |
18 |
C \ev |
C \ev |
81 |
fac = 0. _d 0 |
fac = 0. _d 0 |
82 |
endif |
endif |
83 |
|
|
84 |
write(standardmessageunit,'(21x,a)') |
write(standardmessageunit,'(21x,a)') |
85 |
& 'ctrl_depth_ini: ctrl update R_low,' |
& 'ctrl_depth_ini: ctrl update R_low,' |
86 |
write(standardmessageunit,'(21x,a)') |
write(standardmessageunit,'(21x,a)') |
87 |
& ' adding the control vector.' |
& ' adding the control vector.' |
88 |
|
|
89 |
C Re-initialize hFacC, so that TAMC/TAF can see it |
C Re-initialize hFacC, so that TAMC/TAF can see it |
99 |
enddo |
enddo |
100 |
enddo |
enddo |
101 |
|
|
102 |
C-- |
C-- |
103 |
do bj = jtlo,jthi |
do bj = jtlo,jthi |
104 |
do bi = itlo,ithi |
do bi = itlo,ithi |
105 |
do j = jmin,jmax |
do j = jmin,jmax |
106 |
do i = imin,imax |
do i = imin,imax |
107 |
xx_r_low(i,j,bi,bj) = 0. _d 0 |
xx_r_low(i,j,bi,bj) = 0. _d 0 |
108 |
enddo |
enddo |
109 |
enddo |
enddo |
110 |
enddo |
enddo |
120 |
do bi = itlo,ithi |
do bi = itlo,ithi |
121 |
do j = jmin,jmax |
do j = jmin,jmax |
122 |
do i = imin,imax |
do i = imin,imax |
123 |
xx_r_low(i,j,bi,bj) = r_low(i,j,bi,bj) |
xx_r_low(i,j,bi,bj) = r_low(i,j,bi,bj) |
124 |
& + fac*tmpfld2d(i,j,bi,bj) |
& + fac*tmpfld2d(i,j,bi,bj) |
125 |
enddo |
enddo |
126 |
enddo |
enddo |