Changes between Version 3 and Version 4 of ticket/0664_TRA_TRP

Timestamp:

2010-04-27T12:41:50+02:00 (14 years ago)

Author:

cetlod

Comment:

--

ticket/0664_TRA_TRP

@@ -20 +20 @@
 [[BR]]
 
-'''(1)''' Creation of advection/diffusion modules with 4D tracer array as input argument[[BR]]
+'''(1)''' Creation of advection/diffusion modules with 4D tracer array as input argument[[BR]][[BR]]
+
 
 {{{
@@ -37 +38 @@
       REAL(wp)        , INTENT(inout), DIMENSION(:,:,:,:)     ::   ptraa           ! tracer trend
 }}}
+[[BR]]
 
-'''(2)''' Test for T & S 3D arrays by using the fortran ''RESHAPE'' function [[BR]]
+'''(2)''' Test for T & S 3D arrays by using the fortran '''''RESHAPE''''' function [[BR]][[BR]]
+
 
 
@@ -81 +84 @@
 }}}
 
+[[BR]]
+
+'''(3)''' Apply to the passive tracers [[BR]][[BR]]
+
+
+
+{{{
+   SUBROUTINE trc_adv( kt )
+      !!----------------------------------------------------------------------
+      !!                  ***  ROUTINE trc_adv  ***
+      !!
+      !! ** Purpose :  
+      !!
+      !! ** Method  : 
+      !!----------------------------------------------------------------------
+      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
+      !!
+      INTEGER ::   jk   ! dummy loop index
+      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zun, zvn, zwn   ! effective transports
+      !!----------------------------------------------------------------------
+
+      !                                              ! effective transport
+      DO jk = 1, jpkm1
+#if defined key_trcbbl_adv                                      
+         !                                                ! eulerian + bbl transport
+         zun(:,:,jk) = e2u(:,:) * fse3u(:,:,jk) * ( un(:,:,jk) - u_trc_bbl(:,:,jk) )
+         zvn(:,:,jk) = e1v(:,:) * fse3v(:,:,jk) * ( vn(:,:,jk) - v_trc_bbl(:,:,jk) )
+         zwn(:,:,jk) = e1t(:,:) * e2t(:,:)      * ( wn(:,:,jk) + w_trc_bbl(:,:,jk) )
+#else
+         !                                                ! eulerian transport only
+         zun(:,:,jk) = e2u(:,:) * fse3u(:,:,jk) *   un(:,:,jk)
+         zvn(:,:,jk) = e1v(:,:) * fse3v(:,:,jk) *   vn(:,:,jk)
+         zwn(:,:,jk) = e1t(:,:) * e2t(:,:)      *   wn(:,:,jk)
+#endif
+      END DO
+      zwn(:,:,jpk) = 0.e0                                 ! no transport trough the bottom
+
+      !                                                   ! add the eiv transport (if necessary)
+      IF( lk_trcldf_eiv )   CALL trc_adv_eiv( kt, zun, zvn, zwn )
+
+      CALL tra_adv_cen2( kt, 'TRC', zun, zvn, zwn, trb, trn, tra )    ! 2nd order centered
+}}}
 
 ----