- Timestamp:
- 2010-11-15T22:20:05+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90
r2380 r2392 5 5 !!====================================================================== 6 6 !! History : OPA ! 1987-07 (G. Madec) Original code 7 !! -! 1993-03 (M. Guyon) symetrical conditions (M. Guyon)8 !! -! 1996-01 (G. Madec) suppression of common work arrays7 !! 6.0 ! 1993-03 (M. Guyon) symetrical conditions (M. Guyon) 8 !! 7.0 ! 1996-01 (G. Madec) suppression of common work arrays 9 9 !! - ! 1996-05 (G. Madec) mask computed from tmask and sup- 10 10 !! ! pression of the double computation of bmask 11 !! -! 1997-02 (G. Madec) mesh information put in domhgr.F12 !! -! 1997-07 (G. Madec) modification of mbathy and fmask11 !! 8.0 ! 1997-02 (G. Madec) mesh information put in domhgr.F 12 !! 8.1 ! 1997-07 (G. Madec) modification of mbathy and fmask 13 13 !! - ! 1998-05 (G. Roullet) free surface 14 !! -! 2000-03 (G. Madec) no slip accurate14 !! 8.2 ! 2000-03 (G. Madec) no slip accurate 15 15 !! - ! 2001-09 (J.-M. Molines) Open boundaries 16 16 !! NEMO 1.0 ! 2002-08 (G. Madec) F90: Free form and module … … 44 44 !! NEMO/OPA 3.2 , LODYC-IPSL (2009) 45 45 !! $Id$ 46 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 47 !!---------------------------------------------------------------------- 48 46 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 47 !!---------------------------------------------------------------------- 49 48 CONTAINS 50 49 … … 132 131 ENDIF 133 132 134 IF ( rn_shlat == 0. ) THEN ; IF(lwp) WRITE(numout,*) ' ocean lateral free-slip '133 IF ( rn_shlat == 0. ) THEN ; IF(lwp) WRITE(numout,*) ' ocean lateral free-slip ' 135 134 ELSEIF ( rn_shlat == 2. ) THEN ; IF(lwp) WRITE(numout,*) ' ocean lateral no-slip ' 136 135 ELSEIF ( 0. < rn_shlat .AND. rn_shlat < 2. ) THEN ; IF(lwp) WRITE(numout,*) ' ocean lateral partial-slip ' … … 308 307 IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA_R2 configuration 309 308 ! ! Increased lateral friction near of some straits 310 IF( n _cla == 0 ) THEN309 IF( nn_cla == 0 ) THEN 311 310 ! ! Gibraltar strait : partial slip (fmask=0.5) 312 311 ij0 = 101 ; ij1 = 101 … … 322 321 ! 323 322 ENDIF 324 325 323 ! ! Danish straits : strong slip (fmask > 2) 326 324 ! We keep this as an example but it is instable in this case … … 331 329 ! 332 330 ENDIF 333 ! ! ===================== 334 IF( cp_cfg == "orca" .AND. jp_cfg .eq. 1 ) THEN ! ORCA R1 configuration 335 ! ! ===================== 336 337 ii0 = 283 ; ii1 = 284 ! Gibraltar Strait 338 ij0 = 200 ; ij1 = 200 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 339 IF(lwp) WRITE(numout,*) 340 IF(lwp) WRITE(numout,*) ' orca_r1: Gibraltar : ' 341 342 ii0 = 314 ; ii1 = 315 ! Bhosporus Strait 343 ij0 = 208 ; ij1 = 208 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 344 IF(lwp) WRITE(numout,*) 345 IF(lwp) WRITE(numout,*) ' orca_r1: Bhosporus : ' 346 347 ii0 = 48 ; ii1 = 48 ! Makassar Strait (Top) 348 ij0 = 149 ; ij1 = 150 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3.0 349 IF(lwp) WRITE(numout,*) 350 IF(lwp) WRITE(numout,*) ' orca_r1: Makassar (Top) : ' 351 352 ii0 = 44 ; ii1 = 44 ! Lombok Strait 353 ij0 = 124 ; ij1 = 125 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 354 IF(lwp) WRITE(numout,*) 355 IF(lwp) WRITE(numout,*) ' orca_r1: Lombok : ' 356 357 ii0 = 53 ; ii1 = 53 ! Ombai Strait 358 ij0 = 124 ; ij1 = 125 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 359 IF(lwp) WRITE(numout,*) 360 IF(lwp) WRITE(numout,*) ' orca_r1: Ombai : ' 361 362 ii0 = 56 ; ii1 = 56 ! Timor Passage 363 ij0 = 124 ; ij1 = 125 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 364 IF(lwp) WRITE(numout,*) 365 IF(lwp) WRITE(numout,*) ' orca_r1: ' 366 367 ii0 = 58 ; ii1 = 58 ! West Halmahera Strait 368 ij0 = 141 ; ij1 = 142 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3.0 369 IF(lwp) WRITE(numout,*) 370 IF(lwp) WRITE(numout,*) ' orca_r1: West Halmahera : ' 371 372 ii0 = 55 ; ii1 = 55 ! East Halmahera Strait 373 ij0 = 141 ; ij1 = 142 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3.0 374 IF(lwp) WRITE(numout,*) 375 IF(lwp) WRITE(numout,*) ' orca_r1: East Halmahera : ' 376 377 ! 378 ! 379 ENDIF 331 ! 332 IF( cp_cfg == "orca" .AND. jp_cfg == 1 ) THEN ! ORCA R1 configuration 333 ! ! Increased lateral friction near of some straits 334 IF(lwp) WRITE(numout,*) 335 IF(lwp) WRITE(numout,*) ' orca_r1: increase friction near the following straits : ' 336 IF(lwp) WRITE(numout,*) ' Gibraltar ' 337 ii0 = 283 ; ii1 = 284 ! Gibraltar Strait 338 ij0 = 200 ; ij1 = 200 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 339 340 IF(lwp) WRITE(numout,*) ' Bhosporus ' 341 ii0 = 314 ; ii1 = 315 ! Bhosporus Strait 342 ij0 = 208 ; ij1 = 208 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 343 344 IF(lwp) WRITE(numout,*) ' Makassar (Top) ' 345 ii0 = 48 ; ii1 = 48 ! Makassar Strait (Top) 346 ij0 = 149 ; ij1 = 150 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3.0 347 348 IF(lwp) WRITE(numout,*) ' Lombok ' 349 ii0 = 44 ; ii1 = 44 ! Lombok Strait 350 ij0 = 124 ; ij1 = 125 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 351 352 IF(lwp) WRITE(numout,*) ' Ombai ' 353 ii0 = 53 ; ii1 = 53 ! Ombai Strait 354 ij0 = 124 ; ij1 = 125 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 355 356 IF(lwp) WRITE(numout,*) ' Timor Passage ' 357 ii0 = 56 ; ii1 = 56 ! Timor Passage 358 ij0 = 124 ; ij1 = 125 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2.0 359 360 IF(lwp) WRITE(numout,*) ' West Halmahera ' 361 ii0 = 58 ; ii1 = 58 ! West Halmahera Strait 362 ij0 = 141 ; ij1 = 142 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3.0 363 364 IF(lwp) WRITE(numout,*) ' East Halmahera ' 365 ii0 = 55 ; ii1 = 55 ! East Halmahera Strait 366 ij0 = 141 ; ij1 = 142 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3.0 367 ! 368 ENDIF 380 369 ! 381 370 CALL lbc_lnk( fmask, 'F', 1. ) ! Lateral boundary conditions on fmask
Note: See TracChangeset
for help on using the changeset viewer.
