Changeset 1756 for trunk/NEMO/OPA_SRC/DYN/sshwzv.F90

Timestamp:

2009-11-25T15:15:20+01:00 (14 years ago)

Author:

smasson

Message:

implement AR5 diagnostics, see ticket:610

File:

• trunk/NEMO/OPA_SRC/DYN/sshwzv.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/DYN/sshwzv.F90

@@ -25 +25 @@
    USE obc_par         ! open boundary cond. parameter
    USE obc_oce
+   USE diaar5, ONLY :   lk_diaar5
    USE iom
 
@@ -65 +66 @@
       !!                hu, hv, hur, hvr : ocean depth and its inverse at u-,v-points
       !!----------------------------------------------------------------------
+      USE oce, ONLY :   z3d => ta   ! use ta as 3D workspace
+      !!
       INTEGER, INTENT(in) ::   kt   ! time step
       !!
@@ -71 +74 @@
       REAL(wp) ::   z2dt, zraur     ! temporary scalars
       REAL(wp), DIMENSION(jpi,jpj) ::   zhdiv       ! 2D workspace
+      REAL(wp), DIMENSION(jpi,jpj) ::   z2d         ! 2D workspace
       !!----------------------------------------------------------------------
 
@@ -193 +197 @@
       END DO
       !
-      CALL iom_put( "woce", wn   )                     ! vertical velocity
-      CALL iom_put( "ssh" , sshn )                     ! sea surface height
+      CALL iom_put( "woce", wn                    )   ! vertical velocity
+      CALL iom_put( "ssh" , sshn                  )   ! sea surface height
+      CALL iom_put( "ssh2", sshn(:,:) * sshn(:,:) )   ! square of sea surface height
+      IF( lk_diaar5 ) THEN
+         z2d(:,:) = rau0 * e1t(:,:) * e2t(:,:)
+         DO jk = 1, jpk
+            z3d(:,:,jk) = wn(:,:,jk) * z2d(:,:)
+         END DO
+         CALL iom_put( "w_masstr" , z3d                     )   !           vertical mass transport
+         CALL iom_put( "w_masstr2", z3d(:,:,:) * z3d(:,:,:) )   ! square of vertical mass transport
+      ENDIF
       !
    END SUBROUTINE ssh_wzv