3 |
|
|
4 |
#include "PACKAGES_CONFIG.h" |
#include "PACKAGES_CONFIG.h" |
5 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
6 |
|
#ifdef ALLOW_EXF |
7 |
|
# include EXF_OPTIONS.h" |
8 |
|
#endif |
9 |
|
|
10 |
CBOP |
CBOP |
11 |
C !ROUTINE: EXTERNAL_FORCING_U |
C !ROUTINE: EXTERNAL_FORCING_U |
314 |
c ENDIF |
c ENDIF |
315 |
#endif |
#endif |
316 |
|
|
317 |
|
#ifdef ALLOW_CLIMTEMP_RELAXATION |
318 |
|
IF ( tauThetaClimRelax3Dim .NE. 0. ) THEN |
319 |
|
DO j=1,sNy |
320 |
|
DO i=1,sNx |
321 |
|
gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj) |
322 |
|
& -1./tauThetaClimRelax3Dim |
323 |
|
& *(theta(i,j,klev,bi,bj)-thetaStar(i,j,klev,bi,bj)) |
324 |
|
& *hFacC(i,j,klev,bi,bj)*recip_hFacC(i,j,kLev,bi,bj) |
325 |
|
ENDDO |
326 |
|
ENDDO |
327 |
|
ENDIF |
328 |
|
#endif |
329 |
|
|
330 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
331 |
IF (useOBCS) THEN |
IF (useOBCS) THEN |
332 |
CALL OBCS_SPONGE_T( |
CALL OBCS_SPONGE_T( |
418 |
ENDDO |
ENDDO |
419 |
ENDIF |
ENDIF |
420 |
|
|
421 |
|
#ifdef ALLOW_CLIMSALT_RELAXATION |
422 |
|
IF ( tauSaltClimRelax3Dim .NE. 0. ) THEN |
423 |
|
DO j=1,sNy |
424 |
|
DO i=1,sNx |
425 |
|
gS(i,j,klev,bi,bj) = gS(i,j,klev,bi,bj) |
426 |
|
& -1./tauSaltClimRelax3Dim |
427 |
|
& *(salt(i,j,klev,bi,bj)-saltStar(i,j,klev,bi,bj)) |
428 |
|
& *hFacC(i,j,klev,bi,bj)*recip_hFacC(i,j,kLev,bi,bj) |
429 |
|
ENDDO |
430 |
|
ENDDO |
431 |
|
ENDIF |
432 |
|
#endif |
433 |
|
|
434 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
435 |
IF (useOBCS) THEN |
IF (useOBCS) THEN |
436 |
CALL OBCS_SPONGE_S( |
CALL OBCS_SPONGE_S( |