!-------------------------------------- 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 afcartauv - Adjoint of cartauv. *subroutine afcartauv(pcart,pu,pv,pcoslonob,psinlonob,prsinlatob,knk) 1 ! IMPLICIT NONE integer knk real*8 pcoslonob,psinlonob,prsinlatob real*8 pcart(3,knk),pu(knk),pv(knk) ! !author Luc Fillion - ARMA/MSC - 7 Feb 05: ! from ez_CARTAUV written by Michel roch - april 90 !Revision: L. Fillion - 25 Feb 2005 - Fast version. !arguments ! #include "taglam4d.cdk"
integer jk real*8 zfac1,zfac2 ! !! do jk=1,knk pcart(1,jk)=0.0 pcart(2,jk)=0.0 pcart(3,jk)=0.0 enddo ! zfac1=prsinlatob*pcoslonob zfac2=prsinlatob*psinlonob ! do jk=1,knk pcart(1,jk) = pcart(1,jk)-pv(jk)*zfac1 pcart(2,jk) = pcart(2,jk)-pv(jk)*zfac2 pv(jk) = 0.0 pcart(2,jk) = pcart(2,jk)+pu(jk)*pcoslonob pcart(1,jk) = pcart(1,jk)-pu(jk)*psinlonob pu(jk) = 0.0 enddo ! return end