*IF DEF,C92_1A,OR,DEF,MAKEBC UIE3F404.67
C ******************************COPYRIGHT****************************** GTS2F400.12025
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.12026
C GTS2F400.12027
C Use, duplication or disclosure of this code is subject to the GTS2F400.12028
C restrictions as set forth in the contract. GTS2F400.12029
C GTS2F400.12030
C Meteorological Office GTS2F400.12031
C London Road GTS2F400.12032
C BRACKNELL GTS2F400.12033
C Berkshire UK GTS2F400.12034
C RG12 2SZ GTS2F400.12035
C GTS2F400.12036
C If no contract has been raised with this copy of the code, the use, GTS2F400.12037
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.12038
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.12039
C Modelling at the above address. GTS2F400.12040
C ******************************COPYRIGHT****************************** GTS2F400.12041
C GTS2F400.12042
CLL Subroutine W_LLTOEQ---------------------------------------------- WLLTOE1A.3
CLL WLLTOE1A.4
CLL Purpose: Calculates u and v components of wind on equatorial WLLTOE1A.5
CLL (eq) latitude longitude grid by rotating wind WLLTOE1A.6
CLL components on standard latitude-longitude (eq) WLLTOE1A.7
CLL grid. WLLTOE1A.8
CLL WLLTOE1A.9
CLL Written by A. Dickinson WLLTOE1A.10
CLL WLLTOE1A.11
CLL Model Modification history from model version 3.0: WLLTOE1A.12
CLL version date WLLTOE1A.13
CLL 4.1 31/05/96 The number of v points to be processed on a UIE2F401.386
CLL C grid differs by row_length. u,v therefore UIE2F401.387
CLL treated separately. UIE2F401.388
CLL Author I.Edmond Reviewer D. Goddard UIE2F401.389
CLL WLLTOE1A.14
CLL Programming standard : WLLTOE1A.15
CLL WLLTOE1A.16
CLL Logical components covered : S134 WLLTOE1A.17
CLL WLLTOE1A.18
CLL Project task : WLLTOE1A.19
CLL WLLTOE1A.20
CLL External documentation: WLLTOE1A.21
CLL WLLTOE1A.22
CLL Documentation: The transformation formulae are described in WLLTOE1A.23
CLL unified model on-line documentation paper S1. WLLTOE1A.24
CLL WLLTOE1A.25
CLLEND ----------------------------------------------------------------- WLLTOE1A.26
C WLLTOE1A.27
C*L Arguments:-------------------------------------------------------- WLLTOE1A.28
SUBROUTINE W_LLTOEQ(COEFF1,COEFF2,U,V,U_EQ,V_EQ,POINTS,POINTS2) 7UIE2F401.390
WLLTOE1A.30
IMPLICIT NONE WLLTOE1A.31
WLLTOE1A.32
INTEGER WLLTOE1A.33
* POINTS !IN Number of points to be processed WLLTOE1A.34
*,POINTS2 ! IN Number of v points to be processed UIE2F401.391
WLLTOE1A.35
REAL WLLTOE1A.36
* COEFF1(POINTS) !IN Coefficient of rotation no 1 WLLTOE1A.37
*,COEFF2(POINTS) !IN Coefficient of rotation no 2 WLLTOE1A.38
*,U_EQ(POINTS) !OUT u component of wind on equatorial grid WLLTOE1A.39
*,V_EQ(POINTS) !OUT v component of wind on equatorial grid WLLTOE1A.40
*,U(POINTS) !IN u component of wind on lat-lon grid WLLTOE1A.41
*,V(POINTS) !IN v component of wind on lat-lon grid WLLTOE1A.42
C Workspace usage:----------------------------------------------------- WLLTOE1A.43
C None WLLTOE1A.44
C---------------------------------------------------------------------- WLLTOE1A.45
C External subroutines called:----------------------------------------- WLLTOE1A.46
C None WLLTOE1A.47
C*--------------------------------------------------------------------- WLLTOE1A.48
C Define local varables:----------------------------------------------- WLLTOE1A.49
INTEGER I WLLTOE1A.50
C---------------------------------------------------------------------- WLLTOE1A.51
WLLTOE1A.52
CL 1. Transform wind components WLLTOE1A.53
C WLLTOE1A.54
C Formulae used are from eq (4.13) WLLTOE1A.55
WLLTOE1A.56
DO 100 I = 1,POINTS WLLTOE1A.57
U_EQ(I)=COEFF1(I)*U(I)-COEFF2(I)*V(I) WLLTOE1A.58
100 CONTINUE WLLTOE1A.60
! On a C grid number of u,v points processed differs by row_length UIE2F401.392
DO I = 1,POINTS2 UIE2F401.393
V_EQ(I)=COEFF1(I)*V(I)+COEFF2(I)*U(I) UIE2F401.394
END DO UIE2F401.395
WLLTOE1A.61
RETURN WLLTOE1A.62
END WLLTOE1A.63
*ENDIF WLLTOE1A.64