/[MITgcm]/MITgcm/verification/aim.5l_cs/code/mom_vecinv.F
ViewVC logotype

Diff of /MITgcm/verification/aim.5l_cs/code/mom_vecinv.F

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

revision 1.3 by jmc, Fri Apr 11 13:42:41 2003 UTC revision 1.4 by jmc, Thu Apr 17 13:49:41 2003 UTC
# Line 312  C-    No-slip BCs impose a drag at botto Line 312  C-    No-slip BCs impose a drag at botto
312         ENDDO         ENDDO
313        ENDIF        ENDIF
314    
315  C--   Forcing term  C--   Forcing term (moved to timestep.F)
316        IF (momForcing)  c     IF (momForcing)
317       &  CALL EXTERNAL_FORCING_U(  c    &  CALL EXTERNAL_FORCING_U(
318       I     iMin,iMax,jMin,jMax,bi,bj,k,  c    I     iMin,iMax,jMin,jMax,bi,bj,k,
319       I     myCurrentTime,myThid)  c    I     myCurrentTime,myThid)
320    
321  C--   Metric terms for curvilinear grid systems  C--   Metric terms for curvilinear grid systems
322  c     IF (usingSphericalPolarMTerms) THEN  c     IF (usingSphericalPolarMTerms) THEN
# Line 378  C-    No-slip BCs impose a drag at botto Line 378  C-    No-slip BCs impose a drag at botto
378         ENDDO         ENDDO
379        ENDIF        ENDIF
380    
381  C--   Forcing term  C--   Forcing term (moved to timestep.F)
382        IF (momForcing)  c     IF (momForcing)
383       & CALL EXTERNAL_FORCING_V(  c    & CALL EXTERNAL_FORCING_V(
384       I     iMin,iMax,jMin,jMax,bi,bj,k,  c    I     iMin,iMax,jMin,jMax,bi,bj,k,
385       I     myCurrentTime,myThid)  c    I     myCurrentTime,myThid)
386    
387  C--   Metric terms for curvilinear grid systems  C--   Metric terms for curvilinear grid systems
388  c     IF (usingSphericalPolarMTerms) THEN  c     IF (usingSphericalPolarMTerms) THEN
# Line 396  c      ENDDO Line 396  c      ENDDO
396  c     ENDIF  c     ENDIF
397    
398  C--   Horizontal Coriolis terms  C--   Horizontal Coriolis terms
399        CALL MOM_VI_MASK_VORT3(bi,bj,k,hFacZ,r_hFacZ,vort3,myThid)        IF (useCoriolis .AND. .NOT.useCDscheme) THEN
 c     CALL MOM_VI_CALC_ABSVORT3(bi,bj,k,vort3,omega3,myThid)  
       IF (useCoriolis) THEN  
400         CALL MOM_VI_CORIOLIS(bi,bj,k,uFld,vFld,omega3,hFacZ,r_hFacZ,         CALL MOM_VI_CORIOLIS(bi,bj,k,uFld,vFld,omega3,hFacZ,r_hFacZ,
401       &                      uCf,vCf,myThid)       &                      uCf,vCf,myThid)
402         DO j=jMin,jMax         DO j=jMin,jMax
# Line 410  c     CALL MOM_VI_CALC_ABSVORT3(bi,bj,k, Line 408  c     CALL MOM_VI_CALC_ABSVORT3(bi,bj,k,
408        ENDIF        ENDIF
409    
410        IF (momAdvection) THEN        IF (momAdvection) THEN
411           CALL MOM_VI_MASK_VORT3(bi,bj,k,hFacZ,r_hFacZ,vort3,myThid)
412  C--   Horizontal advection of relative vorticity  C--   Horizontal advection of relative vorticity
413  c      CALL MOM_VI_U_CORIOLIS(bi,bj,K,vFld,omega3,r_hFacZ,uCf,myThid)  c      CALL MOM_VI_U_CORIOLIS(bi,bj,K,vFld,omega3,r_hFacZ,uCf,myThid)
414         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,
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22