Changeset 2450 for branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DTA/dtasal.F90

Timestamp:

2010-12-04T16:20:50+01:00 (14 years ago)

Author:

gm

Message:

v3.3beta: #766 share the deepest ocean level indices continuaton

File:

• branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DTA/dtasal.F90 (modified) (1 diff)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DTA/dtasal.F90

@@ -190 +190 @@
             DO jj = 1, jpj               ! interpolation of salinity at the last ocean level (i.e. the partial step)
                DO ji = 1, jpi
-                  ik = mbathy(ji,jj) - 1
-                  IF( ik > 2 ) THEN
+                  ik = mbkt(ji,jj)
+                  IF( ik > 1 ) THEN
                      zl = ( gdept_0(ik) - fsdept_0(ji,jj,ik) ) / ( gdept_0(ik) - gdept_0(ik-1) )
                      s_dta(ji,jj,ik) = (1.-zl) * s_dta(ji,jj,ik) + zl * s_dta(ji,jj,ik-1)