!-------------------------------------- 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 v4d_uv2img_tr - Equivalent to v4d_uv2img for TRAJECTORY * #include "model_macros_f.h"*
subroutine v4d_uv2img_tr() 2 * #include "impnone.cdk"
* *author * M.Tanguay * *revision * v2_10 - Tanguay M. - initial MPI version * v3_03 - Tanguay M. - use Ind comdeck * *object * see id section * *arguments * none * *implicits #include "glb_ld.cdk"
#include "geomg.cdk"
#include "lun.cdk"
#include "dcst.cdk"
#include "indm.cdk"
* integer i,j,k real*8 pdsc1_8 * if (Lun_out.gt.0) write(Lun_out,1000) * * Calculation of wind images from components * ------------------------------------------ do j= 1, l_nj pdsc1_8 = Geomg_cy_8(j) / Dcst_rayt_8 do k=1,G_nk do i= 1, l_niu Indm_um(i,j,k) = pdsc1_8 * Indm_um(i,j,k) end do end do end do * do j= 1, l_njv pdsc1_8 = Geomg_cyv_8(j) / Dcst_rayt_8 do k=1,G_nk do i= 1, l_ni Indm_vm(i,j,k) = pdsc1_8 * Indm_vm(i,j,k) end do end do end do * 1000 format( +/,'TRAJ OF CONVERT COMPONENTS TO WIND IMAGES: (S/R V4D_UV2IMG_TR)', +/,'==============================================================', +/) * return end