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

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

@@ -74 +74 @@
       & ahtw,                &
       & aht0
+#if defined key_ldfslp
    USE ldfslp        , ONLY: & ! lateral physics: slope of diffusion
       & wslpi,               & !: i_slope at W-points
       & wslpj                  !: j-slope at W-points
+#endif 
 #if defined key_zdfddm
    USE zdfddm        , ONLY: &
       & avs
 #endif 
+   USE traldf_tam
    USE in_out_manager, ONLY: & ! I/O manager 
       & lwp,          &
@@ -106 +109 @@
    PUBLIC tra_zdf_imp_adj       !  routine called by tra_zdf_adj.F90
    PUBLIC tra_zdf_imp_adj_tst   !  routine called by tst.F90
+#if defined key_tst_tlm
    PUBLIC tra_zdf_imp_tlm_tst   !  routine called by tamtst.F90 
+#endif
 
    !! * Substitutions
@@ -483 +488 @@
       !!---------------------------------------------------------------------
 
-      IF( kt == nit000 ) THEN
+      IF( kt == nitend ) THEN
          IF(lwp)WRITE(numout,*)
          IF(lwp)WRITE(numout,*) 'tra_zdf_imp_adj : implicit vertical mixing'
          IF(lwp)WRITE(numout,*) '~~~~~~~~~~~~~~~ '
+         CALL ldf_ctl_tam  ! init of l_traldf_rot
          zavi = 0._wp      ! avoid warning at compilation phase when lk_ldfslp=F
       ENDIF
@@ -985 +991 @@
 
    END SUBROUTINE tra_zdf_imp_adj_tst
-
+#if defined key_tst_tlm
    SUBROUTINE tra_zdf_imp_tlm_tst( kumadt )
       !!-----------------------------------------------------------------------
@@ -1019 +1025 @@
       USE tamtrj              ! writing out state trajectory
       USE par_tlm,    ONLY: &
+        & tlm_bch,          &
         & cur_loop,         &
         & h_ratio
@@ -1077 +1084 @@
       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
@@ -1126 +1133 @@
       ! Output filename Xn=F(X0)
       !--------------------------------------------------------------------
-      file_wop='trj_wop_trazdf_imp'
       CALL tlm_namrd
       gamma = h_ratio
+      file_wop='trj_wop_trazdf_imp'
+      file_xdx='trj_xdx_trazdf_imp'
       !--------------------------------------------------------------------
       ! Initialize the tangent input with random noise: dx
@@ -1169 +1177 @@
       ! Complete Init for Direct
       !-------------------------------------------------------------------
-      CALL istate_p  
+      IF ( tlm_bch /= 2 )      CALL istate_p  
 
       ! *** initialize the reference trajectory
@@ -1192 +1200 @@
       !  Compute the direct model F(X0,t=n) = Xn
       !--------------------------------------------------------------------
-      CALL tra_zdf_imp(nit000, rdttra)
-
-      IF ( cur_loop .EQ. 0) CALL trj_wri_spl(file_wop)
-
+      IF ( tlm_bch /= 2 )  CALL tra_zdf_imp(nit000, rdttra)
+      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  
@@ -1235 +1237 @@
          zta_wop  (:,:,:) = ta  (:,:,:)
          zsa_wop  (:,:,:) = sa  (:,:,:)
-
+         CALL trj_rd_spl(file_xdx) 
+         zta_out  (:,:,:) = ta  (:,:,:)
+         zsa_out  (:,:,:) = sa  (:,:,:)
          !--------------------------------------------------------------------
          ! Compute the Linearization Error 
@@ -1378 +1382 @@
    END SUBROUTINE tra_zdf_imp_tlm_tst
 #endif
+#endif
 END MODULE trazdf_imp_tam