47 |
C !LOCAL VARIABLES: ==================================================== |
C !LOCAL VARIABLES: ==================================================== |
48 |
C i,j :: loop indices |
C i,j :: loop indices |
49 |
C scalingFactor :: scaling factor (from acceleration to density) |
C scalingFactor :: scaling factor (from acceleration to density) |
50 |
|
INTEGER iMin,iMax,jMin,jMax |
51 |
|
PARAMETER( iMin = 0 , iMax = sNx+1 ) |
52 |
|
PARAMETER( jMin = 0 , jMax = sNy+1 ) |
53 |
INTEGER i,j |
INTEGER i,j |
54 |
_RL scalingFactor |
_RL scalingFactor |
55 |
CEOP |
CEOP |
59 |
c scalingFactor=gravitySign*recip_gravity |
c scalingFactor=gravitySign*recip_gravity |
60 |
|
|
61 |
IF ( use3dCoriolis ) THEN |
IF ( use3dCoriolis ) THEN |
62 |
DO j=1-Oly,sNy+Oly-1 |
DO j=jMin,jMax |
63 |
DO i=1-Olx,sNx+Olx-1 |
DO i=iMin,iMax |
64 |
effectiveBuoy(i,j)=effectiveBuoy(i,j) |
effectiveBuoy(i,j)=effectiveBuoy(i,j) |
65 |
& +scalingFactor* |
& +scalingFactor* |
66 |
& fCoriCos(i,j,bi,bj)* |
& fCoriCos(i,j,bi,bj)* |
74 |
ENDIF |
ENDIF |
75 |
|
|
76 |
IF ( useNHMTerms ) THEN |
IF ( useNHMTerms ) THEN |
77 |
DO j=1-Oly,sNy+Oly-1 |
DO j=jMin,jMax |
78 |
DO i=1-Olx,sNx+Olx-1 |
DO i=iMin,iMax |
79 |
effectiveBuoy(i,j)=effectiveBuoy(i,j) |
effectiveBuoy(i,j)=effectiveBuoy(i,j) |
80 |
& +scalingFactor* |
& +scalingFactor* |
81 |
& ( (uFld( i ,j,k,bi,bj)*uFld( i ,j,k,bi,bj) |
& ( (uFld( i ,j,k,bi,bj)*uFld( i ,j,k,bi,bj) |