Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/DYN/dynadv_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/DYN/dynadv_tam.F90 (modified) (8 diffs)

branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/DYN/dynadv_tam.F90

@@ -547 +547 @@
       END IF
    END SUBROUTINE dyn_adv_ctl_tam
-
+#if defined key_tst_tlm
    SUBROUTINE dyn_adv_tlm_tst( kumadt )
       !!-----------------------------------------------------------------------
@@ -580 +580 @@
       USE tamtrj              ! writing out state trajectory
       USE par_tlm,    ONLY: &
+        & tlm_bch,          &
         & cur_loop,         &
         & h_ratio
@@ -641 +642 @@
 
       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):: &
@@ -696 +697 @@
       ! Output filename Xn=F(X0)
       !--------------------------------------------------------------------
-      file_wop='trj_wop_dynadv'
-
       CALL tlm_namrd
       gamma = h_ratio
-
+      file_wop='trj_wop_dynadv'
+      file_xdx='trj_xdx_dynadv'
       !--------------------------------------------------------------------
       ! Initialize the tangent input with random noise: dx
@@ -741 +741 @@
       ! Complete Init for Direct
       !-------------------------------------------------------------------
-      CALL istate_p  
+      IF ( tlm_bch /= 2 )      CALL istate_p  
 
       ! *** initialize the reference trajectory
@@ -769 +769 @@
       !  Compute the direct model F(X0,t=n) = Xn
       !--------------------------------------------------------------------
-      CALL dyn_adv(nit000)
-      IF ( cur_loop .EQ. 0) CALL trj_wri_spl(file_wop)
+      IF ( tlm_bch /= 2 ) CALL dyn_adv(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
-         !--------------------------------------------------------------------  
-         zua_out  (:,:,:) = ua   (:,:,:)
-         zva_out  (:,:,:) = va   (:,:,:)          
-
+      IF ( tlm_bch == 2 ) THEN        
          !--------------------------------------------------------------------
          ! Initialize the tangent variables 
@@ -812 +807 @@
          zua_wop  (:,:,:) = ua  (:,:,:)
          zva_wop  (:,:,:) = va  (:,:,:)
+         CALL trj_rd_spl(file_xdx) 
+         zua_out  (:,:,:) = ua  (:,:,:)
+         zva_out  (:,:,:) = va  (:,:,:)
          !--------------------------------------------------------------------
          ! Compute the Linearization Error 
@@ -944 +942 @@
    END SUBROUTINE dyn_adv_tlm_tst
 #endif
+#endif
 END MODULE dynadv_tam