11 |
I bi,bj,k, |
I bi,bj,k, |
12 |
I uFld, |
I uFld, |
13 |
O vMetricTerms, |
O vMetricTerms, |
14 |
I myThid) |
I myThid ) |
15 |
|
|
16 |
C !DESCRIPTION: |
C !DESCRIPTION: |
17 |
C Calculates the meridional metric term due to sphericity: |
C Calculates the meridional metric term due to sphericity: |
29 |
C !INPUT PARAMETERS: =================================================== |
C !INPUT PARAMETERS: =================================================== |
30 |
C bi,bj :: tile indices |
C bi,bj :: tile indices |
31 |
C k :: vertical level |
C k :: vertical level |
32 |
C uFld :: zonal flow |
C uFld :: zonal flow |
33 |
C myThid :: thread number |
C myThid :: thread number |
34 |
INTEGER bi,bj,k |
INTEGER bi,bj,k |
35 |
_RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
46 |
|
|
47 |
DO j=1-Olx+1,sNy+Oly |
DO j=1-Olx+1,sNy+Oly |
48 |
DO i=1-Olx,sNx+Olx-1 |
DO i=1-Olx,sNx+Olx-1 |
49 |
vMetricTerms(i,j) = -recip_RSphere |
vMetricTerms(i,j) = -recip_rSphere*recip_deepFacC(k) |
50 |
& *0.25*( uFld(i,j )+uFld(i+1,j ) |
& *0.25*( uFld(i,j )+uFld(i+1,j ) |
51 |
& +uFld(i,j-1)+uFld(i+1,j-1) |
& +uFld(i,j-1)+uFld(i+1,j-1) |
52 |
& ) |
& ) |