Changeset 413 for trunk/NEMO/OPA_SRC/SOL/solmat.F90

Timestamp:

2006-03-20T18:27:15+01:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_bugfix_032 : CT : use dt (linked to Euler time scheme) instead of 2*dt for the construction of the matrix used by the elliptic solvers

File:

• trunk/NEMO/OPA_SRC/SOL/solmat.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/SOL/solmat.F90

@@ -19 +19 @@
    USE lbclnk          ! lateral boudary conditions
    USE lib_mpp         ! distributed memory computing
+   USE in_out_manager  ! I/O manager
 
    IMPLICIT NONE
@@ -33 +34 @@
 CONTAINS
 
-   SUBROUTINE sol_mat
+   SUBROUTINE sol_mat( kt )
       !!----------------------------------------------------------------------
       !!                  ***  ROUTINE sol_mat  ***
@@ -67 +68 @@
       !!   9.0  !  05-11  (V. Garnier) Surface pressure gradient organization
       !!----------------------------------------------------------------------
+      !! * Arguments
+      INTEGER, INTENT(in) :: kt
+
       !! * Local declarations
       INTEGER ::   ji, jj                    ! dummy loop indices
@@ -93 +97 @@
       gcdmat(:,:) = 0.e0
       
-      z2dt = 2. * rdt
+      IF( neuler == 0 .AND. kt == nit000 ) THEN
+         z2dt = rdt
+      ELSE
+         z2dt = 2. * rdt
+      ENDIF
 
 #if defined key_dynspg_flt && ! defined key_obc