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

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

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

revision 1.1 by molod, Fri Jul 16 19:37:04 2004 UTC revision 1.2 by molod, Fri Jul 16 19:52:21 2004 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4        SUBROUTINE TILE (        SUBROUTINE TILE (
5       I                 NCH, DTSTEP, ITYP, TRAINL,TRAINC, TSNOW,  UM,       I                 NCH, DTSTEP, ITYP, TRAINL,TRAINC, TSNOW,  UM,
6       I                 ETURB,  DEDQA,  DEDTC,  HSTURB, DHSDQA, DHSDTC,       I                 ETURB,  DEDQA,  DEDTC,  HSTURB, DHSDQA, DHSDTC,
# Line 23  C**** Line 26  C****
26        IMPLICIT NONE        IMPLICIT NONE
27  C  C
28  C****  C****
29        include 'sibber.h'  #include "sibber.h"
30  C****  C****
31        INTEGER NCH        INTEGER NCH
32        INTEGER ITYP(NCH)        INTEGER ITYP(NCH)
# Line 51  C**** Line 54  C****
54       &                    DELEA,    VGPH1(NTYPS),     VGPH2(NTYPS),       &                    DELEA,    VGPH1(NTYPS),     VGPH2(NTYPS),
55       &            VGRPLN(NTYPS),    CSOIL0(NTYPS),          WSOI12,       &            VGRPLN(NTYPS),    CSOIL0(NTYPS),          WSOI12,
56       &             VGBEE(NTYPS),    DELZ12(NTYPS),       &             VGBEE(NTYPS),    DELZ12(NTYPS),
57       &            DELZ23(NTYPS),    SNWMID(NTYPS)       &            DELZ23(NTYPS)
58  C****  C****
59    
60        REAL PHLAY(nch,NLAY), AKLAY(nch,NLAY), SWET12(nch),        REAL PHLAY(nch,NLAY), AKLAY(nch,NLAY), SWET12(nch),
# Line 154  c       to obtain improved albedo, radsw Line 157  c       to obtain improved albedo, radsw
157  c ------------------------------------------------------------------------------------  c ------------------------------------------------------------------------------------
158        DATA CSOIL0 /175000.,175000.,175000.,175000.,175000.,        DATA CSOIL0 /175000.,175000.,175000.,175000.,175000.,
159       .             175000.,175000.,120000.,175000., 70000./       .             175000.,175000.,120000.,175000., 70000./
160        include 'snwmid.h'  #include "snwmid.h"
161    
162  C****  C****
163  C**** ---------------------------------------------------------------------  C**** ---------------------------------------------------------------------
# Line 366  C****         EIRFRC=fraction of total e Line 369  C****         EIRFRC=fraction of total e
369       I              ETURB,  DEDEA,  DEDTC, HSTURB, DHSDEA, DHSDTC,       I              ETURB,  DEDEA,  DEDTC, HSTURB, DHSDEA, DHSDTC,
370       I                 RC, DRCDEA, DRCDTC,       I                 RC, DRCDEA, DRCDTC,
371       I              SWNET, HLWDWN, ALWRAD, BLWRAD, ESNFRC,       I              SWNET, HLWDWN, ALWRAD, BLWRAD, ESNFRC,
372       I                 TM,     EM,  CSOIL,   PSUR, EMAXRT, VGWMAX,                     I                 TM,     EM,  CSOIL,   PSUR, EMAXRT, VGWMAX,
373       U                 TC,     TD,     EA,   SWET,   SNOW,       U                 TC,     TD,     EA,   SWET,   SNOW,
374       O             RUNOFF,   EVAP, SHFLUX,  SMELT,  HLWUP, BOMB,STRDG1,       O             RUNOFF,   EVAP, SHFLUX,  SMELT,  HLWUP, BOMB,STRDG1,
375       O             STRDG2, STRDG3, STRDG4, STRDG5, STRDG6, STRDG7,       O             STRDG2, STRDG3, STRDG4, STRDG5, STRDG6, STRDG7,
# Line 507  C**** runs-off, rather than freezes. Line 510  C**** runs-off, rather than freezes.
510  C****  C****
511        IMPLICIT NONE        IMPLICIT NONE
512  C****  C****
513        include 'sibber.h'  #include "sibber.h"
514  C****  C****
515        INTEGER NCH        INTEGER NCH
516        INTEGER ITYP(NCH), ChNo        INTEGER ITYP(NCH), ChNo
# Line 552  C**** (Time scale TIMFRL for large scale Line 555  C**** (Time scale TIMFRL for large scale
555  C**** to reflect the effective loss of "position memory" when storm  C**** to reflect the effective loss of "position memory" when storm
556  C**** covers entire grid square.)  C**** covers entire grid square.)
557    
558        XTCORR= (1.-TIMFRL) * MIN( 1.,(CAPAC(CHNO)/SATCAP(CHNO))/FWETL )            XTCORR= (1.-TIMFRL) * MIN( 1.,(CAPAC(CHNO)/SATCAP(CHNO))/FWETL )
559    
560  C****  C****
561  C**** Fill interception reservoir with precipitation.  C**** Fill interception reservoir with precipitation.
# Line 586  C**** Line 589  C****
589  C**** Determine XTCORR, the fraction of a storm that falls on a previously  C**** Determine XTCORR, the fraction of a storm that falls on a previously
590  C**** wet surface due to the time correlation of precipitation position.  C**** wet surface due to the time correlation of precipitation position.
591    
592        XTCORR= (1.-TIMFRC) * MIN( 1.,(CAPAC(CHNO)/SATCAP(CHNO))/FWETC )            XTCORR= (1.-TIMFRC) * MIN( 1.,(CAPAC(CHNO)/SATCAP(CHNO))/FWETC )
593    
594  C****  C****
595  C**** Fill interception reservoir with precipitation.  C**** Fill interception reservoir with precipitation.
# Line 683  C**** (p. 1353, Sellers 1985.)  Extincti Line 686  C**** (p. 1353, Sellers 1985.)  Extincti
686  C****  C****
687        IMPLICIT NONE        IMPLICIT NONE
688  C****  C****
689        include 'sibber.h'  #include "sibber.h"
690  C****  C****
691        INTEGER NCH        INTEGER NCH
692        INTEGER ITYP(NCH), ChNo        INTEGER ITYP(NCH), ChNo
# Line 781  C**** This subroutine returns soil moist Line 784  C**** This subroutine returns soil moist
784    
785        IMPLICIT NONE        IMPLICIT NONE
786  C****  C****
787        include 'sibber.h'  #include "sibber.h"
788  C****  C****
789        INTEGER NCH        INTEGER NCH
790        INTEGER ITYP(NCH), ChNo        INTEGER ITYP(NCH), ChNo
# Line 831  C**** from the surface through an energy Line 834  C**** from the surface through an energy
834  C****  C****
835        IMPLICIT NONE        IMPLICIT NONE
836  C****  C****
837        include 'sibber.h'  #include "sibber.h"
838  C****  C****
839        INTEGER NCH        INTEGER NCH
840        INTEGER ITYP(NCH), ChNo        INTEGER ITYP(NCH), ChNo
# Line 1001  C**** Line 1004  C****
1004          CHOKE = .FALSE.          CHOKE = .FALSE.
1005          DEA = EM(CHNO) - EA(CHNO)          DEA = EM(CHNO) - EA(CHNO)
1006          DTC =          DTC =
1007       &   (Q0 + ALHX(CHNO)*(ETURB(ChNo)-EMAXRT(CHNO)) - DHSDEA(CHNO)*DEA)       &   (Q0 + ALHX(CHNO)*(ETURB(ChNo)-EMAXRT(CHNO)) - DHSDEA(CHNO)*DEA)
1008       &          /  ( A11 - ALHX(CHNO)*DEDTC(ChNo) )       &          /  ( A11 - ALHX(CHNO)*DEDTC(ChNo) )
1009    
1010            STRDG1(ChNo)=DTC            STRDG1(ChNo)=DTC
# Line 1131  C**** This subroutine computes the vapor Line 1134  C**** This subroutine computes the vapor
1134  C****  C****
1135        IMPLICIT NONE        IMPLICIT NONE
1136  C****  C****
1137        include 'sibber.h'  #include "sibber.h"
1138  C****  C****
1139        INTEGER NCH        INTEGER NCH
1140        INTEGER ITYP(NCH), ChNo        INTEGER ITYP(NCH), ChNo
# Line 1172  C**** Compute temperature stress factor. Line 1175  C**** Compute temperature stress factor.
1175  C****  C****
1176        IMPLICIT NONE        IMPLICIT NONE
1177  C****  C****
1178        include 'sibber.h'  #include "sibber.h"
1179  C****  C****
1180        INTEGER NCH        INTEGER NCH
1181        INTEGER ITYP(NCH), ChNo, TypPtr        INTEGER ITYP(NCH), ChNo, TypPtr
# Line 1242  C**** leaf water potential stress. Line 1245  C**** leaf water potential stress.
1245  C****  C****
1246        IMPLICIT NONE        IMPLICIT NONE
1247  C****  C****
1248        include 'sibber.h'  #include "sibber.h"
1249  C****  C****
1250        INTEGER NCH        INTEGER NCH
1251        INTEGER ITYP(NCH), ChNo        INTEGER ITYP(NCH), ChNo
# Line 1412  C**** contents of the interception reser Line 1415  C**** contents of the interception reser
1415  C****  C****
1416        IMPLICIT NONE        IMPLICIT NONE
1417  C****  C****
1418        include 'sibber.h'  #include "sibber.h"
1419  C****  C****
1420        INTEGER NCH        INTEGER NCH
1421        INTEGER ITYP(NCH), ChNo        INTEGER ITYP(NCH), ChNo
# Line 1562  C**** This subroutine computes diffusion Line 1565  C**** This subroutine computes diffusion
1565  C****  C****
1566        IMPLICIT NONE        IMPLICIT NONE
1567  C****  C****
1568        include 'sibber.h'  #include "sibber.h"
1569  C****  C****
1570        INTEGER NCH        INTEGER NCH
1571        INTEGER ITYP(NCH), ChNo        INTEGER ITYP(NCH), ChNo
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22