Changeset 15176 for NEMO/branches/UKMO/NEMO_4.0.4_momentum_trends/src/OCE/TRD/trddyn.F90

Timestamp:

2021-08-09T16:08:00+02:00 (3 years ago)

Author:

davestorkey

Message:

UKMO/NEMO_4.0.4_momentum_trends:

1. Output complete budget of 2D trends.
2. Separate out ice-ocean stress from wind stress (for 2D budget).
3. Allow thickness weighting of 2D trends.

File:

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

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

@@ -45 +45 @@
    REAL(wp), ALLOCATABLE, DIMENSION(:,:,:), SAVE :: zutrd_bfr, zvtrd_bfr
    REAL(wp), ALLOCATABLE, DIMENSION(:,:)  , SAVE :: zutrd_tau2d, zvtrd_tau2d
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:)  , SAVE :: zutrd_tau2d_ice, zvtrd_tau2d_ice
    REAL(wp), ALLOCATABLE, DIMENSION(:,:)  , SAVE :: zutrd_tfr2d, zvtrd_tfr2d
    REAL(wp), ALLOCATABLE, DIMENSION(:,:)  , SAVE :: zutrd_bfr2d, zvtrd_bfr2d
@@ -334 +335 @@
           DEALLOCATE( zutrd_pvo, zvtrd_pvo )
 
-      CASE( jpdyn_tau )
+      CASE( jpdyn_tau2d_ice )
+          !
+          ! Save 2D ice-ocean stress trend locally to be subtracted from
+          ! 2D wind stress trend and added to 2D top friction trend. 
+          IF( .NOT. ALLOCATED(zutrd_tau2d_ice) ) ALLOCATE( zutrd_tau2d_ice(jpi,jpj), zvtrd_tau2d_ice(jpi,jpj) )
+          zutrd_tau2d_ice(:,:) = putrd(:,:)
+          zvtrd_tau2d_ice(:,:) = pvtrd(:,:)
+
+      CASE( jpdyn_tau2d )
           !
           ! Save 2D wind forcing trend locally (to be added to ZDF trend)
@@ -341 +350 @@
           zutrd_tau2d(:,:) = putrd(:,:)
           zvtrd_tau2d(:,:) = pvtrd(:,:)
+          IF( ALLOCATED(zutrd_tau2d_ice) ) THEN
+             zutrd_tau2d(:,:) = zutrd_tau2d(:,:) - zutrd_tau2d_ice(:,:) 
+             zvtrd_tau2d(:,:) = zvtrd_tau2d(:,:) - zvtrd_tau2d_ice(:,:) 
+          ENDIF
 
       CASE( jpdyn_tfr )
@@ -348 +361 @@
           zutrd_tfr2d(:,:) = putrd(:,:)
           zvtrd_tfr2d(:,:) = pvtrd(:,:)
+          IF( ALLOCATED(zutrd_tau2d_ice) ) THEN
+             zutrd_tfr2d(:,:) = zutrd_tfr2d(:,:) + zutrd_tau2d_ice(:,:) 
+             zvtrd_tfr2d(:,:) = zvtrd_tfr2d(:,:) + zvtrd_tau2d_ice(:,:) 
+          ENDIF
 
       CASE( jpdyn_bfr )
@@ -480 +497 @@
       CASE( jpdyn_pvo_corr ) ;   CALL iom_put( "utrd_pvo_corr", putrd )   ! planetary vorticity correction
                                  CALL iom_put( "vtrd_pvo_corr", pvtrd )
-      CASE( jpdyn_tau )      ;   CALL iom_put( "utrd_tau2d", putrd )      ! wind stress depth-mean trend
+      CASE( jpdyn_frc2d )    ;   CALL iom_put( "utrd_frc2d", putrd )      ! constant forcing term from barotropic calcn.
+                                 CALL iom_put( "vtrd_frc2d", pvtrd ) 
+      CASE( jpdyn_tau2d )    ;   CALL iom_put( "utrd_tau2d", putrd )      ! wind stress depth-mean trend
                                  CALL iom_put( "vtrd_tau2d", pvtrd )
       CASE( jpdyn_bfr )      ;   CALL iom_put( "utrd_bfr2d", putrd )      ! bottom friction depth-mean trend
@@ -486 +505 @@
       CASE( jpdyn_tfr )      ;   CALL iom_put( "utrd_tfr2d", putrd )      ! top friction depth-mean trend
                                  CALL iom_put( "vtrd_tfr2d", pvtrd )
+      CASE( jpdyn_tot )      ;   CALL iom_put( "utrd_tot2d", putrd )      ! total 2D trend, excluding time filter
+                                 CALL iom_put( "vtrd_tot2d", pvtrd )
       END SELECT
       !