Changeset 14717 for NEMO/releases/r4.0/r4.0-HEAD/src/OCE/SBC/sbcssr.F90

Timestamp:

2021-04-16T11:42:56+02:00 (3 years ago)

Author:

clem

Message:

4.0-HEAD: correctly handle diagnostics of mass, salt and heat budgets (see ticket #2652). And fix Pierre ticket #2642

File:

• NEMO/releases/r4.0/r4.0-HEAD/src/OCE/SBC/sbcssr.F90 (modified) (2 diffs)

NEMO/releases/r4.0/r4.0-HEAD/src/OCE/SBC/sbcssr.F90

@@ -91 +91 @@
          IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN      !    Add restoring term     !
             !                                      ! ========================= !
+            qrp(:,:) = 0._wp ! necessary init
+            erp(:,:) = 0._wp
             !
             IF( nn_sstr == 1 ) THEN                                   !* Temperature restoring term
@@ -140 +142 @@
                      qns(ji,jj) = qns(ji,jj) - zerp * rcp * sst_m(ji,jj)
                      erp(ji,jj) = zerp
+                     qrp(ji,jj) = qrp(ji,jj) - zerp * rcp * sst_m(ji,jj)
                   END DO
                END DO
             ENDIF
+            !
+            ! outputs
+            CALL iom_put( 'hflx_ssr_cea', qrp(:,:) )
+            IF( nn_sssr == 1 )   CALL iom_put( 'sflx_ssr_cea',  erp(:,:) * sss_m(:,:) )
+            IF( nn_sssr == 2 )   CALL iom_put( 'vflx_ssr_cea', -erp(:,:) )
             !
          ENDIF