!-------------------------------------- LICENCE BEGIN ------------------------------------ !Environment Canada - Atmospheric Science and Technology License/Disclaimer, ! version 3; Last Modified: May 7, 2008. !This is free but copyrighted software; you can use/redistribute/modify it under the terms !of the Environment Canada - Atmospheric Science and Technology License/Disclaimer !version 3 or (at your option) any later version that should be found at: !http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html ! !This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; !without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. !See the above mentioned License/Disclaimer for more details. !You should have received a copy of the License/Disclaimer along with this software; !if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec), !CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca !-------------------------------------- LICENCE END -------------------------------------- ***s/r tab_vec - *subroutine tab_vec (tab,Minx,Maxx,Miny,Maxy,nk, 9 $ vec,i0,in,j0,jn,idir) implicit none * integer Minx,Maxx,Miny,Maxy,nk,i0,in,j0,jn,idir real*8 tab(Minx:Maxx,Miny:Maxy,nk),vec(*) * * author Abdessamad Qaddouri - December 2006 * *revision * v3_30 - Qaddouri A. - initial version * integer iloc,i,j,k * * --------------------------------------------------------------- * iloc=0 * if (idir.eq.1) then do k = 1 , nk do j = j0 , jn do i = i0 , in iloc=iloc+1 vec(iloc) =tab(i,j,k) enddo enddo enddo else if (idir.eq.(-1)) then do k = 1 , nk do j = j0 , jn do i = i0 , in iloc=iloc+1 tab(i,j,k) =vec(iloc) enddo enddo enddo endif * * --------------------------------------------------------------- * return end