MODULE HCON
 
      USE MACHINE
      IMPLICIT NONE
      SAVE
 
 
      real (kind=kind_rad) AMOLWT, CSUBP,    DIFFCTR,  G, GINV
     &,                     GRAVDR, O3DIFCTR, P0,       P0INV
     &,                     GP0INV, P0XZP2,   P0XZP8,   P0X2
     &,                     RADCON, RADCON1,  RATCO2MW, RATH2OMW
     &,                     RGASK,  RGASSP,   SECPDA
!
!***THE FOLLOWING DATA ARE LEVEL-INDEPENDENT
      DATA G/980.665/
!
      real (kind=kind_rad) HUNDRED, HNINETY, SIXTY, FIFTY, TEN, EIGHT
     &,                     FIVE,    FOUR,    THREE, TWO,   ONE, HAF
     &,                     QUARTR,  ZERO
!
      real (kind=kind_rad) H83E26,   H71E26,  H1E15 ,  H1E13,   H1E11
     &,                     H1E8,     H2E6,    H1E6,    H69766E5,H4E5
     &,                     H165E5,   H5725E4, H488E4,  H1E4,    H24E3
     &,                     H20788E3, H2075E3, H18E3,   H1224E3
     &,                     H67390E2, H5E2,    H3082E2, H3E2,    H2945E2
     &,                     H29316E2, H26E2,   H25E2,   H23E2,   H2E2
     &,                     H15E2,    H1386E2, H1036E2, H8121E1, H35E1
     &,                     H3116E1,  H28E1,   H181E1,  H18E1,   H161E1
     &,                     H16E1,    H1226E1, H9P94,   H6P08108,H3P6
     &,                     H3P5,     H2P9,    H2P8,    H2P5,    H1P8
     &,                     H1P4387,  H1P41819,H1P4,    H1P25892,H1P082
     &,                     HP816,    HP805,   HP8,     HP60241
     &,                     HP602409, HP6,     HP526315,HP518,   HP5048
     &,                     HP3795,   HP369,   HP26,    HP228,   HP219
     &,                     HP166666, HP144,   HP118666,HP1
       real (kind=kind_rad) H658M2,  H625M2,   H44871M2, H44194M2
     &,                     H42M2,   H41666M2, H28571M2, H2118M2
     &,                     H129M2,  H1M2,     H559M3,   H3M3
     &,                     H235M3,  H1M3,     H987M4,   H323M4
     &,                     H3M4,    H285M4,   H1M4,     H75826M4
     &,                     H6938M5, H394M5,   H37412M5, H15M5
     &,                     H1439M5, H128M5,   H102M5,   H1M5
     &,                     H7M6,    H4999M6,  H451M6,   H25452M6
     &,                     H1M6,    H391M7,   H1174M7,  H8725M8
     &,                     H327M8,  H257M8,   H1M8,     H23M10
     &,                     H14M10,  H11M10,   H1M10,    H83M11
     &,                     H82M11,  H8M11,    H77M11,   H72M11
     &,                     H53M11 , H48M11,   H44M11,   H42M11
     &,                     H37M11,  H35M11,   H32M11,   H3M11
     &,                     H28M11,  H24M11,   H23M11,   H2M11
     &,                     H18M11,  H15M11,   H14M11,   H114M11
     &,                     H11M11,  H1M11,    H96M12,   H93M12
     &,                     H77M12,  H74M12,   H65M12,   H62M12
     &,                     H6M12,   H45M12,   H44M12,   H4M12
     &,                     H38M12,  H37M12,   H3M12,    H29M12
     &,                     H28M12,  H24M12,   H21M12,   H16M12
     &,                     H14M12,  H12M12,   H8M13,    H46M13
     &,                     H36M13,  H135M13,  H12M13,   H1M13
     &,                     H3M14,   H15M14,   H14M14,   H101M16
     &,                     H1M16,   H1M17,    H1M18,    H1M19
     &,                     H1M20,   H1M21,    H1M22,    H1M23
     &,                     H1M24,   H26M30,   H14M30,   H25M31
     &,                     H21M31,  H12M31,   H9M32,    H55M32
     &,                     H45M32,  H4M33,    H62M34
