Changeset 2257 for branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/TRP

Timestamp:

2010-10-13T17:58:28+02:00 (14 years ago)

Author:

cetlod

Message:

Apply the modified leap-frog scheme on runoff & correct a bug on time stepping for hpg implicit case

File:

• branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90 (modified) (1 diff)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90

@@ -120 +120 @@
          ztrdt(:,:,:,:)  = trn(:,:,:,:)
       ENDIF
-
       ! Leap-Frog + Asselin filter time stepping
-      IF( lk_vvl ) THEN   ;   CALL tra_nxt_vvl( kt, 'TRC', trb, trn, tra, jptra )      ! variable volume level (vvl) 
-      ELSE                ;   CALL tra_nxt_fix( kt, 'TRC', trb, trn, tra, jptra )      ! fixed    volume level 
+      IF( neuler == 0 .AND. kt == nit000 ) THEN        ! Euler time-stepping at first time-step
+         !                                             ! (only swap)
+         DO jn = 1, jptra
+            DO jk = 1, jpkm1
+               trn(:,:,jk,jn) = tra(:,:,jk,jn)
+            END DO
+         END DO
+         !                                              
+      ELSE
+         ! Leap-Frog + Asselin filter time stepping
+         IF( lk_vvl ) THEN   ;   CALL tra_nxt_vvl( kt, 'TRC', trb, trn, tra, jptra )      ! variable volume level (vvl) 
+         ELSE                ;   CALL tra_nxt_fix( kt, 'TRC', trb, trn, tra, jptra )      ! fixed    volume level 
+         ENDIF
       ENDIF