New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/TRA/traldf_tam.F90 – NEMO

Ignore:
Timestamp:
2011-02-15T12:58:59+01:00 (13 years ago)
Author:
vidard
Message:

refer to ticket #798

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/TRA/traldf_tam.F90

    r1885 r2587  
    1111   !!          9.0  ! 08-06 (A. Vidard) Skeleton 
    1212   !!          9.0  ! 09-03 (F. Vigilant) adding tra_ldf_lap option 
     13   !!          9.0  ! 10-06 (P.A. Bouttier) adding tra_ldf_bilap option 
    1314   !!---------------------------------------------------------------------- 
    1415 
     
    2526      & tra_ldf_lap_tan,     & 
    2627      & tra_ldf_lap_adj,     & 
    27       & tra_ldf_lap_adj_tst, & 
    28       & tra_ldf_lap_tlm_tst 
     28#if defined key_tst_tlm 
     29      & tra_ldf_lap_tlm_tst, & 
     30#endif 
     31      & tra_ldf_lap_adj_tst 
     32   USE traldf_bilap_tam, ONLY: & !lateral mixing                (tra_ldf_bilap routine) 
     33      & tra_ldf_bilap_tan,   & 
     34      & tra_ldf_bilap_adj 
    2935   USE in_out_manager, ONLY: & ! I/O manager 
    3036      & ctl_stop, nit000, lwp, numout, nitend 
     
    5056   PUBLIC   tra_ldf_adj     ! called by step_tam.F90  
    5157   PUBLIC   tra_ldf_adj_tst ! called by tamtst.F90  
     58#if defined key_tst_tlm 
    5259   PUBLIC   tra_ldf_tlm_tst ! called by tamtst.F90  
     60#endif 
     61   PUBLIC   ldf_ctl_tam     ! called by trazdf_imp (init of l_traldf_rot) 
    5362 
    5463   INTEGER ::  nldf 
     
    7685      CASE ( 0 )   ;   CALL tra_ldf_lap_tan   ( kt )      ! iso-level laplacian 
    7786      CASE ( 1 )   ;   CALL tra_ldf_iso_tan   ( kt )      ! rotated laplacian (except dk[ dk[.] ] part) 
     87      CASE ( 2 )   ;   CALL tra_ldf_bilap_tan ( kt )      ! iso-level bilaplacian 
    7888      END SELECT 
    7989   END SUBROUTINE tra_ldf_tan 
     
    94104      CASE ( 0 )   ;   CALL tra_ldf_lap_adj   ( kt )      ! rotated laplacian (except dk[ dk[.] ] part) 
    95105      CASE ( 1 )   ;   CALL tra_ldf_iso_adj   ( kt )      ! rotated laplacian (except dk[ dk[.] ] part) 
     106      CASE ( 2 )   ;   CALL tra_ldf_bilap_adj ( kt )      ! iso-level bilaplacian 
    96107      END SELECT 
    97108      ! 
     
    206217 
    207218      IF( ln_traldf_bilap ) THEN      ! bilaplacian operator 
    208            CALL ctl_stop( '          You shouldn t have seen this error message, ln_trad_bilap option not impemented yet for tam' ) 
     219         IF ( ln_zco ) THEN                ! z-coordinate 
     220            IF ( ln_traldf_level )   nldf = 2      ! iso-level  (no rotation) 
     221            IF ( ln_traldf_hor   )   nldf = 2      ! horizontal (no rotation) 
     222            IF ( ln_traldf_iso   )   ierr = 2      ! isoneutral (   rotation) 
     223      ENDIF 
     224         IF ( ln_zps ) THEN             ! z-coordinate 
     225            IF ( ln_traldf_level )   ierr = 1      ! iso-level not allowed  
     226            IF ( ln_traldf_hor   )   nldf = 2      ! horizontal (no rotation) 
     227            IF ( ln_traldf_iso   )   ierr = 2      ! isoneutral (   rotation) 
     228         ENDIF 
    209229      ENDIF 
    210230 
    211231      IF( ierr == 1 )   CALL ctl_stop( ' iso-level in z-coordinate - partial step, not allowed' ) 
     232      IF( ierr == 2 )   CALL ctl_stop( ' isoneutral bilaplacian operator does not exist' ) 
    212233      IF( lk_traldf_eiv .AND. .NOT.ln_traldf_iso )   & 
    213234           CALL ctl_stop( '          eddy induced velocity on tracers',   & 
     
    227248   END SUBROUTINE ldf_ctl_tam 
    228249 
    229  
     250#if defined key_tst_tlm 
    230251   SUBROUTINE tra_ldf_tlm_tst( kumadt ) 
    231252      !!----------------------------------------------------------------------- 
     
    266287   !!====================================================================== 
    267288#endif 
     289#endif 
    268290END MODULE traldf_tam 
Note: See TracChangeset for help on using the changeset viewer.