!
      real (kind=kind_rad) HM2M2,    HM6666M2, HMP5,     HMP575
     &,                     HMP66667, HMP805,   HM1EZ,    HM13EZ
     &,                     HM19EZ,   HM1E1,    HM1597E1, HM161E1
     &,                     HM1797E1, HM181E1,  HM8E1,    HM1E2
!
      END MODULE HCON
 
      MODULE RNDDTA
 
      USE MACHINE
      IMPLICIT NONE
 
       SAVE
!
      integer nblw, nblx, nbly, nblm
!
      PARAMETER (NBLW=163,NBLX=47,NBLY=15)
      PARAMETER (NBLM=NBLY-1)
!
!    COMMON BLOCK BANDTA CONTAINS RANDOM BAND PARAMETERS FOR THE LW
!    CALCULATIONS USING 10 CM-1 WIDE BANDS.THE 15 UM CO2 COMPLEX
!    IS 2 BANDS,560-670 AND 670-800 CM-1. OZONE COEFFICIENTS ARE
!    IN 3 BANDS,670-800 (14.1 UM),990-1070 AND 1070-1200 (9.6 UM).
!    THE  (NBLW) BANDS NOW INCLUDE:
!                56 BANDS, 10  CM-1 WIDE    0  -   560  CM-1
!                 2 BANDS, 15 UM COMPLEX  560  -   670  CM-1
!                                         670  -   800  CM-1
!                 3 "CONTINUUM" BANDS     800  -   900  CM-1
!                                         900  -   990  CM-1
!                                        1070  -   1200 CM-1
!                 1 BAND FOR 9.6 UM BAND  990  -   1070 CM-1
!               100 BANDS, 10 CM-1 WIDE  1200  -   2200 CM-1
!                 1 BAND FOR 4.3 UM SRC  2270  -   2380 CM-1
!    THUS NBLW PRESENTLY EQUALS    163
!    ALL BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER
!
!        ARNDM   =   RANDOM "A" PARAMETER FOR (NBLW) BANDS
!        BRNDM   =   RANDOM "B" PARAMETER FOR (NBLW) BANDS
!        BETAD   =   CONTINUUM COEFFICIENTS FOR (NBLW) BANDS
!        AP,BP   =   CAPPHI COEFFICIENTS FOR (NBLW) BANDS
!        ATP,BTP =   CAPPSI COEFFICIENTS FOR (NBLW) BANDS
!        BANDLO  =   LOWEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS
!        BANDHI  =   HIGHEST FREQUENCY IN EACH OF (NBLW) FREQ. BANDS
!        AO3RND  =   RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE
!                    BANDS
!        BO3RND  =   RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE
!                    BANDS
!        AB15    =   THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS
!                    REPRESENTING THE 15 UM BAND COMPLEX OF CO2
!     DATA FOR ARNDM,BRNDM,AP,BP,ATP,BTP,AO3RND,BO3RND ARE OBTAINED BY
!     USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS ARE FROM
!     ROBERTS (1976).
      real (kind=kind_rad) ARNDM(NBLW),BRNDM(NBLW),BETAD(NBLW),AP(NBLW)
     &,                 BP(NBLW),ATP(NBLW),BTP(NBLW),BANDLO(NBLW)
     &,                 BANDHI(NBLW),AO3RND(3),BO3RND(3),AB15(2)
      DATA AO3RND /
     &   0.543368E+02,  0.234676E+04,  0.384881E+02/
      DATA BO3RND /
     &   0.526064E+01,  0.922424E+01,  0.496515E+01/
