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

Timestamp:

2005-03-17T15:02:38+01:00 (19 years ago)

Author:

opalod

Message:

CT : UPDATE151 : New trends organization

File:

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

trunk/NEMO/OPA_SRC/TRA/trasbc.F90

@@ -13 +13 @@
    USE ocesbc          ! surface thermohaline fluxes
    USE phycst          ! physical constant
-   USE traqsr          ! ???
-   USE trdtra_oce      ! ocean active tracer trends 
+   USE traqsr          ! solar radiation penetration
+   USE trdmod          ! ocean trends 
+   USE trdmod_oce      ! ocean variables trends
    USE in_out_manager  ! I/O manager
 
@@ -75 +76 @@
       !!   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
       !!----------------------------------------------------------------------
+      !! * 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
@@ -96 +102 @@
 #endif
 
+      ! Save ta and sa trends
+      IF( l_trdtra )   THEN
+         ztdta(:,:,:) = ta(:,:,:) 
+         ztdsa(:,:,:) = sa(:,:,:) 
+      ENDIF
+
       IF( .NOT.ln_traqsr )   qsr(:,:) = 0.e0   ! no solar radiation penetration
 
@@ -113 +125 @@
             ta(ji,jj,1) = ta(ji,jj,1) + zta
             sa(ji,jj,1) = sa(ji,jj,1) + zsa
-#if defined key_trdtra || defined key_trdmld
-            ! save the trend associated to (t,s) sea surface boundary condition
-            flxtrd(ji,jj,1) = zta
-            flxtrd(ji,jj,2) = zsa
-#endif
          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)
+      ENDIF
       
       IF(l_ctl) THEN         ! print mean trends (used for debugging)