| File |
Line number |
Procedure |
Code |
|
./pkg/fizhi/fizhi_lsm.F |
1008 |
FLUXES |
& (Q0 + ALHX(CHNO)*(ETURB(ChNo)-EMAXRT(CHNO)) - DHSDEA(CHNO)*DEA) |
|
./pkg/fizhi/fizhi_lsm.F |
1009 |
FLUXES |
& / ( A11 - ALHX(CHNO)*DEDTC(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1011 |
FLUXES |
STRDG1(ChNo)=DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1012 |
FLUXES |
STRDG2(ChNo)=DEA |
|
./pkg/fizhi/fizhi_lsm.F |
1013 |
FLUXES |
STRDG8(ChNo)=0. |
|
./pkg/fizhi/fizhi_lsm.F |
1016 |
FLUXES |
SHFLUX(ChNo) = HSTURB(ChNo) + DHSDEA(ChNo)*DEA + |
|
./pkg/fizhi/fizhi_lsm.F |
1017 |
FLUXES |
& DHSDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1041 |
FLUXES |
DTC = ( Q0 + ALHX(CHNO)*ETURB(ChNo) - DHSDEA(CHNO)*DEA ) / |
|
./pkg/fizhi/fizhi_lsm.F |
1042 |
FLUXES |
& ( A11 - ALHX(CHNO)*DEDTC(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1044 |
FLUXES |
STRDG1(ChNo)=DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1045 |
FLUXES |
STRDG2(ChNo)=DEA |
|
./pkg/fizhi/fizhi_lsm.F |
1046 |
FLUXES |
STRDG8(ChNo)=0. |
|
./pkg/fizhi/fizhi_lsm.F |
1049 |
FLUXES |
SHFLUX(ChNo) = HSTURB(ChNo) + DHSDEA(ChNo)*DEA + |
|
./pkg/fizhi/fizhi_lsm.F |
1050 |
FLUXES |
& DHSDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1055 |
FLUXES |
DTC = ( Q0X + ALHX(CHNO)*ETURB(ChNo) - DHSDEA(CHNO)*DEA ) / |
|
./pkg/fizhi/fizhi_lsm.F |
1056 |
FLUXES |
& ( A11 - ALHX(CHNO)*DEDTC(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1062 |
FLUXES |
SHFLUX(ChNo) = HSTURB(ChNo) + DHSDEA(ChNo)*DEA + |
|
./pkg/fizhi/fizhi_lsm.F |
1063 |
FLUXES |
& DHSDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1075 |
FLUXES |
STRDG1(ChNo)=DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1076 |
FLUXES |
STRDG2(ChNo)=DEA |
|
./pkg/fizhi/fizhi_lsm.F |
1077 |
FLUXES |
STRDG8(ChNo)=-A12*DEA/A11 |
|
./pkg/fizhi/fizhi_lsm.F |
1079 |
FLUXES |
EVAP(ChNo) = ETURB(ChNo) + DEDEA(ChNo)*DEA + DEDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1080 |
FLUXES |
SHFLUX(ChNo) = HSTURB(ChNo) + DHSDEA(ChNo)*DEA + |
|
./pkg/fizhi/fizhi_lsm.F |
1081 |
FLUXES |
& DHSDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1092 |
FLUXES |
EVAP(ChNo) = ETURB(ChNo) + DEDEA(ChNo)*DEA + DEDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1093 |
FLUXES |
SHFLUX(ChNo) = HSTURB(ChNo) + DHSDEA(ChNo)*DEA + |
|
./pkg/fizhi/fizhi_lsm.F |
1094 |
FLUXES |
& DHSDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1100 |
FLUXES |
TC(ChNo) = TC(ChNo) + DTC |
|
./pkg/fizhi/fizhi_lsm.F |
1101 |
FLUXES |
EA(ChNo) = EA(ChNo) + DEA |
|
./pkg/fizhi/fizhi_lsm.F |
1102 |
FLUXES |
TD(ChNo) = TD(ChNo) + |
|
./pkg/fizhi/fizhi_lsm.F |
1104 |
FLUXES |
& DTSTEP * CDEEPS * (TC(ChNo) - TD(ChNo)) / |
|
./pkg/fizhi/fizhi_lsm.F |
1105 |
FLUXES |
. (CSOIL(ChNo)*67.73*deepfac(ityp(ChNo))) |
|
./pkg/fizhi/fizhi_lsm.F |
1141 |
VPDFAC |
INTEGER ITYP(NCH), ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
1151 |
VPDFAC |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
1182 |
TMPFAC |
INTEGER ITYP(NCH), ChNo, TypPtr |
|
./pkg/fizhi/fizhi_lsm.F |
1212 |
TMPFAC |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
1214 |
TMPFAC |
TypPtr = MOD(ChNo,MemFac) + (ITYP(ChNo)-1)*MemFac + 1 |
|
./pkg/fizhi/fizhi_lsm.F |
1215 |
TMPFAC |
FTEMP(ChNo) = (TC(ChNo) - VGTLL(TypPtr)) * |
|
./pkg/fizhi/fizhi_lsm.F |
1216 |
TMPFAC |
& (TC(ChNo) - VGTU(TypPtr)) * |
|
./pkg/fizhi/fizhi_lsm.F |
1217 |
TMPFAC |
& ( VGTCF1(TypPtr)*TC(ChNo)*TC(ChNo) + |
|
./pkg/fizhi/fizhi_lsm.F |
1218 |
TMPFAC |
& VGTCF2(TypPtr)*TC(ChNo) + |
|
./pkg/fizhi/fizhi_lsm.F |
1220 |
TMPFAC |
IF ( TC(ChNo) .LE. VGTLL(TypPtr) .OR. TC(ChNo) .GE. VGTU(TypPtr) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1221 |
TMPFAC |
& FTEMP (ChNo) = 1.E-10 |
|
./pkg/fizhi/fizhi_lsm.F |
1222 |
TMPFAC |
FTEMP(CHNO) = MIN( 1. _d 0, MAX( FTEMP(ChNo), 1. _d -10 ) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1252 |
SMFAC |
INTEGER ITYP(NCH), ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
1262 |
SMFAC |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
1264 |
SMFAC |
RCUNTD = RCUN(ChNo) / ( VPDSTR(ChNo)*FTEMP(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1265 |
SMFAC |
RHOAIR = PSUR(ChNo) * 100. / ( RGAS*TC(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1268 |
SMFAC |
CONST = RHOAIR * EPSILON / PSUR(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
1269 |
SMFAC |
DEF = ( ESATTC(ChNo) - EA(ChNo) ) * CONST |
|
./pkg/fizhi/fizhi_lsm.F |
1270 |
SMFAC |
D12 = VGPH1X(ChNo) - VGPH2X(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
1271 |
SMFAC |
DR2 = PHR(ChNo) - Z2(ChNo) - VGPH2X(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
1272 |
SMFAC |
RSOIL = RSOIL1(ChNo) + RSOIL2(ChNo)/SOILCO(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
1273 |
SMFAC |
R0 = ( VGRPLX(ChNo) + RSOIL ) / RHOW |
|
./pkg/fizhi/fizhi_lsm.F |
1277 |
SMFAC |
RC(ChNo) = RCUNTD / RSTFAC |
|
./pkg/fizhi/fizhi_lsm.F |
1300 |
RSURFP |
INTEGER ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
1308 |
RSURFP |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
1310 |
RSURFP |
U2 = UM(ChNo) * U2FAC(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
1312 |
RSURFP |
RSURF = RDC(ChNo) / U2 + 26. + 6. / (1.E-20 + WET(ChNo))**2 |
|
./pkg/fizhi/fizhi_lsm.F |
1327 |
RSURFP |
RC(ChNo) = RC(CHNO) * RSURF / ( RC(ChNo) + RSURF ) |
|
./pkg/fizhi/fizhi_lsm.F |
1355 |
RCANOP |
INTEGER ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
1365 |
RCANOP |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
1372 |
RCANOP |
RC(ChNo)=RC(ChNo)*(1.-POTFRC(CHNO))/ |
|
./pkg/fizhi/fizhi_lsm.F |
1373 |
RCANOP |
& ( 1.+POTFRC(CHNO)*RC(ChNo)/RA(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1381 |
RCANOP |
RC(ChNo)=RC(ChNo)*(1.-SNWFRC(CHNO))/ |
|
./pkg/fizhi/fizhi_lsm.F |
1382 |
RCANOP |
& ( 1.+SNWFRC(CHNO)*RC(ChNo)/RA(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1391 |
RCANOP |
IF ( RAMPFC .GT. 1. ) RC(ChNo) = 0. |
|
./pkg/fizhi/fizhi_lsm.F |
1422 |
WUPDAT |
INTEGER ITYP(NCH), ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
1435 |
WUPDAT |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
1439 |
WUPDAT |
WLAY1 = SWET(ChNo,SFCLY) * VGWMAX(SFCLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1440 |
WUPDAT |
WLAY2 = SWET(ChNo,ROOTLY) * VGWMAX(ROOTLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1513 |
WUPDAT |
SNOW(ChNo) = SNOW(ChNo) - ESNO(CHNO) |
|
./pkg/fizhi/fizhi_lsm.F |
1514 |
WUPDAT |
CAPAC(ChNo) = CAPAC(ChNo) - EINT(CHNO) |
|
./pkg/fizhi/fizhi_lsm.F |
1515 |
WUPDAT |
SWET(ChNo,SFCLY) = (WLAY1 - ELAY1) / VGWMAX(SFCLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1516 |
WUPDAT |
SWET(ChNo,ROOTLY) = (WLAY2 - ELAY2) / VGWMAX(ROOTLY,ITYP(CHNO)) |
|
./pkg/fizhi/fizhi_lsm.F |
1519 |
WUPDAT |
SWET(ChNo,SFCLY) = MIN( 1. _d 0, MAX( 0. _d 0, SWET(ChNo,SFCLY) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1521 |
WUPDAT |
SWET(ChNo,ROOTLY) = MIN( 1. _d 0, MAX( 0. _d 0, SWET(ChNo,ROOTLY) |
|
./pkg/fizhi/fizhi_lsm.F |
1532 |
WUPDAT |
IF( CAPAC(ChNo) .GT. SATCAP(ChNo) ) THEN |
|
./pkg/fizhi/fizhi_lsm.F |
1533 |
WUPDAT |
THRU = CAPAC(ChNo) - SATCAP(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
1534 |
WUPDAT |
CAPAC(ChNo) = SATCAP(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
1535 |
WUPDAT |
SWET(ChNo,SFCLY) = SWET(ChNo,SFCLY) + |
|
./pkg/fizhi/fizhi_lsm.F |
1536 |
WUPDAT |
& THRU / VGWMAX(SFCLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1538 |
WUPDAT |
IF ( SWET(ChNo,SFCLY) .GT. 1. ) THEN |
|
./pkg/fizhi/fizhi_lsm.F |
1539 |
WUPDAT |
DEWRUN = ( SWET(ChNo,SFCLY) - 1. ) * VGWMAX(SFCLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1540 |
WUPDAT |
SWET(ChNo,SFCLY) = 1. |
|
./pkg/fizhi/fizhi_lsm.F |
1541 |
WUPDAT |
RUNOFF(ChNo) = RUNOFF(ChNo) + DEWRUN/DTSTEP |
|
./pkg/fizhi/fizhi_lsm.F |
1542 |
WUPDAT |
RUNSRF(ChNo) = RUNSRF(ChNo) + DEWRUN/DTSTEP |
|
./pkg/fizhi/fizhi_lsm.F |
1574 |
GWATER |
INTEGER ITYP(NCH), ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
1588 |
GWATER |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
1591 |
GWATER |
ZDEP12 = VGZDEX(SFCLY,ChNo) + VGZDEX(ROOTLY,ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
1598 |
GWATER |
DHDZ = 1. + 2.*(PHLAY(ChNo,SFCLY)-PHLAY(ChNo,ROOTLY))/ZDEP12 |
|
./pkg/fizhi/fizhi_lsm.F |
1599 |
GWATER |
AKAVE = AKLAY(ChNo,ROOTLY) |
|
./pkg/fizhi/fizhi_lsm.F |
1610 |
GWATER |
& SWET(ChNo,SFCLY) * WSMAX(SFCLY,ITYP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1611 |
GWATER |
GWFLUX = MIN( GWFLUX, WSMAX(ROOTLY,ITYP(ChNo)) - |
|
./pkg/fizhi/fizhi_lsm.F |
1612 |
GWATER |
& SWET(ChNo,ROOTLY) * WSMAX(ROOTLY,ITYP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1616 |
GWATER |
& SWET(ChNo,ROOTLY) * WSMAX(ROOTLY,ITYP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1617 |
GWATER |
GWFLUX = -MIN( ABS(GWFLUX), WSMAX(SFCLY,ITYP(ChNo)) - |
|
./pkg/fizhi/fizhi_lsm.F |
1618 |
GWATER |
& SWET(ChNo,SFCLY) * WSMAX(SFCLY,ITYP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1629 |
GWATER |
SWET(ChNo,SFCLY) = SWET(ChNo,SFCLY) - |
|
./pkg/fizhi/fizhi_lsm.F |
1630 |
GWATER |
& GWFLUX / WSMAX(SFCLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1631 |
GWATER |
SWET(ChNo,ROOTLY) = SWET(ChNo,ROOTLY) + |
|
./pkg/fizhi/fizhi_lsm.F |
1632 |
GWATER |
& GWFLUX / WSMAX(ROOTLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1636 |
GWATER |
ZDEP23 = VGZDEX(ROOTLY,ChNo) + VGZDEX(RECHLY,ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
1637 |
GWATER |
DHDZ=1. + 2. * (PHLAY(ChNo,ROOTLY) - PHLAY(ChNo,RECHLY)) / ZDEP23 |
|
./pkg/fizhi/fizhi_lsm.F |
1639 |
GWATER |
AKAVE = AKLAY(ChNo,ROOTLY) |
|
./pkg/fizhi/fizhi_lsm.F |
1641 |
GWATER |
AKAVE = AKLAY(ChNo,RECHLY) |
|
./pkg/fizhi/fizhi_lsm.F |
1650 |
GWATER |
& SWET(ChNo,ROOTLY) * WSMAX(ROOTLY,ITYP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1651 |
GWATER |
GWFLUX = MIN( GWFLUX, WSMAX(RECHLY,ITYP(ChNo)) - |
|
./pkg/fizhi/fizhi_lsm.F |
1652 |
GWATER |
& SWET(ChNo,RECHLY) * WSMAX(RECHLY,ITYP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1656 |
GWATER |
& SWET(ChNo,RECHLY) * WSMAX(RECHLY,ITYP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1657 |
GWATER |
GWFLUX = -MIN( ABS(GWFLUX), WSMAX(ROOTLY,ITYP(ChNo)) - |
|
./pkg/fizhi/fizhi_lsm.F |
1658 |
GWATER |
& SWET(ChNo,ROOTLY) * WSMAX(ROOTLY,ITYP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1668 |
GWATER |
SWET(ChNo,ROOTLY) = SWET(ChNo,ROOTLY) - |
|
./pkg/fizhi/fizhi_lsm.F |
1669 |
GWATER |
& GWFLUX / WSMAX(ROOTLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1670 |
GWATER |
SWET(ChNo,RECHLY) = SWET(ChNo,RECHLY) + |
|
./pkg/fizhi/fizhi_lsm.F |
1671 |
GWATER |
& GWFLUX / WSMAX(RECHLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1678 |
GWATER |
GWFLUX = VGSLOX(ChNo) * AKLAY(ChNo,RECHLY) * 1000. * DTSTEP |
|
./pkg/fizhi/fizhi_lsm.F |
1686 |
GWATER |
& SWET(ChNo,RECHLY) * WSMAX(RECHLY,ITYP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
1687 |
GWATER |
SWET(ChNo,RECHLY) = SWET(ChNo,RECHLY) - |
|
./pkg/fizhi/fizhi_lsm.F |
1688 |
GWATER |
& GWFLUX / WSMAX(RECHLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
1690 |
GWATER |
RUNOFF (ChNo) = RUNOFF (ChNo) + GWFLUX/DTSTEP |
|
./pkg/fizhi/fizhi_lsm.F |
170 |
TILE |
DO 50 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
172 |
TILE |
VGBEEX(ChNo) = VGBEE(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
173 |
TILE |
VGPSAX(ChNo) = VGPSAT(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
174 |
TILE |
VGCSAX(ChNo) = VGCSAT(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
175 |
TILE |
VGZDEX(SFCLY ,ChNo) = VGZDEP(SFCLY ,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
176 |
TILE |
VGZDEX(ROOTLY,ChNo) = VGZDEP(ROOTLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
177 |
TILE |
VGZDEX(RECHLY,ChNo) = VGZDEP(RECHLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
178 |
TILE |
VGSLOX(ChNo) = VGSLOP(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
179 |
TILE |
VGPH1X(ChNo) = VGPH1(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
180 |
TILE |
VGPH2X(ChNo) = VGPH2(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
181 |
TILE |
VGRPLX(ChNo) = VGRPLN(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
187 |
TILE |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
206 |
TILE |
SATCAP(ChNo) = 0.1 * ZLAI(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
207 |
TILE |
CSOIL(CHNO) = CSOIL0(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
208 |
TILE |
SWET(ChNo,SFCLY ) = max( min(SWET1(ChNo),1. _d 0), 0. _d 0) |
|
./pkg/fizhi/fizhi_lsm.F |
209 |
TILE |
SWET(ChNo,ROOTLY) = max( min(SWET2(ChNo),1. _d 0), 0. _d 0) |
|
./pkg/fizhi/fizhi_lsm.F |
210 |
TILE |
SWET(ChNo,RECHLY) = max( min(SWET3(ChNo),1. _d 0), 0. _d 0) |
|
./pkg/fizhi/fizhi_lsm.F |
211 |
TILE |
CAPAC(CHNO) = max( min(CAPAC(ChNo),SATCAP(CHNO)), 0. _d 0) |
|
./pkg/fizhi/fizhi_lsm.F |
215 |
TILE |
FCAN(CHNO) = MIN( 1. _d 0, MAX(0. _d 0,CAPAC(ChNo)/SATCAP(ChNo)) ) |
|
./pkg/fizhi/fizhi_lsm.F |
219 |
TILE |
WSOI12=SWET(ChNo,SFCLY ) * VGWMAX(SFCLY ,ITYP(ChNo)) + |
|
./pkg/fizhi/fizhi_lsm.F |
220 |
TILE |
& SWET(ChNo,ROOTLY) * VGWMAX(ROOTLY,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
221 |
TILE |
SWET12(ChNo) = WSOI12 / |
|
./pkg/fizhi/fizhi_lsm.F |
222 |
TILE |
& (VGWMAX(SFCLY,ITYP(ChNo)) + VGWMAX(ROOTLY,ITYP(ChNo))) |
|
./pkg/fizhi/fizhi_lsm.F |
296 |
TILE |
DO 120 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
297 |
TILE |
TX(ChNo) = TC(ChNo) + DELTC |
|
./pkg/fizhi/fizhi_lsm.F |
298 |
TILE |
ESATTX(ChNo) = ESATTC(ChNo) + DESDTC(CHNO) * DELTC |
|
./pkg/fizhi/fizhi_lsm.F |
299 |
TILE |
EAX(ChNo) = EA(ChNo) + DELEA |
|
./pkg/fizhi/fizhi_lsm.F |
319 |
TILE |
DO 125 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
320 |
TILE |
DRCDTC(ChNo) = (RCX(ChNo) - RC(ChNo)) / DELTC |
|
./pkg/fizhi/fizhi_lsm.F |
339 |
TILE |
DO 150 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
340 |
TILE |
DRCDEA(ChNo) = (RCX(ChNo) - RC(ChNo)) / DELEA |
|
./pkg/fizhi/fizhi_lsm.F |
455 |
TILE |
DO 2000 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
460 |
TILE |
SWET1(ChNo) = SWET(ChNo,SFCLY) |
|
./pkg/fizhi/fizhi_lsm.F |
461 |
TILE |
SWET2(ChNo) = SWET(ChNo,ROOTLY) |
|
./pkg/fizhi/fizhi_lsm.F |
462 |
TILE |
SWET3(ChNo) = SWET(ChNo,RECHLY) |
|
./pkg/fizhi/fizhi_lsm.F |
464 |
TILE |
IF(SWET1(ChNo).LT.1.E-10) SWET1(CHNO) = 0.0 |
|
./pkg/fizhi/fizhi_lsm.F |
465 |
TILE |
IF(SWET2(ChNo).LT.1.E-10) SWET2(CHNO) = 0.0 |
|
./pkg/fizhi/fizhi_lsm.F |
466 |
TILE |
IF(SWET3(ChNo).LT.1.E-10) SWET3(CHNO) = 0.0 |
|
./pkg/fizhi/fizhi_lsm.F |
467 |
TILE |
IF(CAPAC(ChNo).LT.1.E-10) CAPAC(CHNO) = 0.0 |
|
./pkg/fizhi/fizhi_lsm.F |
468 |
TILE |
IF(SNOW (ChNo).LT.1.E-10) SNOW (CHNO) = 0.0 |
|
./pkg/fizhi/fizhi_lsm.F |
469 |
TILE |
IF(RUNOFF(ChNo).LT.1.E-10) RUNOFF(CHNO) = 0.0 |
|
./pkg/fizhi/fizhi_lsm.F |
515 |
INTERC |
INTEGER ITYP(NCH), ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
529 |
INTERC |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
52 |
TILE |
INTEGER ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
532 |
INTERC |
SNOW(ChNo) = SNOW(ChNo) + TSNOW(ChNo)*DTSTEP |
|
./pkg/fizhi/fizhi_lsm.F |
535 |
INTERC |
SNOWM = MIN( SNOW(ChNo), |
|
./pkg/fizhi/fizhi_lsm.F |
536 |
INTERC |
& MAX( 0. _d 0, (TC(ChNo)-TF)*CSOIL(ChNo)/ALHM ) ) |
|
./pkg/fizhi/fizhi_lsm.F |
538 |
INTERC |
TC(ChNo) = TC(ChNo) - SNOWM * ALHM / CSOIL(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
542 |
INTERC |
SNOW(ChNo) = SNOW(ChNo) - SNOWM |
|
./pkg/fizhi/fizhi_lsm.F |
567 |
INTERC |
WATADD = TRAINL(ChNo)*DTSTEP + SMELT(CHNO)*DTSTEP |
|
./pkg/fizhi/fizhi_lsm.F |
602 |
INTERC |
WATADD = TRAINC(ChNo)*DTSTEP |
|
./pkg/fizhi/fizhi_lsm.F |
633 |
INTERC |
CAVAIL = ( 1.-WETFRC)*FWETL*WMAX(1,ITYP(ChNo))*(1-WETEQ1(CHNO)) |
|
./pkg/fizhi/fizhi_lsm.F |
636 |
INTERC |
SWET1(ChNo) = SWET1(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
637 |
INTERC |
* + THRUL * ( 1. - WETFRC) / WMAX(1,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
640 |
INTERC |
SWET1(CHNO) = SWET1(CHNO) + CAVAIL / WMAX(1,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
652 |
INTERC |
CAVAIL = (1.-WETFRC)*FWETC*WMAX(1,ITYP(ChNo))*(1-WETEQ1(CHNO)) |
|
./pkg/fizhi/fizhi_lsm.F |
655 |
INTERC |
SWET1(ChNo) = SWET1(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
656 |
INTERC |
* + THRUC * ( 1. - WETFRC) / WMAX(1,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
659 |
INTERC |
SWET1(CHNO) = SWET1(CHNO) + CAVAIL / WMAX(1,ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
662 |
INTERC |
RUNOFF(ChNo) = RUNOFF(ChNo) + (WRUNC+WRUNL)/DTSTEP |
|
./pkg/fizhi/fizhi_lsm.F |
663 |
INTERC |
RUNSRF(ChNo) = RUNSRF(ChNo) + (WRUNC+WRUNL)/DTSTEP |
|
./pkg/fizhi/fizhi_lsm.F |
693 |
RCUNST |
INTEGER ITYP(NCH), ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
730 |
RCUNST |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
735 |
RCUNST |
AA = 0.5 - (0.633 + 0.330*VGCHIL(ITYP(ChNo)))*VGCHIL(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
737 |
RCUNST |
CC = ( AA + BB*SUNANG(ChNo) ) / SUNANG(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
739 |
RCUNST |
EXTK1 = CC * SQSCAT(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
740 |
RCUNST |
EXTK2 = (ONE / VGZMEW(ITYP(ChNo))) * SQSCAT(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
742 |
RCUNST |
DUM1 = PDIR(ChNo) * CC |
|
./pkg/fizhi/fizhi_lsm.F |
743 |
RCUNST |
DUM2 = (ONE-PDIR(ChNo)) * ( BB*(ONE/3.+PIE/4.) + AA*1.5 ) |
|
./pkg/fizhi/fizhi_lsm.F |
747 |
RCUNST |
ZK = PDIR(ChNo) *MIN( EXTK1, 50. _d 0/ZLAI(ChNo) ) + |
|
./pkg/fizhi/fizhi_lsm.F |
748 |
RCUNST |
& (ONE-PDIR(ChNo))*MIN( EXTK2, 50. _d 0/ZLAI(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
752 |
RCUNST |
GAMMA = VGRST1(ITYP(ChNo)) / VGRST3(ITYP(ChNo)) + |
|
./pkg/fizhi/fizhi_lsm.F |
753 |
RCUNST |
& VGRST2(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
755 |
RCUNST |
EKAT = EXP( ZK*ZLAI(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
756 |
RCUNST |
RHO4 = GAMMA / (PAR(ChNo) * (DUM1 + DUM2)) |
|
./pkg/fizhi/fizhi_lsm.F |
758 |
RCUNST |
DUM1 = (VGRST2(ITYP(ChNo)) - GAMMA) / (GAMMA + 1.E-20) |
|
./pkg/fizhi/fizhi_lsm.F |
760 |
RCUNST |
DUM3 = ZK * VGRST3(ITYP(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
761 |
RCUNST |
RCINV = ( DUM1*LOG(DUM2) + ZK*ZLAI(ChNo) ) / DUM3 |
|
./pkg/fizhi/fizhi_lsm.F |
763 |
RCUNST |
RCUN(ChNo) = ONE / (RCINV * GREEN(ChNo) + 1.E-10) |
|
./pkg/fizhi/fizhi_lsm.F |
791 |
SOIL |
INTEGER ITYP(NCH), ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
797 |
SOIL |
DO 100 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
800 |
SOIL |
WEXPB = WET0**VGBEEX(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
802 |
SOIL |
PHR(ChNo) = VGPSAX(ChNo) / WEXPB |
|
./pkg/fizhi/fizhi_lsm.F |
803 |
SOIL |
SOILCO(ChNo) = VGCSAX(ChNo) * WEXPB * WEXPB * WET0 * WET0 * WET0 |
|
./pkg/fizhi/fizhi_lsm.F |
841 |
FLUXES |
INTEGER ITYP(NCH), ChNo |
|
./pkg/fizhi/fizhi_lsm.F |
869 |
FLUXES |
DO 200 ChNo = 1, NCH |
|
./pkg/fizhi/fizhi_lsm.F |
872 |
FLUXES |
CDEEPS = PIE * CSOIL(ChNo) * 2. / 86400. |
|
./pkg/fizhi/fizhi_lsm.F |
873 |
FLUXES |
RHOAIR = PSUR(ChNo) * 100. / (RGAS * TC(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
874 |
FLUXES |
CONST = RHOAIR * EPSILON / PSUR(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
880 |
FLUXES |
A11 = CSOIL(ChNo)/DTSTEP + |
|
./pkg/fizhi/fizhi_lsm.F |
882 |
FLUXES |
& DHSDTC(ChNo) + |
|
./pkg/fizhi/fizhi_lsm.F |
883 |
FLUXES |
& ALHX(CHNO)*DEDTC(ChNo) + |
|
./pkg/fizhi/fizhi_lsm.F |
885 |
FLUXES |
A12 = DHSDEA(ChNo) + ALHX(CHNO) * DEDEA(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
886 |
FLUXES |
Q0 = SWNET(ChNo) + |
|
./pkg/fizhi/fizhi_lsm.F |
887 |
FLUXES |
& HLWDWN(ChNo) - |
|
./pkg/fizhi/fizhi_lsm.F |
889 |
FLUXES |
& HSTURB(ChNo) - |
|
./pkg/fizhi/fizhi_lsm.F |
890 |
FLUXES |
& ALHX(CHNO) * ETURB(ChNo) - |
|
./pkg/fizhi/fizhi_lsm.F |
891 |
FLUXES |
& CDEEPS * (TC(ChNo) - TD(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
893 |
FLUXES |
STRDG3(ChNo)=Q0/A11 |
|
./pkg/fizhi/fizhi_lsm.F |
894 |
FLUXES |
STRDG4(ChNo)=(SWNET(ChNo)+HLWDWN(ChNo)-HLWTC)/A11 |
|
./pkg/fizhi/fizhi_lsm.F |
895 |
FLUXES |
STRDG5(ChNo)=HSTURB(ChNo)/A11 |
|
./pkg/fizhi/fizhi_lsm.F |
896 |
FLUXES |
STRDG6(ChNo)=ALHX(CHNO)*ETURB(ChNo)/A11 |
|
./pkg/fizhi/fizhi_lsm.F |
897 |
FLUXES |
STRDG7(ChNo)=CDEEPS*(TC(ChNo) - TD(ChNo))/A11 |
|
./pkg/fizhi/fizhi_lsm.F |
898 |
FLUXES |
STRDG9(ChNo)=A11 |
|
./pkg/fizhi/fizhi_lsm.F |
910 |
FLUXES |
EPLANT = CONST * (ESATTC(ChNo) - EA(ChNo)) / RC(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
911 |
FLUXES |
IF(EPLANT*ETURB(ChNo).GE.0.) THEN |
|
./pkg/fizhi/fizhi_lsm.F |
912 |
FLUXES |
EHARMN = 2.*EPLANT*ETURB(CHNO) / (EPLANT + ETURB(ChNo)) |
|
./pkg/fizhi/fizhi_lsm.F |
923 |
FLUXES |
A21 = -DEDTC(ChNo)*RC(ChNo) + |
|
./pkg/fizhi/fizhi_lsm.F |
924 |
FLUXES |
& max(0. _d 0, CONST*DESDTC(ChNo) - EHARMN*DRCDTC(ChNo) ) |
|
./pkg/fizhi/fizhi_lsm.F |
925 |
FLUXES |
A22 = -( RC(ChNo)*DEDEA(ChNo) + |
|
./pkg/fizhi/fizhi_lsm.F |
926 |
FLUXES |
& max( 0. _d 0, CONST + EHARMN*DRCDEA(ChNo) ) ) |
|
./pkg/fizhi/fizhi_lsm.F |
928 |
FLUXES |
F0 = RC(ChNo) * (ETURB(ChNo) - EPLANT) |
|
./pkg/fizhi/fizhi_lsm.F |
933 |
FLUXES |
STRDG1(ChNo)=DTC |
|
./pkg/fizhi/fizhi_lsm.F |
934 |
FLUXES |
STRDG2(ChNo)=DEA |
|
./pkg/fizhi/fizhi_lsm.F |
935 |
FLUXES |
STRDG8(ChNo)=-A12*DEA/A11 |
|
./pkg/fizhi/fizhi_lsm.F |
937 |
FLUXES |
EVAP(ChNo) = ETURB(ChNo) + DEDEA(ChNo)*DEA + DEDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
938 |
FLUXES |
SHFLUX(ChNo) = HSTURB(ChNo) + DHSDEA(ChNo)*DEA + |
|
./pkg/fizhi/fizhi_lsm.F |
939 |
FLUXES |
& DHSDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
943 |
FLUXES |
A21 = -DESDTC(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
945 |
FLUXES |
F0 = ESATTC(ChNo) - EA(ChNo) |
|
./pkg/fizhi/fizhi_lsm.F |
949 |
FLUXES |
STRDG1(ChNo)=DTC |
|
./pkg/fizhi/fizhi_lsm.F |
950 |
FLUXES |
STRDG2(ChNo)=DEA |
|
./pkg/fizhi/fizhi_lsm.F |
951 |
FLUXES |
STRDG8(ChNo)=-A12*DEA/A11 |
|
./pkg/fizhi/fizhi_lsm.F |
953 |
FLUXES |
EVAP(ChNo) = ETURB(ChNo) + DEDEA(ChNo)*DEA + DEDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
954 |
FLUXES |
SHFLUX(ChNo) = HSTURB(ChNo) + DHSDEA(ChNo)*DEA + |
|
./pkg/fizhi/fizhi_lsm.F |
955 |
FLUXES |
& DHSDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
987 |
FLUXES |
STRDG1(ChNo)=DTC |
|
./pkg/fizhi/fizhi_lsm.F |
988 |
FLUXES |
STRDG2(ChNo)=DEA |
|
./pkg/fizhi/fizhi_lsm.F |
989 |
FLUXES |
STRDG8(ChNo)=-A12*DEA/A11 |
|
./pkg/fizhi/fizhi_lsm.F |
991 |
FLUXES |
EVAP(ChNo) = ETURB(ChNo) + DEDEA(ChNo)*DEA + |
|
./pkg/fizhi/fizhi_lsm.F |
992 |
FLUXES |
& DEDTC(ChNo)*DTC |
|
./pkg/fizhi/fizhi_lsm.F |
993 |
FLUXES |
SHFLUX(ChNo) = HSTURB(ChNo) + DHSDEA(ChNo)*DEA + |
|
./pkg/fizhi/fizhi_lsm.F |
994 |
FLUXES |
& DHSDTC(ChNo)*DTC |