/[MITgcm]/MITgcm/pkg/mom_vecinv/mom_vecinv.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_vecinv/mom_vecinv.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.40 by jmc, Thu Jun 9 15:57:45 2005 UTC revision 1.41 by jmc, Mon Jun 20 23:07:17 2005 UTC
# Line 415  C- jmc: change it to keep the Coriolis t Line 415  C- jmc: change it to keep the Coriolis t
415        ENDIF        ENDIF
416    
417        IF (momAdvection) THEN        IF (momAdvection) THEN
418  C--   Horizontal advection of relative vorticity  C--   Horizontal advection of relative (or absolute) vorticity
419         IF (useAbsVorticity) THEN         IF (highOrderVorticity.AND.useAbsVorticity) THEN
420          CALL MOM_VI_U_CORIOLIS(bi,bj,K,vFld,omega3,hFacZ,r_hFacZ,          CALL MOM_VI_U_CORIOLIS_C4(bi,bj,k,vFld,omega3,r_hFacZ,
421       &                         uCf,myThid)       &                         uCf,myThid)
422         ELSEIF (highOrderVorticity) THEN         ELSEIF (highOrderVorticity) THEN
423          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,
424         &                         uCf,myThid)
425           ELSEIF (useAbsVorticity) THEN
426            CALL MOM_VI_U_CORIOLIS(bi,bj,K,vFld,omega3,hFacZ,r_hFacZ,
427       &                         uCf,myThid)       &                         uCf,myThid)
428         ELSE         ELSE
429          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,
430       &                         uCf,myThid)       &                         uCf,myThid)
431         ENDIF         ENDIF
432         DO j=jMin,jMax         DO j=jMin,jMax
# Line 431  C--   Horizontal advection of relative v Line 434  C--   Horizontal advection of relative v
434           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)
435          ENDDO          ENDDO
436         ENDDO         ENDDO
437         IF (useAbsVorticity) THEN         IF (highOrderVorticity.AND.useAbsVorticity) THEN
438          CALL MOM_VI_V_CORIOLIS(bi,bj,K,uFld,omega3,hFacZ,r_hFacZ,          CALL MOM_VI_V_CORIOLIS_C4(bi,bj,K,uFld,omega3,r_hFacZ,
439       &                         vCf,myThid)       &                         vCf,myThid)
440         ELSEIF (highOrderVorticity) THEN         ELSEIF (highOrderVorticity) THEN
441          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,
442         &                         vCf,myThid)
443           ELSEIF (useAbsVorticity) THEN
444            CALL MOM_VI_V_CORIOLIS(bi,bj,K,uFld,omega3,hFacZ,r_hFacZ,
445       &                         vCf,myThid)       &                         vCf,myThid)
446         ELSE         ELSE
447          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,
448       &                         vCf,myThid)       &                         vCf,myThid)
449         ENDIF         ENDIF
450         DO j=jMin,jMax         DO j=jMin,jMax
# Line 571  C--   Set du/dt & dv/dt on boundaries to Line 577  C--   Set du/dt & dv/dt on boundaries to
577          ENDIF          ENDIF
578  #endif /*  ALLOW_MNC  */  #endif /*  ALLOW_MNC  */
579        ENDIF        ENDIF
580          
581  #endif /* ALLOW_MOM_VECINV */  #endif /* ALLOW_MOM_VECINV */
582    
583        RETURN        RETURN
Legend:
Removed from v.1.40  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.22