New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12319 for NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/ZDF/zdfosm.F90 – NEMO

Ignore:
Timestamp:
2020-01-14T18:28:28+01:00 (4 years ago)
Author:
agn
Message:

FKOSM reproducible after wrapping hbl, hble, dh with ln_lnk_multi & calc ibld over halo points

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/ZDF/zdfosm.F90

    r12317 r12319  
    141141   !!---------------------------------------------------------------------- 
    142142   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    143    !! $Id$ 
     143   !! $Id: zdfosm.F90 12317 2020-01-14 12:40:47Z agn $ 
    144144   !! Software governed by the CeCILL license (see ./LICENSE) 
    145145   !!---------------------------------------------------------------------- 
     
    483483     !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    484484     ! 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. 
    485488      hbl(:,:) = MAX(hbl(:,:), gdepw_n(:,:,4) ) 
    486489      ibld(:,:) = 4 
    487490      DO jk = 5, jpkm1 
    488          DO jj = 2, jpjm1 
    489             DO ji = 2, jpim1 
     491         DO jj = 1, jpj 
     492            DO ji = 1, jpi 
    490493               IF ( hbl(ji,jj) >= gdepw_n(ji,jj,jk) ) THEN 
    491494                  ibld(ji,jj) = MIN(mbkt(ji,jj), jk) 
     
    526529         END DO 
    527530 
    528       imld = ibld           ! use imld to hold previous blayer index 
     531      imld(:,:) = ibld(:,:)           ! use imld to hold previous blayer index 
    529532      ibld(:,:) = 4 
    530533 
     
    11451148           END DO 
    11461149        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. ) 
    11471152        ! Lateral boundary conditions on final outputs for gham[ts],  on W-grid  (sign unchanged) 
    11481153        ! Lateral boundary conditions on final outputs for gham[uv],  on [UV]-grid  (sign unchanged) 
     
    21922197        CALL iom_set_rstw_var_active('wn') 
    21932198        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 
    21952203     ENDIF 
    21962204   END SUBROUTINE zdf_osm_init 
Note: See TracChangeset for help on using the changeset viewer.