! { dg-do compile } | |
! { dg-options "-Ofast" } | |
SUBROUTINE sedi_1D(QX1d, DZ1d,kdir,BX1d,kbot,ktop) | |
real, dimension(:) :: QX1d,DZ1d | |
real, dimension(size(QX1d)) :: VVQ | |
logical BX_present | |
do k= kbot,ktop,kdir | |
VVQ= VV_Q0 | |
enddo | |
Vxmaxx= min0 | |
if (kdir==1) then | |
dzMIN = minval(DZ1d) | |
endif | |
npassx= Vxmaxx/dzMIN | |
DO nnn= 1,npassx | |
if (BX_present) then | |
do k= ktop,kdir | |
BX1d= iDZ1d0 | |
enddo | |
endif | |
ENDDO | |
END |