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

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