!-------------------------------------- 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 MFOTVT_TL *SUBROUTINE MFOTVT_TL(TV,TT,QQ,TT5,QQ5,NI,NK,N) 1 * #include "impnone.cdk"
INTEGER NI, NK, N REAL TV(NI,NK), TT(NI,NK), QQ(NI,NK), TT5(NI,NK), QQ5(NI,NK) * *Author * A. Zadra (May 2002) * *Object * Tangent linear of MFOTVT * *IMPLICITES #include "consphy.cdk"
*MODULES EXTERNAL INCTPHY * ** *-------------------------------------------------------------------- INTEGER K, I #include "dintern.cdk"
#include "fintern.cdk"
*-------------------------------------------------------------------- #include "initcph.cdk"
* DO K=1,NK DO I=1,N TV(I,K) = DBLE(TT(I,K)) * (1.0D0 + DBLE(DELTA)*DBLE(QQ5(I,K))) + + DBLE(TT5(I,K))* DBLE(DELTA)*DBLE(QQ(I,K)) ENDDO ENDDO * RETURN CONTAINS #include "fintern90.cdk"
END