Changeset 5605 for branches/NERC/dev_r5589_marine_glacier_termini/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

Timestamp:

2015-07-16T20:05:07+02:00 (9 years ago)

Author:

mathiot

Message:

Marine glacier termini: initial commit

File:

• branches/NERC/dev_r5589_marine_glacier_termini/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90 (modified) (2 diffs)

branches/NERC/dev_r5589_marine_glacier_termini/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

@@ -26 +26 @@
    USE trdtra          ! trends manager: tracers 
    !
+   USE traicw          ! ice wall parametrization
    USE in_out_manager  ! I/O manager
    USE prtctl          ! Print control
@@ -287 +288 @@
          END DO  
       ENDIF
+
+
+      ! compute trend
+      IF ( ln_traicw ) THEN
+         tsa(:,:,:,1) = tsa(:,:,:,1) + 0.5_wp * ( ricw_tsc(:,:,:,1) + ricw_tsc_b(:,:,:,1) )
+         tsa(:,:,:,2) = tsa(:,:,:,2) + 0.5_wp * ( ricw_tsc(:,:,:,2) + ricw_tsc_b(:,:,:,2) )
+      END IF
  
       IF( l_trdtra )   THEN                      ! send trends for further diagnostics