Changeset 4292 for branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domstp.F90

Timestamp:

2013-11-20T17:28:04+01:00 (10 years ago)

Author:

cetlod

Message:

dev_MERGE_2013 : 1st step of the merge, see ticket #1185

File:

• branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domstp.F90 (modified) (1 diff)

branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domstp.F90

@@ -90 +90 @@
 
             DO jk = 1, jpk
-               IF( gdept_0(jk) <= rdth ) rdttra(jk) = rdtmin
-               IF( gdept_0(jk) >  rdth ) THEN
+               IF( gdept_1d(jk) <= rdth ) rdttra(jk) = rdtmin
+               IF( gdept_1d(jk) >  rdth ) THEN
                   rdttra(jk) = rdtmin + ( rdtmax - rdtmin )   &
-                                      * ( EXP( ( gdept_0(jk ) - rdth ) / rdth ) - 1. )   &
-                                      / ( EXP( ( gdept_0(jpk) - rdth ) / rdth ) - 1. )
+                                      * ( EXP( ( gdept_1d(jk ) - rdth ) / rdth ) - 1. )   &
+                                      / ( EXP( ( gdept_1d(jpk) - rdth ) / rdth ) - 1. )
                ENDIF
                IF(lwp) WRITE(numout,"(36x,f5.2,5x,i3)") rdttra(jk)/3600., jk