*IF DEF,A14_1A ENGSUM1A.2
C ******************************COPYRIGHT****************************** GTS2F400.2395
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.2396
C GTS2F400.2397
C Use, duplication or disclosure of this code is subject to the GTS2F400.2398
C restrictions as set forth in the contract. GTS2F400.2399
C GTS2F400.2400
C Meteorological Office GTS2F400.2401
C London Road GTS2F400.2402
C BRACKNELL GTS2F400.2403
C Berkshire UK GTS2F400.2404
C RG12 2SZ GTS2F400.2405
C GTS2F400.2406
C If no contract has been raised with this copy of the code, the use, GTS2F400.2407
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.2408
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.2409
C Modelling at the above address. GTS2F400.2410
C ******************************COPYRIGHT****************************** GTS2F400.2411
C GTS2F400.2412
CLL SUBROUTINE ENERGY_SUM--------------------------------------------- ENGSUM1A.3
CLL ENGSUM1A.4
CLL PURPOSE : PART OF ENERGY CORRECTION SUITE OF ROUTINES ENGSUM1A.5
CLL - TO SUM A QUANTITY GLOBALLY ON A LEVEL ENGSUM1A.6
CLL ENGSUM1A.7
CLL NOT SUITABLE FOR SINGLE COLUMN MODEL USE ENGSUM1A.8
CLL ENGSUM1A.9
CLL CODE WRITTEN FOR CRAY Y-MP BY D.GREGORY FEBRUARY 1991 ENGSUM1A.10
CLL ENGSUM1A.11
CLL MODEL MODIFICATION HISTORY FROM MODEL VERSION 3.0: ENGSUM1A.12
CLL VERSION DATE ENGSUM1A.13
! 4.1 28/11/95 Changed interface to ENERGY_SUM to make APB5F401.55
! suitable for MPP use. P.Burton APB5F401.56
CLL ENGSUM1A.14
CLL PROGRAMMING STANDARDS : UNIFIED MODEL DOCUMENTATION PAPER NO. 4 ENGSUM1A.15
CLL VERSION NO. 1 ENGSUM1A.16
CLL ENGSUM1A.17
CLL SYSTEM TASK : P## ENGSUM1A.18
CLL ENGSUM1A.19
CLL DOCUMENTATION : UNIFIED MODEL DOCUMENTATION PAPER P### ENGSUM1A.20
CLL ENGSUM1A.21
CLLEND----------------------------------------------------------------- ENGSUM1A.22
C ENGSUM1A.23
C*L ARGUMENTS--------------------------------------------------------- ENGSUM1A.24
C ENGSUM1A.25
SUBROUTINE ENERGY_SUM (VAR,START_POINT,END_POINT,FIELD_SIZE, 3,1APB5F401.57
& AREA,MASS,RS,CONV_FAC,SUM_VAR) APB5F401.58
C ENGSUM1A.27
IMPLICIT NONE ENGSUM1A.28
C ENGSUM1A.29
*CALL C_G 
ENGSUM1A.30
C ENGSUM1A.31
C---------------------------------------------------------------------- ENGSUM1A.32
C VECTOR LENGTHS ENGSUM1A.33
C---------------------------------------------------------------------- ENGSUM1A.34
C ENGSUM1A.35
INTEGER START_POINT, ! IN point to start sum at APB5F401.59
& END_POINT, ! IN point to end sum at APB5F401.60
& FIELD_SIZE ! IN size of field APB5F401.61
C ENGSUM1A.37
C ENGSUM1A.38
C---------------------------------------------------------------------- ENGSUM1A.39
C VARIABLES WHICH ARE INPUT ENGSUM1A.40
C---------------------------------------------------------------------- ENGSUM1A.41
C ENGSUM1A.42
REAL VAR(FIELD_SIZE), ! IN variable to be summed APB5F401.62
& AREA(FIELD_SIZE),! IN area of grid box APB5F401.63
& RS(FIELD_SIZE), ! IN radius at level of the energy sum APB5F401.64
& MASS(FIELD_SIZE) ! IN mass of grid box APB5F401.65
C ENGSUM1A.51
REAL CONV_FAC ! IN CONVERSION FACTOR TO TRANSLATE ENGSUM1A.52
! FLUX INTO ENERGY UNITS ENGSUM1A.53
C ENGSUM1A.54
C ENGSUM1A.55
C---------------------------------------------------------------------- ENGSUM1A.56
C VARIABLES WHICH ARE IN AND OUT ENGSUM1A.57
C---------------------------------------------------------------------- ENGSUM1A.58
C ENGSUM1A.59
REAL SUM_VAR ! INOUT SUM OF ENERGY ENGSUM1A.60
C ENGSUM1A.61
C ENGSUM1A.62
C---------------------------------------------------------------------- ENGSUM1A.63
C VARIABLES WHICH ARE DEFINED LOCALLY - NONE ENGSUM1A.64
C---------------------------------------------------------------------- ENGSUM1A.65
C ENGSUM1A.66
REAL WORK(FIELD_SIZE) ! work space APB5F401.66
C ENGSUM1A.68
C ENGSUM1A.69
C---------------------------------------------------------------------- ENGSUM1A.70
C INTERNAL LOOP COUNTERS ENGSUM1A.71
C---------------------------------------------------------------------- ENGSUM1A.72
C ENGSUM1A.73
INTEGER I ! LOOP COUNTER ENGSUM1A.74
C ENGSUM1A.75
C ENGSUM1A.76
C---------------------------------------------------------------------- ENGSUM1A.77
C EXTERNAL SUBROUTINE CALLS - NONE ENGSUM1A.78
C---------------------------------------------------------------------- ENGSUM1A.79
C ENGSUM1A.80
C*--------------------------------------------------------------------- ENGSUM1A.81
CL ENGSUM1A.82
CL--------------------------------------------------------------------- ENGSUM1A.83
CL SUM VARIABLE OVER LAYER ENGSUM1A.84
CL--------------------------------------------------------------------- ENGSUM1A.85
CL ENGSUM1A.86
C ENGSUM1A.87
C CALCULATE CONTRIBUTION TO TOTAL ENERGY FROM EACH GRID BOX ENGSUM1A.88
C ENGSUM1A.89
DO I=START_POINT,END_POINT APB5F401.67
WORK(I) = RS(I)*RS(I)*AREA(I)*MASS(I)*VAR(I)*CONV_FAC/G ENGSUM1A.91
END DO ENGSUM1A.92
C ENGSUM1A.93
CALL DO_SUMS(WORK,FIELD_SIZE,START_POINT,END_POINT,1,SUM_VAR) APB5F401.68
C ENGSUM1A.97
RETURN ENGSUM1A.98
END ENGSUM1A.99
*ENDIF ENGSUM1A.100