161 |
INTEGER k, km1, kp1, kup, kDown |
INTEGER k, km1, kp1, kup, kDown |
162 |
|
|
163 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
164 |
_RL tmpFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL tmpFac |
|
LOGICAL DIAGNOSTICS_IS_ON |
|
|
EXTERNAL DIAGNOSTICS_IS_ON |
|
165 |
#endif /* ALLOW_DIAGNOSTICS */ |
#endif /* ALLOW_DIAGNOSTICS */ |
166 |
|
|
167 |
|
|
370 |
#ifdef ALLOW_MOM_FLUXFORM |
#ifdef ALLOW_MOM_FLUXFORM |
371 |
IF (.NOT. vectorInvariantMomentum) CALL MOM_FLUXFORM( |
IF (.NOT. vectorInvariantMomentum) CALL MOM_FLUXFORM( |
372 |
I bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown, |
I bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown, |
373 |
I dPhiHydX,dPhiHydY,KappaRU,KappaRV, |
I KappaRU, KappaRV, |
374 |
U fVerU, fVerV, |
U fVerU, fVerV, |
375 |
|
O guDissip, gvDissip, |
376 |
I myTime, myIter, myThid) |
I myTime, myIter, myThid) |
377 |
#endif |
#endif |
378 |
#ifdef ALLOW_MOM_VECINV |
#ifdef ALLOW_MOM_VECINV |
379 |
IF (vectorInvariantMomentum) CALL MOM_VECINV( |
IF (vectorInvariantMomentum) CALL MOM_VECINV( |
380 |
I bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown, |
I bi,bj,iMin,iMax,jMin,jMax,k,kup,kDown, |
381 |
I dPhiHydX,dPhiHydY,KappaRU,KappaRV, |
I KappaRU, KappaRV, |
382 |
U fVerU, fVerV, |
U fVerU, fVerV, |
383 |
O guDissip, gvDissip, |
O guDissip, gvDissip, |
384 |
I myTime, myIter, myThid) |
I myTime, myIter, myThid) |
487 |
IF ( usediagnostics ) THEN |
IF ( usediagnostics ) THEN |
488 |
|
|
489 |
CALL DIAGNOSTICS_FILL(totPhihyd,'PHIHYD ',0,Nr,0,1,1,myThid) |
CALL DIAGNOSTICS_FILL(totPhihyd,'PHIHYD ',0,Nr,0,1,1,myThid) |
490 |
|
CALL DIAGNOSTICS_FILL(phiHydLow,'PHIBOT ',0, 1,0,1,1,myThid) |
491 |
|
|
492 |
IF ( DIAGNOSTICS_IS_ON('PHIHYDSQ',myThid) ) THEN |
tmpFac = 1. _d 0 |
493 |
DO bj = myByLo(myThid), myByHi(myThid) |
CALL DIAGNOSTICS_SCALE_FILL(totPhihyd,tmpFac,2, |
494 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
& 'PHIHYDSQ',0,Nr,0,1,1,myThid) |
|
DO k = 1,Nr |
|
|
DO j = 1,sNy |
|
|
DO i = 1,sNx |
|
|
tmpFld(i,j) = totPhihyd(i,j,k,bi,bj)*totPhihyd(i,j,k,bi,bj) |
|
|
ENDDO |
|
|
ENDDO |
|
|
CALL DIAGNOSTICS_FILL(tmpFld,'PHIHYDSQ',k,1,2,bi,bj,myThid) |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDIF |
|
495 |
|
|
496 |
CALL DIAGNOSTICS_FILL(phiHydLow,'PHIBOT ',0,1,0,1,1,myThid) |
CALL DIAGNOSTICS_SCALE_FILL(phiHydLow,tmpFac,2, |
497 |
|
& 'PHIBOTSQ',0, 1,0,1,1,myThid) |
|
IF ( DIAGNOSTICS_IS_ON('PHIBOTSQ',myThid) ) THEN |
|
|
DO bj = myByLo(myThid), myByHi(myThid) |
|
|
DO bi = myBxLo(myThid), myBxHi(myThid) |
|
|
DO j = 1,sNy |
|
|
DO i = 1,sNx |
|
|
tmpFld(i,j) = phiHydLow(i,j,bi,bj)*phiHydLow(i,j,bi,bj) |
|
|
ENDDO |
|
|
ENDDO |
|
|
CALL DIAGNOSTICS_FILL(tmpFld,'PHIBOTSQ',0,1,2,bi,bj,myThid) |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDIF |
|
498 |
|
|
499 |
ENDIF |
ENDIF |
500 |
#endif /* ALLOW_DIAGNOSTICS */ |
#endif /* ALLOW_DIAGNOSTICS */ |