Changeset 9075 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC

Timestamp:

2017-12-15T14:50:00+01:00 (6 years ago)

Author:

timgraham

Message:

Fix for uninitialised variable

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90 (modified) (2 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90

@@ -523 +523 @@
       ! Note that the "unclipped" bottom friction parameter is used even with explicit drag
       IF( ln_wd_il ) THEN
+         zztmp = -1._wp / rdtbt
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
-               zu_frc(ji,jj) = zu_frc(ji,jj) + MAX(r1_hu_n(ji,jj) * zCdU_u(ji,jj),-1._wp / rdtbt) * zwx(ji,jj) *  wdrampu(ji,jj)
-               zv_frc(ji,jj) = zv_frc(ji,jj) + MAX(r1_hv_n(ji,jj) * zCdU_v(ji,jj),-1._wp / rdtbt) * zwy(ji,jj) *  wdrampv(ji,jj)
+               zu_frc(ji,jj) = zu_frc(ji,jj) + MAX(r1_hu_n(ji,jj) * zCdU_u(ji,jj), zztmp ) * zwx(ji,jj) *  wdrampu(ji,jj)
+               zv_frc(ji,jj) = zv_frc(ji,jj) + MAX(r1_hv_n(ji,jj) * zCdU_v(ji,jj), zztmp ) * zwy(ji,jj) *  wdrampv(ji,jj)
             END DO
          END DO
@@ -540 +541 @@
       !                                         ! Add top stress contribution from baroclinic velocities:      
       IF( ln_bt_fw ) THEN
+         zztmp = -1._wp / rdtbt
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.