- Timestamp:
- 2020-06-02T17:46:26+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/SWE/dynldf_lap_blp.F90
r12822 r13005 20 20 USE in_out_manager ! I/O manager 21 21 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 22 22 USE lib_mpp ! MPP library 23 ! 24 USE usrdef_nam , ONLY : nn_dynldf_lap_typ ! use laplacian parameter 25 ! 23 26 IMPLICIT NONE 24 27 PRIVATE … … 31 34 INTEGER, PUBLIC, PARAMETER :: np_dynldf_lap_symN = 3 ! symmetric laplacian (cartesian) 32 35 33 INTEGER, PUBLIC, PARAMETER :: ln_dynldf_lap_typ = 1 ! choose type of laplacian (ideally from namelist)36 !INTEGER, PUBLIC, PARAMETER :: nn_dynldf_lap_typ = 1 ! choose type of laplacian (ideally from namelist) 34 37 !!anSYM 35 38 !! * Substitutions … … 77 80 WRITE(numout,*) 'dyn_ldf : iso-level harmonic (laplacian) operator, pass=', kpass 78 81 WRITE(numout,*) '~~~~~~~ ' 79 WRITE(numout,*) ' ln_dynldf_lap_typ = ', ln_dynldf_lap_typ80 SELECT CASE( ln_dynldf_lap_typ ) ! print the choice of operator82 WRITE(numout,*) ' nn_dynldf_lap_typ = ', nn_dynldf_lap_typ 83 SELECT CASE( nn_dynldf_lap_typ ) ! print the choice of operator 81 84 CASE( np_dynldf_lap_rot ) ; WRITE(numout,*) ' ==>>> div-rot laplacian' 82 85 CASE( np_dynldf_lap_sym ) ; WRITE(numout,*) ' ==>>> symmetric laplacian (covariant form)' 83 CASE( np_dynldf_lap_symN) ; WRITE(numout,*) ' ==>>> symmetric laplacian ( simpleform)'86 CASE( np_dynldf_lap_symN) ; WRITE(numout,*) ' ==>>> symmetric laplacian (cartesian form)' 84 87 END SELECT 85 88 ENDIF … … 89 92 ENDIF 90 93 ! 91 SELECT CASE( ln_dynldf_lap_typ )94 SELECT CASE( nn_dynldf_lap_typ ) 92 95 ! 93 96 CASE ( np_dynldf_lap_rot ) !== Vorticity-Divergence form ==! … … 157 160 END DO ! End of slab 158 161 ! 159 CASE ( np_dynldf_lap_symN ) !== Symmetric form ==! ( naiveway)162 CASE ( np_dynldf_lap_symN ) !== Symmetric form ==! (cartesian way) 160 163 ! 161 164 DO jk = 1, jpkm1 ! Horizontal slab … … 190 193 ! 191 194 CASE DEFAULT ! error 192 CALL ctl_stop('STOP','dyn_ldf_lap: wrong value for ln_dynldf_lap_typ' )195 CALL ctl_stop('STOP','dyn_ldf_lap: wrong value for nn_dynldf_lap_typ' ) 193 196 END SELECT 194 197 !
Note: See TracChangeset
for help on using the changeset viewer.