*IF DEF,C92_1A,OR,DEF,MAKEBC UIE3F404.66
C ******************************COPYRIGHT****************************** GTS2F400.11917
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.11918
C GTS2F400.11919
C Use, duplication or disclosure of this code is subject to the GTS2F400.11920
C restrictions as set forth in the contract. GTS2F400.11921
C GTS2F400.11922
C Meteorological Office GTS2F400.11923
C London Road GTS2F400.11924
C BRACKNELL GTS2F400.11925
C Berkshire UK GTS2F400.11926
C RG12 2SZ GTS2F400.11927
C GTS2F400.11928
C If no contract has been raised with this copy of the code, the use, GTS2F400.11929
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.11930
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.11931
C Modelling at the above address. GTS2F400.11932
C ******************************COPYRIGHT****************************** GTS2F400.11933
C GTS2F400.11934
CLL Subroutine W_EQTOLL---------------------------------------------- WEQTOL1A.3
CLL WEQTOL1A.4
CLL Purpose: Calculates u and v components of wind on standard WEQTOL1A.5
CLL latitude-longitude grid by rotating wind WEQTOL1A.6
CLL components on equatorial latitude-longitude (eq) WEQTOL1A.7
CLL grid. WEQTOL1A.8
CLL WEQTOL1A.9
CLL Written by A. Dickinson WEQTOL1A.10
CLL WEQTOL1A.11
CLL Model Modification history from model version 3.0: WEQTOL1A.12
CLL version date WEQTOL1A.13
CLL 4.1 31/05/96 The number of v points to be processed on a UIE2F401.376
CLL C grid differs by row_length. u,v therefore UIE2F401.377
CLL treated separately. UIE2F401.378
CLL Author I.Edmond Reviewer D. Goddard UIE2F401.379
CLL WEQTOL1A.14
CLL Programming standard : WEQTOL1A.15
CLL WEQTOL1A.16
CLL Logical components covered : S133 WEQTOL1A.17
CLL WEQTOL1A.18
CLL Project task : WEQTOL1A.19
CLL WEQTOL1A.20
CLL Documentation: The transformation formulae are described in WEQTOL1A.21
CLL unified model on-line documentation paper S1. WEQTOL1A.22
CLL WEQTOL1A.23
CLLEND ----------------------------------------------------------------- WEQTOL1A.24
C WEQTOL1A.25
C*L Arguments:-------------------------------------------------------- WEQTOL1A.26
SUBROUTINE W_EQTOLL(COEFF1,COEFF2,U_EQ,V_EQ,U,V,POINTS,POINTS2) 11UIE2F401.380
WEQTOL1A.28
IMPLICIT NONE WEQTOL1A.29
WEQTOL1A.30
INTEGER WEQTOL1A.31
* POINTS !IN Number of points to be processed WEQTOL1A.32
*,POINTS2 ! IN Number of v points to be processed UIE2F401.381
WEQTOL1A.33
REAL WEQTOL1A.34
* COEFF1(POINTS) !IN Coefficient of rotation no 1 WEQTOL1A.35
*,COEFF2(POINTS) !IN Coefficient of rotation no 2 WEQTOL1A.36
*,U_EQ(POINTS) !IN u component of wind on equatorial grid WEQTOL1A.37
*,V_EQ(POINTS) !IN v component of wind on equatorial grid WEQTOL1A.38
*,U(POINTS) !OUT u component of wind on lat-lon grid WEQTOL1A.39
*,V(POINTS) !OUT v component of wind on lat-lon grid WEQTOL1A.40
C Workspace usage:----------------------------------------------------- WEQTOL1A.41
C None WEQTOL1A.42
C---------------------------------------------------------------------- WEQTOL1A.43
C External subroutines called:----------------------------------------- WEQTOL1A.44
C None WEQTOL1A.45
C*--------------------------------------------------------------------- WEQTOL1A.46
C Define local varables:----------------------------------------------- WEQTOL1A.47
INTEGER I WEQTOL1A.48
C---------------------------------------------------------------------- WEQTOL1A.49
WEQTOL1A.50
CL 1. Transform wind components WEQTOL1A.51
C WEQTOL1A.52
C Formulae used are from eq (4.14) WEQTOL1A.53
WEQTOL1A.54
DO 100 I = 1,POINTS WEQTOL1A.55
U(I)=COEFF1(I)*U_EQ(I)+COEFF2(I)*V_EQ(I) WEQTOL1A.56
100 CONTINUE WEQTOL1A.58
! On a C grid number of u,v points processed differs by row_length UIE2F401.382
DO I = 1,POINTS2 UIE2F401.383
V(I)=COEFF1(I)*V_EQ(I)-COEFF2(I)*U_EQ(I) UIE2F401.384
END DO UIE2F401.385
WEQTOL1A.59
RETURN WEQTOL1A.60
END WEQTOL1A.61
*ENDIF WEQTOL1A.62