13 |
I xA,yA,uTrans,vTrans,rTrans,maskUp, |
I xA,yA,uTrans,vTrans,rTrans,maskUp, |
14 |
I KappaRS, |
I KappaRS, |
15 |
U fVerS, |
U fVerS, |
16 |
I myCurrentTime, myThid ) |
I myCurrentTime,myIter,myThid ) |
17 |
C /==========================================================\ |
C /==========================================================\ |
18 |
C | SUBROUTINE CALC_GS | |
C | SUBROUTINE CALC_GS | |
19 |
C | o Calculate the salt tendency terms. | |
C | o Calculate the salt tendency terms. | |
48 |
#include "PARAMS.h" |
#include "PARAMS.h" |
49 |
#include "GRID.h" |
#include "GRID.h" |
50 |
#include "FFIELDS.h" |
#include "FFIELDS.h" |
51 |
|
#include "GAD.h" |
52 |
|
|
53 |
C == Routine arguments == |
C == Routine arguments == |
54 |
C fVerS - Flux of salt (S) in the vertical |
C fVerS - Flux of salt (S) in the vertical |
73 |
INTEGER k,kUp,kDown,kM1 |
INTEGER k,kUp,kDown,kM1 |
74 |
INTEGER bi,bj,iMin,iMax,jMin,jMax |
INTEGER bi,bj,iMin,iMax,jMin,jMax |
75 |
_RL myCurrentTime |
_RL myCurrentTime |
76 |
|
INTEGER myIter |
77 |
INTEGER myThid |
INTEGER myThid |
78 |
CEndOfInterface |
CEndOfInterface |
79 |
|
|
471 |
I iMin,iMax,jMin,jMax,bi,bj,k, |
I iMin,iMax,jMin,jMax,bi,bj,k, |
472 |
I myCurrentTime,myThid) |
I myCurrentTime,myThid) |
473 |
|
|
474 |
|
IF ( saltAdvScheme.EQ.ENUM_CENTERED_2ND |
475 |
|
& .OR.saltAdvScheme.EQ.ENUM_UPWIND_3RD |
476 |
|
& .OR.saltAdvScheme.EQ.ENUM_CENTERED_4TH ) THEN |
477 |
|
CALL ADAMS_BASHFORTH2( |
478 |
|
I bi, bj, K, |
479 |
|
U gS, gSnm1, |
480 |
|
I myIter, myThid ) |
481 |
|
ENDIF |
482 |
|
|
483 |
|
#ifdef NONLIN_FRSURF |
484 |
|
IF (nonlinFreeSurf.GT.0) THEN |
485 |
|
CALL FREESURF_RESCALE_G( |
486 |
|
I bi, bj, K, |
487 |
|
U gS, |
488 |
|
I myThid ) |
489 |
|
ENDIF |
490 |
|
#endif /* NONLIN_FRSURF */ |
491 |
|
|
492 |
RETURN |
RETURN |
493 |
END |
END |