- Timestamp:
- 2020-01-14T18:28:28+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/ZDF/zdfosm.F90
r12317 r12319 141 141 !!---------------------------------------------------------------------- 142 142 !! NEMO/OCE 4.0 , NEMO Consortium (2018) 143 !! $Id $143 !! $Id: zdfosm.F90 12317 2020-01-14 12:40:47Z agn $ 144 144 !! Software governed by the CeCILL license (see ./LICENSE) 145 145 !!---------------------------------------------------------------------- … … 483 483 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 484 484 ! BL must be always 4 levels deep. 485 ! For calculation of lateral buoyancy gradients for FK in 486 ! zdf_osm_zmld_horizontal_gradients need halo values for ibld, so must 487 ! previously exist for hbl also. 485 488 hbl(:,:) = MAX(hbl(:,:), gdepw_n(:,:,4) ) 486 489 ibld(:,:) = 4 487 490 DO jk = 5, jpkm1 488 DO jj = 2, jpjm1489 DO ji = 2, jpim1491 DO jj = 1, jpj 492 DO ji = 1, jpi 490 493 IF ( hbl(ji,jj) >= gdepw_n(ji,jj,jk) ) THEN 491 494 ibld(ji,jj) = MIN(mbkt(ji,jj), jk) … … 526 529 END DO 527 530 528 imld = ibld! use imld to hold previous blayer index531 imld(:,:) = ibld(:,:) ! use imld to hold previous blayer index 529 532 ibld(:,:) = 4 530 533 … … 1145 1148 END DO 1146 1149 END DO 1150 ! Lateral boundary conditions on final outputs for hbl, on T-grid (sign unchanged) 1151 CALL lbc_lnk_multi( 'zdfosm', hbl, 'T', 1., dh, 'T', 1., hmle, 'T', 1. ) 1147 1152 ! Lateral boundary conditions on final outputs for gham[ts], on W-grid (sign unchanged) 1148 1153 ! Lateral boundary conditions on final outputs for gham[uv], on [UV]-grid (sign unchanged) … … 2192 2197 CALL iom_set_rstw_var_active('wn') 2193 2198 CALL iom_set_rstw_var_active('hbl') 2194 CALL iom_set_rstw_var_active('hbli') 2199 CALL iom_set_rstw_var_active('dh') 2200 IF( ln_osm_mle ) THEN 2201 CALL iom_set_rstw_var_active('hmle') 2202 END IF 2195 2203 ENDIF 2196 2204 END SUBROUTINE zdf_osm_init
Note: See TracChangeset
for help on using the changeset viewer.