- Timestamp:
- 2020-01-27T15:31:53+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/DOM/domzgr.F90
r12150 r12340 45 45 !! * Substitutions 46 46 # include "vectopt_loop_substitute.h90" 47 # include "do_loop_substitute.h90" 47 48 !!---------------------------------------------------------------------- 48 49 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 150 151 ! 151 152 ! ! ice shelf draft and bathymetry 152 DO jj = 1,jpj 153 DO ji = 1,jpi 154 ikt = mikt(ji,jj) 155 ikb = mbkt(ji,jj) 156 bathy (ji,jj) = gdepw_0(ji,jj,ikb+1) 157 risfdep(ji,jj) = gdepw_0(ji,jj,ikt ) 158 END DO 159 END DO 153 DO_2D_11_11 154 ikt = mikt(ji,jj) 155 ikb = mbkt(ji,jj) 156 bathy (ji,jj) = gdepw_0(ji,jj,ikb+1) 157 risfdep(ji,jj) = gdepw_0(ji,jj,ikt ) 158 END_2D 160 159 ! 161 160 ! ! deepest/shallowest W level Above/Below ~10m … … 315 314 ! ! N.B. top k-index of W-level = mikt 316 315 ! ! bottom k-index of W-level = mbkt+1 317 DO jj = 1, jpjm1 318 DO ji = 1, jpim1 319 miku(ji,jj) = MAX( mikt(ji+1,jj ) , mikt(ji,jj) ) 320 mikv(ji,jj) = MAX( mikt(ji ,jj+1) , mikt(ji,jj) ) 321 mikf(ji,jj) = MAX( mikt(ji ,jj+1) , mikt(ji,jj), mikt(ji+1,jj ), mikt(ji+1,jj+1) ) 322 ! 323 mbku(ji,jj) = MIN( mbkt(ji+1,jj ) , mbkt(ji,jj) ) 324 mbkv(ji,jj) = MIN( mbkt(ji ,jj+1) , mbkt(ji,jj) ) 325 END DO 326 END DO 316 DO_2D_10_10 317 miku(ji,jj) = MAX( mikt(ji+1,jj ) , mikt(ji,jj) ) 318 mikv(ji,jj) = MAX( mikt(ji ,jj+1) , mikt(ji,jj) ) 319 mikf(ji,jj) = MAX( mikt(ji ,jj+1) , mikt(ji,jj), mikt(ji+1,jj ), mikt(ji+1,jj+1) ) 320 ! 321 mbku(ji,jj) = MIN( mbkt(ji+1,jj ) , mbkt(ji,jj) ) 322 mbkv(ji,jj) = MIN( mbkt(ji ,jj+1) , mbkt(ji,jj) ) 323 END_2D 327 324 ! converte into REAL to use lbc_lnk ; impose a min value of 1 as a zero can be set in lbclnk 328 325 zk(:,:) = REAL( miku(:,:), wp ) ; CALL lbc_lnk( 'domzgr', zk, 'U', 1. ) ; miku(:,:) = MAX( NINT( zk(:,:) ), 1 )
Note: See TracChangeset
for help on using the changeset viewer.