Changeset 9630 for branches/UKMO/dev_isf_divg_corr_GO6_package_r9385/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

Timestamp:

2018-05-24T20:50:11+02:00 (6 years ago)

Author:

antsia

Message:

Code for divergence correction

File:

• branches/UKMO/dev_isf_divg_corr_GO6_package_r9385/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90 (modified) (1 diff)

branches/UKMO/dev_isf_divg_corr_GO6_package_r9385/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

@@ -220 +220 @@
       !
       !
+
+      !---start the tendency with tracer source/sink from divergence correction
+      IF ( ln_iscpl .AND. kt == nit000) THEN        
+         IF (lwp) WRITE(numout,*) 'add tendency tsa due to re-meshing at kt = nit000'
+         tsa(:,:,:,jp_tem) = tsa(:,:,:,jp_tem) + rhdivdiff_trac(:,:,:,jp_tem)*tmask(:,:,:)
+         tsa(:,:,:,jp_sal) = tsa(:,:,:,jp_sal) + rhdivdiff_trac(:,:,:,jp_sal)*tmask(:,:,:)
+      ENDIF
+
       !----------------------------------------
       !       Ice Shelf effects (ISF)