!
!    COMMON BLOCK BDWIDE CONTAINS RANDOM BAND PARAMETERS FOR SPECIFIC
!    WIDE BANDS. AT PRESENT,THE INFORMATION CONSISTS OF 1) RANDOM
!    MODEL PARAMETERS FOR THE 15 UM BAND,560-800 CM-1; 2) THE
!    CONTINUUM COEFFICIENT FOR THE 800-990,1070-1200 CM-1 BAND
!        SPECIFICALLY:
!        AWIDE       =   RANDOM "A" PARAMETER FOR  BAND
!        BWIDE       =   RANDOM "B" PARAMETER FOR  BAND
!        BETAWD      =   CONTINUUM COEFFICIENTS FOR BAND
!        APWD,BPWD   =   CAPPHI COEFFICIENTS FOR  BAND
!        ATPWD,BTPWD =   CAPPSI COEFFICIENTS FOR BAND
!        BDLOWD      =   LOWEST FREQUENCY IN EACH  FREQ  BAND
!        BDHIWD      =   HIGHEST FREQUENCY IN EACH FREQ  BAND
!        AB15WD      =   THE PRODUCT ARNDM*BRNDM FOR THE ONE BAND
!                        REPRESENTING THE 15 UM BAND COMPLEX OF CO2
!        BETINW      =   CONT.COEFFICIENT FOR A SPECIFIED WIDE
!                        FREQ.BAND (800-990 AND 1070-1200 CM-1).
!        SKO2D       =   1./BETINW, USED IN SPA88 FOR CONT. COEFFS
!        SKC1R       =   BETAWD/BETINW, USED FOR CONT. COEFF. FOR
!                        15 UM BAND IN FST88
!        SKO3R       =   RATIO OF CONT. COEFF. FOR 9.9 UM BAND TO
!                        BETINW, USED FOR 9.6 UM CONT COEFF IN FST88
!     DATA FOR AWIDE,BWIDE,APWD,BPWD,ATPWD,BTPWD,AO3WD,BO3WD ARE
!     OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS
!     ARE FROM ROBERTS (1976).
      real (kind=kind_rad) AWIDE,BWIDE,BETAWD,APWD,BPWD,ATPWD,BTPWD,
     &                  BDLOWD,BDHIWD,BETINW,AB15WD,SKO2D,SKC1R,SKO3R
      DATA AWIDE  /
     &   0.309801E+01/
      DATA BWIDE  /
     &   0.495357E-01/
      DATA APWD   /
     &   0.177115E-01/
      DATA BPWD   /
     &  -0.545226E-04/
      DATA ATPWD  /
     &   0.187967E-01/
      DATA BTPWD  /
     &  -0.567449E-04/
      DATA BETAWD /
     &   0.347839E+02/
      DATA BETINW /
     &   0.766811E+01/
      DATA BDLOWD /
     &   0.560000E+03/
      DATA BDHIWD /
     &   0.800000E+03/
