897 |
DO bi=1,nSx |
DO bi=1,nSx |
898 |
DO j=1,sNy |
DO j=1,sNy |
899 |
DO i=1,sNx |
DO i=1,sNx |
900 |
fv(i,j,bi,bj) = Qsw(i,j,bi,bj) - |
Qnet(i,j,bi,bj) = Qsw(i,j,bi,bj) - |
901 |
& AREA(i,j,bi,bj) * local(i,j,bi,bj) + |
& AREA(i,j,bi,bj) * local(i,j,bi,bj) + |
902 |
& (1.-AREA(i,j,bi,bj)) * Qnet_tmp(i,j,bi,bj) |
& (1.-AREA(i,j,bi,bj)) * Qnet_tmp(i,j,bi,bj) |
903 |
ENDDO |
ENDDO |
930 |
DO bi=1,nSx |
DO bi=1,nSx |
931 |
DO j=1,sNy |
DO j=1,sNy |
932 |
DO i=1,sNx |
DO i=1,sNx |
933 |
EmPmR(i,j,bi,bj) = - rhoConstFresh * |
EmPmR(i,j,bi,bj) = - AREA(i,j,bi,bj) * local (i,j,bi,bj) + |
934 |
& AREA(i,j,bi,bj) * local(i,j,bi,bj) + |
& ( 1. - AREA(i,j,bi,bj)) * EmPmR_tmp(i,j,bi,bj) |
|
& (1.-AREA(i,j,bi,bj)) * EmPmR_tmp(i,j,bi,bj) |
|
935 |
ENDDO |
ENDDO |
936 |
ENDDO |
ENDDO |
937 |
ENDDO |
ENDDO |