Newer
Older
END DO
CALL iom_rstput( 0, 0, inum, 'flagv_'//cgrid(igrd), ztmp(:,:), ktype = jp_i4 )
ztmp(:,:) = zmask(:,:)
DO ib_bdy = 1, nb_bdy
DO ib = 1, idx_bdy(ib_bdy)%nblenrim(igrd) ! ntreat defined only for rims 0 and 1
ii = idx_bdy(ib_bdy)%nbi(ib,igrd)
ij = idx_bdy(ib_bdy)%nbj(ib,igrd)
ztmp(ii,ij) = REAL(idx_bdy(ib_bdy)%ntreat(ib,igrd), wp) + 10.
IF( zmask(ii,ij) == 0. ) ztmp(ii,ij) = - ztmp(ii,ij)
END DO
END DO
CALL iom_rstput( 0, 0, inum, 'ntreat_'//cgrid(igrd), ztmp(:,:), ktype = jp_i4 )
END DO
CALL iom_close( inum )
END SUBROUTINE bdy_meshwri
!!=================================================================================
END MODULE bdyini