/[MITgcm]/MITgcm/pkg/kpp/kpp_forcing_surf.F
ViewVC logotype

Diff of /MITgcm/pkg/kpp/kpp_forcing_surf.F

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

revision 1.2 by dimitri, Sat Sep 22 17:55:32 2007 UTC revision 1.3 by atn, Fri Dec 21 02:54:34 2007 UTC
# Line 200  cph) Line 200  cph)
200         jp1 = j + 1         jp1 = j + 1
201         DO i = imin, imax         DO i = imin, imax
202          ip1 = i+1          ip1 = i+1
   
203          bo(I,J) = - gravity *          bo(I,J) = - gravity *
204       &       ( TTALPHA(I,J,1) * (surfForcT(i,j,bi,bj)+       &       ( TTALPHA(I,J,1) * (surfForcT(i,j,bi,bj)+
205       &       surfForcTice(i,j,bi,bj)) +       &       surfForcTice(i,j,bi,bj)) +
206       &       SSBETA(I,J,1) * surfForcS(i,j,bi,bj) )       &       SSBETA(I,J,1) * surfForcS(i,j,bi,bj) )
207       &       / rhoSurf(I,J)       &       / rhoSurf(I,J)
   
208          bosol(I,J) = gravity * TTALPHA(I,J,1) * Qsw(i,j,bi,bj) *          bosol(I,J) = gravity * TTALPHA(I,J,1) * Qsw(i,j,bi,bj) *
209       &       recip_Cp*recip_rhoConst       &       recip_Cp*recip_rhoConst
210       &       / rhoSurf(I,J)       &       / rhoSurf(I,J)
 #ifdef ALLOW_SALT_PLUME  
         boplume(I,J) = - gravity * SSBETA(I,J,1)  
      &       * saltPlumeFlux(i,j,bi,bj)  
      &       * recip_rhoConst / rhoSurf(I,J)  
 #endif /* ALLOW_SALT_PLUME */  
211         END DO         END DO
212        END DO        END DO
213    
214    #ifdef ALLOW_SALT_PLUME
215          IF ( useSALT_PLUME ) THEN
216             DO j = jmin, jmax
217                jp1 = j + 1
218                DO i = imin, imax
219                   ip1 = i+1
220                   boplume(I,J) = - gravity * SSBETA(I,J,1)
221         &              * saltPlumeFlux(i,j,bi,bj)
222         &              * recip_rhoConst / rhoSurf(I,J)
223                END DO
224             END DO
225          ENDIF
226    #endif /* ALLOW_SALT_PLUME */
227    
228  cph(  cph(
229  CADJ store ustar = comlev1_kpp, key = ikppkey  CADJ store ustar = comlev1_kpp, key = ikppkey
230  cph)  cph)
231    
232    #ifdef ALLOW_DIAGNOSTICS
233          IF ( useDiagnostics ) THEN
234             CALL DIAGNOSTICS_FILL(bo     ,'KPPbo   ',0,1,1,bi,bj,myThid)
235             CALL DIAGNOSTICS_FILL(bosol  ,'KPPbosol',0,1,1,bi,bj,myThid)
236    #ifdef ALLOW_SALT_PLUME
237             CALL DIAGNOSTICS_FILL(boplume,'KPPboplm',0,1,1,bi,bj,myThid)
238    #endif /* ALLOW_SALT_PLUME */
239          ENDIF
240    #endif /* ALLOW_DIAGNOSTICS */
241    
242  c------------------------------------------------------------------------  c------------------------------------------------------------------------
243  c     velocity shear  c     velocity shear
244  c     --------------  c     --------------
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22