Changeset 1609

Timestamp:

2009-08-12T12:38:09+02:00 (15 years ago)

Author:

ctlod

Message:

Bug in dtatem and dtasal with key_vvl, see ticket: #528

Location:

trunk/NEMO/OPA_SRC/DTA

Files:

• dtasal.F90 (modified) (2 diffs)
• dtatem.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/DTA/dtasal.F90

@@ -238 +238 @@
                  DO ji = 1, jpi
                     DO jk = 1, jpk
-                       zl=fsdept(ji,jj,jk)
+                       zl=fsdept_0(ji,jj,jk)
                        IF(zl <  gdept_0(1)) zsaldta(jk,jl) =  saldta(ji,jj,1,jl)
                        IF(zl >  gdept_0(jpk)) zsaldta(jk,jl) =  saldta(ji,jj,jpkm1,jl)
@@ -271 +271 @@
                        ik = mbathy(ji,jj) - 1
                        IF( ik > 2 ) THEN
-                          zl = ( gdept_0(ik) - fsdept(ji,jj,ik) ) / ( gdept_0(ik) - gdept_0(ik-1) )
+                          zl = ( gdept_0(ik) - fsdept_0(ji,jj,ik) ) / ( gdept_0(ik) - gdept_0(ik-1) )
                           saldta(ji,jj,ik,jl) = (1.-zl) * saldta(ji,jj,ik,jl) +zl * saldta(ji,jj,ik-1,jl)
                        ENDIF

trunk/NEMO/OPA_SRC/DTA/dtatem.F90

@@ -252 +252 @@
                   DO ji = 1, jpi
                      DO jk = 1, jpk
-                        zl=fsdept(ji,jj,jk)
+                        zl=fsdept_0(ji,jj,jk)
                         IF(zl < gdept_0(1)) ztemdta(jk,jl) =  temdta(ji,jj,1,jl)
                         IF(zl > gdept_0(jpk)) ztemdta(jk,jl) =  temdta(ji,jj,jpkm1,jl)
@@ -286 +286 @@
                         ik = mbathy(ji,jj) - 1
                         IF( ik > 2 ) THEN
-                           zl = ( gdept_0(ik) - fsdept(ji,jj,ik) ) / ( gdept_0(ik) - gdept_0(ik-1) )
+                           zl = ( gdept_0(ik) - fsdept_0(ji,jj,ik) ) / ( gdept_0(ik) - gdept_0(ik-1) )
                            temdta(ji,jj,ik,jl) = (1.-zl) * temdta(ji,jj,ik,jl) + zl * temdta(ji,jj,ik-1,jl)
                         ENDIF