Changeset 789 for trunk/NEMO/OPA_SRC/DYN/dynhpg.F90
- Timestamp:
- 2008-01-11T19:04:56+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynhpg.F90
r719 r789 18 18 !! dyn_hpg : update the momentum trend with the now horizontal 19 19 !! gradient of the hydrostatic pressure 20 !! default case : k-j-i loops (vector opt. available)21 20 !! hpg_ctl : initialisation and control of options 22 21 !! hpg_zco : z-coordinate scheme … … 30 29 USE oce ! ocean dynamics and tracers 31 30 USE dom_oce ! ocean space and time domain 32 USE dynhpg_jki !33 31 USE phycst ! physical constants 34 32 USE in_out_manager ! I/O manager … … 42 40 43 41 PUBLIC dyn_hpg ! routine called by step module 44 45 #if defined key_mpp_omp46 !!----------------------------------------------------------------------47 !! 'key_mpp_omp' : j-k-i loop (j-slab)48 !!----------------------------------------------------------------------49 LOGICAL, PUBLIC, PARAMETER :: lk_dynhpg_jki = .TRUE. !: OpenMP hpg flag50 LOGICAL, PUBLIC, PARAMETER :: lk_dynhpg = .FALSE. !: vector hpg flag51 #else52 !!----------------------------------------------------------------------53 !! default case : k-j-i loop (vector opt.)54 !!----------------------------------------------------------------------55 LOGICAL, PUBLIC, PARAMETER :: lk_dynhpg_jki = .FALSE. !: OpenMP hpg flag56 LOGICAL, PUBLIC, PARAMETER :: lk_dynhpg = .TRUE. !: vector hpg flag57 #endif58 42 59 43 !!* Namelist nam_dynhpg : Choice of horizontal pressure gradient computation … … 111 95 CASE ( 5 ) ; CALL hpg_djc ( kt ) ! s-coordinate (Density Jacobian with Cubic polynomial) 112 96 CASE ( 6 ) ; CALL hpg_rot ( kt ) ! s-coordinate (ROTated axes scheme) 113 CASE ( 10 ) ; CALL hpg_zco_jki( kt ) ! z-coordinate (k-j-i)114 CASE ( 11 ) ; CALL hpg_zps_jki( kt ) ! z-coordinate plus partial steps (interpolation) (k-j-i)115 CASE ( 12 ) ; CALL hpg_sco_jki( kt ) ! s-coordinate (standard jacobian formulation) (k-j-i)116 97 END SELECT 117 98 … … 186 167 IF ( ioptio /= 1 ) CALL ctl_stop( ' NO or several hydrostatic pressure gradient options used' ) 187 168 188 IF( lk_dynhpg_jki ) THEN189 nhpg = nhpg + 10190 IF(lwp) WRITE(numout,*)191 IF(lwp) WRITE(numout,*) ' Autotasking or OPENMP: use j-k-i loops (i.e. _jki routines)'192 ENDIF193 169 ! 194 170 END SUBROUTINE hpg_ctl
Note: See TracChangeset
for help on using the changeset viewer.