Changeset 78 for trunk/NEMO/OPA_SRC/DYN
- Timestamp:
- 2004-04-22T14:54:25+02:00 (20 years ago)
- Location:
- trunk/NEMO/OPA_SRC/DYN
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/divcur.F90
r3 r78 123 123 ! open boundaries (div must be zero behind the open boundary) 124 124 ! mpp remark: The zeroing of hdivn can probably be extended to 1->jpi/jpj for the correct row/column 125 IF( lp eastobc) hdivn(nie0p1:nie1p1,nje0 :nje1 ,jk) = 0.e0 ! east126 IF( lp westobc) hdivn(niw0 :niw1 ,njw0 :njw1 ,jk) = 0.e0 ! west127 IF( lp northobc) hdivn(nin0 :nin1 ,njn0p1:njn1p1,jk) = 0.e0 ! north128 IF( lp southobc) hdivn(nis0 :nis1 ,njs0 :njs1 ,jk) = 0.e0 ! south125 IF( lp_obc_east ) hdivn(nie0p1:nie1p1,nje0 :nje1 ,jk) = 0.e0 ! east 126 IF( lp_obc_west ) hdivn(niw0 :niw1 ,njw0 :njw1 ,jk) = 0.e0 ! west 127 IF( lp_obc_north ) hdivn(nin0 :nin1 ,njn0p1:njn1p1,jk) = 0.e0 ! north 128 IF( lp_obc_south ) hdivn(nis0 :nis1 ,njs0 :njs1 ,jk) = 0.e0 ! south 129 129 #endif 130 130 … … 319 319 ! open boundaries (div must be zero behind the open boundary) 320 320 ! mpp remark: The zeroing of hdivn can probably be extended to 1->jpi/jpj for the correct row/column 321 IF( lp eastobc) hdivn(nie0p1:nie1p1,nje0 :nje1 ,jk) = 0.e0 ! east322 IF( lp westobc) hdivn(niw0 :niw1 ,njw0 :njw1 ,jk) = 0.e0 ! west323 IF( lp northobc) hdivn(nin0 :nin1 ,njn0p1:njn1p1,jk) = 0.e0 ! north324 IF( lp southobc) hdivn(nis0 :nis1 ,njs0 :njs1 ,jk) = 0.e0 ! south321 IF( lp_obc_east ) hdivn(nie0p1:nie1p1,nje0 :nje1 ,jk) = 0.e0 ! east 322 IF( lp_obc_west ) hdivn(niw0 :niw1 ,njw0 :njw1 ,jk) = 0.e0 ! west 323 IF( lp_obc_north ) hdivn(nin0 :nin1 ,njn0p1:njn1p1,jk) = 0.e0 ! north 324 IF( lp_obc_south ) hdivn(nis0 :nis1 ,njs0 :njs1 ,jk) = 0.e0 ! south 325 325 #endif 326 326 ! ! -------- -
trunk/NEMO/OPA_SRC/DYN/dynspg_rl.F90
r31 r78 69 69 !! gradient or a successive-over-relaxation method (depending 70 70 !! on nsolv, a namelist parameter). 71 !! -3- add to bsfd the island trends if l _isl=T71 !! -3- add to bsfd the island trends if lk_isl=T 72 72 !! -4- compute the after streamfunction is for further diagnos- 73 73 !! tics using a leap-frog scheme. … … 247 247 ! ----------------------------------- 248 248 249 IF( l _isl ) CALL isl_dyn_spg ! update bsfd249 IF( lk_isl ) CALL isl_dyn_spg ! update bsfd 250 250 251 251 … … 253 253 ! Compute bsf trend for OBC points (not masked) 254 254 255 IF( lp eastobc) THEN255 IF( lp_obc_east ) THEN 256 256 ! compute bsf trend at the boundary from bsfeob, computed in obc_spg 257 257 IF( neuler == 0 .AND. kt == nit000 ) THEN … … 268 268 ENDIF 269 269 270 IF( lp westobc) THEN270 IF( lp_obc_west ) THEN 271 271 ! compute bsf trend at the boundary from bsfwob, computed in obc_spg 272 272 IF( neuler == 0 .AND. kt == nit000 ) THEN … … 283 283 ENDIF 284 284 285 IF( lp northobc) THEN285 IF( lp_obc_north ) THEN 286 286 ! compute bsf trend at the boundary from bsfnob, computed in obc_spg 287 287 IF( neuler == 0 .AND. kt == nit000 ) THEN … … 298 298 ENDIF 299 299 300 IF( lp southobc) THEN300 IF( lp_obc_south ) THEN 301 301 ! compute bsf trend at the boundary from bsfsob, computed in obc_spg 302 302 IF( neuler == 0 .AND. kt == nit000 ) THEN … … 366 366 # if defined key_obc 367 367 ! Swap of boundary arrays 368 IF( lp eastobc) THEN368 IF( lp_obc_east ) THEN 369 369 ! (jped,jpef),nieob 370 370 IF( kt < nit000+3 .AND. .NOT.ln_rstart ) THEN … … 400 400 ENDIF 401 401 402 IF( lp westobc) THEN402 IF( lp_obc_west ) THEN 403 403 ! (jpwd,jpwf),niwob 404 404 IF( kt < nit000+3 .AND. .NOT.ln_rstart ) THEN … … 434 434 ENDIF 435 435 436 IF( lp northobc) THEN436 IF( lp_obc_north ) THEN 437 437 ! njnob,(jpnd,jpnf) 438 438 IF( kt < nit000 + 3 .AND. .NOT.ln_rstart ) THEN … … 468 468 ENDIF 469 469 470 IF( lp southobc) THEN470 IF( lp_obc_south ) THEN 471 471 ! njsob,(jpsd,jpsf) 472 472 IF( kt < nit000+3 .AND. .NOT.ln_rstart ) THEN
Note: See TracChangeset
for help on using the changeset viewer.