73 |
K=1 |
K=1 |
74 |
|
|
75 |
C calculate SCHMIDT NO. for O2 |
C calculate SCHMIDT NO. for O2 |
76 |
DO j=1-OLy,sNy+OLy |
DO j=jmin,jmax |
77 |
DO i=1-OLx,sNx+OLx |
DO i=imin,imax |
78 |
IF (hFacC(i,j,k,bi,bj).NE.0.) THEN |
IF (hFacC(i,j,k,bi,bj).NE.0.) THEN |
79 |
SchmidtNoO2(i,j) = |
SchmidtNoO2(i,j) = |
80 |
& sox1 |
& sox1 |
127 |
END DO |
END DO |
128 |
|
|
129 |
C update surface tendencies |
C update surface tendencies |
130 |
DO j=1-OLy,sNy+OLy |
DO j=jmin,jmax |
131 |
DO i=1-OLx,sNx+OLx |
DO i=imin,imax |
132 |
SGO2(i,j)= maskC(i,j,1,bi,bj)*FluxO2(i,j) |
SGO2(i,j)= maskC(i,j,1,bi,bj)*FluxO2(i,j) |
133 |
& *recip_drF(1) * recip_hFacC(i,j,1,bi,bj) |
& *recip_drF(1) * recip_hFacC(i,j,1,bi,bj) |
134 |
ENDDO |
ENDDO |