Changeset 8542

Timestamp:

2017-09-19T12:19:04+02:00 (7 years ago)

Author:

cetlod

Message:

trunk : minor bugfixes on initialisation sequence of passive tracers

Location:

trunk/NEMOGCM/NEMO/TOP_SRC

Files:

• trcdta.F90 (modified) (1 diff)
• trcini.F90 (modified) (2 diffs)

trunk/NEMOGCM/NEMO/TOP_SRC/trcdta.F90

@@ -227 +227 @@
                      ik = mbkt(ji,jj) 
                      IF( ik > 1 ) THEN
-                        zl = ( gdept_1d(ik) - gdept_n(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) )
+                        zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) )
                         ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik-1)
                      ENDIF
                      ik = mikt(ji,jj)
                      IF( ik > 1 ) THEN
-                        zl = ( gdept_n(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) )
+                        zl = ( gdept_0(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) )
                         ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik+1)
                      ENDIF

trunk/NEMOGCM/NEMO/TOP_SRC/trcini.F90

@@ -69 +69 @@
       !
       CALL trc_ini_sms   ! SMS
-      CALL trc_ini_inv   ! Inventories
       CALL trc_ini_trp   ! passive tracers transport
       CALL trc_ice_ini   ! Tracers in sea ice
@@ -78 +77 @@
       IF( nn_dttrc /= 1 ) &
       CALL trc_sub_ini    ! Initialize variables for substepping passive tracers
+      !
+      CALL trc_ini_inv   ! Inventories
       !
       IF( nn_timing == 1 )   CALL timing_stop('trc_init')