!
!    COMMON BLOCK BDCOMB CONTAINS RANDOM BAND PARAMETERS FOR THE LW
!    CALCULATIONS USING COMBINED WIDE FREQUENCY BANDS BETWEEN 160 AND
!    1200 CM-1,AS WELL AS THE 2270-2380 BAND FOR SOURCE CALC.
!        BANDS 1-8: COMBINED WIDE FREQUENCY BANDS FOR 160-560 CM-1
!        BANDS 9-14: FREQUENCY BANDS,AS IN BANDTA (NARROW BANDS)
!                    FOR 560-1200 CM-1
!        BAND  15:  FREQUENCY BAND 2270-2380 CM-1,USED FOR SOURCE
!                   CALCULATION ONLY
!        THUS NBLY PRESENTLY EQUALS   15
!
!        BANDS ARE ARRANGED IN ORDER OF INCREASING WAVENUMBER
!        ACOMB       =   RANDOM "A" PARAMETER FOR (NBLY) BANDS
!        BCOMB       =   RANDOM "B" PARAMETER FOR (NBLY) BANDS
!        BETACM      =   CONTINUUM COEFFICIENTS FOR (NBLY) BANDS
!        APCM,BPCM   =   CAPPHI COEFFICIENTS FOR (NBLY) BANDS
!        ATPCM,BTPCM =   CAPPSI COEFFICIENTS FOR (NBLY) BANDS
!        BDLOCM      =   LOWEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS
!        BDHICM      =   HIGHEST FREQUENCY IN EACH OF (NBLY) FREQ. BANDS
!        AO3CM       =   RANDOM "A" PARAMETER FOR OZONE IN (3) OZONE
!                        BANDS
!        BO3CM       =   RANDOM "B" PARAMETER FOR OZONE IN (3) OZONE
!                        BANDS
!        AB15CM      =   THE PRODUCT ARNDM*BRNDM FOR THE TWO BANDS
!                        REPRESENTING THE 15 UM BAND COMPLEX OF CO2
!        BETINC      =   CONT.COEFFICIENT FOR A SPECIFIED WIDE
!                        FREQ.BAND (800-990 AND 1070-1200 CM-1).
!        IBAND       =   INDEX NO OF THE 40 WIDE BANDS USED IN
!                        COMBINED WIDE BAND CALCULATIONS. IN OTHER
!                        WORDS,INDEX TELLING WHICH OF THE 40 WIDE
!                        BANDS BETWEEN 160-560 CM-1 ARE INCLUDED IN
!                        EACH OF THE FIRST 8 COMBINED WIDE BANDS
!     DATA FOR ACOMB,BCOMB,APCM,BPCM,ATPCM,BTPCM,AO3CM,BO3CM ARE
!     OBTAINED BY USING THE AFGL 1982 CATALOG. CONTINUUM COEFFICIENTS
!     ARE FROM ROBERTS (1976). IBAND INDEX VALUES ARE OBTAINED BY
!     EXPERIMENTATION.
      integer  IBAND(40)
      real (kind=kind_rad) ACOMB(NBLY),BCOMB(NBLY),
     &                  BETACM(NBLY),APCM(NBLY),BPCM(NBLY),ATPCM(NBLY),
     &                  BTPCM(NBLY),BDLOCM(NBLY),BDHICM(NBLY),BETINC,
     &                  AO3CM(3),BO3CM(3),AB15CM(2)
      DATA ACOMB  /
     &   0.152070E+05,  0.332194E+04,  0.527177E+03,  0.163124E+03,
     &   0.268808E+03,  0.534591E+02,  0.268071E+02,  0.123133E+02,
     &   0.600199E+01,  0.640803E+00,  0.501549E-01,  0.167961E-01,
     &   0.178110E-01,  0.170166E+00,  0.537083E-02/
      DATA BCOMB  /
     &   0.152538E+00,  0.118677E+00,  0.103660E+00,  0.100119E+00,
     &   0.127518E+00,  0.118409E+00,  0.904061E-01,  0.642011E-01,
     &   0.629660E-01,  0.643346E-01,  0.717082E-01,  0.629730E-01,
     &   0.875182E-01,  0.857907E-01,  0.214005E+00/
      DATA APCM   /
     &  -0.671879E-03,  0.654345E-02,  0.143657E-01,  0.923593E-02,
     &   0.117022E-01,  0.159596E-01,  0.181600E-01,  0.145013E-01,
     &   0.170062E-01,  0.233303E-01,  0.256735E-01,  0.274745E-01,
     &   0.279259E-01,  0.197002E-01,  0.349782E-01/
      DATA BPCM   /
     &  -0.113520E-04, -0.323965E-04, -0.448417E-04, -0.230779E-04,
     &  -0.361981E-04, -0.145117E-04,  0.198349E-04, -0.486529E-04,
     &  -0.550050E-04, -0.684057E-04, -0.447093E-04, -0.778390E-04,
     &  -0.982953E-04, -0.772497E-04, -0.748263E-04/
      DATA ATPCM  /
     &  -0.106346E-02,  0.641531E-02,  0.137362E-01,  0.922513E-02,
     &   0.136162E-01,  0.169791E-01,  0.206959E-01,  0.166223E-01,
     &   0.171776E-01,  0.229724E-01,  0.275530E-01,  0.302731E-01,
     &   0.281662E-01,  0.199525E-01,  0.370962E-01/
      DATA BTPCM  /
     &  -0.735731E-05, -0.294149E-04, -0.505592E-04, -0.280894E-04,
     &  -0.492972E-04, -0.341508E-04, -0.362947E-04, -0.250487E-04,
     &  -0.521369E-04, -0.746260E-04, -0.744124E-04, -0.881905E-04,
     &  -0.933645E-04, -0.664045E-04, -0.115290E-03/
      DATA BETACM /
     &   0.000000E+00,  0.000000E+00,  0.000000E+00,  0.000000E+00,
     &   0.188625E+03,  0.144293E+03,  0.174098E+03,  0.909366E+02,
     &   0.497489E+02,  0.221212E+02,  0.113124E+02,  0.754174E+01,
     &   0.589554E+01,  0.495227E+01,  0.000000E+00/
      DATA IBAND  /
     &    2,   1,   2,   2,   1,   2,   1,   3,   2,   2,
     &    3,   2,   2,   4,   2,   4,   2,   3,   3,   2,
     &    4,   3,   4,   3,   7,   5,   6,   7,   6,   5,
     &    7,   6,   7,   8,   6,   6,   8,   8,   8,   8/
