Changeset 789 for trunk/NEMO/OPA_SRC/DYN/dynzdf.F90
- Timestamp:
- 2008-01-11T19:04:56+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynzdf.F90
r719 r789 17 17 USE dynzdf_exp ! vertical diffusion: explicit (dyn_zdf_exp routine) 18 18 USE dynzdf_imp ! vertical diffusion: implicit (dyn_zdf_imp routine) 19 USE dynzdf_imp_jki ! vertical diffusion implicit (dyn_zdf_imp_jki routine)20 19 21 20 USE ldfdyn_oce ! ocean dynamics: lateral physics … … 74 73 ! 75 74 CASE ( 0 ) ; CALL dyn_zdf_exp ( kt, r2dt ) ! explicit scheme 76 CASE ( 1 ) ; CALL dyn_zdf_imp ( kt, r2dt ) ! implicit scheme (k-j-i loop) 77 CASE ( 2 ) ; CALL dyn_zdf_imp_jki( kt, r2dt ) ! implicit scheme (j-k-i loop) 75 CASE ( 1 ) ; CALL dyn_zdf_imp ( kt, r2dt ) ! implicit scheme 78 76 ! 79 77 CASE ( -1 ) ! esopa: test all possibility with control print … … 83 81 CALL dyn_zdf_imp ( kt, r2dt ) 84 82 CALL prt_ctl( tab3d_1=ua, clinfo1=' zdf1 - Ua: ', mask1=umask, & 85 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' )86 CALL dyn_zdf_imp_jki( kt, r2dt )87 CALL prt_ctl( tab3d_1=ua, clinfo1=' zdf2 - Ua: ', mask1=umask, &88 83 & tab3d_2=va, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' ) 89 84 END SELECT … … 109 104 !! ** Method : implicit (euler backward) scheme (default) 110 105 !! explicit (time-splitting) scheme if ln_zdfexp=T 111 !! OpenMP / NEC autotasking: use j-k-i loops112 106 !!---------------------------------------------------------------------- 113 107 USE zdftke … … 125 119 IF( ln_dynldf_hor .AND. ln_sco ) nzdf = 1 ! horizontal lateral physics in s-coordinate 126 120 127 ! OpenMP / NEC autotasking128 #if defined key_mpp_omp129 IF( nzdf == 1 ) nzdf = 2 ! j-k-i loop130 #endif131 132 121 IF( lk_esopa ) nzdf = -1 ! Esopa key: All schemes used 133 122 … … 139 128 IF( nzdf == 0 ) WRITE(numout,*) ' Explicit time-splitting scheme' 140 129 IF( nzdf == 1 ) WRITE(numout,*) ' Implicit (euler backward) scheme' 141 IF( nzdf == 2 ) WRITE(numout,*) ' Implicit (euler backward) scheme with j-k-i loops'142 130 ENDIF 143 131 !
Note: See TracChangeset
for help on using the changeset viewer.