537 |
|
|
538 |
IF (momAdvection) THEN |
IF (momAdvection) THEN |
539 |
C-- Horizontal advection of relative (or absolute) vorticity |
C-- Horizontal advection of relative (or absolute) vorticity |
540 |
IF (highOrderVorticity.AND.useAbsVorticity) THEN |
IF ( (highOrderVorticity.OR.upwindVorticity) |
541 |
|
& .AND.useAbsVorticity ) THEN |
542 |
CALL MOM_VI_U_CORIOLIS_C4(bi,bj,k,vFld,omega3,r_hFacZ, |
CALL MOM_VI_U_CORIOLIS_C4(bi,bj,k,vFld,omega3,r_hFacZ, |
543 |
& uCf,myThid) |
& uCf,myThid) |
544 |
ELSEIF (highOrderVorticity) THEN |
ELSEIF ( (highOrderVorticity.OR.upwindVorticity) ) THEN |
545 |
CALL MOM_VI_U_CORIOLIS_C4(bi,bj,k,vFld,vort3, r_hFacZ, |
CALL MOM_VI_U_CORIOLIS_C4(bi,bj,k,vFld,vort3, r_hFacZ, |
546 |
& uCf,myThid) |
& uCf,myThid) |
547 |
ELSEIF (useAbsVorticity) THEN |
ELSEIF ( useAbsVorticity ) THEN |
548 |
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, |
549 |
& uCf,myThid) |
& uCf,myThid) |
550 |
ELSE |
ELSE |
556 |
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) |
557 |
ENDDO |
ENDDO |
558 |
ENDDO |
ENDDO |
559 |
IF (highOrderVorticity.AND.useAbsVorticity) THEN |
IF ( (highOrderVorticity.OR.upwindVorticity) |
560 |
|
& .AND.useAbsVorticity ) THEN |
561 |
CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K,uFld,omega3,r_hFacZ, |
CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K,uFld,omega3,r_hFacZ, |
562 |
& vCf,myThid) |
& vCf,myThid) |
563 |
ELSEIF (highOrderVorticity) THEN |
ELSEIF ( (highOrderVorticity.OR.upwindVorticity) ) THEN |
564 |
CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K,uFld,vort3, r_hFacZ, |
CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K,uFld,vort3, r_hFacZ, |
565 |
& vCf,myThid) |
& vCf,myThid) |
566 |
ELSEIF (useAbsVorticity) THEN |
ELSEIF ( useAbsVorticity ) THEN |
567 |
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, |
568 |
& vCf,myThid) |
& vCf,myThid) |
569 |
ELSE |
ELSE |