Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/SOL/sol_oce_tam.F90

Timestamp:

2011-02-15T12:58:59+01:00 (13 years ago)

Author:

vidard

Message:

refer to ticket #798

File:

• branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/SOL/sol_oce_tam.F90 (modified) (5 diffs)

branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/SOL/sol_oce_tam.F90

@@ -31 +31 @@
       & nmax_fs
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+   USE sol_oce       , ONLY: &
+      & nmin
+#if defined key_dynspg_flt
    USE solver        , ONLY: & ! Solver
       & solver_init
-
+#endif
    !! * Routine accessibility
    
@@ -41 +44 @@
    PUBLIC &
       & sol_oce_tam_init, & !: routine called by nemovar.F90
+      & sol_oce_tam_deallocate, &
                     !
       & gcx_tl,   & !: Tangent of now solution of the elliptic equation
@@ -48 +52 @@
       & gcx_ad,   & !: Adjoint of solution of the elliptic equation
       & gcxb_ad,  & !: Adjoint of before solution of the elliptic equation
-      & gcb_ad      !: Adjoint of 2nd member of barotropic linear system
-
+      & gcb_ad ,  & !: Adjoint of 2nd member of barotropic linear system
+      & nitsor
    !! * Module variables
    INTEGER, DIMENSION(:), ALLOCATABLE :: &
@@ -89 +93 @@
          & kindic        ! indicate which variables to allocate/initialize
 
+#if defined key_dynspg_flt
       IF ( kindic == 0 .OR. kindic == 1 ) THEN
 
+         IF ( kindic == 0 ) CALL solver_init( nit000 )
+
          IF ( .NOT. ALLOCATED(nitsor) ) THEN
 
             ALLOCATE( nitsor( nitend - nit000 + 1 ) )
-
-         ENDIF
-
-         nitsor(:) = nmax_fs
-
-      IF ( kindic == 0 ) CALL solver_init( nit000 )
+            nitsor(:) = nmin
+
+         ENDIF
+
 
       ENDIF
@@ -162 +167 @@
 
       ENDIF
-
+#endif
    END SUBROUTINE sol_oce_tam_init
-
+   SUBROUTINE sol_oce_tam_deallocate(kindic)
+      !!-----------------------------------------------------------------------
+      !!
+      !!                  ***  ROUTINE sol_oce_tam_deallocate  ***
+      !!
+      !! ** Purpose : 
+      !!
+      !! ** Method  : kindic = 0  deallocate both tl and ad variables
+      !!              kindic = 1  deallocate only tl variables
+      !!              kindic = 2  deallocate only ad variables
+      !!
+      !! ** Action  :
+      !!                   
+      !! References : 
+      !!
+      !! History :
+      !!         ! 2010-06 (A. Vidard) Initial version
+      !!-----------------------------------------------------------------------
+      !! * Arguments
+      INTEGER, INTENT(IN) :: &
+         & kindic        ! indicate which variables to allocate/initialize
+#if defined key_dynspg_flt
+      IF ( kindic == 0 ) THEN
+         IF ( ALLOCATED(nitsor) )  DEALLOCATE( nitsor )
+      END IF
+
+      IF ( kindic == 0 .OR. kindic == 1 ) THEN
+
+         IF ( ALLOCATED(gcx_tl) )  DEALLOCATE( gcx_tl  )
+
+         IF ( ALLOCATED(gcxb_tl) ) DEALLOCATE( gcxb_tl )
+
+         IF ( ALLOCATED(gcb_tl) )  DEALLOCATE( gcb_tl  )
+
+         IF ( ALLOCATED(gcr_tl) )  DEALLOCATE( gcr_tl  )
+
+      ENDIF
+
+      IF ( kindic == 0 .OR. kindic == 2 ) THEN
+
+         IF ( ALLOCATED(gcx_ad) )  DEALLOCATE( gcx_ad  )
+
+         IF ( ALLOCATED(gcxb_ad) ) DEALLOCATE( gcxb_ad )
+
+         IF ( ALLOCATED(gcb_ad) )  DEALLOCATE( gcb_ad  )
+
+      ENDIF
+      
+#endif
+   END SUBROUTINE sol_oce_tam_deallocate
 END MODULE sol_oce_tam