Changeset 3805

Timestamp:

2013-02-12T18:53:37+01:00 (11 years ago)

Author:

acc

Message:

Branch dev_MERGE_2012. #1043 also correct sign of diagnostic poleward diffusive transport

Location:

branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/TRA

Files:

• traldf_bilapg.F90 (modified) (1 diff)
• traldf_iso.F90 (modified) (1 diff)

branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilapg.F90

@@ -248 +248 @@
          ! "Poleward" diffusive heat or salt transport
          IF( cdtype == 'TRA' .AND. ln_diaptr .AND. ( kaht == 2 ) .AND. ( MOD( kt, nn_fptr ) == 0 ) ) THEN
-            IF( jn == jp_tem)   htr_ldf(:) = ptr_vj( zftv(:,:,:) )
-            IF( jn == jp_sal)   str_ldf(:) = ptr_vj( zftv(:,:,:) )
+            ! note sign is reversed to give down-gradient diffusive transports (#1043)
+            IF( jn == jp_tem)   htr_ldf(:) = ptr_vj( -zftv(:,:,:) )
+            IF( jn == jp_sal)   str_ldf(:) = ptr_vj( -zftv(:,:,:) )
          ENDIF
 

branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso.F90

@@ -212 +212 @@
          ! "Poleward" diffusive heat or salt transports (T-S case only)
          IF( cdtype == 'TRA' .AND. ln_diaptr .AND. ( MOD( kt, nn_fptr ) == 0 ) ) THEN
-            IF( jn == jp_tem)   htr_ldf(:) = ptr_vj( zftv(:,:,:) )
-            IF( jn == jp_sal)   str_ldf(:) = ptr_vj( zftv(:,:,:) )
+            ! note sign is reversed to give down-gradient diffusive transports (#1043)
+            IF( jn == jp_tem)   htr_ldf(:) = ptr_vj( -zftv(:,:,:) )
+            IF( jn == jp_sal)   str_ldf(:) = ptr_vj( -zftv(:,:,:) )
          ENDIF