id summary reporter owner description type status priority milestone component version severity resolution keywords cc 1854 Use of jpdom_data in mppini_2.h90 is invalid? jamesharle jamesharle "= Context = When running nemo with jpni x jpnj < jpnij there are iom_get calls in mppini_2.h90 that use jpdom_data: i.e. when reading ztop and zbot (defined as array jpiglo x jpjglo) = Analysis = As far as I can see jpdom_data uses istart(1:2) = (/ mig(1), mjg(1) + jstartrow - 1 /). This is a problem because: 1) mig and mjg are yet to be defined (this happens later in the call to dom_ini) and 2) this references istart to the global position of the local processor - and ztop/zbot require the global data not just local. This is also an issue if applying an additional correction if using the bdy_msk array. = Fix = As in v3.6_stable when bathy_meter.nc was read in - jpdom_unkown can be used. {{{ 69c69 < REAL(wp), DIMENSION(jpiglo,jpjglo) :: zbot, ztop ! global workspace --- > REAL(wp), DIMENSION(jpiglo,jpjglo) :: zbot, ztop, zbdy ! global workspace 96,97c96,97 < CALL iom_get( inum, jpdom_data, 'bottom_level' , zbot ) ! nb of ocean T-points < CALL iom_get( inum, jpdom_data, 'top_level' , ztop ) ! nb of ocean T-points (ISF) --- > CALL iom_get( inum, jpdom_unknown, 'bottom_level' , zbot ) ! nb of ocean T-points > CALL iom_get( inum, jpdom_unknown, 'top_level' , ztop ) ! nb of ocean T-points (ISF) 118c118 < CALL iom_get ( inum, jpdom_data, 'bdy_msk', bdytmask(:,:) ) --- > CALL iom_get ( inum, jpdom_unknown, 'bdy_msk', zbdy ) 120c120 < WHERE ( bdytmask(:,:) <= 0. ) imask = 0 --- > WHERE ( zbdy(:,:) <= 0. ) imask = 0 }}}" Bug closed low OCE trunk fixed MPP, BDY