7 |
C !ROUTINE: GAD_DEL2 |
C !ROUTINE: GAD_DEL2 |
8 |
|
|
9 |
C !INTERFACE: ========================================================== |
C !INTERFACE: ========================================================== |
10 |
SUBROUTINE GAD_DEL2( |
SUBROUTINE GAD_DEL2( |
11 |
I bi,bj,k, |
I bi,bj,k, |
12 |
I dTdx,dTdy, |
I dTdx,dTdy, |
13 |
O del2, |
O del2, |
42 |
C !LOCAL VARIABLES: ==================================================== |
C !LOCAL VARIABLES: ==================================================== |
43 |
C i,j :: loop indices |
C i,j :: loop indices |
44 |
INTEGER i,j |
INTEGER i,j |
|
c _RL dfx(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
|
|
c _RL dfy(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
|
45 |
CEOP |
CEOP |
46 |
|
|
47 |
C-- Initialise to zero. |
C-- Initialise to zero. |
55 |
C and divide by grid-cell volume |
C and divide by grid-cell volume |
56 |
DO j=1-Oly,sNy+Oly-1 |
DO j=1-Oly,sNy+Oly-1 |
57 |
DO i=1-Olx,sNx+Olx-1 |
DO i=1-Olx,sNx+Olx-1 |
58 |
del2(i,j)=recip_rA(i,j,bi,bj) |
del2(i,j)=recip_rA(i,j,bi,bj)*recip_deepFacC(k) |
59 |
& *recip_drF(k) |
& *recip_drF(k)*_recip_hFacC(i,j,k,bi,bj) |
|
& *_recip_hFacC(i,j,k,bi,bj) |
|
60 |
& *( ( dTdx(i+1,j)-dTdx(i,j) ) |
& *( ( dTdx(i+1,j)-dTdx(i,j) ) |
61 |
& +( dTdy(i,j+1)-dTdy(i,j) ) |
& +( dTdy(i,j+1)-dTdy(i,j) ) |
62 |
& ) |
& ) |