Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/DYN/dynhpg_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/dynhpg_tam.F90 (modified) (10 diffs)

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

@@ -109 +109 @@
    PUBLIC   dyn_hpg_adj    ! routine called by step_tam module
    PUBLIC   dyn_hpg_adj_tst! routine called by test module
+#if defined key_tst_tlm
    PUBLIC   dyn_hpg_tlm_tst! routine called by test module
+#endif
 
    !!* Namelist nam_dynhpg : Choice of horizontal pressure gradient computation
@@ -1059 +1061 @@
          &              )
    END SUBROUTINE dyn_hpg_adj_tst
-
+#if defined key_tst_tlm
    SUBROUTINE dyn_hpg_tlm_tst( kumadt )
       !!-----------------------------------------------------------------------
@@ -1089 +1091 @@
       USE tamtrj              ! writing out state trajectory
       USE par_tlm,    ONLY: &
+        & tlm_bch,          &
         & cur_loop,         &
         & h_ratio
@@ -1160 +1163 @@
          & jk
       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):: &
@@ -1228 +1231 @@
       ! Output filename Xn=F(X0)
       !--------------------------------------------------------------------
-      file_wop='trj_wop_dynhpg'
-
       CALL tlm_namrd
       gamma = h_ratio  
-
+      file_wop='trj_wop_dynhpg'
+      file_xdx='trj_xdx_dynhpg'
        !--------------------------------------------------------------------
       ! Initialize the tangent input with random noise: dx
@@ -1262 +1264 @@
       ! Complete Init for Direct
       !-------------------------------------------------------------------
-      CALL istate_p  
+      IF ( tlm_bch /= 2 )      CALL istate_p  
 
       ! *** initialize the reference trajectory
@@ -1296 +1298 @@
       !  Compute the direct model F(X0,t=n) = Xn
       !--------------------------------------------------------------------
-      CALL dyn_hpg(nit000)
-      IF ( cur_loop .EQ. 0) CALL trj_wri_spl(file_wop)
+      IF ( tlm_bch /= 2 ) CALL dyn_hpg(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 
@@ -1324 +1321 @@
          ! Compute the scalar product: ( L(t0,tn) gamma dx0 ) )
          !--------------------------------------------------------------------
-
          zsp2_1    = DOT_PRODUCT( ua_tl, ua_tl  )
          zsp2_2    = DOT_PRODUCT( va_tl, va_tl  )
@@ -1334 +1330 @@
          zua_wop  (:,:,:) = ua  (:,:,:)
          zva_wop  (:,:,:) = va  (:,:,:)
+         CALL trj_rd_spl(file_xdx) 
+         zua_out  (:,:,:) = ua  (:,:,:)
+         zva_out  (:,:,:) = va  (:,:,:)
          !--------------------------------------------------------------------
          ! Compute the Linearization Error 
@@ -1474 +1473 @@
    !!======================================================================
 #endif
+#endif
 END MODULE dynhpg_tam