/[MITgcm]/MITgcm/pkg/fizhi/fizhi_lsm.F

Diff of /MITgcm/pkg/fizhi/fizhi_lsm.F

Parent Directory | Revision Log | View Revision Graph Revision Graph | View Patch Patch

-revision 1.2 by molod,
Fri Jul 16 19:52:21 2004 UTC
+revision 1.6 by ce107,
Thu Jun 16 16:46:12 2005 UTC
 Line 1
  C $Header$
  C $Name$
+ #include "FIZHI_OPTIONS.h"
        SUBROUTINE TILE (
       I                 NCH, DTSTEP, ITYP, TRAINL,TRAINC, TSNOW,  UM,
       I                 ETURB,  DEDQA,  DEDTC,  HSTURB, DHSDQA, DHSDTC,
-Line 30 
 C****
+Line 31 
 C****
  C****
        INTEGER NCH
        INTEGER ITYP(NCH)
-       REAL  DTSTEP,  TRAINL(NCH), TRAINC(NCH), TSNOW(NCH), UM(NCH),
+       _RL  DTSTEP,  TRAINL(NCH), TRAINC(NCH), TSNOW(NCH), UM(NCH),
       &       ETURB(NCH),  DEDQA(NCH), HSTURB(NCH), DHSDTC(NCH),
       &         TM (NCH),     CD(NCH), SUNANG(NCH), DHSDQA(NCH),
       &         QM (NCH), PARDIR(NCH), PARDIF(NCH),  SWNET(NCH),
       &      HLWDWN(NCH),   PSUR(NCH),   ZLAI(NCH),  GREEN(NCH),
       &          Z2(NCH), SQSCAT(NCH),  DEDTC(NCH)
-       REAL  RSOIL1(NCH), RSOIL2(NCH),    RDC(NCH),  U2FAC(NCH),
+       _RL  RSOIL1(NCH), RSOIL2(NCH),    RDC(NCH),  U2FAC(NCH),
       &     QSATTC(NCH), DQSDTC(NCH),  ALWRAD(NCH), BLWRAD(NCH),
       &          TC(NCH),     TD(NCH),     QA(NCH),   BOMB(NCH),
       &       SWET1(NCH),  SWET2(NCH),  SWET3(NCH),  CAPAC(NCH),
       &        SNOW(NCH),   EVAP(NCH), SHFLUX(NCH), RUNOFF(NCH)
-       REAL   EINT(NCH),    ESOI(NCH),   EVEG(NCH),   ESNO(NCH),
+       _RL   EINT(NCH),    ESOI(NCH),   EVEG(NCH),   ESNO(NCH),
       &     STRDG1(NCH),  STRDG2(NCH), STRDG3(NCH), STRDG4(NCH),
       &     STRDG5(NCH),  STRDG6(NCH), STRDG7(NCH), STRDG8(NCH),
       &     STRDG9(NCH),
-Line 49 
 C****
+Line 50 
 C****
       &     RUNSRF(NCH),  FWSOIL(NCH)
  C****
        INTEGER ChNo
-       REAL     SWET(nch,NLAY),   VGPSAT(NTYPS),   VGCSAT (NTYPS),
+       _RL     SWET(nch,NLAY),   VGPSAT(NTYPS),   VGCSAT (NTYPS),
       &       VGZDEP(NLAY,NTYPS),   VGSLOP(NTYPS),            DELTC,
       &                    DELEA,    VGPH1(NTYPS),     VGPH2(NTYPS),
       &            VGRPLN(NTYPS),    CSOIL0(NTYPS),          WSOI12,
-Line 57 
 C****
+Line 58 
 C****
       &            DELZ23(NTYPS)
  C****
-       REAL PHLAY(nch,NLAY), AKLAY(nch,NLAY), SWET12(nch),
+       _RL PHLAY(nch,NLAY), AKLAY(nch,NLAY), SWET12(nch),
       &      CSOIL(nch),
       &       RCUN(nch),     VPDSTR(nch),   ESATTX(nch),
       &       VPDSTX(nch),   VGBEEX(nch)
-       REAL EMAXRT(nch), VGWMAX(NLAY,NTYPS), FTEMP(nch),
+       _RL EMAXRT(nch), VGWMAX(NLAY,NTYPS), FTEMP(nch),
       &        PHR(nch),      SOILCO(nch),    RC(nch),
       &        EAX(nch),          TX(nch),   RCX(nch),
       &     DRCDTC(nch),      SATCAP(nch),   PAR(nch),
       &       PDIR(nch),       DUMMY(nch)
-       REAL FTEMPX(nch),  DRCDEA(nch), VGPSAX(nch), VGCSAX(nch),
+       _RL FTEMPX(nch),  DRCDEA(nch), VGPSAX(nch), VGCSAX(nch),
       &                VGZDEX(NLAY,nch), VGSLOX(nch), VGPH1X(nch),
       &     VGPH2X(nch),  VGRPLX(nch)
-       REAL  DEDEA(nch),  DHSDEA(nch),     EM(nch), ESATTC(nch),
+       _RL  DEDEA(nch),  DHSDEA(nch),     EM(nch), ESATTC(nch),
       &     DESDTC(nch),      EA(nch),     RA(nch),   ALHX(nch),
       &     WETEQ1(nch), WETEQ2(nch),
       &        RX1(nch),     RX2(nch), SNWFRC(nch), POTFRC(nch),
       &     ESNFRC(nch),  EIRFRC(nch),   FCAN(nch),  EPFRC,
       &     DEFCIT, EADJST, RTBS
-       real cmpbug
+       _RL cmpbug
  C****
        DATA VGWMAX   /8.4,    621.6,      840.0,
-Line 185 
 C**** Pre-process input arrays as necess
+Line 186 
 C**** Pre-process input arrays as necess
        DO 100 ChNo = 1, NCH
-       DEDQA(CHNO)  = MAX(  DEDQA(CHNO), 500./ALHE )
+       DEDQA(CHNO)  = MAX(  DEDQA(CHNO), 500. _d 0/ALHE )
-       DEDTC(CHNO)  = MAX(  DEDTC(CHNO),   0. )
+       DEDTC(CHNO)  = MAX(  DEDTC(CHNO),   0. _d 0)
-       DHSDQA(CHNO) = MAX( DHSDQA(CHNO),   0. )
+       DHSDQA(CHNO) = MAX( DHSDQA(CHNO),   0. _d 0)
-       DHSDTC(CHNO) = MAX( DHSDTC(CHNO), -10. )
+       DHSDTC(CHNO) = MAX( DHSDTC(CHNO), -10. _d 0)
        EM(CHNO)     = QM(CHNO) * PSUR(CHNO) / EPSILON
        EA(CHNO)     = QA(CHNO) * PSUR(CHNO) / EPSILON
-Line 204 
 C      ETURB(CHNO)  = ETURB(CHNO) / ALHX
+Line 205 
 C      ETURB(CHNO)  = ETURB(CHNO) / ALHX
        RA(CHNO)     = ONE / ( CD(CHNO) * UM(CHNO) )
        SATCAP(ChNo) = 0.1 * ZLAI(ChNo)
        CSOIL(CHNO)  = CSOIL0(ITYP(ChNo))
-       SWET(ChNo,SFCLY ) = max( min(SWET1(ChNo),1.), 0.)
+       SWET(ChNo,SFCLY ) = max( min(SWET1(ChNo),1. _d 0), 0. _d 0)
-       SWET(ChNo,ROOTLY) = max( min(SWET2(ChNo),1.), 0.)
+       SWET(ChNo,ROOTLY) = max( min(SWET2(ChNo),1. _d 0), 0. _d 0)
-       SWET(ChNo,RECHLY) = max( min(SWET3(ChNo),1.), 0.)
+       SWET(ChNo,RECHLY) = max( min(SWET3(ChNo),1. _d 0), 0. _d 0)
-       CAPAC(CHNO)       = max( min(CAPAC(ChNo),SATCAP(CHNO)), 0.)
+       CAPAC(CHNO)       = max( min(CAPAC(ChNo),SATCAP(CHNO)), 0. _d 0)
  C****
        SNWFRC(CHNO) = SNOW(CHNO) / ( SNOW(CHNO) + SNWMID(ITYP(CHNO)) )
-       FCAN(CHNO) = MIN( 1., MAX(0.,CAPAC(ChNo)/SATCAP(ChNo)) )
+       FCAN(CHNO) = MIN( 1. _d 0, MAX(0. _d 0,CAPAC(ChNo)/SATCAP(ChNo)) )
        POTFRC(CHNO)=1.-(1.-SNWFRC(CHNO))*(1.-FCAN(CHNO))
-Line 230 
 C      SMELT(CHNO)  = 0.
+Line 231 
 C      SMELT(CHNO)  = 0.
        FWSOIL(CHNO) = 0.
  CONTINUE
  C****
  C**** - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
  C****       STEP 1: COMPUTE EFFECTIVE RESISTANCE RC FOR ENERGY BALANCE.
 Line 363 
 C****         EIRFRC=fraction of total e
        ALHX(CHNO) = (1.-ESNFRC(CHNO))*ALHE + ESNFRC(CHNO)*ALHS
 CONTINUE
        CALL FLUXES (
       I                NCH,   ITYP, DTSTEP, ESATTC, DESDTC,   ALHX,
       I              ETURB,  DEDEA,  DEDTC, HSTURB, DHSDEA, DHSDTC,
-Line 514 
 C****
+Line 513 
 C****
  C****
        INTEGER NCH
        INTEGER ITYP(NCH), ChNo
-       REAL  TRAINL(NCH), TRAINC(NCH), TSNOW(NCH),  SATCAP(NCH),
+       _RL  TRAINL(NCH), TRAINC(NCH), TSNOW(NCH),  SATCAP(NCH),
       &      WMAX(NLAY,NTYPS), TC(NCH),   CSOIL(NCH),   CAPAC(NCH),
       &      SNOW(NCH), SWET1(NCH),  RUNOFF(NCH),  SMELT(NCH),
       &     RUNSRF(NCH), FWSOIL(NCH), WETEQ1(NCH), WETINT
-       REAL DTSTEP, SNOWM, WATADD, CAVAIL, THRUC, WRUNC, WRUNL,
+       _RL DTSTEP, SNOWM, WATADD, CAVAIL, THRUC, WRUNC, WRUNL,
       &     TIMFRL, TIMFRC, FWETL, FWETC, THRU1, THRU2, THRUL, XTCORR,
       &     WETFRC
  C****
-Line 534 
 C**** Add to snow cover.  Melt snow if n
+Line 533 
 C**** Add to snow cover.  Melt snow if n
        SNOWM = 0.
        IF( SNOW(CHNO).GT.0. .AND. TC(CHNO).GT.TF ) THEN
          SNOWM = MIN( SNOW(ChNo),
-      &       MAX( 0., (TC(ChNo)-TF)*CSOIL(ChNo)/ALHM )  )
+      &       MAX( 0. _d 0, (TC(ChNo)-TF)*CSOIL(ChNo)/ALHM )  )
          IF( SNOWM .EQ. SNOW(CHNO) ) THEN
              TC(ChNo) = TC(ChNo) - SNOWM * ALHM / CSOIL(ChNo)
              SNOW(CHNO)=0.
-Line 555 
 C**** (Time scale TIMFRL for large scale
+Line 554 
 C**** (Time scale TIMFRL for large scale
  C**** to reflect the effective loss of "position memory" when storm
  C**** covers entire grid square.)
-       XTCORR= (1.-TIMFRL) * MIN( 1.,(CAPAC(CHNO)/SATCAP(CHNO))/FWETL )
+       XTCORR= (1.-TIMFRL) * MIN( 1. _d 0,(CAPAC(CHNO)/SATCAP(CHNO))/FWETL )
  C****
  C**** Fill interception reservoir with precipitation.
-Line 589 
 C****
+Line 588 
 C****
  C**** Determine XTCORR, the fraction of a storm that falls on a previously
  C**** wet surface due to the time correlation of precipitation position.
-       XTCORR= (1.-TIMFRC) * MIN( 1.,(CAPAC(CHNO)/SATCAP(CHNO))/FWETC )
+       XTCORR= (1.-TIMFRC) * MIN( 1. _d 0,(CAPAC(CHNO)/SATCAP(CHNO))/FWETC )
  C****
  C**** Fill interception reservoir with precipitation.
-Line 691 
 C****
+Line 690 
 C****
        INTEGER NCH
        INTEGER ITYP(NCH), ChNo
-       REAL  SUNANG(NCH),   PDIR(NCH),   PAR(NCH),   ZLAI(NCH),
+       _RL  SUNANG(NCH),   PDIR(NCH),   PAR(NCH),   ZLAI(NCH),
       &      SQSCAT(NCH),  GREEN(NCH),  RCUN(NCH)
-       REAL   VGCHIL(NTYPS), VGZMEW(NTYPS),
+       _RL   VGCHIL(NTYPS), VGZMEW(NTYPS),
       &       VGRST1(NTYPS), VGRST2(NTYPS), VGRST3(NTYPS)
-       REAL            RHO4,         EXTK1,         EXTK2,
+       _RL            RHO4,         EXTK1,         EXTK2,
       &               RCINV,         GAMMA,          EKAT,    DUM1,
       &                DUM2,          DUM3,            AA,      BB,
       &                  ZK,            CC
-Line 743 
 C**** (Note: CHIL is constrained to be >
+Line 742 
 C**** (Note: CHIL is constrained to be >
  C**** Bound extinction coefficient by 50./ZLAI:
-       ZK =     PDIR(ChNo) *MIN( EXTK1, 50./ZLAI(ChNo) ) +
+       ZK =     PDIR(ChNo) *MIN( EXTK1, 50. _d 0/ZLAI(ChNo) ) +
-      &    (ONE-PDIR(ChNo))*MIN( EXTK2, 50./ZLAI(ChNo) )
+      &    (ONE-PDIR(ChNo))*MIN( EXTK2, 50. _d 0/ZLAI(ChNo) )
  C**** Now compute unstressed canopy resistance:
-Line 789 
 C****
+Line 788 
 C****
        INTEGER NCH
        INTEGER ITYP(NCH), ChNo
-       REAL    WET(NCH),    PHR(NCH), SOILCO(NCH), VGPSAX(NCH),
+       _RL    WET(NCH),    PHR(NCH), SOILCO(NCH), VGPSAX(NCH),
       &     VGCSAX(NCH), VGBEEX(NCH),  DELZ(NTYPS),  WETEQ(NCH),
       &     WEXPB,  WET0,  PHEQ
        DO 100 ChNo = 1, NCH
-       WET0  = MAX(WET(CHNO),0.01)
+       WET0  = MAX(WET(CHNO),0.01 _d 0)
        WEXPB = WET0**VGBEEX(ChNo)
        PHR(ChNo)    = VGPSAX(ChNo) / WEXPB
-Line 839 
 C****
+Line 838 
 C****
        INTEGER NCH
        INTEGER ITYP(NCH), ChNo
-       REAL       DTSTEP, ESATTC(NCH), DESDTC(NCH),   ALHX(NCH),
+       _RL       DTSTEP, ESATTC(NCH), DESDTC(NCH),   ALHX(NCH),
       &       ETURB(NCH),  DEDEA(NCH),  DEDTC(NCH),
       &      HSTURB(NCH), DHSDEA(NCH), DHSDTC(NCH),
       &          RC(NCH), DRCDEA(NCH), DRCDTC(NCH),
-Line 850 
 C****
+Line 849 
 C****
       &               SWET1(NCH,NLAY),   SNOW(NCH),
       &      RUNOFF(NCH),   EVAP(NCH), SHFLUX(NCH),  SMELT(NCH),
       &       HLWUP(NCH),   BOMB(NCH)
-       REAL STRDG1(NCH),STRDG2(NCH),STRDG3(NCH),STRDG4(NCH)
+       _RL STRDG1(NCH),STRDG2(NCH),STRDG3(NCH),STRDG4(NCH)
-       REAL STRDG5(NCH),STRDG6(NCH),STRDG7(NCH),STRDG8(NCH)
+       _RL STRDG5(NCH),STRDG6(NCH),STRDG7(NCH),STRDG8(NCH)
-       REAL STRDG9(NCH)
+       _RL STRDG9(NCH)
-       REAL     HLWTC,  CDEEPS,      Q0,  RHOAIR,   CONST,  DHLWTC,
+       _RL     HLWTC,  CDEEPS,      Q0,  RHOAIR,   CONST,  DHLWTC,
       &        EPLANT,     A11,     A12,     A21,     A22,      F0,
       &           DEA,     DTC,  SNLEFT,     Q0X,  Q0SNOW,
       &         EANEW,  ESATNW,  EHARMN, DETERM, DENOM
-       LOGICAL DEBUG, CHOKE
+       LOGICAL CHOKE
-       DATA DEBUG /.FALSE./
+       _RL deepfac(ntyps)
-       real deepfac(ntyps)
        DATA deepfac /1.,1.,1.,1.,1.,1.,1.,1.,1.,1./
  C****
  C**** -------------------------------------------------------------------
-Line 876 
 C****
+Line 874 
 C****
  C****
  C**** Compute matrix elements A11, A22, AND Q0 (energy balance equation).
  C****
        A11 = CSOIL(ChNo)/DTSTEP +
       &        DHLWTC +
       &        DHSDTC(ChNo) +
-Line 920 
 C****            or outweighs any decrea
+Line 919 
 C****            or outweighs any decrea
  C****
            A21 =  -DEDTC(ChNo)*RC(ChNo) +
-      &      max(0., CONST*DESDTC(ChNo) - EHARMN*DRCDTC(ChNo) )
+      &      max(0. _d 0, CONST*DESDTC(ChNo) - EHARMN*DRCDTC(ChNo) )
            A22 = -( RC(ChNo)*DEDEA(ChNo) +
-      &               max( 0., CONST + EHARMN*DRCDEA(ChNo) )   )
+      &               max( 0. _d 0, CONST + EHARMN*DRCDEA(ChNo) )   )
            F0 = RC(ChNo) * (ETURB(ChNo) - EPLANT)
-           DETERM = MIN( A12*A21/(A11*A22) - 1., -0.1 )
+           DETERM = MIN( A12*A21/(A11*A22) - 1., -0.1 _d 0)
            DEA = ( Q0*A21 - A11*F0 ) / ( DETERM * A11*A22 )
            DTC = ( Q0 - A12*DEA ) / A11
-Line 1108 
 c  CHANGED THIS: deep layer 2 times deep
+Line 1107 
 c  CHANGED THIS: deep layer 2 times deep
  C**** Make sure EA remains positive
-       EA(CHNO) = MAX(EA(CHNO), 0.0)
+       EA(CHNO) = MAX(EA(CHNO), 0.0 _d 0)
 CONTINUE
-Line 1138 
 C****
+Line 1137 
 C****
  C****
        INTEGER NCH
        INTEGER ITYP(NCH), ChNo
-       REAL    ESATTC(NCH), EA(NCH),  VPDSTR(NCH)
+       _RL    ESATTC(NCH), EA(NCH),  VPDSTR(NCH)
-       REAL  VGDFAC(NTYPS)
+ c     _RL  VGDFAC(NTYPS)
  C****
-       DATA VGDFAC /   .0273,    .0357,    .0310,    .0238,
+ c     DATA VGDFAC /   .0273,    .0357,    .0310,    .0238,
-               .0275,    .0275,       0.,       0.,
+ c    5                .0275,    .0275,       0.,       0.,
-                  0.,      0. /
+ c    9                   0.,      0. /
  C****
  C**** -----------------------------------------------------------------
        DO 100 ChNo = 1, NCH
  C****
  c      VPDSTR(ChNo) = 1. - (ESATTC(ChNo)-EA(ChNo)) * VGDFAC(ITYP(ChNo))
- c      VPDSTR (ChNo) = MIN( 1., MAX( VPDSTR(ChNo), 1.E-10 ) )
+ c      VPDSTR (ChNo) = MIN( 1. _d 0, MAX( VPDSTR(ChNo), 1. _d -10 ) )
        VPDSTR(CHNO) = 1.
  C****
  CONTINUE
-Line 1179 
 C****
+Line 1178 
 C****
  C****
        INTEGER NCH
        INTEGER ITYP(NCH), ChNo, TypPtr
-       REAL TC(NCH), FTEMP(NCH)
+       _RL TC(NCH), FTEMP(NCH)
-       REAL    VGTLL(MemFac*NTYPS),    VGTU(MemFac*NTYPS),
+       _RL    VGTLL(MemFac*NTYPS),    VGTU(MemFac*NTYPS),
       &       VGTCF1(MemFac*NTYPS),  VGTCF2(MemFac*NTYPS),
       &       VGTCF3(MemFac*NTYPS)
  C****
-Line 1218 
 C****
+Line 1217 
 C****
       &                      VGTCF3(TypPtr) )
        IF ( TC(ChNo) .LE. VGTLL(TypPtr) .OR. TC(ChNo) .GE. VGTU(TypPtr) )
       &      FTEMP (ChNo) = 1.E-10
-       FTEMP(CHNO) = MIN( 1., MAX( FTEMP(ChNo), 1.E-10 ) )
+       FTEMP(CHNO) = MIN( 1. _d 0, MAX( FTEMP(ChNo), 1. _d -10 ) )
  C****
  CONTINUE
  C****
-Line 1249 
 C****
+Line 1248 
 C****
  C****
        INTEGER NCH
        INTEGER ITYP(NCH), ChNo
-       REAL  ESATTC(NCH),      EA(NCH),     PHR(NCH),  SOILCO(NCH),
+       _RL  ESATTC(NCH),      EA(NCH),     PHR(NCH),  SOILCO(NCH),
       &        RCUN(NCH),  VPDSTR(NCH),   FTEMP(NCH),      TC(NCH),
       &        PSUR(NCH),      Z2(NCH),  RSOIL1(NCH),  RSOIL2(NCH),
       &      VGPH1X(NCH),  VGPH2X(NCH),  VGRPLX(NCH),      RC(NCH)
-       REAL   RCUNTD,  RHOAIR,   CONST,     DEF,     D12,     DR2,
+       _RL   RCUNTD,  RHOAIR,   CONST,     DEF,     D12,     DR2,
       &        RSOIL,      R0,    EEST,  RSTFAC
  C****
  C**** -----------------------------------------------------------------
-Line 1272 
 C****
+Line 1271 
 C****
        R0 = ( VGRPLX(ChNo) + RSOIL ) / RHOW
        EEST = DEF*DR2 / ( RCUNTD*D12 + DEF*R0 )
        RSTFAC = ( DR2 - R0*EEST ) / D12
-       RSTFAC = MIN( 1., MAX( 0.001, RSTFAC ) )
+       RSTFAC = MIN( 1. _d 0, MAX( 0.001 _d 0, RSTFAC ) )
        RC(ChNo) = RCUNTD / RSTFAC
  C****
  CONTINUE
-Line 1297 
 C****
+Line 1296 
 C****
        IMPLICIT NONE
        INTEGER NCH
        INTEGER ChNo
-       REAL     UM(NCH),   U2FAC(NCH),     Z2(NCH),    RDC(NCH),
+       _RL     UM(NCH),   U2FAC(NCH),     Z2(NCH),    RDC(NCH),
       &        WET(NCH),  ESATTC(NCH),     EA(NCH),
       &         RC(NCH),     RX1(NCH),    RX2(NCH)
-       REAL  U2, RSURF, HESAT
+       _RL  U2, RSURF, HESAT
  C****
  C**** -----------------------------------------------------------------
-Line 1312 
 c      RSURF = RDC(ChNo) / U2 + 30. / (1
+Line 1311 
 c      RSURF = RDC(ChNo) / U2 + 30. / (1
  C**** Account for subsaturated humidity at soil surface:
  C****
-       HESAT = ESATTC(CHNO) * MIN( 1., WET(CHNO)*2. )
+       HESAT = ESATTC(CHNO) * MIN( 1. _d 0, WET(CHNO)*2. _d 0)
        IF( EA(CHNO) .LT. HESAT ) THEN
            RSURF=RSURF*( 1. + (ESATTC(CHNO)-HESAT)/(HESAT-EA(CHNO)) )
          ELSE
-Line 1352 
 C****
+Line 1351 
 C****
        IMPLICIT NONE
        INTEGER NCH
        INTEGER ChNo
-       REAL CAPAC(NCH), SNOW(NCH), SATCAP(NCH), RA(NCH), ETURB(NCH),
+       _RL CAPAC(NCH), SNOW(NCH), SATCAP(NCH), RA(NCH), ETURB(NCH),
       &          RC(NCH), SNWFRC(NCH), POTFRC(NCH)
-       REAL ETCRIT,RAMPFC
+       _RL ETCRIT,RAMPFC
  C**** (Note: ETCRIT arbitrarily set to ~-5 W/m2, or -2.e-6 mm/sec.)
        DATA ETCRIT/ -2.E-6 /
-Line 1419 
 C****
+Line 1418 
 C****
  C****
        INTEGER NCH
        INTEGER ITYP(NCH), ChNo
-       REAL    EVAP(NCH),  SATCAP(NCH),  VGWMAX(NLAY,NTYPS),
+       _RL    EVAP(NCH),  SATCAP(NCH),  VGWMAX(NLAY,NTYPS),
       &          TC(NCH),      RA(NCH),             RC(NCH),
       &       CAPAC(NCH),    SNOW(NCH),    SWET(nch,NLAY),
       &      RUNOFF(NCH),     RX1(NCH),
       &         RX2(NCH),  RUNSRF(NCH),  FWSOIL(NCH),
       &      ESNFRC(NCH),  EIRFRC(NCH)
-       REAL EINT(NCH), ESOI(NCH), EVEG(NCH), ESNO(NCH)
+       _RL EINT(NCH), ESOI(NCH), EVEG(NCH), ESNO(NCH)
-       REAL  DTSTEP, EGRO, FWS, THRU, DEWRUN,
+       _RL  DTSTEP, EGRO, FWS, THRU, DEWRUN,
       &     WTOTAL,WLAY1,WLAY2,ELAY1,ELAY2,EGROI
  C****
  C**** -----------------------------------------------------------------
-Line 1515 
 C****
+Line 1514 
 C****
        SWET(ChNo,ROOTLY) = (WLAY2 - ELAY2) / VGWMAX(ROOTLY,ITYP(CHNO))
  C****
  C**** Ensure against numerical precision problems:
-       SWET(ChNo,SFCLY)  = MIN( 1., MAX( 0., SWET(ChNo,SFCLY) ) )
+       SWET(ChNo,SFCLY)  = MIN( 1. _d 0, MAX( 0. _d 0, SWET(ChNo,SFCLY) ) )
-       SWET(ChNo,ROOTLY) = MIN( 1., MAX( 0., SWET(ChNo,ROOTLY) ) )
+       SWET(ChNo,ROOTLY) = MIN( 1. _d 0, MAX( 0. _d 0, SWET(ChNo,ROOTLY) ) )
  C****
  C****
  C**** -------------------------------------------------
-Line 1569 
 C****
+Line 1568 
 C****
  C****
        INTEGER NCH
        INTEGER ITYP(NCH), ChNo
-       REAL  VGSLOX(NCH), RUNOFF(NCH), GDRAIN(NCH)
+       _RL  VGSLOX(NCH), RUNOFF(NCH), GDRAIN(NCH)
-       REAL  ZDEP12,   AKAVE,  GWFLUX,  ZDEP23,  HALFMX,  DHDZ,
+       _RL  ZDEP12,   AKAVE,  GWFLUX,  ZDEP23,  HALFMX,  DHDZ,
       &    FAREA, TFM2, FRAMP
-       REAL   WSMAX(NLAY,NTYPS),   PHLAY(nch,NLAY),
+       _RL   WSMAX(NLAY,NTYPS),   PHLAY(nch,NLAY),
       &       AKLAY(nch,NLAY),             TC(NCH),
       &                  DTSTEP,    SWET(nch,NLAY),
       &      VGZDEX(NLAY,nch),  WETEQ1(NCH),  WETEQ2(NCH),
-Line 1619 
 C****
+Line 1618 
 C****
  C**** Prevent diffusion when ground is frozen (INCLUDE RAMPING):
        TFM2=TF-2.
        FRAMP=(TC(CHNO)-TFM2)/2.
-       FRAMP=MIN(1., MAX(0.,FRAMP) )
+       FRAMP=MIN(1. _d 0, MAX(0. _d 0,FRAMP) )
        GWFLUX=GWFLUX*FRAMP
  C****
  C**** Update water contents
-Line 1658 
 C**** Test for limits on water holding c
+Line 1657 
 C**** Test for limits on water holding c
  C****
  C**** Prevent diffusion when ground is frozen (INCLUDE RAMPING):
        FRAMP=(TC(CHNO)-TFM2)/2.
-       FRAMP=MIN(1., MAX(0.,FRAMP) )
+       FRAMP=MIN(1. _d 0, MAX(0. _d 0,FRAMP) )
        GWFLUX=GWFLUX*FRAMP
  C****
  C**** Update water contents
-Line 1676 
 C****
+Line 1675 
 C****
  C**** Prevent diffusion when ground is frozen (INCLUDE RAMPING):
        FRAMP=(TC(CHNO)-TFM2)/2.
-       FRAMP=MIN(1., MAX(0.,FRAMP) )
+       FRAMP=MIN(1. _d 0, MAX(0. _d 0,FRAMP) )
        GWFLUX=GWFLUX*FRAMP
        GWFLUX = MIN( GWFLUX,

 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.6
 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.6
-Removed from v.1.2
+Added in v.1.6

	ViewVC Help
Powered by ViewVC 1.1.22