/[MITgcm]/MITgcm/model/src/thermodynamics.F
ViewVC logotype

Diff of /MITgcm/model/src/thermodynamics.F

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

revision 1.16 by jmc, Fri Feb 8 22:15:33 2002 UTC revision 1.17 by adcroft, Mon Mar 4 17:26:41 2002 UTC
# Line 472  C disable this section of code. Line 472  C disable this section of code.
472       I                      myTime,myIter,myThid)       I                      myTime,myIter,myThid)
473           ENDIF           ENDIF
474          ENDIF          ENDIF
475    C Since passive tracers are configurable separately from T,S we
476    C call the multi-dimensional method for PTRACERS regardless
477    C of whether multiDimAdvection is set or not.
478    #ifdef ALLOW_PTRACERS
479            IF ( usePTRACERS ) THEN
480             CALL PTRACERS_ADVECTION( bi,bj,myIter,myTime,myThid )
481            ENDIF
482    #endif /* ALLOW_PTRACERS */
483  #endif /* DISABLE_MULTIDIM_ADVECTION */  #endif /* DISABLE_MULTIDIM_ADVECTION */
484    
485  C--     Start of thermodynamics loop  C--     Start of thermodynamics loop
# Line 561  C        and step forward storing result Line 569  C        and step forward storing result
569       I         myIter,myThid)       I         myIter,myThid)
570           ENDIF           ENDIF
571  #endif  #endif
572    #ifdef ALLOW_PTRACERS
573             IF ( usePTRACERS ) THEN
574               CALL PTRACERS_INTEGERATE(
575         I         bi,bj,k,
576         I         xA,yA,uTrans,vTrans,rTrans,maskUp,
577         X         KappaRS,
578         I         myIter,myTime,myThid)
579             ENDIF
580    #endif /* ALLOW_PTRACERS */
581    
582  #ifdef   ALLOW_OBCS  #ifdef   ALLOW_OBCS
583  C--      Apply open boundary conditions  C--      Apply open boundary conditions
# Line 632  CADJ STORE gTr1(:,:,:,bi,bj) = comlev1_b Line 649  CADJ STORE gTr1(:,:,:,bi,bj) = comlev1_b
649           ENDIF           ENDIF
650  #endif  #endif
651    
652    #ifdef ALLOW_PTRACERS
653    C Vertical diffusion (implicit) for passive tracers
654             IF ( usePTRACERS ) THEN
655               CALL PTRACERS_IMPLDIFF( bi,bj,KappaRS,myThid )
656             ENDIF
657    #endif /* ALLOW_PTRACERS */
658    
659  #ifdef   ALLOW_OBCS  #ifdef   ALLOW_OBCS
660  C--      Apply open boundary conditions  C--      Apply open boundary conditions
661           IF (useOBCS) THEN           IF (useOBCS) THEN
# Line 661  Ccs- Line 685  Ccs-
685        ENDIF        ENDIF
686  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
687    
688    #ifndef DISABLE_DEBUGMODE
689          If (debugMode) THEN
690           CALL DEBUG_STATS_RL(Nr,uVel,'Uvel (THERMODYNAMICS)',myThid)
691           CALL DEBUG_STATS_RL(Nr,vVel,'Vvel (THERMODYNAMICS)',myThid)
692           CALL DEBUG_STATS_RL(Nr,wVel,'Wvel (THERMODYNAMICS)',myThid)
693           CALL DEBUG_STATS_RL(Nr,theta,'Theta (THERMODYNAMICS)',myThid)
694           CALL DEBUG_STATS_RL(Nr,salt,'Salt (THERMODYNAMICS)',myThid)
695           CALL DEBUG_STATS_RL(Nr,Gt,'Gt (THERMODYNAMICS)',myThid)
696           CALL DEBUG_STATS_RL(Nr,Gs,'Gs (THERMODYNAMICS)',myThid)
697           CALL DEBUG_STATS_RL(Nr,GtNm1,'GtNm1 (THERMODYNAMICS)',myThid)
698           CALL DEBUG_STATS_RL(Nr,GsNm1,'GsNm1 (THERMODYNAMICS)',myThid)
699           CALL PTRACERS_DEBUG(myThid)
700          ENDIF
701    #endif
702    
703        RETURN        RETURN
704        END        END
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22