Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/TRA/trasbc_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/TRA/trasbc_tam.F90 (modified) (10 diffs)

branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/TRA/trasbc_tam.F90

@@ -96 +96 @@
    PUBLIC   tra_sbc_adj     ! routine called by step_tam.F90
    PUBLIC   tra_sbc_adj_tst ! routine called by tst.F90
+#if defined key_tst_tlm
    PUBLIC   tra_sbc_tlm_tst ! routine calle  by tamtst.F90
+#endif
 
    !! * Substitutions
@@ -579 +581 @@
    END SUBROUTINE tra_sbc_adj_tst
 
-
+#if defined key_tst_tlm
    SUBROUTINE tra_sbc_tlm_tst ( kumadt ) 
       !!-----------------------------------------------------------------------
@@ -612 +614 @@
       USE tamtrj              ! writing out state trajectory
       USE par_tlm,    ONLY: &
+        & tlm_bch,          &
         & cur_loop,         &
         & h_ratio
@@ -676 +679 @@
          & zgsp7 
       CHARACTER (LEN=14)  :: cl_name
-      CHARACTER (LEN=128) :: file_out, file_wop
+      CHARACTER (LEN=128) :: file_out, file_wop, file_xdx
       CHARACTER (LEN=90)  ::  FMT
       REAL(KIND=wp), DIMENSION(100):: &
@@ -731 +734 @@
       ! Output filename Xn=F(X0)
       !--------------------------------------------------------------------
-      file_wop='trj_wop_trasbc'
       CALL tlm_namrd
       gamma = h_ratio
+      file_wop='trj_wop_trasbc'
+      file_xdx='trj_xdx_trasbc'
       !--------------------------------------------------------------------
       ! Initialize the tangent input with random noise: dx
@@ -778 +782 @@
       ! Complete Init for Direct
       !-------------------------------------------------------------------
-      CALL istate_p  
+      IF ( tlm_bch /= 2 )      CALL istate_p  
 
       ! *** initialize the reference trajectory
@@ -804 +808 @@
       !  Compute the direct model F(X0,t=n) = Xn
       !--------------------------------------------------------------------
-      CALL tra_sbc(nit000)
-
-      IF ( cur_loop .EQ. 0) CALL trj_wri_spl(file_wop)
-
+      IF ( tlm_bch /= 2 ) CALL tra_sbc(nit000)
+      IF ( tlm_bch == 0 ) CALL trj_wri_spl(file_wop)
+      IF ( tlm_bch == 1 ) CALL trj_wri_spl(file_xdx)
       !--------------------------------------------------------------------
       !  Compute the Tangent 
       !--------------------------------------------------------------------
-      IF ( cur_loop .NE. 0) THEN
-         !--------------------------------------------------------------------
-         !  Storing data
-         !--------------------------------------------------------------------  
-         zta_out  (:,:,:) = ta   (:,:,:)
-         zsa_out  (:,:,:) = sa   (:,:,:)          
-
+      IF ( tlm_bch == 2 ) THEN         
          !--------------------------------------------------------------------
          ! Initialize the tangent variables: dy^* = W dy  
@@ -836 +833 @@
          ! Compute the scalar product: ( L(t0,tn) gamma dx0 ) )
          !--------------------------------------------------------------------
-
          zsp2_1    = DOT_PRODUCT( ta_tl, ta_tl  )
          zsp2_2    = DOT_PRODUCT( sa_tl, sa_tl  )
@@ -847 +843 @@
          zta_wop  (:,:,:) = ta  (:,:,:)
          zsa_wop  (:,:,:) = sa  (:,:,:)
+         CALL trj_rd_spl(file_xdx) 
+         zta_out  (:,:,:) = ta  (:,:,:)
+         zsa_out  (:,:,:) = sa  (:,:,:)
          !--------------------------------------------------------------------
          ! Compute the Linearization Error 
@@ -981 +980 @@
 !!======================================================================
 #endif
+#endif
 END MODULE trasbc_tam