- Timestamp:
- 2019-05-16T15:23:56+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/LDF/ldfdyn.F90
r10968 r10986 151 151 WRITE(numout,*) ' lower limit (default 1.0) rn_minfac = ', rn_minfac 152 152 WRITE(numout,*) ' upper limit (default 1.0) rn_maxfac = ', rn_maxfac 153 IF(lflush) CALL FLUSH(numout) 153 154 ENDIF 154 155 … … 224 225 END SELECT 225 226 WRITE(numout,*) 227 IF(lflush) CALL FLUSH(numout) 226 228 ENDIF 227 229 … … 233 235 ! 234 236 IF( ln_dynldf_OFF ) THEN 235 IF(lwp) WRITE(numout,*) ' ==>>> No viscous operator selected. ahmt and ahmf are not allocated' 237 IF(lwp) THEN 238 WRITE(numout,*) ' ==>>> No viscous operator selected. ahmt and ahmf are not allocated' 239 IF(lflush) CALL FLUSH(numout) 240 ENDIF 236 241 RETURN 237 242 ! … … 255 260 ! 256 261 CASE( 0 ) !== constant ==! 257 IF(lwp) WRITE(numout,*) ' ==>>> eddy viscosity. = constant = ', zah0, cl_Units 262 IF(lwp) THEN 263 WRITE(numout,*) ' ==>>> eddy viscosity. = constant = ', zah0, cl_Units 264 IF(lflush) CALL FLUSH(numout) 265 ENDIF 258 266 ahmt(:,:,1:jpkm1) = zah0 259 267 ahmf(:,:,1:jpkm1) = zah0 260 268 ! 261 269 CASE( 10 ) !== fixed profile ==! 262 IF(lwp) WRITE(numout,*) ' ==>>> eddy viscosity = F( depth )' 263 IF(lwp) WRITE(numout,*) ' surface viscous coef. = constant = ', zah0, cl_Units 270 IF(lwp) THEN 271 WRITE(numout,*) ' ==>>> eddy viscosity = F( depth )' 272 WRITE(numout,*) ' surface viscous coef. = constant = ', zah0, cl_Units 273 IF(lflush) CALL FLUSH(numout) 274 ENDIF 264 275 ahmt(:,:,1) = zah0 ! constant surface value 265 276 ahmf(:,:,1) = zah0 … … 267 278 ! 268 279 CASE ( -20 ) !== fixed horizontal shape read in file ==! 269 IF(lwp) WRITE(numout,*) ' ==>>> eddy viscosity = F(i,j) read in eddy_viscosity.nc file' 280 IF(lwp) THEN 281 WRITE(numout,*) ' ==>>> eddy viscosity = F(i,j) read in eddy_viscosity.nc file' 282 IF(lflush) CALL FLUSH(numout) 283 ENDIF 270 284 CALL iom_open( 'eddy_viscosity_2D.nc', inum ) 271 285 CALL iom_get ( inum, jpdom_data, 'ahmt_2d', ahmt(:,:,1) ) … … 278 292 ! 279 293 CASE( 20 ) !== fixed horizontal shape ==! 280 IF(lwp) WRITE(numout,*) ' ==>>> eddy viscosity = F( e1, e2 ) or F( e1^3, e2^3 ) (lap. or blp. case)' 281 IF(lwp) WRITE(numout,*) ' using a fixed viscous velocity = ', rn_Uv ,' m/s and Lv = Max(e1,e2)' 282 IF(lwp) WRITE(numout,*) ' maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, ' for e1=1°)' 294 IF(lwp) THEN 295 WRITE(numout,*) ' ==>>> eddy viscosity = F( e1, e2 ) or F( e1^3, e2^3 ) (lap. or blp. case)' 296 WRITE(numout,*) ' using a fixed viscous velocity = ', rn_Uv ,' m/s and Lv = Max(e1,e2)' 297 WRITE(numout,*) ' maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, ' for e1=1°)' 298 IF(lflush) CALL FLUSH(numout) 299 ENDIF 283 300 CALL ldf_c2d( 'DYN', zUfac , inn , ahmt, ahmf ) ! surface value proportional to scale factor^inn 284 301 ! 285 302 CASE( -30 ) !== fixed 3D shape read in file ==! 286 IF(lwp) WRITE(numout,*) ' ==>>> eddy viscosity = F(i,j,k) read in eddy_viscosity_3D.nc file' 303 IF(lwp) THEN 304 WRITE(numout,*) ' ==>>> eddy viscosity = F(i,j,k) read in eddy_viscosity_3D.nc file' 305 IF(lflush) CALL FLUSH(numout) 306 ENDIF 287 307 CALL iom_open( 'eddy_viscosity_3D.nc', inum ) 288 308 CALL iom_get ( inum, jpdom_data, 'ahmt_3d', ahmt ) … … 291 311 ! 292 312 CASE( 30 ) !== fixed 3D shape ==! 293 IF(lwp) WRITE(numout,*) ' ==>>> eddy viscosity = F( latitude, longitude, depth )' 294 IF(lwp) WRITE(numout,*) ' using a fixed viscous velocity = ', rn_Uv ,' m/s and Ld = Max(e1,e2)' 295 IF(lwp) WRITE(numout,*) ' maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, ' for e1=1°)' 313 IF(lwp) THEN 314 WRITE(numout,*) ' ==>>> eddy viscosity = F( latitude, longitude, depth )' 315 WRITE(numout,*) ' using a fixed viscous velocity = ', rn_Uv ,' m/s and Ld = Max(e1,e2)' 316 WRITE(numout,*) ' maximum reachable coefficient (at the Equator) = ', zah_max, cl_Units, ' for e1=1°)' 317 IF(lflush) CALL FLUSH(numout) 318 ENDIF 296 319 CALL ldf_c2d( 'DYN', zUfac , inn , ahmt, ahmf ) ! surface value proportional to scale factor^inn 297 320 CALL ldf_c1d( 'DYN', ahmt(:,:,1), ahmf(:,:,1), ahmt, ahmf ) ! reduction with depth 298 321 ! 299 322 CASE( 31 ) !== time varying 3D field ==! 300 IF(lwp) WRITE(numout,*) ' ==>>> eddy viscosity = F( latitude, longitude, depth , time )' 301 IF(lwp) WRITE(numout,*) ' proportional to the local velocity : 1/2 |u|e (lap) or 1/12 |u|e^3 (blp)' 323 IF(lwp) THEN 324 WRITE(numout,*) ' ==>>> eddy viscosity = F( latitude, longitude, depth , time )' 325 WRITE(numout,*) ' proportional to the local velocity : 1/2 |u|e (lap) or 1/12 |u|e^3 (blp)' 326 IF(lflush) CALL FLUSH(numout) 327 ENDIF 302 328 ! 303 329 l_ldfdyn_time = .TRUE. ! will be calculated by call to ldf_dyn routine in step.F90 304 330 ! 305 331 CASE( 32 ) !== time varying 3D field ==! 306 IF(lwp) WRITE(numout,*) ' ==>>> eddy viscosity = F( latitude, longitude, depth , time )' 307 IF(lwp) WRITE(numout,*) ' proportional to the local deformation rate and gridscale (Smagorinsky)' 332 IF(lwp) THEN 333 WRITE(numout,*) ' ==>>> eddy viscosity = F( latitude, longitude, depth , time )' 334 WRITE(numout,*) ' proportional to the local deformation rate and gridscale (Smagorinsky)' 335 IF(lflush) CALL FLUSH(numout) 336 ENDIF 308 337 ! 309 338 l_ldfdyn_time = .TRUE. ! will be calculated by call to ldf_dyn routine in step.F90
Note: See TracChangeset
for help on using the changeset viewer.