!
!***THE FOLLOWING DATA ARE LEVEL-INDEPENDENT
!CCCC DATA RCO2/3.3E-4/
!     DATA G/980.665/
!CCCC DATA CTAUDA/.5/
!CCCC DATA CSOLAR/1.96/
!CCCC DATA CCOSZ/.5/
!   B0,B1,B2,B3 ARE COEFFICIENTS USED TO CORRECT FOR THE USE OF 250K IN
!   THE PLANCK FUNCTION USED IN EVALUATING PLANCK-WEIGHTED CO2
!   TRANSMISSION FUNCTIONS. (SEE REF. 4)
!     DATA B0,B1,B2,B3/-.51926410E-4,-.18113332E-3,
!    1 -.10680132E-5,-.67303519E-7/
!    *******************************************************************
!    *                                                                 *
!    *                           B L C K F S  FROM G F D L             *
!    *  UNUSED DATA CLEANED OUT - NOV 86 AND MAR 89 ..K.A.CAMPANA....  *
!    &                                                                 *
!    *******************************************************************
!
!     FOR SEASONAL VARIATION
!        SEASON=1,2,3,4 FOR WINTER,SPRING,SUMMER,FALL ONLY (NOT ACTIVE)
!        SEASON=5 - SEASONAL VARIATION(I.E.INTERPOLATE TO DAY OF FCST)
!
!     INTEGER SEASON
!     COMMON/DIUCON/SEASON,FCSTDA,JTIME(5),DAZ(12),JDNMC,
!    &              FJDNMC,TSLAG,RLAG,TIMIN,TPI,HPI,YEAR,DAY,DHR,IXXXX
!     DATA SEASON/5/
!     DATA TSLAG/45.25/,  RLAG/14.8125/
!     DATA DAY/86400./,  YEAR/365.25/
!     DATA TPI/6.283185308/,  HPI/1.570796327/
!     DATA DAY/86400./,  YEAR/365.25/
!     DATA TPI/6.283185308/,  HPI/1.570796327/
!     DATA JTIME/0,1,0,0,0/
!     DATA DHR/2./
!     DATA DAZ/0.,31.,59.,90.,120.,151.,181.,212.,243.,273.,304.,334./
!
!     SEA SURFACE ALBEDO DATA
!
      real (kind=kind_rad) ALBD(21,20),ZA(20),TRN(21),DZA(19)
      DATA ZA/90.,88.,86.,84.,82.,80.,78.,76.,74.,70.,66.,62.,58.,54.,
     &        50.,40.,30.,20.,10.,0.0/
      DATA TRN/.00,.05,.10,.15,.20,.25,.30,.35,.40,.45,.50,.55,.60,.65,
     &         .70,.75,.80,.85,.90,.95,1.00/
      DATA DZA/8*2.0,6*4.0,5*10.0/
!
      real (kind=kind_rad)
     &      EM1V(5040),EM1VW(5040),T1(5041),
     &      T2(5040),T4(5040),EM3V(5040)
cjfe &      EM1(28,180),EM1WDE(28,180),TABLE1(28,180),
cjfe &      TABLE2(28,180),TABLE3(28,180),EM3(28,180)
!
      real (kind=kind_rad) DELCM(NBLY)
      DATA DELCM  /
     &   0.300000E+02,  0.110000E+03,  0.600000E+02,  0.400000E+02,
     &   0.200000E+02,  0.500000E+02,  0.400000E+02,  0.500000E+02,
     &   0.110000E+03,  0.130000E+03,  0.100000E+03,  0.900000E+02,
     &   0.800000E+02,  0.130000E+03,  0.110000E+03/
!
!      real (kind=kind_rad) SC
!      DATA SC/2.0/
!
      END MODULE RNDDTA
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      MODULE COMCD1
 
      USE MACHINE
      IMPLICIT NONE
 
      SAVE
 
cmy delete rocp alread in module physcons
      REAL (KIND=kind_rad) PTOPC(4,2),CVTOP,VVCLD(2),CLAPSE
      REAL (KIND=kind_rad) CRHRH,PSTRT,CLAPKC,DCLPS,CLPSE
      INTEGER KLOWT,KLOWB,LLYR
 
      END MODULE COMCD1