| 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( |