3 |
|
|
4 |
#include "PACKAGES_CONFIG.h" |
#include "PACKAGES_CONFIG.h" |
5 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
|
#ifdef ALLOW_EXF |
|
|
# include "EXF_OPTIONS.h" |
|
|
#endif |
|
6 |
|
|
7 |
CBOP |
CBOP |
8 |
C !ROUTINE: EXTERNAL_FORCING_U |
C !ROUTINE: EXTERNAL_FORCING_U |
318 |
endif |
endif |
319 |
#endif |
#endif |
320 |
|
|
|
#ifdef ALLOW_CLIMTEMP_RELAXATION |
|
|
IF ( tauThetaClimRelax3Dim .NE. 0. ) THEN |
|
|
DO j=1,sNy |
|
|
DO i=1,sNx |
|
|
gT(i,j,klev,bi,bj) = gT(i,j,klev,bi,bj) |
|
|
& -1./tauThetaClimRelax3Dim |
|
|
& *(theta(i,j,klev,bi,bj)-thetaStar(i,j,klev,bi,bj)) |
|
|
& *hFacC(i,j,klev,bi,bj)*recip_hFacC(i,j,kLev,bi,bj) |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDIF |
|
|
#endif |
|
|
|
|
321 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
322 |
IF (useOBCS) THEN |
IF (useOBCS) THEN |
323 |
CALL OBCS_SPONGE_T( |
CALL OBCS_SPONGE_T( |
416 |
endif |
endif |
417 |
#endif |
#endif |
418 |
|
|
|
#ifdef ALLOW_CLIMSALT_RELAXATION |
|
|
IF ( tauSaltClimRelax3Dim .NE. 0. ) THEN |
|
|
DO j=1,sNy |
|
|
DO i=1,sNx |
|
|
gS(i,j,klev,bi,bj) = gS(i,j,klev,bi,bj) |
|
|
& -1./tauSaltClimRelax3Dim |
|
|
& *(salt(i,j,klev,bi,bj)-saltStar(i,j,klev,bi,bj)) |
|
|
& *hFacC(i,j,klev,bi,bj)*recip_hFacC(i,j,kLev,bi,bj) |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDIF |
|
|
#endif |
|
|
|
|
419 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
420 |
IF (useOBCS) THEN |
IF (useOBCS) THEN |
421 |
CALL OBCS_SPONGE_S( |
CALL OBCS_SPONGE_S( |