Changeset 503 for trunk/NEMO/OPA_SRC/TRA/trasbc.F90

Timestamp:

2006-09-27T10:52:29+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_064 : CT : general trends update including the addition of mean windows analysis possibility in the mixed layer

File:

• trunk/NEMO/OPA_SRC/TRA/trasbc.F90 (modified) (6 diffs)

trunk/NEMO/OPA_SRC/TRA/trasbc.F90

@@ -4 +4 @@
    !! Ocean active tracers:  surface boundary condition
    !!==============================================================================
-   !! History :
-   !!   8.2  !  98-10  (G. Madec, G. Roullet, M. Imbard)  Original code
-   !!   8.2  !  01-02  (D. Ludicone)  sea ice and free surface
-   !!   8.5  !  02-06  (G. Madec)  F90: Free form and module
-   !!   9.0  !  04-08  (C. Talandier) New trends organization
+   !! History :  8.2  !  98-10  (G. Madec, G. Roullet, M. Imbard)  Original code
+   !!            8.2  !  01-02  (D. Ludicone)  sea ice and free surface
+   !!            8.5  !  02-06  (G. Madec)  F90: Free form and module
+   !!----------------------------------------------------------------------
+
    !!----------------------------------------------------------------------
    !!   tra_sbc      : update the tracer trend at ocean surface
    !!----------------------------------------------------------------------
-   !! * Modules used
    USE oce             ! ocean dynamics and active tracers
    USE dom_oce         ! ocean space domain variables
@@ -26 +25 @@
    PRIVATE
 
-   !! * Routine accessibility
-   PUBLIC tra_sbc              ! routine called by step.F90
+   PUBLIC   tra_sbc    ! routine called by step.F90
 
    !! * Substitutions
@@ -35 +33 @@
    !!   OPA 9.0 , LOCEAN-IPSL (2005) 
    !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 
@@ -78 +76 @@
       !!                with the tracer surface boundary condition 
       !!              - save the trend it in ttrd ('key_trdtra')
+      !!----------------------------------------------------------------------
+      USE oce, ONLY :   ztrdt => ua   ! use ua as 3D workspace   
+      USE oce, ONLY :   ztrds => va   ! use va as 3D workspace   
       !!
-      !!----------------------------------------------------------------------
-      !! * Modules used     
-      USE oce, ONLY :    ztdta => ua,      & ! use ua as 3D workspace   
-                         ztdsa => va         ! use va as 3D workspace   
-
-      !! * Arguments
-      INTEGER, INTENT( in ) ::   kt          ! ocean time-step index
-
-      !! * Local declarations
+      INTEGER, INTENT(in) ::   kt     ! ocean time-step index
+      !!
       INTEGER  ::   ji, jj                   ! dummy loop indices
       REAL(wp) ::   zta, zsa, zsrau, zse3t   ! temporary scalars
@@ -98 +92 @@
       ENDIF
 
-      ! 0. initialization
-      zsrau = 1. / rauw
+      zsrau = 1. / rauw             ! initialization
 #if defined key_zco
       zse3t = 1. / e3t_0(1)
 #endif
 
-      ! Save ta and sa trends
-      IF( l_trdtra )   THEN
-         ztdta(:,:,:) = ta(:,:,:) 
-         ztdsa(:,:,:) = sa(:,:,:) 
+      IF( l_trdtra ) THEN           ! Save ta and sa trends
+         ztrdt(:,:,:) = ta(:,:,:) 
+         ztrds(:,:,:) = sa(:,:,:) 
       ENDIF
 
       IF( .NOT.ln_traqsr )   qsr(:,:) = 0.e0   ! no solar radiation penetration
 
-      ! 1. Concentration dillution effect on (t,s)
+      ! Concentration dillution effect on (t,s)
       DO jj = 2, jpj
          DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -118 +110 @@
             zse3t = 1. / fse3t(ji,jj,1)
 #endif
-            ! temperature : heat flux
-            zta = ro0cpr * ( qt(ji,jj) - qsr(ji,jj) ) * zse3t
-
-            ! salinity :  concent./dilut. effect
-            zsa = emps(ji,jj) * zsrau * sn(ji,jj,1) * zse3t
-            
-            ! add the trend to the general tracer trend
-            ta(ji,jj,1) = ta(ji,jj,1) + zta
+            zta = ro0cpr * ( qt(ji,jj) - qsr(ji,jj) ) * zse3t      ! temperature : heat flux
+            zsa = emps(ji,jj) * zsrau * sn(ji,jj,1)   * zse3t      ! salinity :  concent./dilut. effect
+            ta(ji,jj,1) = ta(ji,jj,1) + zta                        ! add the trend to the general tracer trend
             sa(ji,jj,1) = sa(ji,jj,1) + zsa
          END DO
       END DO
 
-      ! save the trends for diagnostic
-      ! sea surface boundary condition tracers trends
-      IF( l_trdtra )   THEN
-         ztdta(:,:,:) = ta(:,:,:) - ztdta(:,:,:)
-         ztdsa(:,:,:) = sa(:,:,:) - ztdsa(:,:,:)
-         CALL trd_mod(ztdta, ztdsa, jpttdnsr, 'TRA', kt)
+      IF( l_trdtra ) THEN      ! save the sbc trends for diagnostic
+         ztrdt(:,:,:) = ta(:,:,:) - ztrdt(:,:,:)
+         ztrds(:,:,:) = sa(:,:,:) - ztrds(:,:,:)
+         CALL trd_mod(ztrdt, ztrds, jptra_trd_nsr, 'TRA', kt)
       ENDIF
-      
-      IF(ln_ctl) THEN         ! print mean trends (used for debugging)
-         CALL prt_ctl(tab3d_1=ta, clinfo1=' sbc  - Ta: ', mask1=tmask, &
-            &         tab3d_2=sa, clinfo2=' Sa: ', mask2=tmask, clinfo3='tra')
-      ENDIF
-
+      !
+      IF(ln_ctl)   CALL prt_ctl( tab3d_1=ta, clinfo1=' sbc  - Ta: ', mask1=tmask,   &
+         &                       tab3d_2=sa, clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' )
+      !
    END SUBROUTINE tra_sbc