419 |
IF (useAbsVorticity) THEN |
IF (useAbsVorticity) THEN |
420 |
CALL MOM_VI_U_CORIOLIS(bi,bj,K,vFld,omega3,hFacZ,r_hFacZ, |
CALL MOM_VI_U_CORIOLIS(bi,bj,K,vFld,omega3,hFacZ,r_hFacZ, |
421 |
& uCf,myThid) |
& uCf,myThid) |
422 |
|
ELSEIF (highOrderVorticity) THEN |
423 |
|
CALL MOM_VI_U_CORIOLIS_C4(bi,bj,k,vFld,vort3,r_hFacZ, |
424 |
|
& uCf,myThid) |
425 |
ELSE |
ELSE |
426 |
CALL MOM_VI_U_CORIOLIS(bi,bj,k,vFld,vort3,hFacZ,r_hFacZ, |
CALL MOM_VI_U_CORIOLIS(bi,bj,k,vFld,vort3,hFacZ,r_hFacZ, |
427 |
& uCf,myThid) |
& uCf,myThid) |
428 |
ENDIF |
ENDIF |
|
c CALL MOM_VI_U_CORIOLIS_C4(bi,bj,K,vFld,vort3,r_hFacZ,uCf,myThid) |
|
429 |
DO j=jMin,jMax |
DO j=jMin,jMax |
430 |
DO i=iMin,iMax |
DO i=iMin,iMax |
431 |
gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+uCf(i,j) |
gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)+uCf(i,j) |
434 |
IF (useAbsVorticity) THEN |
IF (useAbsVorticity) THEN |
435 |
CALL MOM_VI_V_CORIOLIS(bi,bj,K,uFld,omega3,hFacZ,r_hFacZ, |
CALL MOM_VI_V_CORIOLIS(bi,bj,K,uFld,omega3,hFacZ,r_hFacZ, |
436 |
& vCf,myThid) |
& vCf,myThid) |
437 |
|
ELSEIF (highOrderVorticity) THEN |
438 |
|
CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K,uFld,vort3,r_hFacZ, |
439 |
|
& vCf,myThid) |
440 |
ELSE |
ELSE |
441 |
CALL MOM_VI_V_CORIOLIS(bi,bj,k,uFld,vort3,hFacZ,r_hFacZ, |
CALL MOM_VI_V_CORIOLIS(bi,bj,k,uFld,vort3,hFacZ,r_hFacZ, |
442 |
& vCf,myThid) |
& vCf,myThid) |
443 |
ENDIF |
ENDIF |
|
c CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K,uFld,vort3,r_hFacZ,vCf,myThid) |
|
444 |
DO j=jMin,jMax |
DO j=jMin,jMax |
445 |
DO i=iMin,iMax |
DO i=iMin,iMax |
446 |
gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+vCf(i,j) |
gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)+vCf(i,j) |