23 |
#include "PARAMS.h" |
#include "PARAMS.h" |
24 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
25 |
#include "GRID.h" |
#include "GRID.h" |
26 |
|
#include "FFIELDS.h" |
27 |
#ifdef ALLOW_EXF |
#ifdef ALLOW_EXF |
28 |
# include "EXF_OPTIONS.h" |
# include "EXF_OPTIONS.h" |
29 |
# include "EXF_FIELDS.h" |
# include "EXF_FIELDS.h" |
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 |