Changeset 4328 for branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/FLO
- Timestamp:
- 2013-12-06T11:25:13+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/FLO/floblk.F90
r3294 r4328 14 14 USE dom_oce ! ocean space and time domain 15 15 USE phycst ! physical constants 16 USE obc_par ! open boundary condition parameters17 16 USE in_out_manager ! I/O manager 18 17 USE lib_mpp ! distribued memory computing library … … 345 344 IF( lk_mpp ) CALL mpp_sum( ijl , jpnfl ) 346 345 347 ! in the case of open boundaries we need to test if the floats don't348 ! go out of the domain. If it goes out, the float is put at the349 ! middle of the mesh in the domain but the trajectory isn't compute350 ! more time.351 # if defined key_obc352 DO jfl = 1, jpnfl353 IF( lp_obc_east ) THEN354 IF( jped <= zgjfl(jfl) .AND. zgjfl(jfl) <= jpef .AND. nieob-1 <= zgifl(jfl) ) THEN355 zgifl (jfl) = INT(zgifl(jfl)) + 0.5356 zgjfl (jfl) = INT(zgjfl(jfl)) + 0.5357 zagefl(jfl) = rdt358 END IF359 END IF360 IF( lp_obc_west ) THEN361 IF( jpwd <= zgjfl(jfl) .AND. zgjfl(jfl) <= jpwf .AND. niwob >= zgifl(jfl) ) THEN362 zgifl (jfl) = INT(zgifl(jfl)) + 0.5363 zgjfl (jfl) = INT(zgjfl(jfl)) + 0.5364 zagefl(jfl) = rdt365 END IF366 END IF367 IF( lp_obc_north ) THEN368 IF( jpnd <= zgifl(jfl) .AND. zgifl(jfl) <= jpnf .AND. njnob-1 >= zgjfl(jfl) ) THEN369 zgifl (jfl) = INT(zgifl(jfl)) + 0.5370 zgjfl (jfl) = INT(zgjfl(jfl)) + 0.5371 zagefl(jfl) = rdt372 END IF373 END IF374 IF( lp_obc_south ) THEN375 IF( jpsd <= zgifl(jfl) .AND. zgifl(jfl) <= jpsf .AND. njsob >= zgjfl(jfl) ) THEN376 zgifl (jfl) = INT(zgifl(jfl)) + 0.5377 zgjfl (jfl) = INT(zgjfl(jfl)) + 0.5378 zagefl(jfl) = rdt379 END IF380 END IF381 END DO382 #endif383 384 346 ! Test to know if a float hasn't integrated enought time 385 347 IF( ln_argo ) THEN
Note: See TracChangeset
for help on using the changeset viewer.