- Timestamp:
- 2011-10-14T17:18:45+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2802_NOCL_Smagorinsky/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90
r2887 r2926 38 38 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 39 39 !! $Id$ 40 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 41 !!---------------------------------------------------------------------- 42 40 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 41 !!---------------------------------------------------------------------- 43 42 CONTAINS 44 43 … … 63 62 !!---------------------------------------------------------------------- 64 63 INTEGER :: ioptio ! ??? 65 LOGICAL :: ll_print = .FALSE. ! Logical flag for printing viscosity coef.64 LOGICAL :: ll_print = .FALSE. ! Logical flag for printing viscosity coef. 66 65 !! 67 66 NAMELIST/namdyn_ldf/ ln_dynldf_lap , ln_dynldf_bilap, & 68 67 & ln_dynldf_level, ln_dynldf_hor , ln_dynldf_iso, & 69 & rn_ahm_0_lap , rn_ahmb_0 , rn_ahm_0_blp ,cmsmag1,cmsmag268 & rn_ahm_0_lap , rn_ahmb_0 , rn_ahm_0_blp 70 69 !!---------------------------------------------------------------------- 71 70 … … 86 85 WRITE(numout,*) ' background viscosity rn_ahmb_0 = ', rn_ahmb_0 87 86 WRITE(numout,*) ' horizontal bilaplacian eddy viscosity rn_ahm_0_blp = ', rn_ahm_0_blp 88 WRITE(numout,*) ' smagorinsky coefficient for laplacian cmsmag1 = ' , cmsmag189 WRITE(numout,*) ' smagorinsky coefficient for bilaplacian cmsmag2 = ' , cmsmag290 87 ENDIF 91 88 … … 131 128 ! Lateral eddy viscosity 132 129 ! ====================== 133 #if defined key_dynldf_smag && ! defined key_dynldf_c3d134 IF(lwp) WRITE(numout,*) ' key_dynldf_smag can not be used without key_dynldf_c3d'135 #endif136 137 130 #if defined key_dynldf_c3d 138 131 CALL ldf_dyn_c3d( ll_print ) ! ahm = 3D coef. = F( longitude, latitude, depth ) … … 213 206 REAL(wp), INTENT(in ) :: pwam ! width of inflection 214 207 REAL(wp), INTENT(in ) :: pbot ! bottom value (0<pbot<= 1) 215 REAL(wp), INTENT(in ), DIMENSION (jpk) :: pdep ! depth of the gridpoint (T, U, V, F)216 REAL(wp), INTENT(inout), DIMENSION (jpi,jpj,jpk) :: pah ! adimensional vertical profile208 REAL(wp), INTENT(in ), DIMENSION (:) :: pdep ! depth of the gridpoint (T, U, V, F) 209 REAL(wp), INTENT(inout), DIMENSION (:,:,:) :: pah ! adimensional vertical profile 217 210 !! 218 211 INTEGER :: jk ! dummy loop indices … … 255 248 REAL(wp), INTENT(in ) :: pwam ! width of inflection 256 249 REAL(wp), INTENT(in ) :: pbot ! bottom value (0<pbot<= 1) 257 REAL(wp), INTENT(in ), DIMENSION (jpi,jpj,jpk) :: pdep ! dep of the gridpoint (T, U, V, F)258 REAL(wp), INTENT(inout), DIMENSION (jpi,jpj,jpk) :: pah ! adimensional vertical profile250 REAL(wp), INTENT(in ), DIMENSION (:,:,:) :: pdep ! dep of the gridpoint (T, U, V, F) 251 REAL(wp), INTENT(inout), DIMENSION (:,:,:) :: pah ! adimensional vertical profile 259 252 !! 260 253 INTEGER :: jk ! dummy loop indices
Note: See TracChangeset
for help on using the changeset viewer.