Changeset 11903 for NEMO/branches

Timestamp:

2019-11-14T11:34:15+01:00 (4 years ago)

Author:

davestorkey

Message:

UKMO/NEMO_4.0_momentum_trends : modification to bottom friction trend diagnostic.

File:

• NEMO/branches/UKMO/NEMO_4.0_momentum_trends/src/OCE/TRD/trddyn.F90 (modified) (4 diffs)

NEMO/branches/UKMO/NEMO_4.0_momentum_trends/src/OCE/TRD/trddyn.F90

@@ -40 +40 @@
    REAL(wp), ALLOCATABLE, DIMENSION(:,:,:), SAVE :: zutrd_hpg, zvtrd_hpg
    REAL(wp), ALLOCATABLE, DIMENSION(:,:,:), SAVE :: zutrd_pvo, zvtrd_pvo
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:,:), SAVE :: zutrd_bfr, zvtrd_bfr
 
    !! * Substitutions
@@ -61 +62 @@
       INTEGER                   , INTENT(in   ) ::   ktrd           ! trend index
       INTEGER                   , INTENT(in   ) ::   kt             ! time step
+      REAL(wp), ALLOCATABLE, DIMENSION(:,:)     ::   zue, zve       ! temporary 2D arrays
+      INTEGER                                   ::   jk
       !!----------------------------------------------------------------------
       !
@@ -93 +96 @@
             pvtrd(:,:,:) = pvtrd(:,:,:) + zvtrd_hpg(:,:,:) 
             DEALLOCATE( zutrd_hpg, zvtrd_hpg )
+         ENDIF
+
+         IF( ktrd == jpdyn_bfr ) THEN
+            !
+            ! Add 3D part of BFR trend minus its depth-mean part to depth-mean trend already saved.
+            ALLOCATE( zue(jpi,jpj), zve(jpi,jpj) )
+            zue(:,:) = e3u_a(:,:,1) * putrd(:,:,1) * umask(:,:,1)
+            zve(:,:) = e3v_a(:,:,1) * pvtrd(:,:,1) * vmask(:,:,1)
+            DO jk = 2, jpkm1
+               zue(:,:) = zue(:,:) + e3u_a(:,:,jk) * putrd(:,:,jk) * umask(:,:,jk)
+               zve(:,:) = zve(:,:) + e3v_a(:,:,jk) * pvtrd(:,:,jk) * vmask(:,:,jk)
+            END DO
+            DO jk = 1, jpkm1
+               putrd(:,:,jk) = zutrd_bfr(:,:,jk) + putrd(:,:,jk) - zue(:,:) * r1_hu_a(:,:)
+               pvtrd(:,:,jk) = zvtrd_bfr(:,:,jk) + pvtrd(:,:,jk) - zve(:,:) * r1_hv_a(:,:)
+            END DO
+            DEALLOCATE( zue, zve, zutrd_bfr, zvtrd_bfr )
          ENDIF
 
@@ -187 +207 @@
           CALL trd_dyn_3d( zutrd_pvo, zvtrd_pvo, jpdyn_pvo, kt )
           DEALLOCATE( zutrd_pvo, zvtrd_pvo )
+
+      CASE( jpdyn_bfr )
+          !
+          ! Create 3D BFR trend from 2D field and also output 2D field as diagnostic in own right.
+          ALLOCATE( zutrd_bfr(jpi,jpj,jpk), zvtrd_bfr(jpi,jpj,jpk) )
+          zutrd_bfr(:,:,:) = 0.0
+          zvtrd_bfr(:,:,:) = 0.0
+          DO jk = 1, jpkm1
+             zutrd_bfr(:,:,jk) = putrd(:,:)
+             zvtrd_bfr(:,:,jk) = pvtrd(:,:)
+          ENDDO
+          CALL trd_dyn_iom_2d( putrd, pvtrd, jpdyn_bfr, kt ) 
 
       CASE DEFAULT