8 |
I xA,yA,uTrans,vTrans,rTrans,maskup,maskC, |
I xA,yA,uTrans,vTrans,rTrans,maskup,maskC, |
9 |
I K13,K23,KappaRS,KapGM, |
I K13,K23,KappaRS,KapGM, |
10 |
U af,df,fZon,fMer,fVerS, |
U af,df,fZon,fMer,fVerS, |
11 |
I myThid ) |
I myCurrentTime, myThid ) |
12 |
C /==========================================================\ |
C /==========================================================\ |
13 |
C | SUBROUTINE CALC_GS | |
C | SUBROUTINE CALC_GS | |
14 |
C | o Calculate the salt tendency terms. | |
C | o Calculate the salt tendency terms. | |
82 |
INTEGER k,kUp,kDown,kM1 |
INTEGER k,kUp,kDown,kM1 |
83 |
INTEGER bi,bj,iMin,iMax,jMin,jMax |
INTEGER bi,bj,iMin,iMax,jMin,jMax |
84 |
INTEGER myThid |
INTEGER myThid |
85 |
|
_RL myCurrentTime |
86 |
CEndOfInterface |
CEndOfInterface |
87 |
|
|
88 |
C == Local variables == |
C == Local variables == |
262 |
ENDDO |
ENDDO |
263 |
ENDDO |
ENDDO |
264 |
|
|
265 |
C-- External P-E forcing term(s) |
C-- External forcing term(s) |
266 |
|
CALL EXTERNAL_FORCING_S( |
267 |
|
I iMin,iMax,jMin,jMax,bi,bj,k, |
268 |
|
I myCurrentTime,myThid) |
269 |
C o Surface relaxation term |
C o Surface relaxation term |
270 |
IF ( TOP_LAYER ) THEN |
C IF ( TOP_LAYER ) THEN |
271 |
DO j=jMin,jMax |
C DO j=jMin,jMax |
272 |
DO i=iMin,iMax |
C DO i=iMin,iMax |
273 |
gS(i,j,k,bi,bj)=gS(i,j,k,bi,bj) |
C gS(i,j,k,bi,bj)=gS(i,j,k,bi,bj) |
274 |
& +maskC(i,j)*( |
C & +maskC(i,j)*( |
275 |
& -lambdaSaltClimRelax*(salt(i,j,k,bi,bj)-SSS(i,j,bi,bj)) |
C & -lambdaSaltClimRelax*(salt(i,j,k,bi,bj)-SSS(i,j,bi,bj)) |
276 |
& +EmPmR(i,j,bi,bj) ) |
C & +EmPmR(i,j,bi,bj) ) |
277 |
ENDDO |
C ENDDO |
278 |
ENDDO |
C ENDDO |
279 |
ENDIF |
C ENDIF |
280 |
|
|
281 |
|
#ifdef ALLOW_LATITUDE_CIRCLE_FFT_FILTER |
282 |
|
C-- |
283 |
|
CALL FILTER_LATCIRCS_FFT_APPLY( gS, 1, sNy, k, k, bi, bj, 1, myThid) |
284 |
|
#endif |
285 |
|
|
286 |
RETURN |
RETURN |
287 |
END |
END |