!-------------------------------------- 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 fcartauv - compute the components of the winds in the rotated system of * coordinates from the winds in the rotated cartesian space ! Fast version of cartauv used by Variational analysis. *subroutine fcartauv(pu,pv,pcart,pcoslonob,psinlonob,prsinlatob,knk) 1 ! IMPLICIT NONE integer knk real*8 pcart(3,knk),pu(knk),pv(knk) real*8 pcoslonob,psinlonob,prsinlatob ! !author Luc Fillion - ARMA/MSC - 7 Feb 05: ! from ez_CARTAUV written by Michel roch - april 90 ! !Revision: L. Fillion - ARMA/MSC - 25 Feb 05: Fast version from cartauv. ! !arguments integer jk real*8 zfac1,zfac2 #include "taglam4d.cdk"
! !! zfac1=-prsinlatob*pcoslonob zfac2=-prsinlatob*psinlonob ! do jk=1,knk pu(jk) = -psinlonob*pcart(1,jk) + pcoslonob*pcart(2,jk) pv(jk) = zfac1*pcart(1,jk) + zfac2*pcart(2,jk) enddo ! RETURN END