/[MITgcm]/MITgcm/pkg/rbcs/rbcs_add_tendency.F
ViewVC logotype

Diff of /MITgcm/pkg/rbcs/rbcs_add_tendency.F

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

revision 1.8 by gforget, Sat May 21 00:44:53 2011 UTC revision 1.9 by gforget, Tue Jun 26 22:19:17 2012 UTC
# Line 45  c     INTEGER myIter Line 45  c     INTEGER myIter
45  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
46  C  i,j            :: loop indices  C  i,j            :: loop indices
47        INTEGER i,j        INTEGER i,j
48          _RL rbcsVanishingFac
49  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
50        INTEGER iTracer        INTEGER iTracer
51        INTEGER irbc        INTEGER irbc
# Line 53  CEOP Line 54  CEOP
54    
55  #ifdef ALLOW_RBCS  #ifdef ALLOW_RBCS
56    
57          if (rbcsVanishingTime.GT.0. _d 0) then
58            rbcsVanishingFac =
59         &  MAX( 0. _d 0 , 1. _d 0 - mytime / rbcsVanishingTime )
60          else
61            rbcsVanishingFac = 1. _d 0
62          endif
63    
64  #ifndef DISABLE_RBCS_MOM  #ifndef DISABLE_RBCS_MOM
65        IF ( tracerNum.EQ.-1 .AND. useRBCuVel ) THEN        IF ( tracerNum.EQ.-1 .AND. useRBCuVel ) THEN
66          DO j=0,sNy+1          DO j=0,sNy+1
# Line 61  CEOP Line 69  CEOP
69       &       - RBC_maskU(i,j,k,bi,bj)/tauRelaxU       &       - RBC_maskU(i,j,k,bi,bj)/tauRelaxU
70       &        *( uVel(i,j,k,bi,bj)- RBCuVel(i,j,k,bi,bj) )       &        *( uVel(i,j,k,bi,bj)- RBCuVel(i,j,k,bi,bj) )
71  c    &        *maskW(i,j,k,bi,bj)  c    &        *maskW(i,j,k,bi,bj)
72         &        *rbcsVanishingFac        
73           ENDDO           ENDDO
74          ENDDO          ENDDO
75        ENDIF        ENDIF
# Line 71  c    &        *maskW(i,j,k,bi,bj) Line 80  c    &        *maskW(i,j,k,bi,bj)
80       &       - RBC_maskV(i,j,k,bi,bj)/tauRelaxV       &       - RBC_maskV(i,j,k,bi,bj)/tauRelaxV
81       &        *( vVel(i,j,k,bi,bj)- RBCvVel(i,j,k,bi,bj) )       &        *( vVel(i,j,k,bi,bj)- RBCvVel(i,j,k,bi,bj) )
82  c    &        *maskS(i,j,k,bi,bj)  c    &        *maskS(i,j,k,bi,bj)
83         &        *rbcsVanishingFac        
84           ENDDO           ENDDO
85          ENDDO          ENDDO
86        ENDIF        ENDIF
# Line 83  c    &        *maskS(i,j,k,bi,bj) Line 93  c    &        *maskS(i,j,k,bi,bj)
93       &       - RBC_mask(i,j,k,bi,bj,1)/tauRelaxT       &       - RBC_mask(i,j,k,bi,bj,1)/tauRelaxT
94       &        *( theta(i,j,k,bi,bj)- RBCtemp(i,j,k,bi,bj) )       &        *( theta(i,j,k,bi,bj)- RBCtemp(i,j,k,bi,bj) )
95  c    &        *maskC(i,j,k,bi,bj)  c    &        *maskC(i,j,k,bi,bj)
96  #ifdef ALLOW_RBCS_SPIN       &        *rbcsVanishingFac
97       &        *MAX( 0. _d 0 , 1. _d 0 - mytime/( 94608000. _d 0 ) )  
 #endif  
98           ENDDO           ENDDO
99          ENDDO          ENDDO
100        ENDIF        ENDIF
# Line 97  c    &        *maskC(i,j,k,bi,bj) Line 106  c    &        *maskC(i,j,k,bi,bj)
106       &       - RBC_mask(i,j,k,bi,bj,2)/tauRelaxS       &       - RBC_mask(i,j,k,bi,bj,2)/tauRelaxS
107       &        *( salt(i,j,k,bi,bj)- RBCsalt(i,j,k,bi,bj) )       &        *( salt(i,j,k,bi,bj)- RBCsalt(i,j,k,bi,bj) )
108  c    &        *maskC(i,j,k,bi,bj)  c    &        *maskC(i,j,k,bi,bj)
109  #ifdef ALLOW_RBCS_SPIN       &        *rbcsVanishingFac
      &        *MAX( 0. _d 0 , 1. _d 0 - mytime/( 94608000. _d 0 ) )  
 #endif  
110           ENDDO           ENDDO
111          ENDDO          ENDDO
112        ENDIF        ENDIF
# Line 116  c    &        *maskC(i,j,k,bi,bj) Line 123  c    &        *maskC(i,j,k,bi,bj)
123       &        *( pTracer(i,j,k,bi,bj,iTracer)       &        *( pTracer(i,j,k,bi,bj,iTracer)
124       &           - RBC_ptracers(i,j,k,bi,bj,iTracer) )       &           - RBC_ptracers(i,j,k,bi,bj,iTracer) )
125  c    &        *maskC(i,j,k,bi,bj)  c    &        *maskC(i,j,k,bi,bj)
126         &        *rbcsVanishingFac
127           ENDDO           ENDDO
128          ENDDO          ENDDO
129         ENDIF         ENDIF
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22