Changeset 2493
- Timestamp:
- 2010-12-20T14:47:12+01:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r2482 r2493 805 805 INTEGER :: ji, jj, jk ! dummy loop indices 806 806 INTEGER :: ik, it ! temporary integers 807 INTEGER :: nlbelow ! temporary integer808 807 LOGICAL :: ll_print ! Allow control print for debugging 809 808 REAL(wp) :: ze3tp , ze3wp ! Last ocean level thickness at T- and W-points … … 833 832 ! =================== 834 833 zmax = gdepw_0(jpk) + e3t_0(jpk) ! maximum depth (i.e. the last ocean level thickness <= 2*e3t_0(jpkm1) ) 835 mbathy(:,:) = jpkm1 ! initialize mbathy to the maximum ocean level available836 834 bathy(:,:) = MIN( zmax , bathy(:,:) ) ! bounded value of bathy (min already set at the end of zgr_bat) 835 WHERE( bathy(:,:) == 0._wp ) ; mbathy(:,:) = 0 ! land : set mbathy to 0 836 ELSE WHERE ; mbathy(:,:) = jpkm1 ! ocean : initialize mbathy to the max ocean level 837 END WHERE 837 838 838 839 ! Compute mbathy for ocean points (i.e. the number of ocean levels)
Note: See TracChangeset
for help on using the changeset viewer.