- Timestamp:
- 2015-02-11T11:50:34+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO7_STARTHOUR/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90
r4811 r5075 111 111 END DO 112 112 ! ! Inverse of the local depth 113 hur(:,:) = 1._wp / ( hu(:,:) + 1._wp - umask (:,:,1) ) * umask(:,:,1)114 hvr(:,:) = 1._wp / ( hv(:,:) + 1._wp - vmask (:,:,1) ) * vmask(:,:,1)113 hur(:,:) = 1._wp / ( hu(:,:) + 1._wp - umask_i(:,:) ) * umask_i(:,:) 114 hvr(:,:) = 1._wp / ( hv(:,:) + 1._wp - vmask_i(:,:) ) * vmask_i(:,:) 115 115 116 116 CALL dom_stp ! time step … … 367 367 ! 368 368 IF(lk_mpp) THEN 369 CALL mpp_minloc( e1t(:,:), tmask (:,:,1), ze1min, iimi1,ijmi1 )370 CALL mpp_minloc( e2t(:,:), tmask (:,:,1), ze2min, iimi2,ijmi2 )371 CALL mpp_maxloc( e1t(:,:), tmask (:,:,1), ze1max, iima1,ijma1 )372 CALL mpp_maxloc( e2t(:,:), tmask (:,:,1), ze2max, iima2,ijma2 )369 CALL mpp_minloc( e1t(:,:), tmask_i(:,:), ze1min, iimi1,ijmi1 ) 370 CALL mpp_minloc( e2t(:,:), tmask_i(:,:), ze2min, iimi2,ijmi2 ) 371 CALL mpp_maxloc( e1t(:,:), tmask_i(:,:), ze1max, iima1,ijma1 ) 372 CALL mpp_maxloc( e2t(:,:), tmask_i(:,:), ze2max, iima2,ijma2 ) 373 373 ELSE 374 ze1min = MINVAL( e1t(:,:), mask = tmask (:,:,1) == 1._wp )375 ze2min = MINVAL( e2t(:,:), mask = tmask (:,:,1) == 1._wp )376 ze1max = MAXVAL( e1t(:,:), mask = tmask (:,:,1) == 1._wp )377 ze2max = MAXVAL( e2t(:,:), mask = tmask (:,:,1) == 1._wp )378 379 iloc = MINLOC( e1t(:,:), mask = tmask (:,:,1) == 1._wp )374 ze1min = MINVAL( e1t(:,:), mask = tmask_i(:,:) == 1._wp ) 375 ze2min = MINVAL( e2t(:,:), mask = tmask_i(:,:) == 1._wp ) 376 ze1max = MAXVAL( e1t(:,:), mask = tmask_i(:,:) == 1._wp ) 377 ze2max = MAXVAL( e2t(:,:), mask = tmask_i(:,:) == 1._wp ) 378 379 iloc = MINLOC( e1t(:,:), mask = tmask_i(:,:) == 1._wp ) 380 380 iimi1 = iloc(1) + nimpp - 1 381 381 ijmi1 = iloc(2) + njmpp - 1 382 iloc = MINLOC( e2t(:,:), mask = tmask (:,:,1) == 1._wp )382 iloc = MINLOC( e2t(:,:), mask = tmask_i(:,:) == 1._wp ) 383 383 iimi2 = iloc(1) + nimpp - 1 384 384 ijmi2 = iloc(2) + njmpp - 1 385 iloc = MAXLOC( e1t(:,:), mask = tmask (:,:,1) == 1._wp )385 iloc = MAXLOC( e1t(:,:), mask = tmask_i(:,:) == 1._wp ) 386 386 iima1 = iloc(1) + nimpp - 1 387 387 ijma1 = iloc(2) + njmpp - 1 388 iloc = MAXLOC( e2t(:,:), mask = tmask (:,:,1) == 1._wp )388 iloc = MAXLOC( e2t(:,:), mask = tmask_i(:,:) == 1._wp ) 389 389 iima2 = iloc(1) + nimpp - 1 390 390 ijma2 = iloc(2) + njmpp - 1
Note: See TracChangeset
for help on using the changeset viewer.