--- trunk/Sources/filtrez/filtreg_v.f 2015/06/05 19:21:08 141 +++ trunk/filtrez/filtreg_v.f 2018/03/20 09:35:59 265 @@ -12,7 +12,7 @@ ! Matrix filter on longitudes. Matrices have already been ! computed. On v grid, direct filter only. - USE dimens_m, ONLY: iim, jjm + USE dimensions, ONLY: iim, jjm use filtreg_hemisph_m, only: filtreg_hemisph USE inifgn_m, ONLY: sddu, unsddu use inifilr_m, only: jfiltnv, jfiltsv, matricevn, matricevs @@ -32,11 +32,7 @@ call assert(size(champ, 1) == iim + 1, "filtreg_v iim + 1") call assert(size(champ, 2) == jjm, "filtreg_v jjm") - IF (intensive) THEN - sdd = sddu - ELSE - sdd = unsddu - END IF + sdd = merge(sddu, unsddu, intensive) call filtreg_hemisph(champ(:, :jfiltnv, :), sdd, matricevn) call filtreg_hemisph(champ(:, jfiltsv:jjm, :), sdd, matricevs)