Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/TRA/traadv_cen2_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/traadv_cen2_tam.F90 (modified) (9 diffs)

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

@@ -127 +127 @@
    PUBLIC   tra_adv_cen2_adj    ! routine called by traadv_tam.F90
    PUBLIC   tra_adv_cen2_adj_tst! routine called by tst.F90
+#if defined key_tst_tlm
    PUBLIC   tra_adv_cen2_tlm_tst! routine called by tamtst.F90
+#endif
 
    REAL(wp), DIMENSION(jpi,jpj) ::  &
@@ -935 +937 @@
 
    END SUBROUTINE tra_adv_cen2_adj_tst
-
+#if defined key_tst_tlm
 SUBROUTINE tra_adv_cen2_tlm_tst( kumadt )
       !!-----------------------------------------------------------------------
@@ -968 +970 @@
       USE tamtrj              ! writing out state trajectory
       USE par_tlm,    ONLY: &
+        & tlm_bch,          &
         & cur_loop,         &
         & h_ratio
@@ -1027 +1030 @@
          & z3r
       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):: &
@@ -1084 +1087 @@
       ! Output filename Xn=F(X0)
       !--------------------------------------------------------------------
-      file_wop='trj_wop_tradv_cen2'
       CALL tlm_namrd
       gamma = h_ratio     
+      file_wop='trj_wop_tradv_cen2'
+      file_xdx='trj_xdx_tradv_cen2'   
       !--------------------------------------------------------------------
       ! Initialize the tangent input with random noise: dx
@@ -1151 +1155 @@
       ! Complete Init for Direct
       !-------------------------------------------------------------------
-      CALL istate_p  
+      IF ( tlm_bch /= 2 ) CALL istate_p  
 
       ! *** initialize the reference trajectory
@@ -1184 +1188 @@
       !  Compute the direct model F(X0,t=n) = Xn
       !--------------------------------------------------------------------
-      CALL tra_adv_cen2(nit000, un, vn, wn)
-      IF ( cur_loop .EQ. 0) CALL trj_wri_spl(file_wop)
+      IF ( tlm_bch /= 2 ) CALL tra_adv_cen2(nit000, un, vn, wn)
+      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  
@@ -1214 +1213 @@
          ! Compute the scalar product: ( L(t0,tn) gamma dx0 ) )
          !--------------------------------------------------------------------
-
          zsp2_Ta    = DOT_PRODUCT( ta_tl, ta_tl  )
          zsp2_Sa    = DOT_PRODUCT( sa_tl, sa_tl  )
 
          zsp2      = zsp2_Ta + zsp2_Sa
-  
          !--------------------------------------------------------------------
          !  Storing data
          !--------------------------------------------------------------------
          CALL trj_rd_spl(file_wop) 
-
          zta_wop  (:,:,:) = ta  (:,:,:)
          zsa_wop  (:,:,:) = sa  (:,:,:)
-
+         CALL trj_rd_spl(file_xdx) 
+         zta_out  (:,:,:) = ta  (:,:,:)
+         zsa_out  (:,:,:) = sa  (:,:,:)
          !--------------------------------------------------------------------
          ! Compute the Linearization Error 
@@ -1372 +1370 @@
   END SUBROUTINE  tra_adv_cen2_tlm_tst
 #endif
-
+#endif
    !!======================================================================
 END MODULE traadv_cen2_tam