Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/TRA/traldf_tam.F90
- Timestamp:
- 2011-02-15T12:58:59+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/TRA/traldf_tam.F90
r1885 r2587 11 11 !! 9.0 ! 08-06 (A. Vidard) Skeleton 12 12 !! 9.0 ! 09-03 (F. Vigilant) adding tra_ldf_lap option 13 !! 9.0 ! 10-06 (P.A. Bouttier) adding tra_ldf_bilap option 13 14 !!---------------------------------------------------------------------- 14 15 … … 25 26 & tra_ldf_lap_tan, & 26 27 & 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 29 35 USE in_out_manager, ONLY: & ! I/O manager 30 36 & ctl_stop, nit000, lwp, numout, nitend … … 50 56 PUBLIC tra_ldf_adj ! called by step_tam.F90 51 57 PUBLIC tra_ldf_adj_tst ! called by tamtst.F90 58 #if defined key_tst_tlm 52 59 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) 53 62 54 63 INTEGER :: nldf … … 76 85 CASE ( 0 ) ; CALL tra_ldf_lap_tan ( kt ) ! iso-level laplacian 77 86 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 78 88 END SELECT 79 89 END SUBROUTINE tra_ldf_tan … … 94 104 CASE ( 0 ) ; CALL tra_ldf_lap_adj ( kt ) ! rotated laplacian (except dk[ dk[.] ] part) 95 105 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 96 107 END SELECT 97 108 ! … … 206 217 207 218 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 209 229 ENDIF 210 230 211 231 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' ) 212 233 IF( lk_traldf_eiv .AND. .NOT.ln_traldf_iso ) & 213 234 CALL ctl_stop( ' eddy induced velocity on tracers', & … … 227 248 END SUBROUTINE ldf_ctl_tam 228 249 229 250 #if defined key_tst_tlm 230 251 SUBROUTINE tra_ldf_tlm_tst( kumadt ) 231 252 !!----------------------------------------------------------------------- … … 266 287 !!====================================================================== 267 288 #endif 289 #endif 268 290 END MODULE traldf_tam
Note: See TracChangeset
for help on using the changeset viewer.