*IF DEF,OCEAN                                                              DRODT.2      
C ******************************COPYRIGHT******************************    DRODT.3      
C (c) CROWN COPYRIGHT 1998, METEOROLOGICAL OFFICE, All Rights Reserved.    DRODT.4      
C                                                                          DRODT.5      
C Use, duplication or disclosure of this code is subject to the            DRODT.6      
C restrictions as set forth in the contract.                               DRODT.7      
C                                                                          DRODT.8      
C                Meteorological Office                                     DRODT.9      
C                London Road                                               DRODT.10     
C                BRACKNELL                                                 DRODT.11     
C                Berkshire UK                                              DRODT.12     
C                RG12 2SZ                                                  DRODT.13     
C                                                                          DRODT.14     
C If no contract has been raised with this copy of the code, the use,      DRODT.15     
C duplication or disclosure of it is strictly prohibited.  Permission      DRODT.16     
C to do so must first be obtained in writing from the Head of Numerical    DRODT.17     
C Modelling at the above address.                                          DRODT.18     
C ******************************COPYRIGHT******************************    DRODT.19     
C                                                                          DRODT.20     
CLL   Subroutine DRODT                                                     DRODT.21     
CLL                                                                        DRODT.22     
CLL   Author: M J Roberts                                                  DRODT.23     
CLL                                                                        DRODT.24     
CLL   Description: Calculate expansion coefficient for temperature         DRODT.25     
CLL                                                                        DRODT.26     
CLL   External documentation:                                              DRODT.27     
CLL       Modular Ocean Model 2 manual and code                            DRODT.28     
CLL                                                                        DRODT.29     
CLL   Implemented at UM vn 4.5 26 June 1998                                DRODT.30     
CLL                                                                        DRODT.31     
CLL   Modification History  :                                              DRODT.32     
CLL   Version   Date     Comment & Name                                    DRODT.33     
CLL   ------- --------  --------------------------------------------       DRODT.34     
CLL                                                                        DRODT.35     
CLL   Subroutine dependencies.                                             DRODT.36     
CLL                                                                        DRODT.37     
CLLEND------------------------------------------------------------------   DRODT.38     
C                                                                          DRODT.39     

      SUBROUTINE DRODT(tx,sx,alpha,imt,km)                                  11DRODT.40     
                                                                           DRODT.41     
      IMPLICIT NONE                                                        DRODT.42     
                                                                           DRODT.43     
*CALL OARRYSIZ                                                             DRODT.44     
*CALL COCSTATE                                                             DRODT.45     
*CALL CNTLOCN                                                              DRODT.46     
                                                                           DRODT.47     
C calculate the expansion coefficient for temperature                      DRODT.48     
                                                                           DRODT.49     
       INTEGER k,imt,km,i                                                  DRODT.50     
       REAL tq,sq,tx(imt,km),sx(imt,km),alpha(imt,km)                      DRODT.51     
                                                                           DRODT.52     
       do k=1,km                                                           DRODT.53     
         do i=1,imt                                                        DRODT.54     
           tq=tx(i,k)-to(k)                                                DRODT.55     
           sq=sx(i,k)-so(k)                                                DRODT.56     
           alpha(i,k) = c(k,1) + (c(k,4) + c(k,7)*sq)*sq                   DRODT.57     
     &                + (2.0*c(k,3) + 2.0*c(k,8)*sq + 3.0*c(k,6)*tq)*tq    DRODT.58     
         enddo                                                             DRODT.59     
       enddo                                                               DRODT.60     
                                                                           DRODT.61     
       RETURN                                                              DRODT.62     
       END                                                                 DRODT.63     
*ENDIF                                                                     DRODT.64     
                                                                           DRODT.65