42 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
43 |
#include "PARAMS.h" |
#include "PARAMS.h" |
44 |
#include "GRID.h" |
#include "GRID.h" |
45 |
|
#include "FFIELDS.h" |
46 |
|
|
47 |
C == Routine arguments == |
C == Routine arguments == |
48 |
C fZon - Work array for flux of temperature in the east-west |
C fZon - Work array for flux of temperature in the east-west |
253 |
ENDDO |
ENDDO |
254 |
|
|
255 |
C-- External thermal forcing term(s) |
C-- External thermal forcing term(s) |
256 |
|
C o Surface relaxation term |
257 |
|
IF ( TOP_LAYER ) THEN |
258 |
|
DO j=jMin,jMax |
259 |
|
DO i=iMin,iMax |
260 |
|
gT(i,j,k,bi,bj)=gT(i,j,k,bi,bj) |
261 |
|
& -lambdaThetaClimRelax*(theta(i,j,k,bi,bj)-thetaClim(i,j,k,bi,bj)) |
262 |
|
ENDDO |
263 |
|
ENDDO |
264 |
|
ENDIF |
265 |
|
|
266 |
RETURN |
RETURN |
267 |
END |
END |