Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/TRA/zpshde_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/zpshde_tam.F90 (modified) (11 diffs)

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

@@ -78 +78 @@
    PUBLIC zps_hde_adj      ! routine called by step_tam.F90
    PUBLIC zps_hde_adj_tst  ! routine called by tst.F90
+#if defined key_tst_tlm
    PUBLIC  zps_hde_tlm_tst ! routine called by tamtst.F90
+#endif
 
    !! * module variables
@@ -946 +948 @@
       
    END SUBROUTINE zps_hde_adj_tst
-
+#if defined key_tst_tlm
    SUBROUTINE zps_hde_tlm_tst( kumadt )
       !!-----------------------------------------------------------------------
@@ -984 +986 @@
       USE tamtrj              ! writing out state trajectory
       USE par_tlm,    ONLY: &
+        & tlm_bch,          &
         & cur_loop,         &
         & h_ratio
@@ -1049 +1052 @@
          & 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):: &
@@ -1104 +1107 @@
       zgru_out(:,:)   = 0.0_wp    
       zgrv_out(:,:)   = 0.0_wp
+      IF ( tlm_bch == 2 ) THEN
       gtu_tl(:,:)     = 0.0_wp
       gtv_tl(:,:)     = 0.0_wp    
@@ -1110 +1114 @@
       gru_tl(:,:)     = 0.0_wp    
       grv_tl(:,:)     = 0.0_wp 
-
+      ENDIF
       zscgtu(:)         = 0.0_wp
       zscgtv(:)         = 0.0_wp
@@ -1132 +1136 @@
       ! Output filename Xn=F(X0)
       !--------------------------------------------------------------------
-      file_wop='trj_wop_zps'
       CALL tlm_namrd
       gamma = h_ratio
+      file_wop='trj_wop_zps'
+      file_xdx='trj_xdx_zps'
       !--------------------------------------------------------------------
       ! Initialize the tangent input with random noise: dx
@@ -1167 +1172 @@
       ! Complete Init for Direct
       !-------------------------------------------------------------------
-      CALL istate_p  
-
+      IF ( tlm_bch /= 2 )  CALL istate_p  
       ! *** initialize the reference trajectory
       ! ------------
@@ -1187 +1191 @@
       !  Compute the direct model F(X0,t=n) = Xn
       !--------------------------------------------------------------------
-      CALL zps_hde(nit000, tn, sn, rhd, gtu, gsu, gru, gtv, gsv, grv)
-
-      IF ( cur_loop .EQ. 0) CALL trj_wri_spl(file_wop)
-
+      IF ( tlm_bch /= 2 ) CALL zps_hde(nit000, tn, sn, rhd, gtu, gsu, gru, gtv, gsv, grv)
+      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
-         !--------------------------------------------------------------------  
-         zgtu_out  (:,:) = gtu   (:,:)
-         zgtv_out  (:,:) = gtv   (:,:)  
-         zgsu_out  (:,:) = gsu   (:,:)
-         zgsv_out  (:,:) = gsv   (:,:) 
-         zgru_out  (:,:) = gru   (:,:)
-         zgrv_out  (:,:) = grv   (:,:)         
-
+      IF ( tlm_bch == 2 ) THEN
          !--------------------------------------------------------------------
          ! Initialize the tangent variables:
@@ -1241 +1234 @@
          zgru_wop  (:,:) = gru  (:,:)
          zgrv_wop  (:,:) = grv  (:,:)
+         CALL trj_rd_spl(file_xdx) 
+         zgtu_out  (:,:) = gtu  (:,:)
+         zgtv_out  (:,:) = gtv  (:,:)
+         zgsu_out  (:,:) = gsu  (:,:)
+         zgsv_out  (:,:) = gsv  (:,:)
+         zgru_out  (:,:) = gru  (:,:)
+         zgrv_out  (:,:) = grv  (:,:)
          !--------------------------------------------------------------------
          ! Compute the Linearization Error 
@@ -1471 +1471 @@
    !!======================================================================
 #endif
+#endif
 END MODULE zpshde_tam