Changeset 9746
- Timestamp:
- 2018-06-06T10:57:59+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_isf_remapping_UKESM_GO6package_r9314/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90
r9732 r9746 1263 1263 !! ** Purpose : check the bathymetry in levels 1264 1264 !! 1265 !! ** Method : T He water column have to contained at least 2 cells1265 !! ** Method : The water column have to contained at least 2 cells 1266 1266 !! Bathymetry and isfdraft are modified (dig/close) to respect 1267 1267 !! this criterion. 1268 !! 0.0 = definition of minal ice shelf draft 1268 1269 !! digging and filling base on depth criterion only 1269 1270 !! 1.0 = set iceshelf to the minimum depth allowed … … 1291 1292 ! 1292 1293 REAL(wp) :: zdepth, vmskjp1, vmskjm1, umskip1, umskim1, umskip1_r, umskim1_r, vmskjp1_r, vmskjm1_r 1294 REAL(wp) :: zisfdep_min 1293 1295 REAL(wp), POINTER, DIMENSION(:,:) :: zdummy, zmask, zrisfdep ! 2D workspace (ISH) 1294 1296 ! … … 1334 1336 END IF 1335 1337 ! 1338 ! 0.0 variable definition of minimal ice shelf draft allowed 1339 ! an ice shelf draft have to be larger than the first level thickness 1340 zisfdep_min = MAX(rn_isfdep_min,e3t_1d(1)) 1341 1336 1342 ! 1.0 set iceshelf to the minimum depth allowed 1337 WHERE(risfdep(:,:) > 0.0_wp .AND. risfdep(:,:) < MAX(rn_isfdep_min,e3t_1d(1)))1338 risfdep(:,:)= rn_isfdep_min1343 WHERE(risfdep(:,:) > 0.0_wp .AND. risfdep(:,:) < zisfdep_min) 1344 risfdep(:,:)=zisfdep_min 1339 1345 END WHERE 1340 1346 ! … … 1351 1357 DO ji = 1, jpi 1352 1358 IF (bathy(ji,jj)-risfdep(ji,jj) < rn_isfhw_min) THEN 1353 IF ((bathy(ji,jj) - rn_isfhw_min) > MAX(rn_isfdep_min,e3t_1d(1))) THEN1354 risfdep(ji,jj) = bathy(ji,jj) - rn_isfhw_min1355 ELSE1359 risfdep(ji,jj) = bathy(ji,jj) - rn_isfhw_min 1360 ! sanity check on risfdep (if < zisfdep_min) => we ground it 1361 IF (risfdep(ji,jj) < zisfdep_min) THEN 1356 1362 risfdep(ji,jj)=0.0 ; misfdep(ji,jj) = 1 1357 1363 bathy (ji,jj)=0.0 ; mbathy (ji,jj) = 0
Note: See TracChangeset
for help on using the changeset viewer.