| File |
Line number |
Procedure |
Code |
|
./pkg/exch2/w2_eeboot.F |
40 |
W2_EEBOOT |
INTEGER ii, I, J |
|
./pkg/exch2/w2_eeboot.F |
47 |
W2_EEBOOT |
DO I = 1,W2_maxNbTiles |
|
./pkg/exch2/w2_eeboot.F |
48 |
W2_EEBOOT |
exch2_tNx(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
49 |
W2_EEBOOT |
exch2_tNy(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
50 |
W2_EEBOOT |
exch2_tBasex(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
51 |
W2_EEBOOT |
exch2_tBasey(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
52 |
W2_EEBOOT |
exch2_txGlobalo(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
53 |
W2_EEBOOT |
exch2_tyGlobalo(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
54 |
W2_EEBOOT |
exch2_isWedge(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
55 |
W2_EEBOOT |
exch2_isNedge(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
56 |
W2_EEBOOT |
exch2_isEedge(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
57 |
W2_EEBOOT |
exch2_isSedge(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
58 |
W2_EEBOOT |
exch2_myFace(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
59 |
W2_EEBOOT |
exch2_mydNx(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
60 |
W2_EEBOOT |
exch2_mydNy(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
61 |
W2_EEBOOT |
exch2_nNeighbours(I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
63 |
W2_EEBOOT |
exch2_neighbourId(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
64 |
W2_EEBOOT |
exch2_opposingSend(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
66 |
W2_EEBOOT |
exch2_pij(ii,J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
68 |
W2_EEBOOT |
exch2_oi(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
69 |
W2_EEBOOT |
exch2_oj(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
70 |
W2_EEBOOT |
exch2_iLo(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
71 |
W2_EEBOOT |
exch2_iHi(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
72 |
W2_EEBOOT |
exch2_jLo(J,I) = 0 |
|
./pkg/exch2/w2_eeboot.F |
73 |
W2_EEBOOT |
exch2_jHi(J,I) = 0 |
|
./pkg/exch2/w2_map_procs.F |
45 |
W2_MAP_PROCS |
INTEGER I, J |
|
./pkg/exch2/w2_map_procs.F |
50 |
W2_MAP_PROCS |
DO I = 1,W2_maxNbTiles |
|
./pkg/exch2/w2_map_procs.F |
51 |
W2_MAP_PROCS |
W2_tileProc(I) = 0 |
|
./pkg/exch2/w2_map_procs.F |
52 |
W2_MAP_PROCS |
W2_tileIndex(I) = 0 |
|
./pkg/exch2/w2_map_procs.F |
77 |
W2_MAP_PROCS |
DO I=1,exch2_nTiles |
|
./pkg/exch2/w2_map_procs.F |
78 |
W2_MAP_PROCS |
IF ( exch2_myFace(I) .NE. 0 ) THEN |
|
./pkg/exch2/w2_map_procs.F |
94 |
W2_MAP_PROCS |
W2_tileProc(I) = np |
|
./pkg/exch2/w2_map_procs.F |
95 |
W2_MAP_PROCS |
W2_tileIndex(I)= bi + (bj-1)*nSx |
|
./pkg/exch2/w2_map_procs.F |
96 |
W2_MAP_PROCS |
W2_procTileList(bi,bj,np) = I |
|
./pkg/exch2/w2_map_procs.F |
97 |
W2_MAP_PROCS |
IF ( np.EQ.thisProc ) W2_myTileList(bi,bj) = I |
|
./pkg/exf/exf_mapfields.F |
342 |
EXF_MAPFIELDS |
saltFlux(I,J,bi,bj) = saltflx(I,J,bi,bj) |
|
./pkg/exf/exf_mapfields.F |
352 |
EXF_MAPFIELDS |
exf_iceFraction(I,J,bi,bj) = |
|
./pkg/exf/exf_mapfields.F |
353 |
EXF_MAPFIELDS |
& MIN( MAX(exf_iceFraction(I,J,bi,bj),zeroRS), oneRS ) |
|
./pkg/fizhi/fizhi_clockstuff.F |
1001 |
ASTRO |
HC = COS( FAC+ALON(I)*DG2RD ) |
|
./pkg/fizhi/fizhi_clockstuff.F |
1002 |
ASTRO |
SJ = SIN(ALAT(I)*DG2RD) |
|
./pkg/fizhi/fizhi_clockstuff.F |
1005 |
ASTRO |
COSZ(I) = SJ*ZS + CJ*ZC*HC |
|
./pkg/fizhi/fizhi_clockstuff.F |
1006 |
ASTRO |
IF( COSZ(I).LT.ZERO ) COSZ(I) = ZERO |
|
./pkg/fizhi/fizhi_clockstuff.F |
1138 |
TIME2FREQ2 |
DO I=15,48 |
|
./pkg/fizhi/fizhi_clockstuff.F |
1139 |
TIME2FREQ2 |
MNDY(I) = MNDY(I-12) + 365 |
|
./pkg/fizhi/fizhi_clockstuff.F |
361 |
GET_TIME |
DO I=15,48 |
|
./pkg/fizhi/fizhi_clockstuff.F |
363 |
GET_TIME |
MNDY(I) = MNDY(I-12) + 365 |
|
./pkg/fizhi/fizhi_clockstuff.F |
490 |
INTERP_TIME |
DO I=15,48 |
|
./pkg/fizhi/fizhi_clockstuff.F |
491 |
INTERP_TIME |
MNDY(I) = MNDY(I-12) + 365 |
|
./pkg/fizhi/fizhi_clockstuff.F |
636 |
TIC_TIME |
DO I=15,48 |
|
./pkg/fizhi/fizhi_clockstuff.F |
637 |
TIC_TIME |
MNDY(I) = MNDY(I-12) + 365 |
|
./pkg/fizhi/fizhi_clockstuff.F |
710 |
TIC_TIME |
DO I=15,48 |
|
./pkg/fizhi/fizhi_clockstuff.F |
711 |
TIC_TIME |
MNDY(I) = MNDY(I-12) + 365 |
|
./pkg/fizhi/fizhi_clockstuff.F |
951 |
ASTRO |
DO 100 I=15,48 |
|
./pkg/fizhi/fizhi_clockstuff.F |
952 |
ASTRO |
MNDY(I) = MNDY(I-12) + 365 |
|
./pkg/fizhi/fizhi_clockstuff.F |
999 |
ASTRO |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_lwrad.F |
145 |
LWRIO |
DO I = 1,ISTRIP |
|
./pkg/fizhi/fizhi_lwrad.F |
146 |
LWRIO |
ADELPL(I,L) = convrt / ( ple(I,L+1)-ple(I,L) ) |
|
./pkg/fizhi/fizhi_lwrad.F |
154 |
LWRIO |
DO I = 1,ISTRIP |
|
./pkg/fizhi/fizhi_lwrad.F |
155 |
LWRIO |
CLRO(I,L) = min( 1.0 _d 0,clro(i,L) ) |
|
./pkg/fizhi/fizhi_lwrad.F |
156 |
LWRIO |
CLMO(I,L) = min( 1.0 _d 0,clmo(i,L) ) |
|
./pkg/fizhi/fizhi_lwrad.F |
164 |
LWRIO |
DO I = 1,ISTRIP |
|
./pkg/fizhi/fizhi_lwrad.F |
165 |
LWRIO |
TZL(I,L) = TZL(I,L)*pk(I,L) |
|
./pkg/fizhi/fizhi_lwrad.F |
168 |
LWRIO |
DO I = 1,ISTRIP |
|
./pkg/fizhi/fizhi_lwrad.F |
170 |
LWRIO |
TZL(I,lm) = ( 2*tzl(i,lm)+tsurf(i) )/3.0 |
|
./pkg/fizhi/fizhi_lwrad.F |
236 |
LWRIO |
DO I = 1,ISTRIP |
|
./pkg/fizhi/fizhi_lwrad.F |
263 |
LWRIO |
DO I = 1,ISTRIP |
|
./pkg/fizhi/fizhi_lwrad.F |
264 |
LWRIO |
DTRAD (I,L) = ple(i,lm+1) * DTRAD (I,L)/pk(I,L) |
|
./pkg/fizhi/fizhi_lwrad.F |
265 |
LWRIO |
DTRADC(I,L) = ple(i,lm+1) * DTRADC(I,L)/pk(I,L) |
|
./pkg/fizhi/fizhi_lwrad.F |
266 |
LWRIO |
dtdtg(I,L) = ple(i,lm+1) * dtdtg (I,L)/pk(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
1236 |
RAS |
do I = 1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1271 |
RAS |
DO I=1,ncrnd |
|
./pkg/fizhi/fizhi_moist.F |
1272 |
RAS |
IRND(I) = (RND(I)-0.0005)*(KCR-KRMIN+1) |
|
./pkg/fizhi/fizhi_moist.F |
1273 |
RAS |
IRND(I) = IRND(I) + KRMIN |
|
./pkg/fizhi/fizhi_moist.F |
1340 |
RAS |
DO I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1341 |
RAS |
thbef(I,L) = POI(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
1342 |
RAS |
POI(I,L) = POI(I,L) + TCU(I,L) * DT * rhfrac(i) |
|
./pkg/fizhi/fizhi_moist.F |
1343 |
RAS |
QOI(I,L) = QOI(I,L) + QCU(I,L) * DT * rhfrac(i) |
|
./pkg/fizhi/fizhi_moist.F |
1348 |
RAS |
DO I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1349 |
RAS |
UOI(I,L+nltop-1,NT)=UOI(I,L+nltop-1,NT)+UCU(I,L,NT)*DT*rhfrac(i) |
|
./pkg/fizhi/fizhi_moist.F |
1353 |
RAS |
DO I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1354 |
RAS |
rains(I,ib) = rains(I,ib) + PCU(I)*dt * rhfrac(i) |
|
./pkg/fizhi/fizhi_moist.F |
1365 |
RAS |
temp = TCU(I,L) * DT * rhfrac(i) |
|
./pkg/fizhi/fizhi_moist.F |
1690 |
CLOUD |
DO 2030 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1691 |
CLOUD |
PRH(I,L) = (PRJ(I,L+1)*PRS(I,L+1) - PRJ(I,L)*PRS(I,L)) |
|
./pkg/fizhi/fizhi_moist.F |
1692 |
CLOUD |
* / ((PRS(I,L+1)-PRS(I,L)) * RKAPP1) |
|
./pkg/fizhi/fizhi_moist.F |
1697 |
CLOUD |
DO 2060 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1698 |
CLOUD |
TX5(I) = POI(I,L) * PRH(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
1699 |
CLOUD |
TX1(I) = (PRS(I,L) + PRS(I,L+1)) * 0.5 |
|
./pkg/fizhi/fizhi_moist.F |
1700 |
CLOUD |
TX3(I) = TX5(I) |
|
./pkg/fizhi/fizhi_moist.F |
1701 |
CLOUD |
CALL QSAT(TX3(I), TX1(I), TX2(I), TX4(I), .TRUE.) |
|
./pkg/fizhi/fizhi_moist.F |
1702 |
CLOUD |
ALF(I,L) = TX2(I) - TX4(I) * TX5(I) |
|
./pkg/fizhi/fizhi_moist.F |
1703 |
CLOUD |
BET(I,L) = TX4(I) * PRH(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
1704 |
CLOUD |
GAM(I,L) = 1.0 / ((1.0 + TX4(I)*ALBCP) * PRH(I,L)) |
|
./pkg/fizhi/fizhi_moist.F |
1705 |
CLOUD |
PRI(I,L) = (CP/CMB2PA) / (PRS(I,L+1) - PRS(I,L)) |
|
./pkg/fizhi/fizhi_moist.F |
1713 |
CLOUD |
DO 10 I=1,lng |
|
./pkg/fizhi/fizhi_moist.F |
1714 |
CLOUD |
TCU(I,L) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1715 |
CLOUD |
QCU(I,L) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1716 |
CLOUD |
CMASS(I,L) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1721 |
CLOUD |
do I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1722 |
CLOUD |
ucu(I,L,nt) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1727 |
CLOUD |
DO 30 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1728 |
CLOUD |
TX1(I) = PRJ(I,K+1) * POI(I,K) |
|
./pkg/fizhi/fizhi_moist.F |
1729 |
CLOUD |
QS1(I) = ALF(I,K) + BET(I,K)*POI(I,K) |
|
./pkg/fizhi/fizhi_moist.F |
1730 |
CLOUD |
QOL(I,K) = MIN(QS1(I)*RHMAX,QOI(I,K)) |
|
./pkg/fizhi/fizhi_moist.F |
1732 |
CLOUD |
HOL(I,K) = TX1(I)*CP + QOL(I,K)*ALHL |
|
./pkg/fizhi/fizhi_moist.F |
1733 |
CLOUD |
ETA(I,K) = ZERO |
|
./pkg/fizhi/fizhi_moist.F |
1734 |
CLOUD |
TX2(I) = (PRJ(I,K+1) - PRJ(I,K)) * POI(I,K) * CP |
|
./pkg/fizhi/fizhi_moist.F |
1739 |
CLOUD |
DO 50 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1740 |
CLOUD |
QS1(I) = ALF(I,L) + BET(I,L)*POI(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
1741 |
CLOUD |
QOL(I,L) = MIN(QS1(I)*RHMAX,QOI(I,L)) |
|
./pkg/fizhi/fizhi_moist.F |
1743 |
CLOUD |
TEM1 = TX2(I) + PRJ(I,L+1) * POI(I,L) * CP |
|
./pkg/fizhi/fizhi_moist.F |
1744 |
CLOUD |
HOL(I,L) = TEM1 + QOL(I,L )* ALHL |
|
./pkg/fizhi/fizhi_moist.F |
1745 |
CLOUD |
HST(I,L) = TEM1 + QS1(I) * ALHL |
|
./pkg/fizhi/fizhi_moist.F |
1747 |
CLOUD |
TX1(I) = (PRJ(I,L+1) - PRJ(I,L)) * POI(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
1748 |
CLOUD |
ETA(I,L) = ETA(I,L+1) + TX1(I)*CPBG |
|
./pkg/fizhi/fizhi_moist.F |
1749 |
CLOUD |
TX2(I) = TX2(I) + TX1(I)*CP |
|
./pkg/fizhi/fizhi_moist.F |
1756 |
CLOUD |
DO 70 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1757 |
CLOUD |
HOL(I,IC) = TX2(I) |
|
./pkg/fizhi/fizhi_moist.F |
1758 |
CLOUD |
QS1(I) = ALF(I,IC) + BET(I,IC)*POI(I,IC) |
|
./pkg/fizhi/fizhi_moist.F |
1759 |
CLOUD |
QOL(I,IC) = MIN(QS1(I)*RHMAX,QOI(I,IC)) |
|
./pkg/fizhi/fizhi_moist.F |
1761 |
CLOUD |
TEM1 = TX2(I) + PRJ(I,IC1) * POI(I,IC) * CP |
|
./pkg/fizhi/fizhi_moist.F |
1762 |
CLOUD |
HOL(I,IC) = TEM1 + QOL(I,IC) * ALHL |
|
./pkg/fizhi/fizhi_moist.F |
1763 |
CLOUD |
HST(I,IC) = TEM1 + QS1(I) * ALHL |
|
./pkg/fizhi/fizhi_moist.F |
1765 |
CLOUD |
TX3(I ) = (PRJ(I,IC1) - PRH(I,IC)) * POI(I,IC) |
|
./pkg/fizhi/fizhi_moist.F |
1766 |
CLOUD |
ETA(I,IC) = ETA(I,IC1) + CPBG * TX3(I) |
|
./pkg/fizhi/fizhi_moist.F |
1769 |
CLOUD |
DO 130 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1770 |
CLOUD |
TX2(I) = HOL(I,K) - HST(I,IC) |
|
./pkg/fizhi/fizhi_moist.F |
1771 |
CLOUD |
TX1(I) = ZERO |
|
./pkg/fizhi/fizhi_moist.F |
1778 |
CLOUD |
DO 160 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1779 |
CLOUD |
TX1(I) = TX1(I) + (HST(I,IC) - HOL(I,L)) * (ETA(I,L) - ETA(I,L+1)) |
|
./pkg/fizhi/fizhi_moist.F |
1785 |
CLOUD |
DO 195 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1786 |
CLOUD |
IF (TX1(I) .GT. ZERO .AND. TX2(I) .GT. ZERO |
|
./pkg/fizhi/fizhi_moist.F |
1789 |
CLOUD |
IA(LEN1) = I |
|
./pkg/fizhi/fizhi_moist.F |
1790 |
CLOUD |
ALM(LEN1) = TX2(I) / TX1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1796 |
CLOUD |
DO 196 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1797 |
CLOUD |
IF (TX2(I) .LE. 0.0 .AND. (HOL(I,K) .GT. HST(I,IC1)) |
|
./pkg/fizhi/fizhi_moist.F |
1800 |
CLOUD |
IA(LEN2) = I |
|
./pkg/fizhi/fizhi_moist.F |
1812 |
CLOUD |
DO I = 1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1813 |
CLOUD |
HST(I,L) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1814 |
CLOUD |
QOL(I,L) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1817 |
CLOUD |
DO 5020 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1818 |
CLOUD |
PCU(I) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1826 |
CLOUD |
DO 250 I=1,LEN2 |
|
./pkg/fizhi/fizhi_moist.F |
1827 |
CLOUD |
ETA(I,K) = 1.0 |
|
./pkg/fizhi/fizhi_moist.F |
1828 |
CLOUD |
II = IA(I) |
|
./pkg/fizhi/fizhi_moist.F |
1829 |
CLOUD |
TX2(I) = 0.5 * (PRS(II,IC) + PRS(II,IC1)) |
|
./pkg/fizhi/fizhi_moist.F |
1830 |
CLOUD |
TX4(I) = PRS(II,K) |
|
./pkg/fizhi/fizhi_moist.F |
1833 |
CLOUD |
DO 252 I=LEN11,LEN2 |
|
./pkg/fizhi/fizhi_moist.F |
1834 |
CLOUD |
WFN(I) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1835 |
CLOUD |
II = IA(I) |
|
./pkg/fizhi/fizhi_moist.F |
1837 |
CLOUD |
TX6(I) = (HST(II,IC1)-HOL(II,K))/(HST(II,IC1)-HST(II,IC)) |
|
./pkg/fizhi/fizhi_moist.F |
1839 |
CLOUD |
TX6(I) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1841 |
CLOUD |
TX2(I) = 0.5 * (PRS(II,IC1)+PRS(II,IC1+1)) * (1.0-TX6(I)) |
|
./pkg/fizhi/fizhi_moist.F |
1842 |
CLOUD |
* + TX2(I) * TX6(I) |
|
./pkg/fizhi/fizhi_moist.F |
1848 |
CLOUD |
DO 255 I=1,LEN2 |
|
./pkg/fizhi/fizhi_moist.F |
1849 |
CLOUD |
TX1(I) = ETA(IA(I),L) |
|
./pkg/fizhi/fizhi_moist.F |
1851 |
CLOUD |
DO 260 I=1,LEN2 |
|
./pkg/fizhi/fizhi_moist.F |
1852 |
CLOUD |
ETA(I,L) = 1.0 + ALM(I) * TX1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1858 |
CLOUD |
DO 270 I=1,LEN1 |
|
./pkg/fizhi/fizhi_moist.F |
1859 |
CLOUD |
II = IA(I) |
|
./pkg/fizhi/fizhi_moist.F |
1860 |
CLOUD |
WFN(I) = - GAM(II,IC) * (PRJ(II,IC1) - PRH(II,IC)) |
|
./pkg/fizhi/fizhi_moist.F |
1861 |
CLOUD |
* * HST(II,IC) * ETA(I,IC1) |
|
./pkg/fizhi/fizhi_moist.F |
1865 |
CLOUD |
DO 290 I=1,LEN2 |
|
./pkg/fizhi/fizhi_moist.F |
1866 |
CLOUD |
II = IA(I) |
|
./pkg/fizhi/fizhi_moist.F |
1867 |
CLOUD |
TX1(I) = HOL(II,K) |
|
./pkg/fizhi/fizhi_moist.F |
1873 |
CLOUD |
DO 380 I=1,LEN2 |
|
./pkg/fizhi/fizhi_moist.F |
1874 |
CLOUD |
II = IA(I) |
|
./pkg/fizhi/fizhi_moist.F |
1875 |
CLOUD |
TEM = TX1(I) + (ETA(I,L) - ETA(I,L+1)) * HOL(II,L) |
|
./pkg/fizhi/fizhi_moist.F |
1877 |
CLOUD |
PCU(I) = PRJ(II,L+1) - PRH(II,L) |
|
./pkg/fizhi/fizhi_moist.F |
1878 |
CLOUD |
TEM1 = ETA(I,L+1) * PCU(I) |
|
./pkg/fizhi/fizhi_moist.F |
1879 |
CLOUD |
TX1(I) = TX1(I)*PCU(I) |
|
./pkg/fizhi/fizhi_moist.F |
1881 |
CLOUD |
PCU(I) = PRH(II,L) - PRJ(II,L) |
|
./pkg/fizhi/fizhi_moist.F |
1882 |
CLOUD |
TEM1 = (TEM1 + ETA(I,L) * PCU(I)) * HST(II,L) |
|
./pkg/fizhi/fizhi_moist.F |
1883 |
CLOUD |
TX1(I) = TX1(I) + TEM*PCU(I) |
|
./pkg/fizhi/fizhi_moist.F |
1885 |
CLOUD |
WFN(I) = WFN(I) + (TX1(I) - TEM1) * GAM(II,L) |
|
./pkg/fizhi/fizhi_moist.F |
1886 |
CLOUD |
TX1(I) = TEM |
|
./pkg/fizhi/fizhi_moist.F |
1892 |
CLOUD |
DO 512 I=1,LEN1 |
|
./pkg/fizhi/fizhi_moist.F |
1893 |
CLOUD |
II = IA(I) |
|
./pkg/fizhi/fizhi_moist.F |
1894 |
CLOUD |
WFN(I) = WFN(I) + TX1(I) * GAM(II,IC)*(PRJ(II,IC1)-PRH(II,IC)) |
|
./pkg/fizhi/fizhi_moist.F |
1895 |
CLOUD |
* - TX3(I) |
|
./pkg/fizhi/fizhi_moist.F |
1896 |
CLOUD |
IF (WFN(I) .GT. 0.0) THEN |
|
./pkg/fizhi/fizhi_moist.F |
1898 |
CLOUD |
I1(LENA) = IA(I) |
|
./pkg/fizhi/fizhi_moist.F |
1899 |
CLOUD |
I2(LENA) = I |
|
./pkg/fizhi/fizhi_moist.F |
1900 |
CLOUD |
TX1(LENA) = WFN(I) |
|
./pkg/fizhi/fizhi_moist.F |
1901 |
CLOUD |
TX2(LENA) = QS1(IA(I)) |
|
./pkg/fizhi/fizhi_moist.F |
1907 |
CLOUD |
DO 515 I=LEN11,LEN2 |
|
./pkg/fizhi/fizhi_moist.F |
1908 |
CLOUD |
WFN(I) = WFN(I) - TX3(I) |
|
./pkg/fizhi/fizhi_moist.F |
1909 |
CLOUD |
IF (WFN(I) .GT. 0.0 .AND. TX6(I) .GT. 0.0) THEN |
|
./pkg/fizhi/fizhi_moist.F |
1911 |
CLOUD |
I1(LENB) = IA(I) |
|
./pkg/fizhi/fizhi_moist.F |
1912 |
CLOUD |
I2(LENB) = I |
|
./pkg/fizhi/fizhi_moist.F |
1913 |
CLOUD |
TX1(LENB) = WFN(I) |
|
./pkg/fizhi/fizhi_moist.F |
1914 |
CLOUD |
TX2(LENB) = QS1(IA(I)) |
|
./pkg/fizhi/fizhi_moist.F |
1915 |
CLOUD |
TX4(LENB) = TX6(I) |
|
./pkg/fizhi/fizhi_moist.F |
1925 |
CLOUD |
DO I = 1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1926 |
CLOUD |
HST(I,L) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1927 |
CLOUD |
QOL(I,L) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1930 |
CLOUD |
DO 5040 I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
1931 |
CLOUD |
PCU(I) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
1937 |
CLOUD |
DO 516 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
1938 |
CLOUD |
WFN(I) = TX1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1939 |
CLOUD |
QS1(I) = TX2(I) |
|
./pkg/fizhi/fizhi_moist.F |
1943 |
CLOUD |
DO 517 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
1944 |
CLOUD |
TX1(I) = ETA(I2(I),L) |
|
./pkg/fizhi/fizhi_moist.F |
1946 |
CLOUD |
DO 520 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
1947 |
CLOUD |
ETA(I,L) = TX1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1952 |
CLOUD |
DO 510 I=1,LENA |
|
./pkg/fizhi/fizhi_moist.F |
1953 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1954 |
CLOUD |
TX8(I) = HST(II,IC) - HOL(II,IC) |
|
./pkg/fizhi/fizhi_moist.F |
1956 |
CLOUD |
DO 530 I=LENA1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
1957 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1958 |
CLOUD |
TX6(I) = TX4(I) |
|
./pkg/fizhi/fizhi_moist.F |
1959 |
CLOUD |
TEM = TX6(I) * (HOL(II,IC)-HOL(II,IC1)) + HOL(II,IC1) |
|
./pkg/fizhi/fizhi_moist.F |
1960 |
CLOUD |
TX8(I) = HOL(II,K) - TEM |
|
./pkg/fizhi/fizhi_moist.F |
1962 |
CLOUD |
TEM1 = TX6(I) * (QOL(II,IC)-QOL(II,IC1)) + QOL(II,IC1) |
|
./pkg/fizhi/fizhi_moist.F |
1963 |
CLOUD |
TX5(I) = TEM - TEM1 * ALHL |
|
./pkg/fizhi/fizhi_moist.F |
1964 |
CLOUD |
QS1(I) = TEM1 + TX8(I)*(ONE/ALHL) |
|
./pkg/fizhi/fizhi_moist.F |
1965 |
CLOUD |
TX3(I) = HOL(II,IC) |
|
./pkg/fizhi/fizhi_moist.F |
1969 |
CLOUD |
DO 620 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
1970 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1971 |
CLOUD |
WLQ(I) = QOL(II,K) - QS1(I) * ETA(I,IC) |
|
./pkg/fizhi/fizhi_moist.F |
1972 |
CLOUD |
TX7(I) = HOL(II,K) |
|
./pkg/fizhi/fizhi_moist.F |
1975 |
CLOUD |
DO 621 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
1976 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1977 |
CLOUD |
UHT(I,NT) = UOI(II,K+nltop-1,NT)-UOI(II,IC+nltop-1,NT) * ETA(I,IC) |
|
./pkg/fizhi/fizhi_moist.F |
1982 |
CLOUD |
DO 630 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
1983 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1984 |
CLOUD |
TEM = ETA(I,L) - ETA(I,L+1) |
|
./pkg/fizhi/fizhi_moist.F |
1985 |
CLOUD |
WLQ(I) = WLQ(I) + TEM * QOL(II,L) |
|
./pkg/fizhi/fizhi_moist.F |
1990 |
CLOUD |
DO I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
1991 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
1992 |
CLOUD |
TEM = ETA(I,L) - ETA(I,L+1) |
|
./pkg/fizhi/fizhi_moist.F |
1993 |
CLOUD |
UHT(I,NT) = UHT(I,NT) + TEM * UOI(II,L+nltop-1,NT) |
|
./pkg/fizhi/fizhi_moist.F |
2000 |
CLOUD |
DO 690 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2001 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2004 |
CLOUD |
HOL(I,K) = TEM * (PRJ(II,K)-PRH(II,KM1))*PRH(II,K)*PRI(II,K) |
|
./pkg/fizhi/fizhi_moist.F |
2005 |
CLOUD |
HOL(I,KM1) = TEM * (PRH(II,K)-PRJ(II,K))*PRH(II,KM1)*PRI(II,KM1) |
|
./pkg/fizhi/fizhi_moist.F |
2006 |
CLOUD |
AKM(I) = ZERO |
|
./pkg/fizhi/fizhi_moist.F |
2007 |
CLOUD |
TX2(I) = 0.5 * (PRS(II,IC) + PRS(II,IC1)) |
|
./pkg/fizhi/fizhi_moist.F |
2012 |
CLOUD |
DO 750 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2013 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2014 |
CLOUD |
TEM = (POI(II,L-1) - POI(II,L)) * ETA(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2017 |
CLOUD |
HOL(I,L) = TEM * (PRJ(II,L)-PRH(II,L-1)) * PRH(II,L) |
|
./pkg/fizhi/fizhi_moist.F |
2018 |
CLOUD |
* * PRI(II,L) + HOL(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2019 |
CLOUD |
HOL(I,L-1) = TEM * (PRH(II,L)-PRJ(II,L)) * PRH(II,L-1) |
|
./pkg/fizhi/fizhi_moist.F |
2022 |
CLOUD |
AKM(I) = AKM(I) - HOL(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2023 |
CLOUD |
* * (ETA(I,L) * (PRH(II,L)-PRJ(II,L)) + |
|
./pkg/fizhi/fizhi_moist.F |
2024 |
CLOUD |
* ETA(I,L+1) * (PRJ(II,L+1)-PRH(II,L))) / PRH(II,L) |
|
./pkg/fizhi/fizhi_moist.F |
2031 |
CLOUD |
DO 770 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2032 |
CLOUD |
TX2(I) = (ONE - TX1(I)) * WLQ(I) |
|
./pkg/fizhi/fizhi_moist.F |
2033 |
CLOUD |
WLQ(I) = TX1(I) * WLQ(I) |
|
./pkg/fizhi/fizhi_moist.F |
2035 |
CLOUD |
TX1(I) = HOL(I,IC) |
|
./pkg/fizhi/fizhi_moist.F |
2037 |
CLOUD |
DO 790 I=LENA1, LENB |
|
./pkg/fizhi/fizhi_moist.F |
2038 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2039 |
CLOUD |
TX1(I) = TX1(I) + (TX5(I)-TX3(I)+QOL(II,IC)*ALHL)*(PRI(II,IC)/CP) |
|
./pkg/fizhi/fizhi_moist.F |
2042 |
CLOUD |
DO 800 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2043 |
CLOUD |
HOL(I,IC) = TX1(I) - TX2(I) * ALBCP * PRI(I1(I),IC) |
|
./pkg/fizhi/fizhi_moist.F |
2047 |
CLOUD |
DO 810 I=1,LENA |
|
./pkg/fizhi/fizhi_moist.F |
2048 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2049 |
CLOUD |
AKM(I) = AKM(I) - ETA(I,IC1) * (PRJ(II,IC1) - PRH(II,IC)) |
|
./pkg/fizhi/fizhi_moist.F |
2050 |
CLOUD |
* * TX1(I) / PRH(II,IC) |
|
./pkg/fizhi/fizhi_moist.F |
2056 |
CLOUD |
DO 830 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2057 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2058 |
CLOUD |
TX3(I) = QOL(II,KM1) - QOL(II,K) |
|
./pkg/fizhi/fizhi_moist.F |
2059 |
CLOUD |
GMH(I,K) = HOL(I,K) + TX3(I) * PRI(II,K) * (ALBCP) |
|
./pkg/fizhi/fizhi_moist.F |
2061 |
CLOUD |
AKM(I) = AKM(I) + GAM(II,KM1)*(PRJ(II,K)-PRH(II,KM1)) |
|
./pkg/fizhi/fizhi_moist.F |
2062 |
CLOUD |
* * GMH(I,K) |
|
./pkg/fizhi/fizhi_moist.F |
2063 |
CLOUD |
TX3(I) = zero |
|
./pkg/fizhi/fizhi_moist.F |
2068 |
CLOUD |
DO 840 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2069 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2070 |
CLOUD |
TX2(I) = TX3(I) |
|
./pkg/fizhi/fizhi_moist.F |
2071 |
CLOUD |
TX3(I) = (QOL(II,L-1) - QOL(II,L)) * ETA(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2072 |
CLOUD |
TX2(I) = TX2(I) + TX3(I) |
|
./pkg/fizhi/fizhi_moist.F |
2074 |
CLOUD |
GMH(I,L) = HOL(I,L) + TX2(I) * PRI(II,L) * (ALBCP*HALF) |
|
./pkg/fizhi/fizhi_moist.F |
2079 |
CLOUD |
DO 850 I=LENA1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2080 |
CLOUD |
TX3(I) = TX3(I) + TWOBAL |
|
./pkg/fizhi/fizhi_moist.F |
2081 |
CLOUD |
* * (TX7(I) - TX8(I) - TX5(I) - QOL(I1(I),IC)*ALHL) |
|
./pkg/fizhi/fizhi_moist.F |
2083 |
CLOUD |
DO 860 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2084 |
CLOUD |
GMH(I,IC) = TX1(I) + PRI(I1(I),IC) * ONEBCP |
|
./pkg/fizhi/fizhi_moist.F |
2085 |
CLOUD |
* * (TX3(I)*(ALHL*HALF) + ETA(I,IC) * TX8(I)) |
|
./pkg/fizhi/fizhi_moist.F |
2091 |
CLOUD |
DO 870 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2092 |
CLOUD |
TX1(I) = GMH(I,K) |
|
./pkg/fizhi/fizhi_moist.F |
2095 |
CLOUD |
DO 880 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2096 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2097 |
CLOUD |
TX1(I) = TX1(I) + (ETA(I,L) - ETA(I,L+1)) * GMH(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2098 |
CLOUD |
TX2(I) = GAM(II,L-1) * (PRJ(II,L) - PRH(II,L-1)) |
|
./pkg/fizhi/fizhi_moist.F |
2102 |
CLOUD |
DO 890 I=LENA1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2103 |
CLOUD |
TX2(I) = ZERO |
|
./pkg/fizhi/fizhi_moist.F |
2106 |
CLOUD |
DO 900 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2107 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2108 |
CLOUD |
AKM(I) = AKM(I) + TX1(I) * |
|
./pkg/fizhi/fizhi_moist.F |
2109 |
CLOUD |
* (TX2(I) + GAM(II,L)*(PRH(II,L)-PRJ(II,L))) |
|
./pkg/fizhi/fizhi_moist.F |
2114 |
CLOUD |
DO 920 I=LENA1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2115 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2116 |
CLOUD |
TX2(I) = 0.5 * (PRS(II,IC) + PRS(II,IC1)) |
|
./pkg/fizhi/fizhi_moist.F |
2117 |
CLOUD |
* + 0.5*(PRS(II,IC+2) - PRS(II,IC)) * (ONE-TX6(I)) |
|
./pkg/fizhi/fizhi_moist.F |
2119 |
CLOUD |
TX1(I) = PRS(II,IC1) |
|
./pkg/fizhi/fizhi_moist.F |
2120 |
CLOUD |
TX5(I) = 0.5 * (PRS(II,IC1) + PRS(II,IC+2)) |
|
./pkg/fizhi/fizhi_moist.F |
2122 |
CLOUD |
IF ((TX2(I) .GE. TX1(I)) .AND. (TX2(I) .LT. TX5(I))) THEN |
|
./pkg/fizhi/fizhi_moist.F |
2123 |
CLOUD |
TX6(I) = ONE - (TX2(I) - TX1(I)) / (TX5(I) - TX1(I)) |
|
./pkg/fizhi/fizhi_moist.F |
2126 |
CLOUD |
HOL(I,IC1) = HOL(I,IC1) + HOL(I,IC) * TEM |
|
./pkg/fizhi/fizhi_moist.F |
2127 |
CLOUD |
HOL(I,IC) = ZERO |
|
./pkg/fizhi/fizhi_moist.F |
2129 |
CLOUD |
GMH(I,IC1) = GMH(I,IC1) + GMH(I,IC) * TEM |
|
./pkg/fizhi/fizhi_moist.F |
2130 |
CLOUD |
GMH(I,IC) = ZERO |
|
./pkg/fizhi/fizhi_moist.F |
2131 |
CLOUD |
ELSEIF (TX2(I) .LT. TX1(I)) THEN |
|
./pkg/fizhi/fizhi_moist.F |
2132 |
CLOUD |
TX6(I) = 1.0 |
|
./pkg/fizhi/fizhi_moist.F |
2134 |
CLOUD |
TX6(I) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
2139 |
CLOUD |
DO I=1,LENC |
|
./pkg/fizhi/fizhi_moist.F |
2140 |
CLOUD |
PCU(I) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
2143 |
CLOUD |
DO 970 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2144 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2145 |
CLOUD |
IF (AKM(I) .LT. ZERO .AND. WLQ(I) .GE. 0.0) THEN |
|
./pkg/fizhi/fizhi_moist.F |
2146 |
CLOUD |
WFN(I) = - TX6(I) * WFN(I) * RASALF / AKM(I) |
|
./pkg/fizhi/fizhi_moist.F |
2148 |
CLOUD |
WFN(I) = ZERO |
|
./pkg/fizhi/fizhi_moist.F |
2151 |
CLOUD |
WFN(I) = MIN(WFN(I), TEM) |
|
./pkg/fizhi/fizhi_moist.F |
2161 |
CLOUD |
PCU(II) = WLQ(I) * WFN(I) * ONEBG |
|
./pkg/fizhi/fizhi_moist.F |
2165 |
CLOUD |
TX4(I) = WFN(I) * (1.0/ALHL) |
|
./pkg/fizhi/fizhi_moist.F |
2166 |
CLOUD |
TX5(I) = WFN(I) * ONEBCP |
|
./pkg/fizhi/fizhi_moist.F |
2172 |
CLOUD |
DO I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2173 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2175 |
CLOUD |
CMASS(II,L) = ETA(I,L+1) * WFN(I) * ONEBG |
|
./pkg/fizhi/fizhi_moist.F |
2177 |
CLOUD |
CMASS(II,L) = WFN(I) * ONEBG |
|
./pkg/fizhi/fizhi_moist.F |
2193 |
CLOUD |
DO 980 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2194 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2195 |
CLOUD |
TEM = (GMH(I,L) - HOL(I,L)) * TX4(I) |
|
./pkg/fizhi/fizhi_moist.F |
2196 |
CLOUD |
TEM1 = HOL(I,L) * TX5(I) |
|
./pkg/fizhi/fizhi_moist.F |
2216 |
CLOUD |
DO 995 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2217 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2218 |
CLOUD |
TEM = half*TX5(I) * PRI(II,K) |
|
./pkg/fizhi/fizhi_moist.F |
2219 |
CLOUD |
TX1(I) = ( UOI(II,KM1+nltop-1,nt) - UOI(II,K+nltop-1,nt)) |
|
./pkg/fizhi/fizhi_moist.F |
2220 |
CLOUD |
ucu(II,K,nt) = TEM * TX1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2226 |
CLOUD |
DO 1010 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2227 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2228 |
CLOUD |
TEM = half*TX5(I) * PRI(II,L) |
|
./pkg/fizhi/fizhi_moist.F |
2229 |
CLOUD |
TEM1 = TX1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2230 |
CLOUD |
TX1(I) = (UOI(II,L-1+nltop-1,nt)-UOI(II,L+nltop-1,nt)) * ETA(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2231 |
CLOUD |
TX3(I) = (TX1(I) + TEM1) * TEM |
|
./pkg/fizhi/fizhi_moist.F |
2233 |
CLOUD |
DO 1020 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2234 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2235 |
CLOUD |
ucu(II,L,nt) = TX3(I) |
|
./pkg/fizhi/fizhi_moist.F |
2238 |
CLOUD |
DO 1030 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2239 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2240 |
CLOUD |
IF (TX6(I) .GE. 1.0) THEN |
|
./pkg/fizhi/fizhi_moist.F |
2241 |
CLOUD |
TEM = half*TX5(I) * PRI(II,IC) |
|
./pkg/fizhi/fizhi_moist.F |
2245 |
CLOUD |
TX1(I) = (TX1(I) + UHT(I,nt) + UHT(I,nt)) * TEM |
|
./pkg/fizhi/fizhi_moist.F |
2247 |
CLOUD |
DO 1040 I=1,LENB |
|
./pkg/fizhi/fizhi_moist.F |
2248 |
CLOUD |
II = I1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2249 |
CLOUD |
ucu(II,IC,nt) = TX1(I) |
|
./pkg/fizhi/fizhi_moist.F |
2281 |
RNCL |
DO 10 I=1,lng |
|
./pkg/fizhi/fizhi_moist.F |
2292 |
RNCL |
CLF(I) = CUCLD |
|
./pkg/fizhi/fizhi_moist.F |
2300 |
RNCL |
IF (PL(I) .GE. P4 .AND. PL(I) .LE. P6) THEN |
|
./pkg/fizhi/fizhi_moist.F |
2301 |
RNCL |
CLF(I) = (P6-PL(I))*CFAC |
|
./pkg/fizhi/fizhi_moist.F |
2302 |
RNCL |
ELSEIF (PL(I) .GT. P6 ) THEN |
|
./pkg/fizhi/fizhi_moist.F |
2303 |
RNCL |
CLF(I) = 0.0 |
|
./pkg/fizhi/fizhi_moist.F |
2432 |
RNEVP |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_moist.F |
2433 |
RNEVP |
RLAR(I) = 0. |
|
./pkg/fizhi/fizhi_moist.F |
2452 |
RNEVP |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_moist.F |
2453 |
RNEVP |
TEMP2(I,L) = PL(I,L)*ZP001 |
|
./pkg/fizhi/fizhi_moist.F |
2454 |
RNEVP |
TEMP2(I,L) = SQRT(TEMP2(I,L)) |
|
./pkg/fizhi/fizhi_moist.F |
2461 |
RNEVP |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_moist.F |
2462 |
RNEVP |
TEMP3(I,L) = GRAVITY*ZP01 / DP(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2470 |
RNEVP |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_moist.F |
2471 |
RNEVP |
TEMP1(I,3) = TL(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2472 |
RNEVP |
TEMP1(I,4) = QL(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2479 |
RNEVP |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_moist.F |
2481 |
RNEVP |
TEMP1(I,5)=TEMP1(I,2)-TEMP1(I,4) |
|
./pkg/fizhi/fizhi_moist.F |
2482 |
RNEVP |
TEMP1(I,6)=TEMP1(I,6)*ELOCP |
|
./pkg/fizhi/fizhi_moist.F |
2483 |
RNEVP |
TEMP1(I,5)=TEMP1(I,5)/(ONE+TEMP1(I,6)) |
|
./pkg/fizhi/fizhi_moist.F |
2484 |
RNEVP |
TEMP1(I,4)=TEMP1(I,4)+TEMP1(I,5)*RELAX |
|
./pkg/fizhi/fizhi_moist.F |
2485 |
RNEVP |
TEMP1(I,3)=TEMP1(I,3)-TEMP1(I,5)*ELOCP*RELAX |
|
./pkg/fizhi/fizhi_moist.F |
2489 |
RNEVP |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_moist.F |
2490 |
RNEVP |
EVP9(I,L) = (TEMP1(I,4) - QL(I,L))/TEMP3(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2505 |
RNEVP |
DO 20 I=1,IRUN |
|
./pkg/fizhi/fizhi_moist.F |
2511 |
RNEVP |
iceevap(I) = EVP9(I,L)*fracice |
|
./pkg/fizhi/fizhi_moist.F |
2513 |
RNEVP |
EVP9(I,L)=EVP9(I,L)-iceevap(I) |
|
./pkg/fizhi/fizhi_moist.F |
2519 |
RNEVP |
exparg = ZM1P04*tmscl*((water(i)*RPHF*TEMP2(I,L))**ZP578) |
|
./pkg/fizhi/fizhi_moist.F |
2521 |
RNEVP |
watevap(I) = EVP9(I,L)*AREARAT*fracwat |
|
./pkg/fizhi/fizhi_moist.F |
2522 |
RNEVP |
IF(watevap(I).GE.water(i)) watevap(I) = water(i) |
|
./pkg/fizhi/fizhi_moist.F |
2523 |
RNEVP |
EVP9(I,L)=EVP9(I,L)-watevap(I) |
|
./pkg/fizhi/fizhi_moist.F |
2527 |
RNEVP |
QL(I,L) = QL(I,L)+(iceevap(i)+watevap(i))*TEMP3(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2528 |
RNEVP |
TL(I,L) = TL(I,L)-(iceevap(i)+watevap(i))*TEMP3(I,L)*ELOCP |
|
./pkg/fizhi/fizhi_moist.F |
2543 |
RNEVP |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_moist.F |
2559 |
RNEVP |
if( pl(i,L).ge.10.0 ) CLSBTH(I,L) = fact |
|
./pkg/fizhi/fizhi_moist.F |
2560 |
RNEVP |
RLAR(I) = RLAR(I)-rainmax(I) |
|
./pkg/fizhi/fizhi_moist.F |
2561 |
RNEVP |
QL(I,L) = QL(I,L)+rainmax(I)*TEMP3(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2562 |
RNEVP |
TL(I,L) = TL(I,L)-rainmax(I)*TEMP3(I,L)*ELOCP |
|
./pkg/fizhi/fizhi_moist.F |
2568 |
RNEVP |
DO I=1,IRUN |
|
./pkg/fizhi/fizhi_moist.F |
2569 |
RNEVP |
IF((RLAR(I).GT.0.0).AND.(rainmax(I).GT.0.0))THEN |
|
./pkg/fizhi/fizhi_moist.F |
2570 |
RNEVP |
RPOW=(RLAR(I)*RPHF*TEMP2(I,L))**ZP578 |
|
./pkg/fizhi/fizhi_moist.F |
2573 |
RNEVP |
TEMP1(I,7) = rainmax(I)*AREARAT |
|
./pkg/fizhi/fizhi_moist.F |
2574 |
RNEVP |
IF(TEMP1(I,7).GE.RLAR(I)) TEMP1(I,7) = RLAR(I) |
|
./pkg/fizhi/fizhi_moist.F |
2575 |
RNEVP |
RLAR(I) = RLAR(I)-TEMP1(I,7) |
|
./pkg/fizhi/fizhi_moist.F |
2576 |
RNEVP |
QL(I,L) = QL(I,L)+TEMP1(I,7)*TEMP3(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
2577 |
RNEVP |
TL(I,L) = TL(I,L)-TEMP1(I,7)*TEMP3(I,L)*ELOCP |
|
./pkg/fizhi/fizhi_moist.F |
372 |
MOISTIO |
DO I = 1,ISTRIP |
|
./pkg/fizhi/fizhi_moist.F |
373 |
MOISTIO |
TH(I,L) = TH(I,L) * P0KAPPA |
|
./pkg/fizhi/fizhi_moist.F |
374 |
MOISTIO |
CLMAXO(I,L) = 0. |
|
./pkg/fizhi/fizhi_moist.F |
375 |
MOISTIO |
CLBOTH(I,L) = 0. |
|
./pkg/fizhi/fizhi_moist.F |
376 |
MOISTIO |
cldmas(I,L) = 0. |
|
./pkg/fizhi/fizhi_moist.F |
377 |
MOISTIO |
detrain(I,L) = 0. |
|
./pkg/fizhi/fizhi_moist.F |
409 |
MOISTIO |
DO 600 I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
410 |
MOISTIO |
TMP1(I,2) = 0. |
|
./pkg/fizhi/fizhi_moist.F |
411 |
MOISTIO |
TMP1(I,3) = 0. |
|
./pkg/fizhi/fizhi_moist.F |
416 |
MOISTIO |
DO 602 I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
417 |
MOISTIO |
TMP1(I,2) = TMP1(I,2) + (PLE(I,L+1)-PLE(I,L))*TH (I,L)/sp(i) |
|
./pkg/fizhi/fizhi_moist.F |
418 |
MOISTIO |
TMP1(I,3) = TMP1(I,3) + (PLE(I,L+1)-PLE(I,L))*SHL(I,L)/sp(i) |
|
./pkg/fizhi/fizhi_moist.F |
421 |
MOISTIO |
DO 603 I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
422 |
MOISTIO |
TMP1(I,4) = 1. / ( (PLE(I,lm+1)-PLE(I,NSUBCL))/sp(I) ) |
|
./pkg/fizhi/fizhi_moist.F |
423 |
MOISTIO |
TH(I,NSUBCL) = TMP1(I,2)*TMP1(I,4) |
|
./pkg/fizhi/fizhi_moist.F |
424 |
MOISTIO |
SHL(I,NSUBCL) = TMP1(I,3)*TMP1(I,4) |
|
./pkg/fizhi/fizhi_moist.F |
435 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
436 |
MOISTIO |
TMP1(I,2) = 0. |
|
./pkg/fizhi/fizhi_moist.F |
439 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
440 |
MOISTIO |
TMP1(I,2) = TMP1(I,2)+(PLE(I,L+1)-PLE(I,L))*UL(I,L,NT)/sp(i) |
|
./pkg/fizhi/fizhi_moist.F |
443 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
444 |
MOISTIO |
UL(I,NSUBCL,NT) = TMP1(I,2)*TMP1(I,4) |
|
./pkg/fizhi/fizhi_moist.F |
452 |
MOISTIO |
DO 112 I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
453 |
MOISTIO |
TMP4(I,L) = PLE(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
456 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
457 |
MOISTIO |
TMP5(I,1) = PTOPKAP / P0KAPPA |
|
./pkg/fizhi/fizhi_moist.F |
460 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
461 |
MOISTIO |
TMP5(I,L) = PLKE(I,L)*P0KINV |
|
./pkg/fizhi/fizhi_moist.F |
464 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
465 |
MOISTIO |
TMP4(I,lm+1) = PLE (I,lm+1) |
|
./pkg/fizhi/fizhi_moist.F |
466 |
MOISTIO |
TMP5(I,lm+1) = PLKE(I,lm+1)*P0KINV |
|
./pkg/fizhi/fizhi_moist.F |
468 |
MOISTIO |
DO 113 I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
469 |
MOISTIO |
TMP4(I,NSUBCL+1) = PLE (I,lm+1) |
|
./pkg/fizhi/fizhi_moist.F |
470 |
MOISTIO |
TMP5(I,NSUBCL+1) = PLKE(I,lm+1)*P0KINV |
|
./pkg/fizhi/fizhi_moist.F |
482 |
MOISTIO |
rh = SHL(I,NSUBCL) / tmp2(i,3) |
|
./pkg/fizhi/fizhi_moist.F |
541 |
MOISTIO |
do I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
550 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
551 |
MOISTIO |
CVTH(I,L) = (TH (I,L) - saveth(i,l)) |
|
./pkg/fizhi/fizhi_moist.F |
552 |
MOISTIO |
CVQ (I,L) = (SHL(I,L) - saveq (i,l)) |
|
./pkg/fizhi/fizhi_moist.F |
556 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
557 |
MOISTIO |
CVTH(I,L) = cvth(i,nsubcl) |
|
./pkg/fizhi/fizhi_moist.F |
558 |
MOISTIO |
CVQ (I,L) = cvq (i,nsubcl) |
|
./pkg/fizhi/fizhi_moist.F |
563 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
564 |
MOISTIO |
TH (I,L) = saveth(i,l) + cvth(i,l) |
|
./pkg/fizhi/fizhi_moist.F |
565 |
MOISTIO |
SHL(I,L) = saveq (i,l) + cvq (i,l) |
|
./pkg/fizhi/fizhi_moist.F |
569 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
570 |
MOISTIO |
CVTH(I,L) = CVTH(I,L) *P0KINV*SP(I)*tminv |
|
./pkg/fizhi/fizhi_moist.F |
571 |
MOISTIO |
CVQ (I,L) = CVQ (I,L) *SP(I)*tminv |
|
./pkg/fizhi/fizhi_moist.F |
579 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
580 |
MOISTIO |
CVU(I,L,nt) = ( UL(I,L,nt)-saveu(i,l,nt) )*sp(i)*tminv |
|
./pkg/fizhi/fizhi_moist.F |
584 |
MOISTIO |
DO I=num,num+nindeces(nsubcl)-1 |
|
./pkg/fizhi/fizhi_moist.F |
645 |
MOISTIO |
DO I=1,ISTRIP |
|
./pkg/fizhi/fizhi_moist.F |
646 |
MOISTIO |
TH(I,L) = TH(I,L) + CVTH(I,L)*tmstp/SP(I) |
|
./pkg/fizhi/fizhi_moist.F |
647 |
MOISTIO |
SHL(I,L) = SHL(I,L) + CVQ(I,L)*tmstp/SP(I) |
|
./pkg/fizhi/fizhi_moist.F |
648 |
MOISTIO |
TL(I,L) = TH(I,L)*PLK(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
649 |
MOISTIO |
saveth(I,L) = th(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
650 |
MOISTIO |
saveq (I,L) = SHL(I,L) |
|
./pkg/fizhi/fizhi_moist.F |
664 |
MOISTIO |
DO I =1,ISTRIP |
|
./pkg/fizhi/fizhi_moist.F |
665 |
MOISTIO |
TMP1(I,L) = sp(i) * (SHL(I,L)-saveq(I,L)) * tminv |
|
./pkg/fizhi/fizhi_moist.F |
669 |
MOISTIO |
DO I =1,ISTRIP |
|
./pkg/fizhi/fizhi_moist.F |
670 |
MOISTIO |
TMP1(I,L) = sp(i) * ((TL(I,L)/PLK(I,L))-saveth(i,l)) * tminv |
|
./pkg/fizhi/fizhi_moist.F |
676 |
MOISTIO |
DO I =1,ISTRIP |
|
./pkg/fizhi/fizhi_moist.F |
677 |
MOISTIO |
TMP1(I,L) = ((TL(I,L)/PLK(I,L))-saveth(i,l))*plk(i,l)*sday*tminv |
|
./pkg/fizhi/fizhi_moist.F |
681 |
MOISTIO |
DO I =1,ISTRIP |
|
./pkg/fizhi/fizhi_moist.F |
682 |
MOISTIO |
TMP1(I,L) = (SHL(I,L)-saveq(I,L)) * 1000. * sday * tminv |
|
./pkg/fizhi/fizhi_moist.F |
728 |
MOISTIO |
prep(i) = PRECIP(I) + PCNET(I) |
|
./pkg/fizhi/fizhi_moist.F |
729 |
MOISTIO |
tmp1(i,1) = PRECIP(I) |
|
./pkg/fizhi/fizhi_moist.F |
981 |
MOISTIO |
tmpdiag(i,j) = (lsp_new(I,j) + snow_new(I,j) + conv_new(i,j)) |
|
./pkg/fizhi/fizhi_swrad.F |
152 |
SWRIO |
zenith(I,j) = TEMP1(I,j) + TEMP2(I,j) |
|
./pkg/fizhi/fizhi_swrad.F |
153 |
SWRIO |
IF( zenith(I,j) .GT. 0.0 ) THEN |
|
./pkg/fizhi/fizhi_swrad.F |
154 |
SWRIO |
zenith(I,j) = (2./3.)*( zenith(I,j)-TEMP2(I,j)*TEMP1(I,j) |
|
./pkg/fizhi/fizhi_swrad.F |
155 |
SWRIO |
. / zenith(I,j) ) |
|
./pkg/fizhi/fizhi_turb.F |
1473 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1474 |
TRBFLX |
wu(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1475 |
TRBFLX |
wv(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1476 |
TRBFLX |
eu(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1477 |
TRBFLX |
et(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1482 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1483 |
TRBFLX |
XLDIAG(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1488 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1489 |
TRBFLX |
FREQDG(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1492 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1495 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1498 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1501 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1504 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1510 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1513 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1516 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1519 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1522 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1525 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1528 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1531 |
TRBFLX |
do I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1537 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1538 |
TRBFLX |
XLSAVE(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1539 |
TRBFLX |
KHSAVE(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1542 |
TRBFLX |
DO I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1543 |
TRBFLX |
CTSAVE(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1544 |
TRBFLX |
XXSAVE(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1545 |
TRBFLX |
YYSAVE(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1546 |
TRBFLX |
ZETASAVE(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1553 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1554 |
TRBFLX |
ADZ1(I,L) = (CP/GRAV)*(PLKE(I,L+1)-PLKE(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
1555 |
TRBFLX |
ADZ1(I,L) = THV(I,L) * ADZ1(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1556 |
TRBFLX |
DZ1TMP(I,L) = ADZ1(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1560 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1561 |
TRBFLX |
ADZ2(I,L) = 0.5 * (ADZ1(I,L)+ADZ1(I,L+1)) |
|
./pkg/fizhi/fizhi_turb.F |
1566 |
TRBFLX |
DO 9042 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1567 |
TRBFLX |
HS(I) = 0.5 * ADZ1(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1571 |
TRBFLX |
DO 9044 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1572 |
TRBFLX |
DZ3(I,1) = HALPHA * ADZ1(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
1575 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1576 |
TRBFLX |
DZ3(I,L) = ALPHA * ADZ1(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1579 |
TRBFLX |
DO 9048 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1580 |
TRBFLX |
DZ3(I,NLEV) = ALPHA * HS(I) |
|
./pkg/fizhi/fizhi_turb.F |
1586 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1587 |
TRBFLX |
TEMP(I,L) = VK * ADZ1(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1590 |
TRBFLX |
DO 9052 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1591 |
TRBFLX |
VKZE(I,NLEVM1) = TEMP(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1596 |
TRBFLX |
DO 9054 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1597 |
TRBFLX |
VKZE(I,L) = VKZE(I,LP1) + TEMP(I,LP1) |
|
./pkg/fizhi/fizhi_turb.F |
1601 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1602 |
TRBFLX |
VKZM(I,L) = VKZE(I,L) - 0.5 * TEMP(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
1609 |
TRBFLX |
DO 9058 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1610 |
TRBFLX |
FAC1 = DPSTR(I,L) / ( DPSTR(I,L) + DPSTR(I,LP1) ) |
|
./pkg/fizhi/fizhi_turb.F |
1612 |
TRBFLX |
RHODZ2(I,L) = FAC1 * THV(I,LP1) |
|
./pkg/fizhi/fizhi_turb.F |
1613 |
TRBFLX |
RHODZ2(I,L) = RHODZ2(I,L) + FAC2 * THV(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1617 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1618 |
TRBFLX |
RHODZ2(I,L) = (RGAS*0.01) * RHODZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1619 |
TRBFLX |
TEMP(I,L) = PLKE(I,L+1) * ADZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1620 |
TRBFLX |
RHODZ2(I,L) = TEMP(I,L) * RHODZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1621 |
TRBFLX |
RHODZ2(I,L) = PLE(I,L+1) / RHODZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1622 |
TRBFLX |
RHOZPK(I,L) = RHODZ2(I,L) * PLKE(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
1623 |
TRBFLX |
RHODZ1(I,L) = (RGAS*0.01) * THV(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
1624 |
TRBFLX |
TEMP(I,L) = PLK(I,L+1) * ADZ1(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
1625 |
TRBFLX |
RHODZ1(I,L) = TEMP(I,L) * RHODZ1(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1626 |
TRBFLX |
RHODZ1(I,L) = PL(I,L+1) / RHODZ1(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1633 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1634 |
TRBFLX |
FLXFPK(I,L) = PLE(I,L+1) - PLE(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1635 |
TRBFLX |
FLXFPK(I,L) = FLXFPK(I,L) * PLK(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1636 |
TRBFLX |
FLXFPK(I,L) = (GRAV*DTAU*0.01) / FLXFPK(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1639 |
TRBFLX |
DO 9064 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1640 |
TRBFLX |
FLXFPK(I,NLEVP1) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1642 |
TRBFLX |
DO 9066 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1643 |
TRBFLX |
IF (IWATER(I).EQ.0 ) FLXFPK(I,NLEVP1) = 1. / PLKE(I,NLEVP1) |
|
./pkg/fizhi/fizhi_turb.F |
1646 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1647 |
TRBFLX |
FLXFAC(I,L) = FLXFPK(I,L) * PLK(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1650 |
TRBFLX |
DO 9070 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1651 |
TRBFLX |
FLXFAC(I,NLEVP1) = FLXFPK(I,NLEVP1) |
|
./pkg/fizhi/fizhi_turb.F |
1653 |
TRBFLX |
DO 9074 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1654 |
TRBFLX |
FLXFPK(I,NLEVP1) = CP * FLXFPK(I,NLEVP1) |
|
./pkg/fizhi/fizhi_turb.F |
1657 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1658 |
TRBFLX |
FLXFCE(I,L) = PL(I,L+1) - PL(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1659 |
TRBFLX |
FLXFCE(I,L) = (GRAV*DTAU*0.01) / FLXFCE(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1665 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1666 |
TRBFLX |
ADZ1(I,L) = 1. / ADZ1(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1670 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1671 |
TRBFLX |
ADZ2(I,L) = 1. / ADZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1674 |
TRBFLX |
DO 9088 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1675 |
TRBFLX |
AHS(I) = 1. / HS(I) |
|
./pkg/fizhi/fizhi_turb.F |
1680 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1681 |
TRBFLX |
DPK(I,L) = ( PLK(I,L+1)-PLK(I,L) ) * ADZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1684 |
TRBFLX |
DO 9092 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1685 |
TRBFLX |
DPK(I,NLEV) = GBYCP / THV(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1690 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1691 |
TRBFLX |
Q(I,L) = 2. * QQ(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1692 |
TRBFLX |
Q(I,L) = SQRT( Q(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
1708 |
TRBFLX |
DO I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1709 |
TRBFLX |
CT(I) = CTSAVE(I) |
|
./pkg/fizhi/fizhi_turb.F |
1710 |
TRBFLX |
XX(I) = XXSAVE(I) |
|
./pkg/fizhi/fizhi_turb.F |
1711 |
TRBFLX |
YY(I) = YYSAVE(I) |
|
./pkg/fizhi/fizhi_turb.F |
1712 |
TRBFLX |
ZETA(I) = ZETASAVE(I) |
|
./pkg/fizhi/fizhi_turb.F |
1716 |
TRBFLX |
DO I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1717 |
TRBFLX |
TL(I,NLEV) = TH(I,NLEV)*PLK(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1721 |
TRBFLX |
DO I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1722 |
TRBFLX |
BB(I,NLEV) = FACEPS*SHSAT(I,NLEV)/(TL(I,NLEV)*TL(I,NLEV)) |
|
./pkg/fizhi/fizhi_turb.F |
1723 |
TRBFLX |
AA(I,NLEV) = 1. / (1. + CLH * BB(I,NLEV) ) |
|
./pkg/fizhi/fizhi_turb.F |
1724 |
TRBFLX |
BB(I,NLEV) = BB(I,NLEV) * AA(I,NLEV) * plk(I,nlev) |
|
./pkg/fizhi/fizhi_turb.F |
1725 |
TRBFLX |
DTH(I,NLEV) = TH(I,NLEV)-TH(I,NLEVP1) |
|
./pkg/fizhi/fizhi_turb.F |
1726 |
TRBFLX |
DSH(I,NLEV) = SH(I,NLEV)-SH(I,NLEVP1) |
|
./pkg/fizhi/fizhi_turb.F |
1727 |
TRBFLX |
SBAR(I,NLEV) = AA(I,NLEV) * (SH(I,NLEV) - SHSAT(I,NLEV)) |
|
./pkg/fizhi/fizhi_turb.F |
1728 |
TRBFLX |
SSDEV(I,NLEV)=CT(I)*(AA(I,NLEV)*DSH(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1729 |
TRBFLX |
1 -BB(I,NLEV)*DTH(I,NLEV)) |
|
./pkg/fizhi/fizhi_turb.F |
1730 |
TRBFLX |
XXZETA(I) = XX(I)-ZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
1731 |
TRBFLX |
IF(XXZETA(I).LT.0.1*XX(I)) XXZETA(I)=0.1*XX(I) |
|
./pkg/fizhi/fizhi_turb.F |
1732 |
TRBFLX |
IF(XXZETA(I).LE.0.) XXZETA(I)=0.1 |
|
./pkg/fizhi/fizhi_turb.F |
1733 |
TRBFLX |
QBYU(I) =QBUSTR * XXZETA(I) ** ONETHRD |
|
./pkg/fizhi/fizhi_turb.F |
1734 |
TRBFLX |
SSDEV(I,NLEV) = B2*YY(I)*SSDEV(I,NLEV)*SSDEV(I,NLEV)/QBYU(I) |
|
./pkg/fizhi/fizhi_turb.F |
1735 |
TRBFLX |
SVAR(I,NLEV) = SQRT(SSDEV(I,NLEV)) |
|
./pkg/fizhi/fizhi_turb.F |
1736 |
TRBFLX |
IF ( SVAR(I,NLEV).LT.Z1PEM25) SVAR(I,NLEV) = Z1PEM25 |
|
./pkg/fizhi/fizhi_turb.F |
1737 |
TRBFLX |
Q1M(I,NLEV) = SBAR(I,NLEV) / SVAR(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1738 |
TRBFLX |
FCC(I,NLEV) = (1./2.) * ( 1. + ERRF( P5SR*Q1M(I,NLEV) ) ) |
|
./pkg/fizhi/fizhi_turb.F |
1739 |
TRBFLX |
SHL(I,NLEV) = FCC(I,NLEV) * SBAR(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1740 |
TRBFLX |
ARG(I,NLEV) = (1./2.)*Q1M(I,NLEV)*Q1M(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1741 |
TRBFLX |
IF(ARG(I,NLEV).LE.ARGMAX) |
|
./pkg/fizhi/fizhi_turb.F |
1742 |
TRBFLX |
1 SHL(I,NLEV) = SHL(I,NLEV)+RSQ2PI*SVAR(I,NLEV)*EXP(-ARG(I,NLEV)) |
|
./pkg/fizhi/fizhi_turb.F |
1743 |
TRBFLX |
BETAT(I,NLEV) = 1. + VIRTCON*SH(I,NLEV) - VRT1CON*SHL(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1744 |
TRBFLX |
BETAW(I,NLEV) = VIRTCON * |
|
./pkg/fizhi/fizhi_turb.F |
1745 |
TRBFLX |
1 ( TH(I,NLEV) + CLH * SHL(I,NLEV) * (1./plk(i,nlev)) ) |
|
./pkg/fizhi/fizhi_turb.F |
1746 |
TRBFLX |
BETAL(I,NLEV) = (1.+VIRTCON*SH(I,NLEV)-TWO*VRT1CON*SHL(I,NLEV)) |
|
./pkg/fizhi/fizhi_turb.F |
1747 |
TRBFLX |
1 * (1./plk(i,nlev)) * CLH - VRT1CON * TH(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1748 |
TRBFLX |
BETAT1(I,NLEV) = BETAT(I,NLEV) - BB(I,NLEV)*FCC(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1749 |
TRBFLX |
1 * BETAL(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1750 |
TRBFLX |
BETAW1(I,NLEV) = BETAW(I,NLEV) + AA(I,NLEV) * FCC(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1751 |
TRBFLX |
1 * BETAL(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1752 |
TRBFLX |
DTHV(I,NLEV) = BETAT1(I,NLEV)*DTH(I,NLEV) + |
|
./pkg/fizhi/fizhi_turb.F |
1753 |
TRBFLX |
1 BETAW1(I,NLEV)*DSH(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1754 |
TRBFLX |
THV(I,NLEVP1) = THV(I,NLEV) - DTHV(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1785 |
TRBFLX |
DO 9098 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1786 |
TRBFLX |
Q(I,NLEV) = QBUSTR * USTAR(I) |
|
./pkg/fizhi/fizhi_turb.F |
1787 |
TRBFLX |
QQ(I,NLEV) = 0.5 * Q(I,NLEV) * Q(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1794 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1795 |
TRBFLX |
DU(I,L) = ( U(I,L)- U(I,L+1) ) * ADZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1796 |
TRBFLX |
DV(I,L) = ( V(I,L)- V(I,L+1) ) * ADZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1805 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1806 |
TRBFLX |
XL(I,L) = XLSAVE(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1812 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1813 |
TRBFLX |
DTH(I,L) = ( TH(I,L)-TH(I,L+1) ) * ADZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1814 |
TRBFLX |
DSH(I,L) = ( SH(I,L)-SH(I,L+1) ) * ADZ2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1815 |
TRBFLX |
TL(I,L) = TH(I,L)*PLK(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1819 |
TRBFLX |
DO I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1824 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1825 |
TRBFLX |
BB(I,L) = FACEPS*SHSAT(I,L)/(TL(I,L)*TL(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
1826 |
TRBFLX |
AA(I,L) = 1. / (1. + CLH * BB(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
1828 |
TRBFLX |
BB(I,L) = BB(I,L) * AA(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1829 |
TRBFLX |
SBAR(I,L) = AA(I,L) * (SH(I,L) - SHSAT(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
1832 |
TRBFLX |
DO I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1834 |
TRBFLX |
SSDEV(I,1) = XL(I,1)*(AA(I,1)*DSH(I,1)- |
|
./pkg/fizhi/fizhi_turb.F |
1835 |
TRBFLX |
1 BB(I,1)*plke(I,2)*DTH(I,1)) |
|
./pkg/fizhi/fizhi_turb.F |
1836 |
TRBFLX |
SSDEV(I,1) = B2 * KHSAVE(I,1) * SSDEV(I,1) * SSDEV(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
1837 |
TRBFLX |
SVAR(I,1) = SQRT(SSDEV(I,1)) |
|
./pkg/fizhi/fizhi_turb.F |
1838 |
TRBFLX |
IF ( SVAR(I,1).LT.Z1PEM25) SVAR(I,1) = Z1PEM25 |
|
./pkg/fizhi/fizhi_turb.F |
1841 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1843 |
TRBFLX |
SSDEV(I,L) = XL(I,L-1)*(AA(I,L)*DSH(I,L-1)- |
|
./pkg/fizhi/fizhi_turb.F |
1844 |
TRBFLX |
1 BB(I,L)*plke(I,L)*DTH(I,L-1)) |
|
./pkg/fizhi/fizhi_turb.F |
1845 |
TRBFLX |
SSDEV(I,L) = B2 * KHSAVE(I,L-1) * SSDEV(I,L-1) * SSDEV(I,L-1) |
|
./pkg/fizhi/fizhi_turb.F |
1846 |
TRBFLX |
SVAR(I,L) = SQRT(SSDEV(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
1848 |
TRBFLX |
SSDEV(I,L) = XL(I,L)*(AA(I,L)*DSH(I,L)- |
|
./pkg/fizhi/fizhi_turb.F |
1849 |
TRBFLX |
1 BB(I,L)*plke(I,L+1)*DTH(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
1850 |
TRBFLX |
SSDEV(I,L) = B2 * KHSAVE(I,L) * SSDEV(I,L) * SSDEV(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1851 |
TRBFLX |
TEMP(I,L) = SQRT(SSDEV(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
1852 |
TRBFLX |
SVAR(I,L) = (1./2.) * (SVAR(I,L) + TEMP(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
1853 |
TRBFLX |
IF ( SVAR(I,L).LT.Z1PEM25) SVAR(I,L) = Z1PEM25 |
|
./pkg/fizhi/fizhi_turb.F |
1857 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1858 |
TRBFLX |
Q1M(I,L) = SBAR(I,L) / SVAR(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1859 |
TRBFLX |
FCC(I,L) = (1./2.) * ( 1. + ERRF( P5SR*Q1M(I,L) ) ) |
|
./pkg/fizhi/fizhi_turb.F |
1860 |
TRBFLX |
SHL(I,L) = FCC(I,L) * SBAR(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1861 |
TRBFLX |
ARG(I,L) = (1./2.)*Q1M(I,L)*Q1M(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1862 |
TRBFLX |
IF(ARG(I,L).LE.ARGMAX) |
|
./pkg/fizhi/fizhi_turb.F |
1863 |
TRBFLX |
1 SHL(I,L) = SHL(I,L)+RSQ2PI*SVAR(I,L)*EXP(-ARG(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
1864 |
TRBFLX |
BETAT(I,L) = 1. + VIRTCON * SH(I,L) - VRT1CON * SHL(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1865 |
TRBFLX |
BETAW(I,L) = VIRTCON * |
|
./pkg/fizhi/fizhi_turb.F |
1866 |
TRBFLX |
1 ( TH(I,L) + (CLH/plk(I,L)) * SHL(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
1867 |
TRBFLX |
BETAL(I,L) = ( 1. + VIRTCON*SH(I,L) - TWO*VRT1CON*SHL(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
1868 |
TRBFLX |
1 * (CLH/plke(I,L+1)) - VRT1CON * TH(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1870 |
TRBFLX |
BETAT1(I,L) = BETAT(I,L) - |
|
./pkg/fizhi/fizhi_turb.F |
1871 |
TRBFLX |
1 BB(I,L)*plk(i,L) * FCC(I,L) * BETAL(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1872 |
TRBFLX |
BETAW1(I,L) = BETAW(I,L) + AA(I,L) * FCC(I,L) * BETAL(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1876 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1877 |
TRBFLX |
DTHV(I,L) = (1./2.)*((BETAT1(I,L)+BETAT1(I,L+1))*DTH(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1878 |
TRBFLX |
1 + (BETAW1(I,L)+BETAW1(I,L+1))*DSH(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
1884 |
TRBFLX |
DO 9102 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1885 |
TRBFLX |
DU(I,NLEV) = CU(I)*XX(I)*AHS(I)*RVK |
|
./pkg/fizhi/fizhi_turb.F |
1886 |
TRBFLX |
DV(I,NLEV) = V(I,NLEV) * DU(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1887 |
TRBFLX |
DU(I,NLEV) = U(I,NLEV) * DU(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
1888 |
TRBFLX |
DTHV(I,NLEV) = CT(I) * YY(I) * |
|
./pkg/fizhi/fizhi_turb.F |
1889 |
TRBFLX |
1 ((THV(I,NLEV)-THV(I,NLEVP1)) * RVK)* AHS(I) |
|
./pkg/fizhi/fizhi_turb.F |
1895 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1896 |
TRBFLX |
STRT(I,L) = CP * DTHV(I,L) * DPK(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1897 |
TRBFLX |
DW2(I,L) = DU(I,L) * DU(I,L) + DV(I,L) * DV(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1898 |
TRBFLX |
IF ( DW2(I,L) .LE. 1.e-4 ) DW2(I,L) = 1.e-4 |
|
./pkg/fizhi/fizhi_turb.F |
1899 |
TRBFLX |
RI(I,L) = STRT(I,L) / DW2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1906 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1907 |
TRBFLX |
SRI(I,L) = RI(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1910 |
TRBFLX |
DO 9108 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1911 |
TRBFLX |
SRI(I,NLEV) = RIB(I) |
|
./pkg/fizhi/fizhi_turb.F |
1913 |
TRBFLX |
DO 9110 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1914 |
TRBFLX |
SWINDS(I) = WS(I) |
|
./pkg/fizhi/fizhi_turb.F |
1919 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1920 |
TRBFLX |
KH(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1921 |
TRBFLX |
KM(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1922 |
TRBFLX |
QQE(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1923 |
TRBFLX |
QE(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1924 |
TRBFLX |
P3(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1925 |
TRBFLX |
IBITSTB(I,L) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
1926 |
TRBFLX |
IF ( QQ(I,L) .GT. 1.e-8 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
1927 |
TRBFLX |
INTQ(I,L) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
1929 |
TRBFLX |
INTQ(I,L) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
1931 |
TRBFLX |
IF ( QQ(I,L).LE.1.e-8 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
1932 |
TRBFLX |
QQ(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1933 |
TRBFLX |
Q(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
1940 |
TRBFLX |
DO 9116 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1941 |
TRBFLX |
IF ( QQ(I,LMINQ).GT.1.e-8 ) IBIT = IBIT + 1 |
|
./pkg/fizhi/fizhi_turb.F |
1967 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1968 |
TRBFLX |
QQ(I,L) = QQE(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1969 |
TRBFLX |
Q(I,L) = QE(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1984 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1985 |
TRBFLX |
IF ( (IBITSTB(I,L).EQ.1) .AND. |
|
./pkg/fizhi/fizhi_turb.F |
1986 |
TRBFLX |
1 ( Q(I,L) .LE. QE(I,L) ) ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
1987 |
TRBFLX |
IBITSTB(I,L) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
1989 |
TRBFLX |
IBITSTB(I,L) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
1994 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
1995 |
TRBFLX |
IF(IBITSTB(I,L).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
1996 |
TRBFLX |
TEMP(I,L) = Q(I,L) / QE(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1997 |
TRBFLX |
KH(I,L) = TEMP(I,L) * KH(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
1998 |
TRBFLX |
KM(I,L) = TEMP(I,L) * KM(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2000 |
TRBFLX |
TEMP(I,L) = 0.01 * QQE(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2001 |
TRBFLX |
IF((IBITSTB(I,L).EQ.1) .AND. |
|
./pkg/fizhi/fizhi_turb.F |
2002 |
TRBFLX |
1 ( QQ(I,L) .LE. TEMP(I,L) )) THEN |
|
./pkg/fizhi/fizhi_turb.F |
2003 |
TRBFLX |
QQ(I,L) = TEMP(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2004 |
TRBFLX |
Q(I,L) = 0.1 * QE(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2006 |
TRBFLX |
IF(IBITSTB(I,L).EQ.1 ) P3(I,L) = (2.*B3) * |
|
./pkg/fizhi/fizhi_turb.F |
2007 |
TRBFLX |
1 ( QE(I,L) - Q(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
2025 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2026 |
TRBFLX |
P3(I,L) = P3(I,L) * DTAU / XL(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2027 |
TRBFLX |
TEMP(I,L) = QQE(I,L) * P3(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2028 |
TRBFLX |
XQ(I,L) = QQE(I,L) - QQ(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2032 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2033 |
TRBFLX |
IF( ( (IBITSTB(I,L).EQ.1) .AND. |
|
./pkg/fizhi/fizhi_turb.F |
2034 |
TRBFLX |
1 ( XQ(I,L) .LT. TEMP(I,L) ) ) |
|
./pkg/fizhi/fizhi_turb.F |
2036 |
TRBFLX |
3 ( (IBITSTB(I,L).EQ.0) .AND. |
|
./pkg/fizhi/fizhi_turb.F |
2037 |
TRBFLX |
4 ( XQ(I,L) .GT. TEMP(I,L) ) ) ) |
|
./pkg/fizhi/fizhi_turb.F |
2038 |
TRBFLX |
5 P3(I,L) = XQ(I,L) / QQE(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2046 |
TRBFLX |
DO 9118 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2047 |
TRBFLX |
RIBIN(I) = RIB(I) |
|
./pkg/fizhi/fizhi_turb.F |
2048 |
TRBFLX |
CUIN(I) = CU(I) |
|
./pkg/fizhi/fizhi_turb.F |
2049 |
TRBFLX |
CTIN(I) = CT(I) |
|
./pkg/fizhi/fizhi_turb.F |
2050 |
TRBFLX |
USTARIN(I) = USTAR(I) |
|
./pkg/fizhi/fizhi_turb.F |
2051 |
TRBFLX |
RHOSIN(I) = RHODZ2(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
2052 |
TRBFLX |
Z0IN(I) = Z0(I) |
|
./pkg/fizhi/fizhi_turb.F |
2053 |
TRBFLX |
ZETAIN(I) = ZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
2056 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2057 |
TRBFLX |
RIINIT(I,L) = RI(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2058 |
TRBFLX |
QQINIT(I,L) = QQ(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2080 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2081 |
TRBFLX |
RHOKDZ(I,L) = RHODZ1(I,L) * QXLM(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2090 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2091 |
TRBFLX |
KHSAVE(I,L) = KH(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2101 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2102 |
TRBFLX |
KM(I,L) = KMBG |
|
./pkg/fizhi/fizhi_turb.F |
2103 |
TRBFLX |
KH(I,L) = KHBG |
|
./pkg/fizhi/fizhi_turb.F |
2111 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2112 |
TRBFLX |
Q(I,L) = 2. * QQ(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2113 |
TRBFLX |
Q(I,L) = SQRT(Q(I,L)) |
|
./pkg/fizhi/fizhi_turb.F |
2114 |
TRBFLX |
XQ(I,L)= XL(I,L) * Q(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2115 |
TRBFLX |
KM(I,L)=XQ(I,L)*KM(I,L)+KMBG |
|
./pkg/fizhi/fizhi_turb.F |
2116 |
TRBFLX |
KH(I,L)=XQ(I,L)*KH(I,L)+KHBG |
|
./pkg/fizhi/fizhi_turb.F |
2137 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2138 |
TRBFLX |
TEMP(I,L) = RHOZPK(I,L) * KH(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2139 |
TRBFLX |
DELTH(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
2142 |
TRBFLX |
DO 9132 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2143 |
TRBFLX |
DELTH(I,NLEVP1) = 1. |
|
./pkg/fizhi/fizhi_turb.F |
2162 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2163 |
TRBFLX |
RHOKDZ(I,L) = RHODZ2(I,L) * KH(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2164 |
TRBFLX |
DELSH(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
2167 |
TRBFLX |
DO 9140 I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2168 |
TRBFLX |
DELSH(I,NLEVP1) = 1. |
|
./pkg/fizhi/fizhi_turb.F |
2205 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2206 |
TRBFLX |
RHOKDZ(I,L) = RHODZ2(I,L) * KM(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2212 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2213 |
TRBFLX |
WU(I,L) = WU(I,L) + RHOKDZ(I,L) * ( U(I,L+1) - U(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
2214 |
TRBFLX |
WV(I,L) = WV(I,L) + RHOKDZ(I,L) * ( V(I,L+1) - V(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
2218 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2219 |
TRBFLX |
IF ( QQ(I,L) .GT. QQMIN ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
2220 |
TRBFLX |
IBITSTB(I,L) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
2222 |
TRBFLX |
IBITSTB(I,L) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
2224 |
TRBFLX |
IF( IBITSTB(I,L).EQ.1 ) FREQDG(I,L) = FREQDG(I,L) + aitr |
|
./pkg/fizhi/fizhi_turb.F |
2233 |
TRBFLX |
DO I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2234 |
TRBFLX |
IF ( (qq(i,l).gt.qqcolmax(I)).and.(levpbl(i).eq.nlev))then |
|
./pkg/fizhi/fizhi_turb.F |
2236 |
TRBFLX |
qqcolmin(i) = 0.1*qqcolmax(I) |
|
./pkg/fizhi/fizhi_turb.F |
2245 |
TRBFLX |
pbldpth(I) = pbldpth(I) + ( (PLE(I,nlev+1)-PLE(I,Lp+2)) + |
|
./pkg/fizhi/fizhi_turb.F |
2249 |
TRBFLX |
pbldpth(I) = pbldpth(I) + ( (PLE(I,nlev+1)-PLE(I,2)) + |
|
./pkg/fizhi/fizhi_turb.F |
2261 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2262 |
TRBFLX |
EU(I,L) = EU(I,L) + AITR*KM(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2263 |
TRBFLX |
ET(I,L) = ET(I,L) + AITR*KH(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2266 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2267 |
TRBFLX |
scu(I) = scu(I) + AITR*cu(I) |
|
./pkg/fizhi/fizhi_turb.F |
2269 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2270 |
TRBFLX |
sct(I) = sct(I) + AITR*ct(I) |
|
./pkg/fizhi/fizhi_turb.F |
2274 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2275 |
TRBFLX |
XLDIAG(I,L) = XLDIAG(I,L) + AITR*XL(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2285 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2286 |
TRBFLX |
XLSAVE(I,L) = XL(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
2289 |
TRBFLX |
DO I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2290 |
TRBFLX |
CTSAVE(I) = CT(I) |
|
./pkg/fizhi/fizhi_turb.F |
2291 |
TRBFLX |
XXSAVE(I) = XX(I) |
|
./pkg/fizhi/fizhi_turb.F |
2292 |
TRBFLX |
YYSAVE(I) = YY(I) |
|
./pkg/fizhi/fizhi_turb.F |
2293 |
TRBFLX |
ZETASAVE(I) = ZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
2298 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2299 |
TRBFLX |
turbfcc(I,L) = turbfcc(I,L) + fcc(I,L) * aitr |
|
./pkg/fizhi/fizhi_turb.F |
2300 |
TRBFLX |
qliq(I,L) = qliq(I,L) + shl(I,L) * aitr |
|
./pkg/fizhi/fizhi_turb.F |
2308 |
TRBFLX |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
2309 |
TRBFLX |
WU(I,L) = WU(I,L) * AITR |
|
./pkg/fizhi/fizhi_turb.F |
2310 |
TRBFLX |
WV(I,L) = WV(I,L) * AITR |
|
./pkg/fizhi/fizhi_turb.F |
2426 |
SFCFLX |
DO 9000 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2427 |
SFCFLX |
VWS(I) = VUS(I) * VUS(I) + VVS(I) * VVS(I) |
|
./pkg/fizhi/fizhi_turb.F |
2428 |
SFCFLX |
IF ( VWS(I) .LE. 1.e-4) VWS(I) = 1.e-4 |
|
./pkg/fizhi/fizhi_turb.F |
2429 |
SFCFLX |
VRIB(I) = ( CP * (VPKE(I)-VPK(I)) ) * |
|
./pkg/fizhi/fizhi_turb.F |
2430 |
SFCFLX |
1 (VTHV1(I) - VTHV2(I)) / VWS(I) |
|
./pkg/fizhi/fizhi_turb.F |
2431 |
SFCFLX |
VWS(I) = SQRT( VWS(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2440 |
SFCFLX |
DO 9002 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2441 |
SFCFLX |
IF (IVWATER(I).EQ.1) IWATER = IWATER + 1 |
|
./pkg/fizhi/fizhi_turb.F |
2447 |
SFCFLX |
DO 9004 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2448 |
SFCFLX |
IF (IVWATER(I).EQ.1) VZ0(I) = 0.0003 |
|
./pkg/fizhi/fizhi_turb.F |
2458 |
SFCFLX |
DO 9006 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2459 |
SFCFLX |
VHZ(I) = VHS(I) / VZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
2460 |
SFCFLX |
VPSIM(I) = LOG( VHZ(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2461 |
SFCFLX |
VAPSIM(I) = 1. / VPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
2462 |
SFCFLX |
VCU(I) = VK * VAPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
2463 |
SFCFLX |
VUSTAR(I) = VCU(I) * VWS(I) |
|
./pkg/fizhi/fizhi_turb.F |
2465 |
SFCFLX |
VPSIG(I) = VH0(I) * VUSTAR(I) - USTH0S |
|
./pkg/fizhi/fizhi_turb.F |
2466 |
SFCFLX |
if(VPSIG(I).lt.0.) VPSIG(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
2467 |
SFCFLX |
VPSIG(I) = SQRT( VPSIG(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2468 |
SFCFLX |
VPSIG(I) = BMDL * VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2469 |
SFCFLX |
VPSIHG(I) = VPSIM(I) + VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2475 |
SFCFLX |
DO 9008 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2476 |
SFCFLX |
VTEMP(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
2488 |
SFCFLX |
DO 9010 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2489 |
SFCFLX |
IF ( IVWATER(I).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
2490 |
SFCFLX |
VCU(I) = VCU(I) * (1. - VDPSIM(I)*VAPSIM(I)) |
|
./pkg/fizhi/fizhi_turb.F |
2491 |
SFCFLX |
VZ0(I) = VZ0(I) + VDZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
2493 |
SFCFLX |
IF ( IVWATER(I).EQ.1) THEN |
|
./pkg/fizhi/fizhi_turb.F |
2494 |
SFCFLX |
IF ( VZ0(I) .LE. Z0MIN ) VZ0(I) = Z0MIN |
|
./pkg/fizhi/fizhi_turb.F |
2496 |
SFCFLX |
VPSIG(I) = VH0(I) * VCU(I) * VWS(I) - USTH0S |
|
./pkg/fizhi/fizhi_turb.F |
2497 |
SFCFLX |
if(VPSIG(I).lt.0.) VPSIG(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
2498 |
SFCFLX |
VPSIG(I) = SQRT( VPSIG(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2499 |
SFCFLX |
VPSIG(I) = BMDL * VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2500 |
SFCFLX |
VPSIHG(I) = VPSIM(I) + VDPSIH(I) + VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2507 |
SFCFLX |
DO 9012 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2508 |
SFCFLX |
VZETA(I) = VK2 * VRIB(I) / (VCU(I) * VCU(I) * VPSIHG(I)) |
|
./pkg/fizhi/fizhi_turb.F |
2513 |
SFCFLX |
DO 9014 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2514 |
SFCFLX |
VZH(I) = VZ0(I) * VAHS(I) |
|
./pkg/fizhi/fizhi_turb.F |
2519 |
SFCFLX |
DO 9016 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2520 |
SFCFLX |
VCU(I) = VK / VPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
2521 |
SFCFLX |
VPSIG(I) = VH0(I) * VCU(I) * VWS(I) - USTH0S |
|
./pkg/fizhi/fizhi_turb.F |
2522 |
SFCFLX |
if(VPSIG(I).lt.0.) VPSIG(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
2523 |
SFCFLX |
VPSIG(I) = SQRT(VPSIG(I)) |
|
./pkg/fizhi/fizhi_turb.F |
2524 |
SFCFLX |
VPSIG(I) = BMDL * VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2525 |
SFCFLX |
VPSIHG(I) = VPSIM(I) + VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2526 |
SFCFLX |
VZETA(I) = VK2 * VRIB(I) / (VCU(I) * VCU(I) * VPSIHG(I)) |
|
./pkg/fizhi/fizhi_turb.F |
2531 |
SFCFLX |
DO 9018 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2532 |
SFCFLX |
IF (IVWATER(I).EQ.1) VUSTAR(I) = VCU(I) * VWS(I) |
|
./pkg/fizhi/fizhi_turb.F |
2535 |
SFCFLX |
DO 9020 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2536 |
SFCFLX |
IF (IVWATER(I).EQ.1 ) then |
|
./pkg/fizhi/fizhi_turb.F |
2537 |
SFCFLX |
VZ0(I) = VTEMP(I) |
|
./pkg/fizhi/fizhi_turb.F |
2538 |
SFCFLX |
IF ( VZ0(I) .LE. Z0MIN ) VZ0(I) = Z0MIN |
|
./pkg/fizhi/fizhi_turb.F |
2561 |
SFCFLX |
DO 9022 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2562 |
SFCFLX |
VZETA(I) = VZETA(I) + VZETAL(I) * VDZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
2563 |
SFCFLX |
IF (IVBITRIB(I).EQ.1 )VZETA(I) = |
|
./pkg/fizhi/fizhi_turb.F |
2564 |
SFCFLX |
1 VPSIM(I) * VPSIM(I) * VRIB(I) * VCT(I) * RVK |
|
./pkg/fizhi/fizhi_turb.F |
2568 |
SFCFLX |
DO 9024 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2569 |
SFCFLX |
IF (IVWATER(I).EQ.1 ) then |
|
./pkg/fizhi/fizhi_turb.F |
2570 |
SFCFLX |
VZ0(I) = VZ0(I) + VZ0L(I) * VDZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
2571 |
SFCFLX |
IF (VZ0(I) .LE. Z0MIN ) VZ0(I) = Z0MIN |
|
./pkg/fizhi/fizhi_turb.F |
2583 |
SFCFLX |
DO 9026 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2584 |
SFCFLX |
VZH(I) = VZ0(I) * VAHS(I) |
|
./pkg/fizhi/fizhi_turb.F |
2589 |
SFCFLX |
DO 9028 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2590 |
SFCFLX |
VCU(I) = VK / VPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
2591 |
SFCFLX |
VUSTAR(I) = VCU(I) * VWS(I) |
|
./pkg/fizhi/fizhi_turb.F |
2593 |
SFCFLX |
VPSIG(I) = VH0(I) * VUSTAR(I) - USTH0S |
|
./pkg/fizhi/fizhi_turb.F |
2594 |
SFCFLX |
if(VPSIG(I).lt.0.) VPSIG(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
2595 |
SFCFLX |
VPSIG(I) = SQRT(VPSIG(I)) |
|
./pkg/fizhi/fizhi_turb.F |
2596 |
SFCFLX |
VPSIG(I) = BMDL * VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2597 |
SFCFLX |
VPSIHG(I) = VPSIH(I) + VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2601 |
SFCFLX |
VAPSIM(I) = VCU(I) * RVK |
|
./pkg/fizhi/fizhi_turb.F |
2602 |
SFCFLX |
VAPSIHG(I) = 1. / VPSIHG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2603 |
SFCFLX |
VRIB1(I) = VAPSIM(I) * VAPSIM(I) * VPSIHG(I) * VZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
2624 |
SFCFLX |
DO 9030 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2625 |
SFCFLX |
VZETAL(I) = VZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
2626 |
SFCFLX |
VZ0L(I) = VZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
2630 |
SFCFLX |
DO 9032 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2631 |
SFCFLX |
VZETA(I) = VZETA(I) * ( 1. + VDZETA(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2632 |
SFCFLX |
IF (IVBITRIB(I).EQ.1 ) VZETA(I) = |
|
./pkg/fizhi/fizhi_turb.F |
2633 |
SFCFLX |
1 VPSIM(I) * VPSIM(I) * VRIB(I) * VAPSIHG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2637 |
SFCFLX |
DO 9034 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2638 |
SFCFLX |
IF (IVWATER(I).EQ.1 ) then |
|
./pkg/fizhi/fizhi_turb.F |
2639 |
SFCFLX |
VZ0(I) = VZ0(I) * ( 1. + VDZ0(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2640 |
SFCFLX |
IF (VZ0(I) .LE. Z0MIN ) VZ0(I) = Z0MIN |
|
./pkg/fizhi/fizhi_turb.F |
2647 |
SFCFLX |
DO 9036 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2648 |
SFCFLX |
VPSIM(I) = VPSIM(I) + VDPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
2649 |
SFCFLX |
VCU(I) = VK / VPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
2650 |
SFCFLX |
VUSTAR(I) = VCU(I) * VWS(I) |
|
./pkg/fizhi/fizhi_turb.F |
2652 |
SFCFLX |
VPSIG(I) = VH0(I) * VUSTAR(I) - USTH0S |
|
./pkg/fizhi/fizhi_turb.F |
2653 |
SFCFLX |
if(VPSIG(I).lt.0.) VPSIG(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
2654 |
SFCFLX |
VPSIG(I) = SQRT(VPSIG(I)) |
|
./pkg/fizhi/fizhi_turb.F |
2655 |
SFCFLX |
VPSIG(I) = BMDL * VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2656 |
SFCFLX |
VPSIHG(I) = VPSIH(I) + VDPSIH(I) + VPSIG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2657 |
SFCFLX |
VCT(I) = VK / VPSIHG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2664 |
SFCFLX |
DO 9038 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2665 |
SFCFLX |
VRIB1(I) = VRIB(I) |
|
./pkg/fizhi/fizhi_turb.F |
2666 |
SFCFLX |
VWS1(I) = VWS(I) |
|
./pkg/fizhi/fizhi_turb.F |
2675 |
SFCFLX |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2676 |
SFCFLX |
VTEMP(I) = 10. * VAHS(I) * VZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
2677 |
SFCFLX |
VZH(I) = VZ0(I) * 0.1 |
|
./pkg/fizhi/fizhi_turb.F |
2682 |
SFCFLX |
DO I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2683 |
SFCFLX |
VTEMP(I) = ( VPSIH2(I) + VPSIG(I) ) / VPSIHG(I) |
|
./pkg/fizhi/fizhi_turb.F |
2684 |
SFCFLX |
VRHO(I) = VPKE(I)*( VTH2(I) + VTEMP(I) * (VTH1(I)-VTH2(I)) ) |
|
./pkg/fizhi/fizhi_turb.F |
2685 |
SFCFLX |
VRHO(I) = VPE(I)*100. / ( RGAS * VRHO(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2730 |
SFCFLX |
DO 9044 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
2731 |
SFCFLX |
VRHOZPK(I) = VRHO(I) * VPKE(I) |
|
./pkg/fizhi/fizhi_turb.F |
2732 |
SFCFLX |
VKH(I) = VUSTAR(I) * VCT(I) |
|
./pkg/fizhi/fizhi_turb.F |
2733 |
SFCFLX |
VKM(I) = VUSTAR(I) * VCU(I) |
|
./pkg/fizhi/fizhi_turb.F |
2952 |
PHI |
DO 9000 I = 1,N |
|
./pkg/fizhi/fizhi_turb.F |
2953 |
PHI |
IF(Z(I).GE.0.15)IBIT1 = IBIT1 + 1 |
|
./pkg/fizhi/fizhi_turb.F |
2954 |
PHI |
IF(Z(I).GT.2. )IBIT2 = IBIT2 + 1 |
|
./pkg/fizhi/fizhi_turb.F |
2959 |
PHI |
DO 9002 I = 1,N |
|
./pkg/fizhi/fizhi_turb.F |
2960 |
PHI |
ZSTAR(I) = 100. * Z(I) - 14. |
|
./pkg/fizhi/fizhi_turb.F |
2964 |
PHI |
DO 9004 I = 1,N |
|
./pkg/fizhi/fizhi_turb.F |
2965 |
PHI |
TEMP1(I) = Z(I)*0.5 |
|
./pkg/fizhi/fizhi_turb.F |
2966 |
PHI |
IF( Z(I) .LE. 2. )TEMP1(I) = 1. |
|
./pkg/fizhi/fizhi_turb.F |
2967 |
PHI |
TEMP1(I) = LOG10(TEMP1(I)) |
|
./pkg/fizhi/fizhi_turb.F |
2968 |
PHI |
TEMP1(I) = (TEMP1(I) + 9.3) * 20. |
|
./pkg/fizhi/fizhi_turb.F |
2969 |
PHI |
IF( Z(I) .GT. 2. ) ZSTAR(I) = TEMP1(I) |
|
./pkg/fizhi/fizhi_turb.F |
2970 |
PHI |
IF( Z(I).GT.1.78e10 ) ZSTAR(I) = 384.9999 |
|
./pkg/fizhi/fizhi_turb.F |
2975 |
PHI |
DO 9006 I = 1,N |
|
./pkg/fizhi/fizhi_turb.F |
2976 |
PHI |
I1(I) = ZSTAR(I) |
|
./pkg/fizhi/fizhi_turb.F |
2977 |
PHI |
I2(I) = I1(I) + 1 |
|
./pkg/fizhi/fizhi_turb.F |
2978 |
PHI |
TEMP1(I) = ZSTAR(I) - I1(I) |
|
./pkg/fizhi/fizhi_turb.F |
2983 |
PHI |
DO 9008 I = 1,N |
|
./pkg/fizhi/fizhi_turb.F |
2985 |
PHI |
E1(I) = PHIM0( I1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2986 |
PHI |
E2(I) = PHIM0( I2(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2987 |
PHI |
PHIM(I) = TEMP1(I) * ( E2(I)-E1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2988 |
PHI |
PHIM(I) = PHIM(I) + E1(I) |
|
./pkg/fizhi/fizhi_turb.F |
2995 |
PHI |
DO 9010 I = 1,N |
|
./pkg/fizhi/fizhi_turb.F |
2997 |
PHI |
E1(I) = PHIH0( I1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2998 |
PHI |
E2(I) = PHIH0( I2(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
2999 |
PHI |
PHIH(I) = TEMP1(I) * ( E2(I)-E1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3000 |
PHI |
PHIH(I) = PHIH(I) + E1(I) |
|
./pkg/fizhi/fizhi_turb.F |
3007 |
PHI |
DO 9012 I = 1,N |
|
./pkg/fizhi/fizhi_turb.F |
3008 |
PHI |
ZSTAR(I) = -Z(I) |
|
./pkg/fizhi/fizhi_turb.F |
3012 |
PHI |
DO 9014 I = 1,N |
|
./pkg/fizhi/fizhi_turb.F |
3013 |
PHI |
IF( Z(I) .LT. 0.15 ) PHIM(I) = 1. + ZSTAR(I) |
|
./pkg/fizhi/fizhi_turb.F |
3014 |
PHI |
2 *(0.25+ZSTAR(I)*(0.09375+ZSTAR(I)* |
|
./pkg/fizhi/fizhi_turb.F |
3015 |
PHI |
3 (0.03125+0.00732422 * ZSTAR(I)))) |
|
./pkg/fizhi/fizhi_turb.F |
3020 |
PHI |
DO 9016 I = 1,N |
|
./pkg/fizhi/fizhi_turb.F |
3021 |
PHI |
IF( Z(I) .LT. 0.15 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3022 |
PHI |
PHIH(I) =1.+ Z(I) * (0.5+ZSTAR(I)*(0.375+ZSTAR(I)* |
|
./pkg/fizhi/fizhi_turb.F |
3023 |
PHI |
1 (0.5+ZSTAR(I)*(0.8203125+ZSTAR(I)* |
|
./pkg/fizhi/fizhi_turb.F |
3024 |
PHI |
2 (1.5+2.93262*ZSTAR(I)))))) |
|
./pkg/fizhi/fizhi_turb.F |
3025 |
PHI |
PHIH(I) = 1. / PHIH(I) |
|
./pkg/fizhi/fizhi_turb.F |
3091 |
PSI |
DO 9000 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3092 |
PSI |
VPSIM(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3093 |
PSI |
VPSIH(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3094 |
PSI |
VX(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3095 |
PSI |
VXS(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3096 |
PSI |
VY(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3097 |
PSI |
VYS(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3098 |
PSI |
ZZ0(I) = VZH(I)*VZZ(I) |
|
./pkg/fizhi/fizhi_turb.F |
3101 |
PSI |
DO 9122 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3102 |
PSI |
IF(VZZ(I).LE.-1.e-7)IBIT = IBIT + 1 |
|
./pkg/fizhi/fizhi_turb.F |
3104 |
PSI |
DO 9022 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3105 |
PSI |
IF(VZZ(I).LE.-1.e-7)THEN |
|
./pkg/fizhi/fizhi_turb.F |
3106 |
PSI |
INTSTB(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
3108 |
PSI |
INTSTB(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
3119 |
PSI |
DO 9002 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3120 |
PSI |
IF (INTSTB(I).EQ.1)THEN |
|
./pkg/fizhi/fizhi_turb.F |
3122 |
PSI |
Z(indx) = VZZ(I) |
|
./pkg/fizhi/fizhi_turb.F |
3123 |
PSI |
Z0(indx) = ZZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3127 |
PSI |
DO 9004 I = 1,IBIT |
|
./pkg/fizhi/fizhi_turb.F |
3128 |
PSI |
Z(I) = -18. * Z(I) |
|
./pkg/fizhi/fizhi_turb.F |
3129 |
PSI |
Z0(I) = -18. * Z0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3141 |
PSI |
DO 9006 I = 1,IBIT |
|
./pkg/fizhi/fizhi_turb.F |
3142 |
PSI |
ARG1(I) = 1. - X1(I) |
|
./pkg/fizhi/fizhi_turb.F |
3143 |
PSI |
IF ( Z(I) .LT. 0.013 ) ARG1(I) = |
|
./pkg/fizhi/fizhi_turb.F |
3144 |
PSI |
1 Z(I) * ( 0.25 - 0.09375 * Z(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3146 |
PSI |
ARG0(I) = 1. - X0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3147 |
PSI |
IF ( Z0(I) .LT. 0.013 ) ARG0(I) = |
|
./pkg/fizhi/fizhi_turb.F |
3148 |
PSI |
1 Z0(I) * ( 0.25 - 0.09375 * Z0(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3150 |
PSI |
ARG1(I) = ARG1(I) * ( 1.+X0(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3151 |
PSI |
ARG0(I) = ARG0(I) * ( 1.+X1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3152 |
PSI |
DX(I) = X1(I) - X0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3153 |
PSI |
ARG1(I) = ARG1(I) / ARG0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3154 |
PSI |
ARG0(I) = -DX(I) / ( 1. + X1(I)*X0(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3155 |
PSI |
ARG0(I) = ATAN( ARG0(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3156 |
PSI |
ARG1(I) = LOG( ARG1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3157 |
PSI |
PSI2(I) = 2. * ARG0(I) + ARG1(I) |
|
./pkg/fizhi/fizhi_turb.F |
3158 |
PSI |
PSI2(I) = PSI2(I) + DX(I) |
|
./pkg/fizhi/fizhi_turb.F |
3162 |
PSI |
DO 9008 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3163 |
PSI |
IF( INTSTB(I).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3165 |
PSI |
VPSIM(I) = PSI2(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3166 |
PSI |
VX(I) = X1(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3167 |
PSI |
VXS(I) = X0(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3178 |
PSI |
DO 9010 I = 1,IBIT |
|
./pkg/fizhi/fizhi_turb.F |
3179 |
PSI |
ARG1(I) = 1. - Y1(I) |
|
./pkg/fizhi/fizhi_turb.F |
3180 |
PSI |
IF( Z(I) .LT. 0.0065 ) ARG1(I) = |
|
./pkg/fizhi/fizhi_turb.F |
3181 |
PSI |
1 Z(I) * ( 0.5 - 0.625 * Z(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3183 |
PSI |
ARG0(I) = 1. - Y0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3184 |
PSI |
IF( Z0(I) .LT. 0.0065 ) ARG0(I) = |
|
./pkg/fizhi/fizhi_turb.F |
3185 |
PSI |
1 Z0(I) * ( 0.5 - 0.625 * Z0(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3187 |
PSI |
ARG1(I) = ARG1(I) * ( 1. + Y0(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3188 |
PSI |
ARG0(I) = ARG0(I) * ( 1. + Y1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3189 |
PSI |
ARG1(I) = ARG1(I) / ARG0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3190 |
PSI |
PSI2(I) = LOG( ARG1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3191 |
PSI |
PSI2(I) = PSI2(I) - Y1(I) + Y0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3195 |
PSI |
DO 9012 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3196 |
PSI |
IF( INTSTB(I).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3198 |
PSI |
VPSIH(I) = PSI2(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3199 |
PSI |
VY(I) = Y1(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3200 |
PSI |
VYS(I) = Y0(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3210 |
PSI |
DO 9114 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3211 |
PSI |
IF(VZZ(I).GT.-1.e-7)THEN |
|
./pkg/fizhi/fizhi_turb.F |
3215 |
PSI |
DO 9014 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3216 |
PSI |
IF(VZZ(I).GT.-1.e-7)THEN |
|
./pkg/fizhi/fizhi_turb.F |
3217 |
PSI |
INTSTB(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
3219 |
PSI |
INTSTB(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
3227 |
PSI |
DO 9016 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3228 |
PSI |
IF (INTSTB(I).EQ.1)THEN |
|
./pkg/fizhi/fizhi_turb.F |
3230 |
PSI |
Z(indx) = VZZ(I) |
|
./pkg/fizhi/fizhi_turb.F |
3231 |
PSI |
Z0(indx) = ZZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3232 |
PSI |
ARG1(indx) = VZH(I) |
|
./pkg/fizhi/fizhi_turb.F |
3239 |
PSI |
DO 9018 I = 1,IBIT |
|
./pkg/fizhi/fizhi_turb.F |
3240 |
PSI |
HZ(I) = 1. / ARG1(I) |
|
./pkg/fizhi/fizhi_turb.F |
3241 |
PSI |
Z1(I) = Z(I) |
|
./pkg/fizhi/fizhi_turb.F |
3242 |
PSI |
Z2(I) = ZWM |
|
./pkg/fizhi/fizhi_turb.F |
3244 |
PSI |
IF ( Z(I) .GT. ZWM ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3245 |
PSI |
Z1(I) = ZWM |
|
./pkg/fizhi/fizhi_turb.F |
3246 |
PSI |
Z2(I) = Z(I) |
|
./pkg/fizhi/fizhi_turb.F |
3249 |
PSI |
IF ( Z0(I) .GT. Z0M ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3250 |
PSI |
Z0(I) = Z0M |
|
./pkg/fizhi/fizhi_turb.F |
3251 |
PSI |
INTZ0(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
3253 |
PSI |
INTZ0(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
3256 |
PSI |
X1NUM(I) = 1. + 5. * Z1(I) |
|
./pkg/fizhi/fizhi_turb.F |
3257 |
PSI |
X0NUM(I) = 1. + 5. * Z0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3258 |
PSI |
X1DEN(I) = 1. / (1. + CM1 * (X1NUM(I) * Z1(I)) ) |
|
./pkg/fizhi/fizhi_turb.F |
3259 |
PSI |
X0DEN(I) = 1. + CM1 * (X0NUM(I) * Z0(I)) |
|
./pkg/fizhi/fizhi_turb.F |
3261 |
PSI |
IF ( (INTZ0(I).EQ.1) .OR. (Z(I).GT.ZWM) ) |
|
./pkg/fizhi/fizhi_turb.F |
3262 |
PSI |
1 HZ(I) = Z1(I) / Z0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3263 |
PSI |
ARG1(I) = HZ(I)*HZ(I)*X0DEN(I)*X1DEN(I) |
|
./pkg/fizhi/fizhi_turb.F |
3264 |
PSI |
ARG1(I) = LOG( ARG1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3265 |
PSI |
ARG1(I) = 0.5 * ARG1(I) |
|
./pkg/fizhi/fizhi_turb.F |
3266 |
PSI |
ARG0(I) = (Z1(I) + 0.1) * (Z0(I) + 0.1) |
|
./pkg/fizhi/fizhi_turb.F |
3267 |
PSI |
ARG0(I) = CM3 + ARG0(I) * CM4 |
|
./pkg/fizhi/fizhi_turb.F |
3268 |
PSI |
ARG0(I) = ( Z1(I) - Z0(I) ) / ARG0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3269 |
PSI |
ARG0(I) = ATAN( ARG0(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3270 |
PSI |
TEMP(I) = ARG1(I) + CM5 * ARG0(I) |
|
./pkg/fizhi/fizhi_turb.F |
3272 |
PSI |
X0(I) = X0NUM(I) / X0DEN(I) |
|
./pkg/fizhi/fizhi_turb.F |
3273 |
PSI |
IF ( INTZ0(I).EQ.1 ) X0(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3274 |
PSI |
Z2ZWM(I) = Z2(I) * RZWM |
|
./pkg/fizhi/fizhi_turb.F |
3283 |
PSI |
DO 9020 I = 1,IBIT |
|
./pkg/fizhi/fizhi_turb.F |
3284 |
PSI |
X1(I) = X1NUM(I) * X1DEN(I) |
|
./pkg/fizhi/fizhi_turb.F |
3285 |
PSI |
ARG1(I) = LOG( Z2ZWM(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3286 |
PSI |
PSI2(I) = TEMP(I) + CM6 * ARG1(I) |
|
./pkg/fizhi/fizhi_turb.F |
3290 |
PSI |
DO 9030 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3291 |
PSI |
IF( INTSTB(I).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3293 |
PSI |
VPSIM(I) = PSI2(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3294 |
PSI |
VX(I) = X1(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3295 |
PSI |
VXS(I) = X0(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3306 |
PSI |
DO 9024 I = 1,IBIT |
|
./pkg/fizhi/fizhi_turb.F |
3307 |
PSI |
Y1DEN(I) = 1. + CM1 * ( X1NUM(I) * Z(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3308 |
PSI |
Y1(I) = X1NUM(I) / Y1DEN(I) |
|
./pkg/fizhi/fizhi_turb.F |
3309 |
PSI |
ARG1(I) = CM7 * Z2ZWM(I) / ( CM2 + Z2(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3310 |
PSI |
ARG0(I) = 6. |
|
./pkg/fizhi/fizhi_turb.F |
3311 |
PSI |
IF ( Z2(I) .GT. ZCM ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3312 |
PSI |
Y1(I) = YCM |
|
./pkg/fizhi/fizhi_turb.F |
3313 |
PSI |
ARG1(I) = Z2(I) * RZCM |
|
./pkg/fizhi/fizhi_turb.F |
3314 |
PSI |
ARG0(I) = YCM |
|
./pkg/fizhi/fizhi_turb.F |
3315 |
PSI |
TEMP(I) = TEMP(I) + CM8 |
|
./pkg/fizhi/fizhi_turb.F |
3317 |
PSI |
ARG1(I) = LOG( ARG1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
3318 |
PSI |
PSI2(I) = TEMP(I) + ARG0(I) * ARG1(I) |
|
./pkg/fizhi/fizhi_turb.F |
3322 |
PSI |
DO 9026 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
3323 |
PSI |
IF( INTSTB(I).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3325 |
PSI |
VPSIH(I) = PSI2(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3326 |
PSI |
VY(I) = Y1(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3327 |
PSI |
VYS(I) = X0(indx) |
|
./pkg/fizhi/fizhi_turb.F |
3419 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3420 |
TRBLEN |
STBFCN(I,L) = STRT(I,L) - RICR * DW2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3421 |
TRBLEN |
INT1(I,L) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
3422 |
TRBLEN |
IF( STBFCN(I,L).LE.0. ) INT1(I,L) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
3427 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3428 |
TRBLEN |
INT2(I,L) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
3429 |
TRBLEN |
IF( (INT1(I,L).EQ.1) .NEQV. (INT1(I,L+1).EQ.1) ) INT2(I,L) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
3435 |
TRBLEN |
DO 30 I=1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3436 |
TRBLEN |
IBIT = IBIT + INT1(I,LMIN) |
|
./pkg/fizhi/fizhi_turb.F |
3445 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3446 |
TRBLEN |
XL0(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3449 |
TRBLEN |
DO 70 I=1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3450 |
TRBLEN |
XL0(I,NLEV) = DZ3(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
3465 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3466 |
TRBLEN |
INT2(I,L) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
3467 |
TRBLEN |
IF ( INT1(I,L).EQ.1 .OR. INT1(I,L+1).EQ.1 ) INT2(I,L) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
3468 |
TRBLEN |
IF ( INT2(I,L).EQ.1 ) |
|
./pkg/fizhi/fizhi_turb.F |
3469 |
TRBLEN |
& XL0(I,L) = (0.5+DZITRP(I,L)) * DZ3(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
3472 |
TRBLEN |
DO 90 I=1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3473 |
TRBLEN |
INT2(I,NLEVM1) = INT1(I,NLEV) |
|
./pkg/fizhi/fizhi_turb.F |
3477 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3478 |
TRBLEN |
IF ( INT2(I,L).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3479 |
TRBLEN |
XL0(I,L+1) = XL0(I,L+1) + ( (0.5-DZITRP(I,L)) * DZ3(I,L+1) ) |
|
./pkg/fizhi/fizhi_turb.F |
3484 |
TRBLEN |
DO 110 I=1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3485 |
TRBLEN |
IF( INT1(I,1).EQ.1 ) XL0(I,1) = XL0(I,1) + DZ3(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3493 |
TRBLEN |
DO 120 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3494 |
TRBLEN |
IF( INT1(I,LM1).EQ.1 ) XL0(I,L) = XL0(I,L) + XL0(I,LM1) |
|
./pkg/fizhi/fizhi_turb.F |
3499 |
TRBLEN |
DO 130 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3500 |
TRBLEN |
IF( INT1(I,NLEVM1).EQ.1 .AND. INT1(I,NLEV).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3501 |
TRBLEN |
XL0(I,NLEV) = XL0(I,NLEV) + XL0(I,NLEVM1) |
|
./pkg/fizhi/fizhi_turb.F |
3509 |
TRBLEN |
DO 140 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3510 |
TRBLEN |
IF( INT1(I,LP1).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3511 |
TRBLEN |
IF( INT1(I,L) .EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3512 |
TRBLEN |
XL0(I,L) = XL0(I,LP1) |
|
./pkg/fizhi/fizhi_turb.F |
3514 |
TRBLEN |
XL0(I,L) = XL0(I,L) + XL0(I,LP1) |
|
./pkg/fizhi/fizhi_turb.F |
3522 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3523 |
TRBLEN |
IF ( XL0(I,L).LT.XL0CNV ) XL0(I,L) = XL0CNV |
|
./pkg/fizhi/fizhi_turb.F |
3536 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3537 |
TRBLEN |
INT2(I,L) = 1 - INT1(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3544 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3545 |
TRBLEN |
IF( INT1(I,L).EQ.0 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3546 |
TRBLEN |
XL0(I,L) = Q(I,L) / XL0(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3547 |
TRBLEN |
XL0(I,L) = XL0(I,L) * XL0(I,L) + 1.0E-20 |
|
./pkg/fizhi/fizhi_turb.F |
3548 |
TRBLEN |
XL0(I,L) = STBFCN(I,L) + XL0(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3549 |
TRBLEN |
XL0(I,L) = SQRT( XL0(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
3550 |
TRBLEN |
XL0(I,L) = Q(I,L) / XL0(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3552 |
TRBLEN |
INT2(I,L) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
3553 |
TRBLEN |
IF( XL0(I,L).LT.XL0MIN ) INT2(I,L) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
3562 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3563 |
TRBLEN |
INT2(I,L) = 1 - INT1(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3564 |
TRBLEN |
IF ( INT2(I,L).EQ.1 ) XL0(I,L) = XL0MIN |
|
./pkg/fizhi/fizhi_turb.F |
3569 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3570 |
TRBLEN |
IF ( INT2(I,L).EQ.1 ) XL0(I,L) = XL0MIN |
|
./pkg/fizhi/fizhi_turb.F |
3581 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3582 |
TRBLEN |
XL(I,L) = XL0(I,L) * VKZE(I,L) / ( XL0(I,L)+VKZE(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
3593 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3594 |
TRBLEN |
Q1(I,L) = Q(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3595 |
TRBLEN |
INT1(I,L) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
3596 |
TRBLEN |
IF ( Q(I,L).LE.Q(I,L+1) ) INT1(I,L) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
3597 |
TRBLEN |
IF ( INT1(I,L).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3598 |
TRBLEN |
XL0(I,L) = XL0(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
3599 |
TRBLEN |
Q1(I,L) = Q(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
3605 |
TRBLEN |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3606 |
TRBLEN |
QXLM(I,L) = XL0(I,L)*VKZM(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3607 |
TRBLEN |
& / ( XL0(I,L)+VKZM(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
3608 |
TRBLEN |
QXLM(I,L) = CLMT53 * Q1(I,L)*QXLM(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3662 |
TRBITP |
integer I,L |
|
./pkg/fizhi/fizhi_turb.F |
3671 |
TRBITP |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3672 |
TRBITP |
DZITRP(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3676 |
TRBITP |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3677 |
TRBITP |
IF ( INTCHG(I,L).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3678 |
TRBITP |
DDD(I,L) = ( CP *(DTHV(I,L+1)*DPK(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3679 |
TRBITP |
& + DTHV(I,L)*DPK(I,L+1)) ) |
|
./pkg/fizhi/fizhi_turb.F |
3680 |
TRBITP |
& - ( (2.*RICR) * ( DU(I,L+1)* DU(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3681 |
TRBITP |
& + DV(I,L+1)* DV(I,L)) ) |
|
./pkg/fizhi/fizhi_turb.F |
3682 |
TRBITP |
AAA(I,L) = STBFCN(I,L) + STBFCN(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
3683 |
TRBITP |
BBB(I,L) = STBFCN(I,L) - STBFCN(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
3684 |
TRBITP |
CCC(I,L) = 1. / BBB(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3685 |
TRBITP |
DZITRP(I,L) = AAA(I,L) * CCC(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3686 |
TRBITP |
AAA(I,L) = AAA(I,L) - DDD(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3687 |
TRBITP |
DDD(I,L) = ( DDD(I,L) * DDD(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
3688 |
TRBITP |
& - 4. * (STBFCN(I,L+1) * STBFCN(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
3689 |
TRBITP |
DDD(I,L) = DDD(I,L)*CCC(I,L)*CCC(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3690 |
TRBITP |
DDD(I,L) = SQRT( DDD(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
3693 |
TRBITP |
IF( INTCHG(I,L).EQ.1 .AND. AAA(I,L).NE.0. ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3694 |
TRBITP |
DZITRP(I,L) = ( BBB(I,L)*(1.-DDD(I,L)) ) / AAA(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3697 |
TRBITP |
DZITRP(I,L) = 0.5 * DZITRP(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3767 |
TRBL20 |
DO 10 I=1,ISTNLV |
|
./pkg/fizhi/fizhi_turb.F |
3768 |
TRBL20 |
EE(I,1) = D1HALF * RI(I,1) + D2HALF |
|
./pkg/fizhi/fizhi_turb.F |
3769 |
TRBL20 |
RF(I,1) = EE(I,1)* EE(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3770 |
TRBL20 |
RF(I,1) = RF(I,1)- D3*RI(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3771 |
TRBL20 |
RF(I,1) = SQRT( RF(I,1) ) |
|
./pkg/fizhi/fizhi_turb.F |
3772 |
TRBL20 |
RF(I,1) = EE(I,1) - RF(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3774 |
TRBL20 |
IF( RI(I,1).LE.1.e-4 .AND. RI(I,1).GE.-1.e-4 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3775 |
TRBL20 |
RF(I,1) = D3B2*RI(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3783 |
TRBL20 |
IF( RI(I,1).LT.RICR .AND. RF(I,1).LT.RFC ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3784 |
TRBL20 |
ZKH(I,1) = ( RFC-RF(I,1) ) / (1.-RF(I,1)) |
|
./pkg/fizhi/fizhi_turb.F |
3785 |
TRBL20 |
ZKM(I,1) = CM * (RF1-RF(I,1)) |
|
./pkg/fizhi/fizhi_turb.F |
3786 |
TRBL20 |
ZKM(I,1) = ZKH(I,1)*ZKM(I,1) / (RF2-RF(I,1)) |
|
./pkg/fizhi/fizhi_turb.F |
3787 |
TRBL20 |
ZKH(I,1) = CH *ZKH(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3788 |
TRBL20 |
QE(I,1) = ZKM(I,1)*DW2(I,1) - ZKH(I,1)*STRT(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3791 |
TRBL20 |
IF( QE(I,1).LT.1.e-14 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3792 |
TRBL20 |
INTSTB(I,1) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
3793 |
TRBL20 |
QE(I,1) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3795 |
TRBL20 |
INTSTB(I,1) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
3796 |
TRBL20 |
QE(I,1) = B1*QE(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3797 |
TRBL20 |
QE(I,1) = SQRT( QE(I,1) ) |
|
./pkg/fizhi/fizhi_turb.F |
3798 |
TRBL20 |
QE(I,1) = XL(I,1)*QE(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3800 |
TRBL20 |
QQE(I,1) = 0.5 * QE(I,1) * QE(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3866 |
TRBL25 |
DO 10 I=1,ISTNLV |
|
./pkg/fizhi/fizhi_turb.F |
3867 |
TRBL25 |
IF( INTQ(I,1).EQ.1 .AND. INTSTB(I,1).EQ.0 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
3868 |
TRBL25 |
XQ(I,1) = XL(I,1) / Q(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3869 |
TRBL25 |
XQ(I,1) = XQ(I,1) * XQ(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3870 |
TRBL25 |
STRT(I,1) = XQ(I,1) * STRT(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3871 |
TRBL25 |
DW2(I,1) = XQ(I,1) * DW2(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3872 |
TRBL25 |
F2(I,1) = 1.+FF2 * STRT(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3873 |
TRBL25 |
F3(I,1) = 1.+FF3 * STRT(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3874 |
TRBL25 |
F4(I,1) = 1.+FF4 * STRT(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3875 |
TRBL25 |
ZKH(I,1) = (F4(I,1) * F2(I,1)) |
|
./pkg/fizhi/fizhi_turb.F |
3876 |
TRBL25 |
& + A4 * (F3(I,1) * DW2(I,1)) |
|
./pkg/fizhi/fizhi_turb.F |
3877 |
TRBL25 |
ZKH(I,1) = (F2(I,1) + A3*DW2(I,1)) |
|
./pkg/fizhi/fizhi_turb.F |
3878 |
TRBL25 |
& / ZKH(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3879 |
TRBL25 |
ZKM(I,1) = A1 * (F3(I,1)*ZKH(I,1)+A5) |
|
./pkg/fizhi/fizhi_turb.F |
3880 |
TRBL25 |
& / F2(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3881 |
TRBL25 |
ZKH(I,1) = A2 * ZKH(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3882 |
TRBL25 |
P3(I,1) = ZKH(I,1)*STRT(I,1) + B3 |
|
./pkg/fizhi/fizhi_turb.F |
3883 |
TRBL25 |
P3(I,1) = 2. * ( ZKM(I,1)*DW2(I,1) - P3(I,1) ) |
|
./pkg/fizhi/fizhi_turb.F |
3884 |
TRBL25 |
P3(I,1) = P3(I,1)*Q(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
3938 |
TRBDIF |
integer I,L |
|
./pkg/fizhi/fizhi_turb.F |
3946 |
TRBDIF |
DO 10 I=1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3947 |
TRBDIF |
CC(I,1) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
3950 |
TRBDIF |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3951 |
TRBDIF |
CC(I,L+1) = RHOKDZ(I,L) * FLXFAC(I,L+1) |
|
./pkg/fizhi/fizhi_turb.F |
3955 |
TRBDIF |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3956 |
TRBDIF |
BB(I,L) = RHOKDZ(I,L) * FLXFAC(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
3957 |
TRBDIF |
AA(I,L) = 1. + CC(I,L) + BB(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
395 |
TURBIO |
do I =1,nchptot |
|
./pkg/fizhi/fizhi_turb.F |
3964 |
TRBDIF |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3965 |
TRBDIF |
AA(I,L) = AA(I,L) - XX2(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
396 |
TURBIO |
utility(I) = pmsc(i) + ptop |
|
./pkg/fizhi/fizhi_turb.F |
3976 |
TRBDIF |
DO 50 I=1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3977 |
TRBDIF |
DXX1G(I) = CC(I,NLEVP1) * ( XX1(I,NLEV)-XX1(I,NLEVP1) ) |
|
./pkg/fizhi/fizhi_turb.F |
3982 |
TRBDIF |
DO 60 I=1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3983 |
TRBDIF |
DXX2G(I) = CC(I,NLEVP1) * ( XX2(I,NLEV)-XX2(I,NLEVP1) ) |
|
./pkg/fizhi/fizhi_turb.F |
3994 |
TRBDIF |
DO I =1,irun |
|
./pkg/fizhi/fizhi_turb.F |
3995 |
TRBDIF |
IF ( XX1(I,L).LT.EPSL ) XX1(I,L) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
4014 |
VTRI0 |
DO 9000 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
4015 |
VTRI0 |
A(I,1) = 1. / A(I,1) |
|
./pkg/fizhi/fizhi_turb.F |
4020 |
VTRI0 |
DO 9002 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
4021 |
VTRI0 |
C(I,L) = C(I,L) * A(I,LM1) |
|
./pkg/fizhi/fizhi_turb.F |
4022 |
VTRI0 |
A(I,L) = 1. / ( A(I,L) - B(I,LM1) * C(I,L) ) |
|
./pkg/fizhi/fizhi_turb.F |
4023 |
VTRI0 |
F(I,L) = F(I,L) + F(I,LM1) * C(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
4028 |
VTRI0 |
DO 9004 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
4029 |
VTRI0 |
Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) * A(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
4047 |
VTRI1 |
DO 9000 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
4048 |
VTRI1 |
Y(I,L) = B(I,L) * Y(I,L+1) * A(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
4067 |
VTRI2 |
DO 9000 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
4068 |
VTRI2 |
F(I,L) = F(I,L) + F(I,L-1) * C(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
4073 |
VTRI2 |
DO 9002 I = 1,irun |
|
./pkg/fizhi/fizhi_turb.F |
4074 |
VTRI2 |
Y(I,L) = (F(I,L) + B(I,L) * Y(I,L+1)) * A(I,L) |
|
./pkg/fizhi/fizhi_turb.F |
4169 |
LINADJ |
DO 9000 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4170 |
LINADJ |
IF (IWATER(I).EQ.1) VX0PSIM(I) = VAPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4174 |
LINADJ |
DO 9002 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4175 |
LINADJ |
VX0PSIM(I) = VX0(I) * VAPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4180 |
LINADJ |
DO 9004 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4181 |
LINADJ |
VDZ0(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
4182 |
LINADJ |
VG(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
4183 |
LINADJ |
VG0(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
4184 |
LINADJ |
VR1MG0(I) = 1. |
|
./pkg/fizhi/fizhi_turb.F |
4190 |
LINADJ |
DO 9006 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4191 |
LINADJ |
IF ( IWATER(I).EQ.1) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4192 |
LINADJ |
VAZ0(I) = 1. / VZ1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4193 |
LINADJ |
VG(I) = VDZSEA(I) * VAZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4194 |
LINADJ |
VG0(I) = VX0PSIM(I) * VG(I) |
|
./pkg/fizhi/fizhi_turb.F |
4195 |
LINADJ |
VR1MG0(I) = 1. / ( 1. - VG0(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
4196 |
LINADJ |
VDZ0(I) = ( VZ2(I) - VZ1(I) ) * VR1MG0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4203 |
LINADJ |
DO 9008 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4204 |
LINADJ |
IF (IWATER(I).EQ.1) VDZ0(I) = VDZ0(I) * VAZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4211 |
LINADJ |
DO 9010 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4212 |
LINADJ |
VXNUM1(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
4213 |
LINADJ |
IF (VRIB1(I).EQ.0.) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4214 |
LINADJ |
INTRIB(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
4216 |
LINADJ |
INTRIB(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
4218 |
LINADJ |
IF ( INTRIB(I).EQ.0 ) VXNUM1(I) = 1. / VRIB1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4219 |
LINADJ |
VPSIGB2(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
4220 |
LINADJ |
if(vpsig(i).gt.0.)VPSIGB2(I) = |
|
./pkg/fizhi/fizhi_turb.F |
4222 |
LINADJ |
VDX(I) = VX(I) - VX0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4223 |
LINADJ |
VDXPSIM(I) = VDX(I) * VAPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4224 |
LINADJ |
VDY(I) = VY(I) - VY0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4225 |
LINADJ |
VXNUM3(I) = - VPSIGB2(I) |
|
./pkg/fizhi/fizhi_turb.F |
4229 |
LINADJ |
IF (IWATER(I).EQ.1) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4230 |
LINADJ |
VDXPSIM(I) = VDXPSIM(I) * VR1MG0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4231 |
LINADJ |
VXNUM3(I) = VXNUM3(I) + VG(I) * ( VY0(I) - VPSIGB2(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
4232 |
LINADJ |
VXNUM2(I) = VY0(I) - VPSIGB2(I) - VX0PSIM(I) * VPSIGB2(I) |
|
./pkg/fizhi/fizhi_turb.F |
4233 |
LINADJ |
VXNUM2(I) = (VXNUM2(I) * VAPSIHG(I)) - 2. * VX0PSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4234 |
LINADJ |
VXNUM2(I) = VXNUM2(I) * VDZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4238 |
LINADJ |
VDEN(I) = VDY(I) + VDXPSIM(I) * VXNUM3(I) |
|
./pkg/fizhi/fizhi_turb.F |
4239 |
LINADJ |
VDEN(I) = ( 1. + VDEN(I) * VAPSIHG(I) ) - 2. * VDXPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4244 |
LINADJ |
DO 9012 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4245 |
LINADJ |
VAWS1(I) = VR1MG0(I) / VWS1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4246 |
LINADJ |
VXNUM3(I) = VXNUM3(I) * VAPSIHG(I) |
|
./pkg/fizhi/fizhi_turb.F |
4250 |
LINADJ |
IF(IWATER(I).EQ.1) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4251 |
LINADJ |
VXNUM3(I) = VXNUM3(I) - 2. * VG0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4252 |
LINADJ |
VXNUM3(I) = VAWS1(I) * VXNUM3(I) |
|
./pkg/fizhi/fizhi_turb.F |
4261 |
LINADJ |
DO 9014 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4262 |
LINADJ |
VXNUM(I) = VRIB2(I) - VRIB1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4263 |
LINADJ |
IF( (VX0(I).GT.XX0MAX).AND.(VXNUM(I).GE.0.) )VXNUM(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
4264 |
LINADJ |
VXNUM(I) = VXNUM1(I) * VXNUM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4269 |
LINADJ |
DO 9016 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4270 |
LINADJ |
VDZETA1(I) = VDZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
4271 |
LINADJ |
VXNUM(I) = VXNUM(I) + VXNUM3(I) * ( VWS2(I) - VWS1(I) ) |
|
./pkg/fizhi/fizhi_turb.F |
4276 |
LINADJ |
DO 9018 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4277 |
LINADJ |
VDZETA1(I) = VXNUM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4278 |
LINADJ |
IF(LWATER.AND.(IWATER(I).EQ.1)) VXNUM(I) = VXNUM(I) + VXNUM2(I) |
|
./pkg/fizhi/fizhi_turb.F |
4279 |
LINADJ |
IF ( VDEN(I) .LT.0.1 ) VDEN(I) = 0.1 |
|
./pkg/fizhi/fizhi_turb.F |
4284 |
LINADJ |
DO 9020 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4285 |
LINADJ |
VDZETA(I) = VXNUM(I) / VDEN(I) |
|
./pkg/fizhi/fizhi_turb.F |
4289 |
LINADJ |
DO 9022 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4290 |
LINADJ |
IF ( (VRIB2(I).EQ.0.) .OR. (VDZETA(I).LE.-1.) ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4291 |
LINADJ |
VINT1(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
4293 |
LINADJ |
VINT1(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
4295 |
LINADJ |
IF ( VINT1(I).EQ.1 ) VDZETA(I) = VDZETA1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4299 |
LINADJ |
DO 9024 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4300 |
LINADJ |
VDZETA2(I) = VDZETA(I) + VDZETA1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4301 |
LINADJ |
IF ( (VRIB2(I).EQ.0.) .OR. (VDZETA2(I).LE.-1.) ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4302 |
LINADJ |
VINT1(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
4304 |
LINADJ |
VINT1(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
4306 |
LINADJ |
IF(VINT1(I).EQ.1)VDZETA(I)=VXNUM1(I)*VRIB2(I) - 1. - VDZETA1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4313 |
LINADJ |
DO 9026 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4314 |
LINADJ |
IF( IWATER(I).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4315 |
LINADJ |
VZCOEF2(I) = VG(I) * VDXPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4316 |
LINADJ |
VDZ0(I) = VDZ0(I) - VZCOEF2(I) * VDZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
4322 |
LINADJ |
DO 9028 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4323 |
LINADJ |
IF(IWATER(I).EQ.1) VZCOEF1(I) = VG(I) * VAWS1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4328 |
LINADJ |
DO 9030 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4329 |
LINADJ |
IF (IWATER(I).EQ.1) VDZ0(I) = |
|
./pkg/fizhi/fizhi_turb.F |
4330 |
LINADJ |
1 VZCOEF1(I) * ( VWS2(I) - VWS1(I) ) - VZCOEF2(I) * VDZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
4337 |
LINADJ |
DO 9032 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4338 |
LINADJ |
IF (IWATER(I).EQ.1) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4339 |
LINADJ |
VDPSIM(I) = - VDZ0(I) * VAZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4340 |
LINADJ |
VDPSIH(I) = VDPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4346 |
LINADJ |
DO 9034 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4347 |
LINADJ |
VDPSIM(I) = VDX(I) * VDZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
4348 |
LINADJ |
VDPSIH(I) = VDY(I) * VDZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
4350 |
LINADJ |
IF (IWATER(I).EQ.1 ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4351 |
LINADJ |
VDPSIM(I) = VDPSIM(I) - VX0(I) * VDZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4352 |
LINADJ |
VDPSIH(I) = VDPSIH(I) - VY0(I) * VDZ0(I) |
|
./pkg/fizhi/fizhi_turb.F |
4361 |
LINADJ |
DO 9036 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4362 |
LINADJ |
VDPSIMC(I) = -0.9 - VDPSIM(I) * VAPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4363 |
LINADJ |
VDPSIHC(I) = -0.9 * VPSIH(I) - VDPSIH(I) |
|
./pkg/fizhi/fizhi_turb.F |
4364 |
LINADJ |
IF ( VDPSIMC(I).GT.0. ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4365 |
LINADJ |
VINT1(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
4367 |
LINADJ |
VINT1(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
4369 |
LINADJ |
IF ( VDPSIHC(I).GT.0. ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4370 |
LINADJ |
VINT2(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
4372 |
LINADJ |
VINT2(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
4374 |
LINADJ |
VDZETA1(I) = 0. |
|
./pkg/fizhi/fizhi_turb.F |
4375 |
LINADJ |
IF(VINT1(I).EQ.1) VDZETA1(I) = VDPSIMC(I) / VDXPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4376 |
LINADJ |
IF((VINT1(I).EQ.1).OR.(VINT2(I).EQ.1)) VTEMPLIN(I) = |
|
./pkg/fizhi/fizhi_turb.F |
4377 |
LINADJ |
1 VDY(I) + VY0(I) * VG(I) * VDXPSIM(I) |
|
./pkg/fizhi/fizhi_turb.F |
4378 |
LINADJ |
IF (VINT2(I).EQ.1) then |
|
./pkg/fizhi/fizhi_turb.F |
4379 |
LINADJ |
VDZETA2(I) = VDPSIHC(I) / VTEMPLIN(I) |
|
./pkg/fizhi/fizhi_turb.F |
4380 |
LINADJ |
IF ( VDZETA2(I).LT.VDZETA1(I) ) VDZETA1(I) = VDZETA2(I) |
|
./pkg/fizhi/fizhi_turb.F |
4382 |
LINADJ |
IF((VINT1(I).EQ.1).OR.(VINT2(I).EQ.1)) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4383 |
LINADJ |
VDZETA(I) = VDZETA1(I) + VDZETA(I) |
|
./pkg/fizhi/fizhi_turb.F |
4384 |
LINADJ |
VDPSIM(I) = VDPSIM(I) + VDX(I) * VR1MG0(I) * VDZETA1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4385 |
LINADJ |
VDPSIH(I) = VDPSIH(I) + VTEMPLIN(I) * VDZETA1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4386 |
LINADJ |
IF ( IWATER(I).EQ.1 ) |
|
./pkg/fizhi/fizhi_turb.F |
4387 |
LINADJ |
1 VDZ0(I) = VDZ0(I) - VG(I) * VDXPSIM(I) * VDZETA1(I) |
|
./pkg/fizhi/fizhi_turb.F |
4449 |
ZCSUB |
DO 9000 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4450 |
ZCSUB |
IF(VUSTAR(I) .LT. 1.e-6)THEN |
|
./pkg/fizhi/fizhi_turb.F |
4451 |
ZCSUB |
INT3(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
4453 |
ZCSUB |
INT3(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
4456 |
ZCSUB |
DO 9002 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4457 |
ZCSUB |
IF(INT3(I).EQ.1) VUSTAR(I) = 1.e-6 |
|
./pkg/fizhi/fizhi_turb.F |
4464 |
ZCSUB |
DO 9004 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4465 |
ZCSUB |
IF( (VUSTAR(I) .GT. USTMX1) .AND. (IWATER(I).EQ.1) ) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4466 |
ZCSUB |
INT4(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
4468 |
ZCSUB |
INT4(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
4471 |
ZCSUB |
DO 9006 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4472 |
ZCSUB |
IF(VUSTAR(I) .GT. USTMX2) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4473 |
ZCSUB |
INT3(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
4475 |
ZCSUB |
INT3(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
4478 |
ZCSUB |
DO 9008 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4479 |
ZCSUB |
IF(VUSTAR(I) .GE. USTMX3) THEN |
|
./pkg/fizhi/fizhi_turb.F |
4480 |
ZCSUB |
INT2(I) = 1 |
|
./pkg/fizhi/fizhi_turb.F |
4482 |
ZCSUB |
INT2(I) = 0 |
|
./pkg/fizhi/fizhi_turb.F |
4487 |
ZCSUB |
DO 9010 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4488 |
ZCSUB |
AA(I,K) = AA1(K) |
|
./pkg/fizhi/fizhi_turb.F |
4489 |
ZCSUB |
IF( INT2(I).EQ.1 ) AA(I,K) = AA2(K) |
|
./pkg/fizhi/fizhi_turb.F |
4490 |
ZCSUB |
IF( INT3(I).EQ.1 ) AA(I,K) = AA3(K) |
|
./pkg/fizhi/fizhi_turb.F |
4491 |
ZCSUB |
IF( INT4(I).EQ.1 ) AA(I,K) = AA4(K) |
|
./pkg/fizhi/fizhi_turb.F |
4499 |
ZCSUB |
DO 9012 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4500 |
ZCSUB |
VDZSEA(I) = ( AA(I,4) + AA(I,5) * VUSTAR(I) ) * VUSTAR(I) |
|
./pkg/fizhi/fizhi_turb.F |
4501 |
ZCSUB |
VZSEA(I) = AA(I,2) + ( AA(I,3) + VDZSEA(I) ) * VUSTAR(I) |
|
./pkg/fizhi/fizhi_turb.F |
4502 |
ZCSUB |
TEMP(I) = AA(I,1) / VUSTAR(I) |
|
./pkg/fizhi/fizhi_turb.F |
4503 |
ZCSUB |
VZSEA(I) = VZSEA(I) + TEMP(I) |
|
./pkg/fizhi/fizhi_turb.F |
4511 |
ZCSUB |
DO 9014 I = 1,IRUN |
|
./pkg/fizhi/fizhi_turb.F |
4512 |
ZCSUB |
VDZSEA(I) = 3. * VDZSEA(I) -(AA(I,4)*VUSTAR(I) - AA(I,3)) |
|
./pkg/fizhi/fizhi_turb.F |
4513 |
ZCSUB |
VDZSEA(I) = VDZSEA(I) * VUSTAR(I) - TEMP(I) |
|
./pkg/fizhi/fizhi_utils.F |
102 |
STRIPINT |
DO 100 I=1,Lena |
|
./pkg/fizhi/fizhi_utils.F |
103 |
STRIPINT |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
107 |
STRIPINT |
DO 300 I=1,Lena |
|
./pkg/fizhi/fizhi_utils.F |
108 |
STRIPINT |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
110 |
STRIPINT |
DO 400 I=1,IB-Lena |
|
./pkg/fizhi/fizhi_utils.F |
111 |
STRIPINT |
B(Lena+I,J) = A(Lena+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
130 |
PASTE |
DO 100 I=1,Lena |
|
./pkg/fizhi/fizhi_utils.F |
131 |
PASTE |
A(I+OFFSET-1,J) = B(I,J) |
|
./pkg/fizhi/fizhi_utils.F |
148 |
PSTBMP |
DO 100 I=1,Lena |
|
./pkg/fizhi/fizhi_utils.F |
149 |
PSTBMP |
A(I+OFFSET-1,J) = A(I+OFFSET-1,J) + B(I,J) |
|
./pkg/fizhi/fizhi_utils.F |
167 |
STRINT |
DO 100 I=1,Lena |
|
./pkg/fizhi/fizhi_utils.F |
168 |
STRINT |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
172 |
STRINT |
DO 300 I=1,Lena |
|
./pkg/fizhi/fizhi_utils.F |
173 |
STRINT |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
175 |
STRINT |
DO 400 I=1,IB-Lena |
|
./pkg/fizhi/fizhi_utils.F |
176 |
STRINT |
B(Lena+I,J) = A(Lena+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
663 |
GRD2MSC |
DO I=1,NCHP |
|
./pkg/fizhi/fizhi_utils.F |
665 |
GRD2MSC |
B(I) = A(IGRD(I)) |
|
./pkg/fizhi/fizhi_utils.F |
690 |
MSC2GRD |
DO I=1,IM*JM |
|
./pkg/fizhi/fizhi_utils.F |
692 |
MSC2GRD |
VT1(I) = ZERO |
|
./pkg/fizhi/fizhi_utils.F |
695 |
MSC2GRD |
DO I=1,NCHP |
|
./pkg/fizhi/fizhi_utils.F |
697 |
MSC2GRD |
VT1(IGRD(I)) = VT1(IGRD(I)) + B(I)*CHFR(I) |
|
./pkg/fizhi/fizhi_utils.F |
700 |
MSC2GRD |
DO I=1,IM*JM |
|
./pkg/fizhi/fizhi_utils.F |
702 |
MSC2GRD |
A(I) = A(I)*(ONE-FRACG(I)) + VT1(I) |
|
./pkg/fizhi/fizhi_utils.F |
73 |
STRIP |
DO 100 I=1,Lena |
|
./pkg/fizhi/fizhi_utils.F |
74 |
STRIP |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
78 |
STRIP |
DO 300 I=1,Lena |
|
./pkg/fizhi/fizhi_utils.F |
79 |
STRIP |
B(I,J) = A(I+OFFSET-1,J) |
|
./pkg/fizhi/fizhi_utils.F |
81 |
STRIP |
DO 400 I=1,IB-Lena |
|
./pkg/fizhi/fizhi_utils.F |
82 |
STRIP |
B(Lena+I,J) = A(Lena+OFFSET-1,J) |
|
./pkg/fizhi/getcon.F |
151 |
NO PROCEDURE |
DO 10 I=1,MAXCON |
|
./pkg/fizhi/getcon.F |
152 |
NO PROCEDURE |
IF(NAME.EQ.ANAME(I)) THEN |
|
./pkg/fizhi/getcon.F |
153 |
NO PROCEDURE |
GETCON = ACON(I) |
|
./pkg/fizhi/getcon.F |
25 |
NO PROCEDURE |
integer I |
|
./pkg/fizhi/slprs.F |
67 |
SLPRS |
TM(I) = TM(I) * (PLE(I,LM+1)/1000.)**AK + BETA*PHIS(I)/(TWO*G) |
|
./pkg/fizhi/slprs.F |
69 |
SLPRS |
TM(I) = THZ(I,LM)*(PLE(I,LM+1)/1000.)**AK + BETA*PHIS(I)/(TWO*G) |
|
./pkg/fizhi/slprs.F |
72 |
SLPRS |
SLP(I) = PHIS(I) / ( R*TM(I) ) |
|
./pkg/fizhi/slprs.F |
73 |
SLPRS |
SLP(I) = PLE(I,LM+1) * EXP( SLP(I) ) |
|
./pkg/fizhi/update_chemistry_exports.F |
322 |
INTERP_OZ |
DO 34 I=1,IRUN |
|
./pkg/fizhi/update_chemistry_exports.F |
333 |
INTERP_OZ |
slope(i) = ( XLAT(I) -ozlats(jlat-1) ) |
|
./pkg/fizhi/update_chemistry_exports.F |
343 |
INTERP_OZ |
QPR1(I) = ozone(jlatm,k) |
|
./pkg/fizhi/update_chemistry_exports.F |
344 |
INTERP_OZ |
QPR2(I) = ozone(jlatp,k) |
|
./pkg/fizhi/update_chemistry_exports.F |
347 |
INTERP_OZ |
DO 38 I=1,IRUN |
|
./pkg/fizhi/update_chemistry_exports.F |
359 |
INTERP_OZ |
DO 44 I = 1,IRUN |
|
./pkg/fizhi/update_chemistry_exports.F |
361 |
INTERP_OZ |
IF( ozlevs(L1).GT.PLEVS(I,L2) ) THEN |
|
./pkg/fizhi/update_chemistry_exports.F |
375 |
INTERP_OZ |
PR1(I) = ozlevs (LM) |
|
./pkg/fizhi/update_chemistry_exports.F |
376 |
INTERP_OZ |
PR2(I) = ozlevs (LP) |
|
./pkg/fizhi/update_chemistry_exports.F |
377 |
INTERP_OZ |
QPR1(I) = O3INT1(I,LM) |
|
./pkg/fizhi/update_chemistry_exports.F |
378 |
INTERP_OZ |
QPR2(I) = O3INT1(I,LP) |
|
./pkg/fizhi/update_chemistry_exports.F |
381 |
INTERP_OZ |
DO 48 I=1,IRUN |
|
./pkg/fizhi/update_chemistry_exports.F |
382 |
INTERP_OZ |
SLOPE(I) = ( QPR1(I)-QPR2(I) ) |
|
./pkg/fizhi/update_chemistry_exports.F |
383 |
INTERP_OZ |
. / ( PR1(I)- PR2(I) ) |
|
./pkg/fizhi/update_chemistry_exports.F |
384 |
INTERP_OZ |
ozrad(I,L2) = QPR2(I) + ( PLEVS(I,L2)-PR2(I) )*SLOPE(I) |
|
./pkg/fizhi/update_chemistry_exports.F |
398 |
INTERP_OZ |
DO 60 I=1,IRUN |
|
./pkg/fizhi/update_chemistry_exports.F |
401 |
INTERP_OZ |
ozrad (I,L2) = ozrad(I,L2) * VOLTOMAS |
|
./pkg/fizhi/update_earth_exports.F |
242 |
SIBALB |
INTEGER I, LAI |
|
./pkg/fizhi/update_earth_exports.F |
257 |
SIBALB |
DATA (ALVDR (I, 1, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
261 |
SIBALB |
DATA (ALVDR (I, 2, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
265 |
SIBALB |
DATA (ALVDR (I, 1, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
269 |
SIBALB |
DATA (ALVDR (I, 2, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
273 |
SIBALB |
DATA (ALVDR (I, 1, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
277 |
SIBALB |
DATA (ALVDR (I, 2, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
281 |
SIBALB |
DATA (ALVDR (I, 1, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
286 |
SIBALB |
DATA (ALVDR (I, 2, 4), I = 1, 14) /14*0.1637/ |
|
./pkg/fizhi/update_earth_exports.F |
289 |
SIBALB |
DATA (ALVDR (I, 1, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
293 |
SIBALB |
DATA (ALVDR (I, 2, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
297 |
SIBALB |
DATA (ALVDR (I, 1, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
301 |
SIBALB |
DATA (ALVDR (I, 2, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
305 |
SIBALB |
DATA (ALVDR (I, 1, 7), I = 1, 14) /14*ALVDRS/ |
|
./pkg/fizhi/update_earth_exports.F |
306 |
SIBALB |
DATA (ALVDR (I, 2, 7), I = 1, 14) /14*ALVDRS/ |
|
./pkg/fizhi/update_earth_exports.F |
309 |
SIBALB |
DATA (ALVDR (I, 1, 8), I = 1, 14) /14*ALVDRDL/ |
|
./pkg/fizhi/update_earth_exports.F |
310 |
SIBALB |
DATA (ALVDR (I, 2, 8), I = 1, 14) /14*ALVDRDL/ |
|
./pkg/fizhi/update_earth_exports.F |
313 |
SIBALB |
DATA (ALVDR (I, 1, 9), I = 1, 14) /14*ALVDRI/ |
|
./pkg/fizhi/update_earth_exports.F |
314 |
SIBALB |
DATA (ALVDR (I, 2, 9), I = 1, 14) /14*ALVDRI/ |
|
./pkg/fizhi/update_earth_exports.F |
317 |
SIBALB |
DATA (ALVDR (I, 1, 10), I = 1, 14) /14*ALVDRDD/ |
|
./pkg/fizhi/update_earth_exports.F |
318 |
SIBALB |
DATA (ALVDR (I, 2, 10), I = 1, 14) /14*ALVDRDD/ |
|
./pkg/fizhi/update_earth_exports.F |
320 |
SIBALB |
DATA (BTVDR (I, 1, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
324 |
SIBALB |
DATA (BTVDR (I, 2, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
328 |
SIBALB |
DATA (BTVDR (I, 1, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
332 |
SIBALB |
DATA (BTVDR (I, 2, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
336 |
SIBALB |
DATA (BTVDR (I, 1, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
340 |
SIBALB |
DATA (BTVDR (I, 2, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
344 |
SIBALB |
DATA (BTVDR (I, 1, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
348 |
SIBALB |
DATA (BTVDR (I, 2, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
352 |
SIBALB |
DATA (BTVDR (I, 1, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
357 |
SIBALB |
DATA (BTVDR (I, 2, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
362 |
SIBALB |
DATA (BTVDR (I, 1, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
367 |
SIBALB |
DATA (BTVDR (I, 2, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
372 |
SIBALB |
DATA (BTVDR (I, 1, 7), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
373 |
SIBALB |
DATA (BTVDR (I, 2, 7), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
375 |
SIBALB |
DATA (BTVDR (I, 1, 8), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
376 |
SIBALB |
DATA (BTVDR (I, 2, 8), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
378 |
SIBALB |
DATA (BTVDR (I, 1, 9), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
379 |
SIBALB |
DATA (BTVDR (I, 2, 9), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
381 |
SIBALB |
DATA (BTVDR (I, 1, 10), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
382 |
SIBALB |
DATA (BTVDR (I, 2, 10), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
385 |
SIBALB |
DATA (GMVDR (I, 1, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
389 |
SIBALB |
DATA (GMVDR (I, 2, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
393 |
SIBALB |
DATA (GMVDR (I, 1, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
397 |
SIBALB |
DATA (GMVDR (I, 2, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
401 |
SIBALB |
DATA (GMVDR (I, 1, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
405 |
SIBALB |
DATA (GMVDR (I, 2, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
409 |
SIBALB |
DATA (GMVDR (I, 1, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
413 |
SIBALB |
DATA (GMVDR (I, 2, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
417 |
SIBALB |
DATA (GMVDR (I, 1, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
422 |
SIBALB |
DATA (GMVDR (I, 2, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
427 |
SIBALB |
DATA (GMVDR (I, 1, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
432 |
SIBALB |
DATA (GMVDR (I, 2, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
437 |
SIBALB |
DATA (GMVDR (I, 1, 7), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
438 |
SIBALB |
DATA (GMVDR (I, 2, 7), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
440 |
SIBALB |
DATA (GMVDR (I, 1, 8), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
441 |
SIBALB |
DATA (GMVDR (I, 2, 8), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
443 |
SIBALB |
DATA (GMVDR (I, 1, 9), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
444 |
SIBALB |
DATA (GMVDR (I, 2, 9), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
446 |
SIBALB |
DATA (GMVDR (I, 1, 10), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
447 |
SIBALB |
DATA (GMVDR (I, 2, 10), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
451 |
SIBALB |
DATA (ALIDR (I, 1, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
455 |
SIBALB |
DATA (ALIDR (I, 2, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
457 |
SIBALB |
DATA (ALIDR (I, 1, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
459 |
SIBALB |
DATA (ALIDR (I, 2, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
461 |
SIBALB |
DATA (ALIDR (I, 1, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
463 |
SIBALB |
DATA (ALIDR (I, 2, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
467 |
SIBALB |
DATA (ALIDR (I, 1, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
471 |
SIBALB |
DATA (ALIDR (I, 2, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
475 |
SIBALB |
DATA (ALIDR (I, 1, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
477 |
SIBALB |
DATA (ALIDR (I, 2, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
481 |
SIBALB |
DATA (ALIDR (I, 1, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
483 |
SIBALB |
DATA (ALIDR (I, 2, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
488 |
SIBALB |
DATA (ALIDR (I, 1, 7), I = 1, 14) /14*ALIDRS/ |
|
./pkg/fizhi/update_earth_exports.F |
489 |
SIBALB |
DATA (ALIDR (I, 2, 7), I = 1, 14) /14*ALIDRS/ |
|
./pkg/fizhi/update_earth_exports.F |
491 |
SIBALB |
DATA (ALIDR (I, 1, 8), I = 1, 14) /14*ALIDRDL/ |
|
./pkg/fizhi/update_earth_exports.F |
492 |
SIBALB |
DATA (ALIDR (I, 2, 8), I = 1, 14) /14*ALIDRDL/ |
|
./pkg/fizhi/update_earth_exports.F |
494 |
SIBALB |
DATA (ALIDR (I, 1, 9), I = 1, 14) /14*ALIDRI/ |
|
./pkg/fizhi/update_earth_exports.F |
495 |
SIBALB |
DATA (ALIDR (I, 2, 9), I = 1, 14) /14*ALIDRI/ |
|
./pkg/fizhi/update_earth_exports.F |
497 |
SIBALB |
DATA (ALIDR (I, 1, 10), I = 1, 14) /14*ALIDRDD/ |
|
./pkg/fizhi/update_earth_exports.F |
498 |
SIBALB |
DATA (ALIDR (I, 2, 10), I = 1, 14) /14*ALIDRDD/ |
|
./pkg/fizhi/update_earth_exports.F |
501 |
SIBALB |
DATA (BTIDR (I, 1, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
505 |
SIBALB |
DATA (BTIDR (I, 2, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
509 |
SIBALB |
DATA (BTIDR (I, 1, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
513 |
SIBALB |
DATA (BTIDR (I, 2, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
517 |
SIBALB |
DATA (BTIDR (I, 1, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
521 |
SIBALB |
DATA (BTIDR (I, 2, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
525 |
SIBALB |
DATA (BTIDR (I, 1, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
529 |
SIBALB |
DATA (BTIDR (I, 2, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
533 |
SIBALB |
DATA (BTIDR (I, 1, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
538 |
SIBALB |
DATA (BTIDR (I, 2, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
543 |
SIBALB |
DATA (BTIDR (I, 1, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
548 |
SIBALB |
DATA (BTIDR (I, 2, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
554 |
SIBALB |
DATA (BTIDR (I, 1, 7), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
555 |
SIBALB |
DATA (BTIDR (I, 2, 7), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
557 |
SIBALB |
DATA (BTIDR (I, 1, 8), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
558 |
SIBALB |
DATA (BTIDR (I, 2, 8), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
560 |
SIBALB |
DATA (BTIDR (I, 1, 9), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
561 |
SIBALB |
DATA (BTIDR (I, 2, 9), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
563 |
SIBALB |
DATA (BTIDR (I, 1, 10), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
564 |
SIBALB |
DATA (BTIDR (I, 2, 10), I = 1, 14) /14*0./ |
|
./pkg/fizhi/update_earth_exports.F |
567 |
SIBALB |
DATA (GMIDR (I, 1, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
571 |
SIBALB |
DATA (GMIDR (I, 2, 1), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
575 |
SIBALB |
DATA (GMIDR (I, 1, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
579 |
SIBALB |
DATA (GMIDR (I, 2, 2), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
583 |
SIBALB |
DATA (GMIDR (I, 1, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
587 |
SIBALB |
DATA (GMIDR (I, 2, 3), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
591 |
SIBALB |
DATA (GMIDR (I, 1, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
595 |
SIBALB |
DATA (GMIDR (I, 2, 4), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
599 |
SIBALB |
DATA (GMIDR (I, 1, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
604 |
SIBALB |
DATA (GMIDR (I, 2, 5), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
609 |
SIBALB |
DATA (GMIDR (I, 1, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
614 |
SIBALB |
DATA (GMIDR (I, 2, 6), I = 1, 14) |
|
./pkg/fizhi/update_earth_exports.F |
619 |
SIBALB |
DATA (GMIDR (I, 1, 7), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
620 |
SIBALB |
DATA (GMIDR (I, 2, 7), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
622 |
SIBALB |
DATA (GMIDR (I, 1, 8), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
623 |
SIBALB |
DATA (GMIDR (I, 2, 8), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
625 |
SIBALB |
DATA (GMIDR (I, 1, 9), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
626 |
SIBALB |
DATA (GMIDR (I, 2, 9), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
628 |
SIBALB |
DATA (GMIDR (I, 1, 10), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
629 |
SIBALB |
DATA (GMIDR (I, 2, 10), I = 1, 14) /14*1./ |
|
./pkg/fizhi/update_earth_exports.F |
654 |
SIBALB |
DO 100 I=1,IRUN |
|
./pkg/fizhi/update_earth_exports.F |
655 |
SIBALB |
ALA = MIN (MAX (ZERO, VLAI(I)), ALATRM) |
|
./pkg/fizhi/update_earth_exports.F |
658 |
SIBALB |
DY = (VGRN(I)- GRN(1)) * (ONE/(GRN(2) - GRN(1))) |
|
./pkg/fizhi/update_earth_exports.F |
660 |
SIBALB |
ALPHA = COEFF (ALVDR (1, 1, ITYP (I)), NLAI, LAI ,DX, DY) |
|
./pkg/fizhi/update_earth_exports.F |
661 |
SIBALB |
BETA = COEFF (BTVDR (1, 1, ITYP (I)), NLAI, LAI ,DX, DY) |
|
./pkg/fizhi/update_earth_exports.F |
662 |
SIBALB |
GAMMA = COEFF (GMVDR (1, 1, ITYP (I)), NLAI, LAI ,DX, DY) |
|
./pkg/fizhi/update_earth_exports.F |
664 |
SIBALB |
AVISDR(I) = ALPHA - ZTH(I)*BETA / (GAMMA+ZTH(I)) |
|
./pkg/fizhi/update_earth_exports.F |
665 |
SIBALB |
AVISDF(I) = ALPHA-BETA |
|
./pkg/fizhi/update_earth_exports.F |
668 |
SIBALB |
ALPHA = COEFF (ALIDR (1, 1, ITYP (I)), NLAI, LAI ,DX, DY) |
|
./pkg/fizhi/update_earth_exports.F |
669 |
SIBALB |
BETA = COEFF (BTIDR (1, 1, ITYP (I)), NLAI, LAI ,DX, DY) |
|
./pkg/fizhi/update_earth_exports.F |
670 |
SIBALB |
GAMMA = COEFF (GMIDR (1, 1, ITYP (I)), NLAI, LAI ,DX, DY) |
|
./pkg/fizhi/update_earth_exports.F |
672 |
SIBALB |
ANIRDR(I) = ALPHA - ZTH(I)*BETA / (GAMMA+ZTH(I)) |
|
./pkg/fizhi/update_earth_exports.F |
673 |
SIBALB |
ANIRDF(I) = ALPHA-BETA |
|
./pkg/fizhi/update_earth_exports.F |
676 |
SIBALB |
IF (SNW (I) .GT. ZERO) THEN |
|
./pkg/fizhi/update_earth_exports.F |
677 |
SIBALB |
FAC = SNW(I) / (SNW(I) + SNWMID(ITYP(I))) |
|
./pkg/fizhi/update_earth_exports.F |
679 |
SIBALB |
AVISDR(I) = AVISDR(I) + (SNWALB(1,ITYP(I)) - AVISDR(I)) * FAC |
|
./pkg/fizhi/update_earth_exports.F |
680 |
SIBALB |
ANIRDR(I) = ANIRDR(I) + (SNWALB(2,ITYP(I)) - ANIRDR(I)) * FAC |
|
./pkg/fizhi/update_earth_exports.F |
681 |
SIBALB |
AVISDF(I) = AVISDF(I) + (SNWALB(3,ITYP(I)) - AVISDF(I)) * FAC |
|
./pkg/fizhi/update_earth_exports.F |
682 |
SIBALB |
ANIRDF(I) = ANIRDF(I) + (SNWALB(4,ITYP(I)) - ANIRDF(I)) * FAC |
|
./pkg/fizhi/update_earth_exports.F |
798 |
GETLGR |
DO 220 I=1,NCHPS |
|
./pkg/fizhi/update_earth_exports.F |
800 |
GETLGR |
IF(ALAT(I,bi,bj).GT.0.) THEN |
|
./pkg/fizhi/update_earth_exports.F |
808 |
GETLGR |
ALAI(I,bi,bj) = VGLA(KK2,ITYP(I,bi,bj))*FAC+ |
|
./pkg/fizhi/update_earth_exports.F |
809 |
GETLGR |
& VGLA(KK1,ITYP(I,bi,bj))*(ONE-FAC) |
|
./pkg/fizhi/update_earth_exports.F |
810 |
GETLGR |
AGRN(I,bi,bj) = VGGR(KK2,ITYP(I,bi,bj))*FAC+ |
|
./pkg/fizhi/update_earth_exports.F |
811 |
GETLGR |
& VGGR(KK1,ITYP(I,bi,bj))*(ONE-FAC) |
|
./pkg/fizhi/update_earth_exports.F |
887 |
GETALB |
DO I=1,IM |
|
./pkg/fizhi/update_earth_exports.F |
889 |
GETALB |
ALBOC(I,J) = A0 + (A1 + (A2 + A3*cosz(I,J))*cosz(I,J))*cosz(I,J) |
|
./pkg/fizhi/update_earth_exports.F |
890 |
GETALB |
ALBVR(I,J,bi,bj) = ALBSI*FRACI(I,J) + ALBOC(I,J)*(ONE-FRACI(I,J)) |
|
./pkg/fizhi/update_earth_exports.F |
891 |
GETALB |
ALBNR(I,J,bi,bj) = ALBVR(I,J,bi,bj) |
|
./pkg/fizhi/update_earth_exports.F |
892 |
GETALB |
ALBVF(I,J,bi,bj) = ALBSI * FRACI(I,J) + OCNALB * (ONE-FRACI(I,J)) |
|
./pkg/fizhi/update_earth_exports.F |
893 |
GETALB |
ALBNF(I,J,bi,bj) = ALBVF(I,J,bi,bj) |
|
./pkg/fizhi/update_earth_exports.F |
908 |
GETALB |
DO I=1,IM |
|
./pkg/fizhi/update_earth_exports.F |
916 |
GETALB |
DO I=1,IM |
|
./pkg/fizhi/update_earth_exports.F |
921 |
GETALB |
DO I=1,IM |
|
./pkg/fizhi/update_earth_exports.F |
928 |
GETALB |
DO I=1,IM |
|
./pkg/fizhi/update_earth_exports.F |
933 |
GETALB |
DO I=1,IM |
|
./pkg/fizhi/update_earth_exports.F |
940 |
GETALB |
DO I=1,IM |
|
./pkg/fizhi/update_earth_exports.F |
945 |
GETALB |
DO I=1,IM |
|
./pkg/fizhi/update_earth_exports.F |
952 |
GETALB |
DO I=1,IM |
|
./pkg/flt/exch2_send_put_vec.F |
127 |
EXCH2_SEND_PUT_VEC_RL |
DO I=1,myd1 |
|
./pkg/flt/exch2_send_put_vec.F |
128 |
EXCH2_SEND_PUT_VEC_RL |
bufRec(I,sBi,sBj,sDir) = array(I,bi,bj,dir) |
|
./pkg/flt/exch2_send_put_vec.F |
67 |
EXCH2_SEND_PUT_VEC_RL |
INTEGER I,J,sBi,sBj, nN, thisTile, tgT, tgN, dir, sDir |
|
./pkg/flt/exch_send_put_vec.F |
143 |
EXCH_SEND_PUT_VEC_X_RL |
DO I=1,myd1 |
|
./pkg/flt/exch_send_put_vec.F |
144 |
EXCH_SEND_PUT_VEC_X_RL |
bufRecE(I,biW,bjW) = arrayW(I,bi,bj) |
|
./pkg/flt/exch_send_put_vec.F |
178 |
EXCH_SEND_PUT_VEC_X_RL |
DO I=1,myd1 |
|
./pkg/flt/exch_send_put_vec.F |
179 |
EXCH_SEND_PUT_VEC_X_RL |
bufRecW(I,biE,bjE) = arrayE(I,bi,bj) |
|
./pkg/flt/exch_send_put_vec.F |
289 |
EXCH_SEND_PUT_VEC_Y_RL |
INTEGER I |
|
./pkg/flt/exch_send_put_vec.F |
360 |
EXCH_SEND_PUT_VEC_Y_RL |
DO I=1,myd1 |
|
./pkg/flt/exch_send_put_vec.F |
361 |
EXCH_SEND_PUT_VEC_Y_RL |
bufRecN(I,biS,bjS) = arrayS(I,bi,bj) |
|
./pkg/flt/exch_send_put_vec.F |
386 |
EXCH_SEND_PUT_VEC_Y_RL |
DO I=1,myd1 |
|
./pkg/flt/exch_send_put_vec.F |
387 |
EXCH_SEND_PUT_VEC_Y_RL |
bufRecS(I,biN,bjN) = arrayN(I,bi,bj) |
|
./pkg/flt/exch_send_put_vec.F |
64 |
EXCH_SEND_PUT_VEC_X_RL |
INTEGER I |
|
./pkg/generic_advdiff/gad_grad_y.F |
60 |
GAD_GRAD_Y |
dTdy(I,1-Oly)=0. |
|
./pkg/ggl90/ggl90_calc.F |
158 |
GGL90_CALC |
recip_hFacI(I,J,K)=0. _d 0 |
|
./pkg/ggl90/ggl90_calc.F |
160 |
GGL90_CALC |
& recip_hFacI(I,J,K)=1. _d 0/hFac |
|
./pkg/ggl90/ggl90_init_varia.F |
105 |
GGL90_INIT_VARIA |
DO I=1-OLx,sNx+OLx |
|
./pkg/ggl90/ggl90_init_varia.F |
124 |
GGL90_INIT_VARIA |
DO I=1-OLx,sNx+OLx |
|
./pkg/ggl90/ggl90_init_varia.F |
149 |
GGL90_INIT_VARIA |
DO I=1-OLx,sNx+OLx |
|
./pkg/ggl90/ggl90_init_varia.F |
150 |
GGL90_INIT_VARIA |
GGL90TKE(I,J,K,bi,bj) = MAX(GGL90TKE(I,J,K,bi,bj), |
|
./pkg/ggl90/ggl90_init_varia.F |
151 |
GGL90_INIT_VARIA |
& GGL90TKEmin)*maskC(I,J,K,bi,bj) |
|
./pkg/ggl90/ggl90_init_varia.F |
37 |
GGL90_INIT_VARIA |
INTEGER I, J, K |
|
./pkg/ggl90/ggl90_init_varia.F |
49 |
GGL90_INIT_VARIA |
DO I=1-OLx,sNx+OLx |
|
./pkg/ggl90/ggl90_init_varia.F |
50 |
GGL90_INIT_VARIA |
GGL90viscArU(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
51 |
GGL90_INIT_VARIA |
GGL90viscArV(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
52 |
GGL90_INIT_VARIA |
GGL90diffKr(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
54 |
GGL90_INIT_VARIA |
GGL90TKE(I,J,K,bi,bj)=GGL90eps*maskC(I,J,K,bi,bj) |
|
./pkg/ggl90/ggl90_init_varia.F |
57 |
GGL90_INIT_VARIA |
GGL90TKE(I,J,K,bi,bj)=GGL90TKEmin*maskC(I,J,K,bi,bj) |
|
./pkg/ggl90/ggl90_init_varia.F |
78 |
GGL90_INIT_VARIA |
DO I=1-OLx,sNx+OLx |
|
./pkg/ggl90/ggl90_init_varia.F |
79 |
GGL90_INIT_VARIA |
IDEMIX_E(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
80 |
GGL90_INIT_VARIA |
IDEMIX_V0(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
81 |
GGL90_INIT_VARIA |
IDEMIX_tau_d(I,J,K,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
87 |
GGL90_INIT_VARIA |
DO I=1-OLx,sNx+OLx |
|
./pkg/ggl90/ggl90_init_varia.F |
88 |
GGL90_INIT_VARIA |
IDEMIX_F_B(I,J,bi,bj) = 0. _d 0 |
|
./pkg/ggl90/ggl90_init_varia.F |
89 |
GGL90_INIT_VARIA |
IDEMIX_F_S(I,J,bi,bj) = 0. _d 0 |
|
./pkg/grdchk/grdchk_get_position.F |
171 |
GRDCHK_GET_POSITION |
& .and. j.eq. OB_Jn(I,bi,bj) ) then |
|
./pkg/grdchk/grdchk_get_position.F |
179 |
GRDCHK_GET_POSITION |
& .and. j.eq. OB_Js(I,bi,bj) ) then |
|
./pkg/grdchk/grdchk_loc.F |
222 |
GRDCHK_LOC |
& .and. j.eq. OB_Jn(I,bi,bj) ) then |
|
./pkg/grdchk/grdchk_loc.F |
225 |
GRDCHK_LOC |
jtmp = OB_Jn(I,bi,bj) |
|
./pkg/grdchk/grdchk_loc.F |
231 |
GRDCHK_LOC |
& .and. j.eq. OB_Js(I,bi,bj) ) then |
|
./pkg/grdchk/grdchk_loc.F |
234 |
GRDCHK_LOC |
jtmp = OB_Js(I,bi,bj) |
|
./pkg/icefront/icefront_init_fixed.F |
100 |
ICEFRONT_INIT_FIXED |
& K_icefront(I,J,bi,bj) = K |
|
./pkg/icefront/icefront_init_fixed.F |
111 |
ICEFRONT_INIT_FIXED |
DO I = 1, sNx |
|
./pkg/icefront/icefront_init_fixed.F |
113 |
ICEFRONT_INIT_FIXED |
IF (ICEFRONTlength(I,J,bi,bj) .GT. 0. _d 0) THEN |
|
./pkg/icefront/icefront_init_fixed.F |
115 |
ICEFRONT_INIT_FIXED |
IsInterface=Isinterface + K_icefront(I,J,bi,bj) |
|
./pkg/icefront/icefront_init_fixed.F |
37 |
ICEFRONT_INIT_FIXED |
INTEGER I, J, K, bi, bj |
|
./pkg/icefront/icefront_init_fixed.F |
86 |
ICEFRONT_INIT_FIXED |
DO I = 1-OLx, sNx+OLx |
|
./pkg/icefront/icefront_init_fixed.F |
87 |
ICEFRONT_INIT_FIXED |
R_icefront(I,J,bi,bj) = ABS(R_icefront(I,J,bi,bj)) |
|
./pkg/icefront/icefront_init_fixed.F |
96 |
ICEFRONT_INIT_FIXED |
DO I = 1-OLx, sNx+OLx |
|
./pkg/icefront/icefront_init_fixed.F |
99 |
ICEFRONT_INIT_FIXED |
IF ( R_icefront(I,J,bi,bj) .GT. ABS(rF(K))) |
|
./pkg/icefront/icefront_init_varia.F |
31 |
ICEFRONT_INIT_VARIA |
INTEGER I, J, K, bi, bj |
|
./pkg/icefront/icefront_init_varia.F |
38 |
ICEFRONT_INIT_VARIA |
DO I = 1, sNx |
|
./pkg/icefront/icefront_thermodynamics.F |
100 |
ICEFRONT_THERMODYNAMICS |
& *ICEFRONTheatTransCoeff(I,J) |
|
./pkg/icefront/icefront_thermodynamics.F |
103 |
ICEFRONT_THERMODYNAMICS |
eps1 = rUnit2mass*HeatCapacity_Cp*ICEFRONTheatTransCoeff(I,J) |
|
./pkg/icefront/icefront_thermodynamics.F |
104 |
ICEFRONT_THERMODYNAMICS |
eps2 = rUnit2mass*ICEFRONTlatentHeat*ICEFRONTsaltTransCoeff(I,J) |
|
./pkg/icefront/icefront_thermodynamics.F |
106 |
ICEFRONT_THERMODYNAMICS |
& *ICEFRONTsaltTransCoeff(I,J) |
|
./pkg/icefront/icefront_thermodynamics.F |
113 |
ICEFRONT_THERMODYNAMICS |
tLoc = theta(I,J,K,bi,bj) |
|
./pkg/icefront/icefront_thermodynamics.F |
114 |
ICEFRONT_THERMODYNAMICS |
sLoc = MAX(salt(I,J,K,bi,bj), 0. _d 0) |
|
./pkg/icefront/icefront_thermodynamics.F |
121 |
ICEFRONT_THERMODYNAMICS |
IF ( K .EQ. K_icefront(I,J,bi,bj)) THEN |
|
./pkg/icefront/icefront_thermodynamics.F |
122 |
ICEFRONT_THERMODYNAMICS |
pLoc = 0.5*(ABS(R_icefront(I,J,bi,bj))+ABS(rF(K))) |
|
./pkg/icefront/icefront_thermodynamics.F |
125 |
ICEFRONT_THERMODYNAMICS |
& (R_icefront(I,J,bi,bj)-pLoc) |
|
./pkg/icefront/icefront_thermodynamics.F |
126 |
ICEFRONT_THERMODYNAMICS |
& / R_icefront(I,J,bi,bj) |
|
./pkg/icefront/icefront_thermodynamics.F |
146 |
ICEFRONT_THERMODYNAMICS |
FreshWaterFlux(I,J) = maskC(I,J,K,bi,bj) * |
|
./pkg/icefront/icefront_thermodynamics.F |
150 |
ICEFRONT_THERMODYNAMICS |
HeatFlux(I,J) = maskC(I,J,K,bi,bj) * HeatCapacity_Cp * |
|
./pkg/icefront/icefront_thermodynamics.F |
151 |
ICEFRONT_THERMODYNAMICS |
& ( -rUnit2mass*ICEFRONTheatTransCoeff(I,J) + |
|
./pkg/icefront/icefront_thermodynamics.F |
152 |
ICEFRONT_THERMODYNAMICS |
& FreshWaterFlux(I,J) ) * ( thetaFreeze - tLoc ) |
|
./pkg/icefront/icefront_thermodynamics.F |
155 |
ICEFRONT_THERMODYNAMICS |
icefront_TendT(i,j,K,bi,bj) = - HeatFlux(I,J)* eps5 |
|
./pkg/icefront/icefront_thermodynamics.F |
156 |
ICEFRONT_THERMODYNAMICS |
icefront_TendS(i,j,K,bi,bj) = FreshWaterFlux(I,J) * |
|
./pkg/icefront/icefront_thermodynamics.F |
171 |
ICEFRONT_THERMODYNAMICS |
& * (ABS(R_icefront(I,J,bi,bj))-ABS(rF(K))) |
|
./pkg/icefront/icefront_thermodynamics.F |
175 |
ICEFRONT_THERMODYNAMICS |
& * (ABS(R_icefront(I,J,bi,bj))-ABS(rF(K))) |
|
./pkg/icefront/icefront_thermodynamics.F |
181 |
ICEFRONT_THERMODYNAMICS |
HeatFlux (I,J) = 0. _d 0 |
|
./pkg/icefront/icefront_thermodynamics.F |
182 |
ICEFRONT_THERMODYNAMICS |
FreshWaterFlux(I,J) = 0. _d 0 |
|
./pkg/icefront/icefront_thermodynamics.F |
60 |
ICEFRONT_THERMODYNAMICS |
INTEGER I,J,K |
|
./pkg/icefront/icefront_thermodynamics.F |
89 |
ICEFRONT_THERMODYNAMICS |
DO I = 1, sNx |
|
./pkg/icefront/icefront_thermodynamics.F |
92 |
ICEFRONT_THERMODYNAMICS |
IF( ICEFRONTlength(I,J,bi,bj) .GT. 0. _d 0 |
|
./pkg/icefront/icefront_thermodynamics.F |
93 |
ICEFRONT_THERMODYNAMICS |
& .AND. K .LE. K_icefront(I,J,bi,bj) ) THEN |
|
./pkg/icefront/icefront_thermodynamics.F |
94 |
ICEFRONT_THERMODYNAMICS |
ICEFRONTheatTransCoeff(I,J) = 1.0 _d -02 |
|
./pkg/icefront/icefront_thermodynamics.F |
95 |
ICEFRONT_THERMODYNAMICS |
& *abs(wVEL(I,J,K,bi,bj)) |
|
./pkg/icefront/icefront_thermodynamics.F |
97 |
ICEFRONT_THERMODYNAMICS |
ICEFRONTheatTransCoeff(I,J) = max |
|
./pkg/icefront/icefront_thermodynamics.F |
98 |
ICEFRONT_THERMODYNAMICS |
& (ICEFRONTheatTransCoeff(I,J),1. _d -04) |
|
./pkg/icefront/icefront_thermodynamics.F |
99 |
ICEFRONT_THERMODYNAMICS |
ICEFRONTsaltTransCoeff(I,J) = 5.05 _d -3 |
|
./pkg/kl10/kl10_calc.F |
105 |
KL10_CALC |
IF ( (rhoS(K).LT.rhoS(K-1)).AND.(maskC(I,J,K,bi |
|
./pkg/kl10/kl10_calc.F |
122 |
KL10_CALC |
KLdiffKr(I,J,1,bi,bj) = MAX(KLviscAr(I,J,1,bi,bj), |
|
./pkg/kl10/kl10_calc.F |
124 |
KL10_CALC |
& diffKr(I,J,1,bi,bj) ) |
|
./pkg/kl10/kl10_calc.F |
165 |
KL10_CALC |
KLviscAr(I,J,K,bi,bj) = MAX(KTemp,viscArNr(K)) |
|
./pkg/kl10/kl10_calc.F |
166 |
KL10_CALC |
KLdiffKr(I,J,K,bi,bj) = MAX(KTemp, |
|
./pkg/kl10/kl10_calc.F |
168 |
KL10_CALC |
& diffKr(I,J,K,bi,bj) ) |
|
./pkg/kl10/kl10_calc.F |
187 |
KL10_CALC |
IF (hfacW(I+di,J-1,K,bi,bj).GT.0.000001) THEN |
|
./pkg/kl10/kl10_calc.F |
188 |
KL10_CALC |
dz = 0.5*(drF(K)*hfacW(I+di,J-1,K,bi,bj) |
|
./pkg/kl10/kl10_calc.F |
189 |
KL10_CALC |
& +drF(Km1)*hfacW(I+di,J-1,Km1,bi,bj)) |
|
./pkg/kl10/kl10_calc.F |
191 |
KL10_CALC |
tempu = (uVel(I+di,J-1,Km1,bi,bj)-uVel(I+di,J |
|
./pkg/kl10/kl10_calc.F |
193 |
KL10_CALC |
epsw(di)=tempu*tempu*KLviscAr(I+di,J-1,K,bi |
|
./pkg/kl10/kl10_calc.F |
201 |
KL10_CALC |
dz=0.5*(drF(Km1)*hfacW(I+di,J-1,Km1,bi ,bj)) |
|
./pkg/kl10/kl10_calc.F |
206 |
KL10_CALC |
IF (hfacS(I-1,J+di,K,bi,bj).GT.0.000001) THEN |
|
./pkg/kl10/kl10_calc.F |
207 |
KL10_CALC |
dz = 0.5*(drF(K)*hfacS(I-1,J+di,K,bi,bj) |
|
./pkg/kl10/kl10_calc.F |
208 |
KL10_CALC |
& +drF(Km1)*hfacS(I-1,J+di,Km1,bi,bj)) |
|
./pkg/kl10/kl10_calc.F |
210 |
KL10_CALC |
tempu = (vVel(I-1,J+di,Km1,bi,bj)-vVel(I-1,J |
|
./pkg/kl10/kl10_calc.F |
212 |
KL10_CALC |
epss(di)=tempu*tempu*KLviscAr(I-1,J+di,K,bi |
|
./pkg/kl10/kl10_calc.F |
220 |
KL10_CALC |
dz=+0.5*(drF(Km1)*hfacS(I-1,J+di,Km1 ,bi,bj)) |
|
./pkg/kl10/kl10_calc.F |
242 |
KL10_CALC |
KLeps(I-1,J-1,Km1,bi,bj) = 0.5*(ep+ec) |
|
./pkg/kl10/kl10_calc.F |
244 |
KL10_CALC |
KLeps(I-1,J-1,Km1,bi,bj) = ec |
|
./pkg/kl10/kl10_calc.F |
259 |
KL10_CALC |
KLeps(I-1,J-1,Nr,bi,bj) =ep |
|
./pkg/kl10/kl10_calc.F |
61 |
KL10_CALC |
INTEGER I, J, K, Km1, JJ |
|
./pkg/kl10/kl10_calc.F |
80 |
KL10_CALC |
DO I=iMin,iMax |
|
./pkg/kl10/kl10_calc.F |
82 |
KL10_CALC |
CALL FIND_RHO_SCALAR(theta(I,J,K,bi,bj), salt(I,J,K,bi,bj), |
|
./pkg/kl10/kl10_calc.F |
83 |
KL10_CALC |
& totPhiHyd(I,J,K,bi,bj),rhot,myThid ) |
|
./pkg/kl10/kl10_calc.F |
87 |
KL10_CALC |
KLeps(I-1,J-1,1,bi,bj)=0.0 |
|
./pkg/kl10/kl10_calc.F |
92 |
KL10_CALC |
KLviscAr(I,J,1,bi,bj) = viscArNr(1) |
|
./pkg/kl10/kl10_calc.F |
93 |
KL10_CALC |
KLviscold = KLviscAr(I,J,1,bi,bj) ! at previous cell center |
|
./pkg/kl10/kl10_calc.F |
97 |
KL10_CALC |
CALL FIND_RHO_SCALAR(theta(I,J,K,bi,bj), salt(I,J,K,bi |
|
./pkg/kl10/kl10_calc.F |
98 |
KL10_CALC |
& ,bj),totPhiHyd(I,J,K,bi,bj),rhot,myThid ) |
|
./pkg/kpp/kpp_forcing_surf.F |
168 |
KPP_FORCING_SURF |
bo (I,J) = p0 |
|
./pkg/kpp/kpp_forcing_surf.F |
169 |
KPP_FORCING_SURF |
bosol(I,J) = p0 |
|
./pkg/kpp/kpp_forcing_surf.F |
172 |
KPP_FORCING_SURF |
boplume(I,J,k) = p0 |
|
./pkg/kpp/kpp_forcing_surf.F |
174 |
KPP_FORCING_SURF |
boplume(I,J,Nrp1) = p0 |
|
./pkg/kpp/kpp_forcing_surf.F |
212 |
KPP_FORCING_SURF |
bo(I,J) = - gravity * |
|
./pkg/kpp/kpp_forcing_surf.F |
213 |
KPP_FORCING_SURF |
& ( TTALPHA(I,J,1) * (surfForcT(i,j,bi,bj)+ |
|
./pkg/kpp/kpp_forcing_surf.F |
215 |
KPP_FORCING_SURF |
& SSBETA(I,J,1) * surfForcS(i,j,bi,bj) ) |
|
./pkg/kpp/kpp_forcing_surf.F |
216 |
KPP_FORCING_SURF |
& / rhoSurf(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
217 |
KPP_FORCING_SURF |
bosol(I,J) = gravity * TTALPHA(I,J,1) * Qsw(i,j,bi,bj) * |
|
./pkg/kpp/kpp_forcing_surf.F |
219 |
KPP_FORCING_SURF |
& / rhoSurf(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
232 |
KPP_FORCING_SURF |
temparray(I,J) = - gravity * |
|
./pkg/kpp/kpp_forcing_surf.F |
233 |
KPP_FORCING_SURF |
& ( SSBETA(I,J,k) * SPforcS(i,j,k) + |
|
./pkg/kpp/kpp_forcing_surf.F |
234 |
KPP_FORCING_SURF |
& TTALPHA(I,J,k)* SPforcT(i,j,k) / HeatCapacity_Cp ) |
|
./pkg/kpp/kpp_forcing_surf.F |
235 |
KPP_FORCING_SURF |
& * recip_rhoConst / rhoSurf(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
236 |
KPP_FORCING_SURF |
boplume(I,J,kp1) = boplume(I,J,k)+temparray(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
244 |
KPP_FORCING_SURF |
boplume(I,J,k) = 0. _d 0 |
|
./pkg/kpp/kpp_forcing_surf.F |
246 |
KPP_FORCING_SURF |
boplume(I,J,1) = - gravity * SSBETA(I,J,1) |
|
./pkg/kpp/kpp_forcing_surf.F |
248 |
KPP_FORCING_SURF |
& * recip_rhoConst / rhoSurf(I,J) |
|
./pkg/kpp/kpp_forcing_surf.F |
307 |
KPP_FORCING_SURF |
& maskC(I,J,k,bi,bj) .GT. 0. .AND. |
|
./pkg/kpp/kpp_init_varia.F |
45 |
KPP_INIT_VARIA |
nzmax(i,j,bi,bj) = kLowC(I,J,bi,bj) |
|
./pkg/kpp/kpp_routines.F |
1698 |
STATEKPP |
INTEGER I, J, K |
|
./pkg/kpp/kpp_routines.F |
1732 |
STATEKPP |
DO I = 1-OLx, sNx+OLx |
|
./pkg/kpp/kpp_routines.F |
1733 |
STATEKPP |
RHO1(I,J) = WORK1(I,J) + rhoConst |
|
./pkg/kpp/kpp_routines.F |
1734 |
STATEKPP |
TTALPHA(I,J,1) = WORK2(I,J) |
|
./pkg/kpp/kpp_routines.F |
1735 |
STATEKPP |
SSBETA(I,J,1) = WORK3(I,J) |
|
./pkg/kpp/kpp_routines.F |
1736 |
STATEKPP |
DBSFC(I,J,1) = 0. |
|
./pkg/kpp/kpp_routines.F |
1801 |
STATEKPP |
DO I = 1-OLx, sNx+OLx |
|
./pkg/kpp/kpp_routines.F |
1802 |
STATEKPP |
TTALPHA(I,J,K) = WORK1 (I,J) |
|
./pkg/kpp/kpp_routines.F |
1803 |
STATEKPP |
SSBETA(I,J,K) = WORK2 (I,J) |
|
./pkg/kpp/kpp_routines.F |
1804 |
STATEKPP |
DBLOC(I,J,K-1) = gravity * (RHOK(I,J) - RHOKM1(I,J)) / |
|
./pkg/kpp/kpp_routines.F |
1805 |
STATEKPP |
& (RHOK(I,J) + rhoConst) |
|
./pkg/kpp/kpp_routines.F |
1806 |
STATEKPP |
DBSFC(I,J,K) = gravity * (RHOK(I,J) - RHO1K (I,J)) / |
|
./pkg/kpp/kpp_routines.F |
1807 |
STATEKPP |
& (RHOK(I,J) + rhoConst) |
|
./pkg/kpp/kpp_routines.F |
1815 |
STATEKPP |
DO I = 1-OLx, sNx+OLx |
|
./pkg/kpp/kpp_routines.F |
1816 |
STATEKPP |
TTALPHA(I,J,Nrp1) = TTALPHA(I,J,Nr) |
|
./pkg/kpp/kpp_routines.F |
1817 |
STATEKPP |
SSBETA(I,J,Nrp1) = SSBETA(I,J,Nr) |
|
./pkg/kpp/kpp_routines.F |
1818 |
STATEKPP |
DBLOC(I,J,Nr) = 0. |
|
./pkg/kpp/kpp_routines.F |
1892 |
KPP_DOUBLEDIFF |
INTEGER I, J, K |
|
./pkg/kpp/kpp_routines.F |
1925 |
KPP_DOUBLEDIFF |
DO I = 1-OLx, sNx+OLx |
|
./pkg/kpp/kpp_routines.F |
1926 |
KPP_DOUBLEDIFF |
alphaDT(I,J) = ( theta(I,J,Km1,bi,bj)-theta(I,J,K,bi,bj) ) |
|
./pkg/kpp/kpp_routines.F |
1927 |
KPP_DOUBLEDIFF |
& * 0.5 _d 0 * ABS( TTALPHA(I,J,Km1) + TTALPHA(I,J,K) ) |
|
./pkg/kpp/kpp_routines.F |
1928 |
KPP_DOUBLEDIFF |
betaDS(I,J) = ( salt(I,J,Km1,bi,bj)-salt(I,J,K,bi,bj) ) |
|
./pkg/kpp/kpp_routines.F |
1929 |
KPP_DOUBLEDIFF |
& * 0.5 _d 0 * ( SSBETA(I,J,Km1) + SSBETA(I,J,K) ) |
|
./pkg/kpp/kpp_routines.F |
1930 |
KPP_DOUBLEDIFF |
nuddt(I,J) = 0. _d 0 |
|
./pkg/kpp/kpp_routines.F |
1931 |
KPP_DOUBLEDIFF |
nudds(I,J) = 0. _d 0 |
|
./pkg/kpp/kpp_routines.F |
1936 |
KPP_DOUBLEDIFF |
DO I = iMin, iMax |
|
./pkg/kpp/kpp_routines.F |
1941 |
KPP_DOUBLEDIFF |
IF ( alphaDT(I,J) .GT. betaDS(I,J) |
|
./pkg/kpp/kpp_routines.F |
1942 |
KPP_DOUBLEDIFF |
& .AND. betaDS(I,J) .GT. 0. _d 0 ) THEN |
|
./pkg/kpp/kpp_routines.F |
1943 |
KPP_DOUBLEDIFF |
Rrho = MIN( alphaDT(I,J)/betaDS(I,J), Rrho0 ) |
|
./pkg/kpp/kpp_routines.F |
1947 |
KPP_DOUBLEDIFF |
nudds(I,J) = dsfmax * nutmp * nutmp * nutmp |
|
./pkg/kpp/kpp_routines.F |
1949 |
KPP_DOUBLEDIFF |
nuddt(I,J) = 0.7 _d 0 * nudds(I,J) |
|
./pkg/kpp/kpp_routines.F |
1950 |
KPP_DOUBLEDIFF |
ELSEIF ( alphaDT(I,J) .LT. 0. _d 0 |
|
./pkg/kpp/kpp_routines.F |
1951 |
KPP_DOUBLEDIFF |
& .AND. betaDS(I,J) .LT. 0. _d 0 |
|
./pkg/kpp/kpp_routines.F |
1952 |
KPP_DOUBLEDIFF |
& .AND.alphaDT(I,J) .GT. betaDS(I,J) ) THEN |
|
./pkg/kpp/kpp_routines.F |
1957 |
KPP_DOUBLEDIFF |
Rrho = alphaDT(I,J)/betaDS(I,J) |
|
./pkg/kpp/kpp_routines.F |
1959 |
KPP_DOUBLEDIFF |
nuddt(I,J) = numol * 0.909 _d 0 |
|
./pkg/kpp/kpp_routines.F |
1966 |
KPP_DOUBLEDIFF |
nudds(I,J) = nuddt(I,J) * MAX( 0.15 _d 0 * Rrho, |
|
./pkg/kpp/kpp_routines.F |
1979 |
KPP_DOUBLEDIFF |
DO I = 1-OLx, sNx+OLx |
|
./pkg/kpp/kpp_routines.F |
1980 |
KPP_DOUBLEDIFF |
kappaRT(I,J,K) = kappaRT(I,J,K) + nuddt(I,J) |
|
./pkg/kpp/kpp_routines.F |
1981 |
KPP_DOUBLEDIFF |
kappaRS(I,J,K) = kappaRS(I,J,K) + nudds(I,J) |
|
./pkg/kpp/kpp_transport_s.F |
63 |
KPP_TRANSPORT_S |
INTEGER I, J |
|
./pkg/land/land_ini_vars.F |
64 |
LAND_INI_VARS |
DO I=1-Olx,sNx+Olx |
|
./pkg/land/land_ini_vars.F |
74 |
LAND_INI_VARS |
DO I=1-Olx,sNx+Olx |
|
./pkg/mdsio/mdsio_gl.F |
1151 |
MDSREADFIELD_2D_GL |
DO I=1,Nx |
|
./pkg/mdsio/mdsio_gl.F |
1152 |
MDSREADFIELD_2D_GL |
global(I,J) = xy_buffer_r4(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
1161 |
MDSREADFIELD_2D_GL |
DO I=1,Nx |
|
./pkg/mdsio/mdsio_gl.F |
1162 |
MDSREADFIELD_2D_GL |
global(I,J) = xy_buffer_r8(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
1178 |
MDSREADFIELD_2D_GL |
DO I=1,sNx |
|
./pkg/mdsio/mdsio_gl.F |
1179 |
MDSREADFIELD_2D_GL |
II=((ip-1)*nSx+(bi-1))*sNx+I |
|
./pkg/mdsio/mdsio_gl.F |
1383 |
MDSWRITEFIELD_2D_GL |
DO I=1,sNx |
|
./pkg/mdsio/mdsio_gl.F |
1384 |
MDSWRITEFIELD_2D_GL |
II=((ip-1)*nSx+(bi-1))*sNx+I |
|
./pkg/mdsio/mdsio_gl.F |
1397 |
MDSWRITEFIELD_2D_GL |
DO I=1,Nx |
|
./pkg/mdsio/mdsio_gl.F |
1398 |
MDSWRITEFIELD_2D_GL |
xy_buffer_r4(I,J) = global(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
1407 |
MDSWRITEFIELD_2D_GL |
DO I=1,Nx |
|
./pkg/mdsio/mdsio_gl.F |
1408 |
MDSWRITEFIELD_2D_GL |
xy_buffer_r8(I,J) = global(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
364 |
MDSREADFIELD_3D_GL |
DO I=1,Nx |
|
./pkg/mdsio/mdsio_gl.F |
365 |
MDSREADFIELD_3D_GL |
global(I,J) = xy_buffer_r4(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
374 |
MDSREADFIELD_3D_GL |
DO I=1,Nx |
|
./pkg/mdsio/mdsio_gl.F |
375 |
MDSREADFIELD_3D_GL |
global(I,J) = xy_buffer_r8(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
391 |
MDSREADFIELD_3D_GL |
DO I=1,sNx |
|
./pkg/mdsio/mdsio_gl.F |
392 |
MDSREADFIELD_3D_GL |
II=((ip-1)*nSx+(bi-1))*sNx+I |
|
./pkg/mdsio/mdsio_gl.F |
594 |
MDSWRITEFIELD_3D_GL |
DO I=1,sNx |
|
./pkg/mdsio/mdsio_gl.F |
595 |
MDSWRITEFIELD_3D_GL |
II=((ip-1)*nSx+(bi-1))*sNx+I |
|
./pkg/mdsio/mdsio_gl.F |
608 |
MDSWRITEFIELD_3D_GL |
DO I=1,Nx |
|
./pkg/mdsio/mdsio_gl.F |
609 |
MDSWRITEFIELD_3D_GL |
xy_buffer_r4(I,J) = global(I,J) |
|
./pkg/mdsio/mdsio_gl.F |
618 |
MDSWRITEFIELD_3D_GL |
DO I=1,Nx |
|
./pkg/mdsio/mdsio_gl.F |
619 |
MDSWRITEFIELD_3D_GL |
xy_buffer_r8(I,J) = global(I,J) |
|
./pkg/mom_common/mom_calc_absvort3.F |
46 |
MOM_CALC_ABSVORT3 |
omega3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
113 |
MOM_CALC_RELVORT3 |
I=1 |
|
./pkg/mom_common/mom_calc_relvort3.F |
117 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
118 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
119 |
MOM_CALC_RELVORT3 |
& ( vFld(I,J)*dyC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
120 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
121 |
MOM_CALC_RELVORT3 |
& + uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
131 |
MOM_CALC_RELVORT3 |
vort3(I-1,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
132 |
MOM_CALC_RELVORT3 |
& recip_rAz(I-1,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
133 |
MOM_CALC_RELVORT3 |
& vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
134 |
MOM_CALC_RELVORT3 |
& -vFld(I-2,J)*dyC(I-2,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
135 |
MOM_CALC_RELVORT3 |
& -uFld(I-1,J)*dxC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
136 |
MOM_CALC_RELVORT3 |
& +vFld(I+0,J-1)*dyC(I+0,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
140 |
MOM_CALC_RELVORT3 |
vort3(I,J-1)=vort3(I-1,J) |
|
./pkg/mom_common/mom_calc_relvort3.F |
154 |
MOM_CALC_RELVORT3 |
I=sNx+1 |
|
./pkg/mom_common/mom_calc_relvort3.F |
159 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
160 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
161 |
MOM_CALC_RELVORT3 |
& (-uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
162 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J)*dyC(I-1,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
163 |
MOM_CALC_RELVORT3 |
& + uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
166 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
167 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
168 |
MOM_CALC_RELVORT3 |
& (-vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
169 |
MOM_CALC_RELVORT3 |
& +uFld(I,J-1)*dxC(I,J-1,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
170 |
MOM_CALC_RELVORT3 |
& - uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
173 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
174 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
175 |
MOM_CALC_RELVORT3 |
& (+uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
176 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
177 |
MOM_CALC_RELVORT3 |
& - vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
188 |
MOM_CALC_RELVORT3 |
vort3(I+1,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
189 |
MOM_CALC_RELVORT3 |
& recip_rAz(I+1,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
190 |
MOM_CALC_RELVORT3 |
& vFld(I+1,J)*dyC(I+1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
191 |
MOM_CALC_RELVORT3 |
& -vFld(I-0,J)*dyC(I-0,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
192 |
MOM_CALC_RELVORT3 |
& -uFld(I+1,J)*dxC(I+1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
193 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J-1)*dyC(I-1,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
197 |
MOM_CALC_RELVORT3 |
vort3(I,J-1)=vort3(I+1,J) |
|
./pkg/mom_common/mom_calc_relvort3.F |
211 |
MOM_CALC_RELVORT3 |
I=1 |
|
./pkg/mom_common/mom_calc_relvort3.F |
216 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
217 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
218 |
MOM_CALC_RELVORT3 |
& (+uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
219 |
MOM_CALC_RELVORT3 |
& +vFld(I,J)*dyC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
220 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
223 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
224 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
225 |
MOM_CALC_RELVORT3 |
& (-uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
226 |
MOM_CALC_RELVORT3 |
& +uFld(I,J-1)*dxC(I,J-1,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
227 |
MOM_CALC_RELVORT3 |
& + vFld(I,J)*dyC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
230 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
231 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
232 |
MOM_CALC_RELVORT3 |
& (+vFld(I,J)*dyC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
233 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
234 |
MOM_CALC_RELVORT3 |
& + uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
245 |
MOM_CALC_RELVORT3 |
vort3(I-1,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
246 |
MOM_CALC_RELVORT3 |
& recip_rAz(I-1,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
247 |
MOM_CALC_RELVORT3 |
& vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
248 |
MOM_CALC_RELVORT3 |
& -vFld(I-2,J)*dyC(I-2,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
249 |
MOM_CALC_RELVORT3 |
& +vFld(I-0,J+1)*dyC(I-0,J+1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
250 |
MOM_CALC_RELVORT3 |
& +uFld(I-1,J-1)*dxC(I-1,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
254 |
MOM_CALC_RELVORT3 |
vort3(I,J+1)=vort3(I-1,J) |
|
./pkg/mom_common/mom_calc_relvort3.F |
268 |
MOM_CALC_RELVORT3 |
I=sNx+1 |
|
./pkg/mom_common/mom_calc_relvort3.F |
273 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
274 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
275 |
MOM_CALC_RELVORT3 |
& (-uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
276 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J)*dyC(I-1,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
277 |
MOM_CALC_RELVORT3 |
& + uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
280 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
281 |
MOM_CALC_RELVORT3 |
& +recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
282 |
MOM_CALC_RELVORT3 |
& (+uFld(I,J-1)*dxC(I,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
283 |
MOM_CALC_RELVORT3 |
& -uFld(I,J)*dxC(I,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
284 |
MOM_CALC_RELVORT3 |
& - vFld(I-1,J)*dyC(I-1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
295 |
MOM_CALC_RELVORT3 |
vort3(I+1,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
296 |
MOM_CALC_RELVORT3 |
& recip_rAz(I+1,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
297 |
MOM_CALC_RELVORT3 |
& vFld(I+1,J)*dyC(I+1,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
298 |
MOM_CALC_RELVORT3 |
& -vFld(I-0,J)*dyC(I-0,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
299 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J+1)*dyC(I-1,J+1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
300 |
MOM_CALC_RELVORT3 |
& +uFld(I+1,J-1)*dxC(I+1,J-1,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
304 |
MOM_CALC_RELVORT3 |
vort3(I,J+1)=vort3(I+1,J) |
|
./pkg/mom_common/mom_calc_relvort3.F |
47 |
MOM_CALC_RELVORT3 |
DO I=1-OLx,sNx+OLx |
|
./pkg/mom_common/mom_calc_relvort3.F |
48 |
MOM_CALC_RELVORT3 |
vort3(I,J) = 0. _d 0 |
|
./pkg/mom_common/mom_calc_relvort3.F |
54 |
MOM_CALC_RELVORT3 |
DO I=2-OLx,sNx+OLx |
|
./pkg/mom_common/mom_calc_relvort3.F |
57 |
MOM_CALC_RELVORT3 |
vort3(I,J)= |
|
./pkg/mom_common/mom_calc_relvort3.F |
58 |
MOM_CALC_RELVORT3 |
& recip_rAz(I,J,bi,bj)*( |
|
./pkg/mom_common/mom_calc_relvort3.F |
59 |
MOM_CALC_RELVORT3 |
& ( vFld(I,J)*dyC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
60 |
MOM_CALC_RELVORT3 |
& -vFld(I-1,J)*dyC(I-1,J,bi,bj) ) |
|
./pkg/mom_common/mom_calc_relvort3.F |
61 |
MOM_CALC_RELVORT3 |
& -( uFld(I,J)*dxC(I,J,bi,bj) |
|
./pkg/mom_common/mom_calc_relvort3.F |
62 |
MOM_CALC_RELVORT3 |
& -uFld(I,J-1)*dxC(I,J-1,bi,bj) ) |
|
./pkg/mom_fluxform/mom_u_del2u.F |
46 |
MOM_U_DEL2U |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_u_xviscflux.F |
50 |
MOM_U_XVISCFLUX |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_u_yviscflux.F |
51 |
MOM_U_YVISCFLUX |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_v_del2v.F |
45 |
MOM_V_DEL2V |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_v_xviscflux.F |
51 |
MOM_V_XVISCFLUX |
INTEGER I,J |
|
./pkg/mom_fluxform/mom_v_yviscflux.F |
50 |
MOM_V_YVISCFLUX |
INTEGER I,J |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
100 |
MOM_VI_CORIOLIS |
& )*uBarXY*recip_dyC(i,j,bi,bj)*_maskS(I,J,K,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
34 |
MOM_VI_CORIOLIS |
INTEGER I,J |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
42 |
MOM_VI_CORIOLIS |
DO I=2-Olx,sNx+Olx |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
52 |
MOM_VI_CORIOLIS |
& )*vBarXY*recip_dxC(i,j,bi,bj)*_maskW(I,J,K,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
58 |
MOM_VI_CORIOLIS |
DO I=2-Olx,sNx+Olx |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
67 |
MOM_VI_CORIOLIS |
& )*vBarXY*recip_dxC(i,j,bi,bj)*_maskW(I,J,K,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
75 |
MOM_VI_CORIOLIS |
DO I=1-Olx,sNx+Olx-1 |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
85 |
MOM_VI_CORIOLIS |
& )*uBarXY*recip_dyC(i,j,bi,bj)*_maskS(I,J,K,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_coriolis.F |
91 |
MOM_VI_CORIOLIS |
DO I=1-Olx,sNx+Olx-1 |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
103 |
MOM_VI_U_VERTSHEAR |
uZp=(mask_Kp1*uFld(I,J,Kp1,bi,bj)-uFld(I,J,K,bi,bj))*rkSign |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
123 |
MOM_VI_U_VERTSHEAR |
uShearTerm(I,J)=-0.5* |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
130 |
MOM_VI_U_VERTSHEAR |
uShearTerm(I,J)=-0.5*(wBarXp*uZp+wBarXm*uZm) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
31 |
MOM_VI_U_VERTSHEAR |
INTEGER I,J,Kp1,Km1 |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
57 |
MOM_VI_U_VERTSHEAR |
DO I=2-OLx,sNx+OLx |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
69 |
MOM_VI_U_VERTSHEAR |
& wFld(I,J,K,bi,bj)*rA(i,j,bi,bj)*maskC(I,J,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
70 |
MOM_VI_U_VERTSHEAR |
& +wFld(I-1,J,K,bi,bj)*rA(i-1,j,bi,bj)*maskC(I-1,J,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
76 |
MOM_VI_U_VERTSHEAR |
& wFld(I,J,Kp1,bi,bj)*rA(i,j,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
77 |
MOM_VI_U_VERTSHEAR |
& +wFld(I-1,J,Kp1,bi,bj)*rA(i-1,j,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
83 |
MOM_VI_U_VERTSHEAR |
& wFld(I,J,K,bi,bj)*maskC(I,J,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
84 |
MOM_VI_U_VERTSHEAR |
& +wFld(I-1,J,K,bi,bj)*maskC(I-1,J,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
89 |
MOM_VI_U_VERTSHEAR |
& wFld(I,J,Kp1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
90 |
MOM_VI_U_VERTSHEAR |
& +wFld(I-1,J,Kp1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_u_vertshear.F |
96 |
MOM_VI_U_VERTSHEAR |
uZm=(uFld(I,J,K,bi,bj)-mask_Km1*uFld(I,J,Km1,bi,bj))*rkSign |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
103 |
MOM_VI_V_VERTSHEAR |
vZp=(mask_Kp1*vFld(I,J,Kp1,bi,bj)-vFld(I,J,K,bi,bj))*rkSign |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
123 |
MOM_VI_V_VERTSHEAR |
vShearTerm(I,J)=-0.5* |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
130 |
MOM_VI_V_VERTSHEAR |
vShearTerm(I,J)=-0.5*(wBarYp*vZp+wBarYm*vZm) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
31 |
MOM_VI_V_VERTSHEAR |
INTEGER I,J,Kp1,Km1 |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
57 |
MOM_VI_V_VERTSHEAR |
DO I=1-OLx,sNx+OLx |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
69 |
MOM_VI_V_VERTSHEAR |
& wFld(I,J,K,bi,bj)*rA(i,j,bi,bj)*maskC(i,j,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
70 |
MOM_VI_V_VERTSHEAR |
& +wFld(I,J-1,K,bi,bj)*rA(i,j-1,bi,bj)*maskC(i,j-1,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
76 |
MOM_VI_V_VERTSHEAR |
& wFld(I,J,Kp1,bi,bj)*rA(i,j,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
77 |
MOM_VI_V_VERTSHEAR |
& +wFld(I,J-1,Kp1,bi,bj)*rA(i,j-1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
83 |
MOM_VI_V_VERTSHEAR |
& wFld(I,J,K,bi,bj)*maskC(i,j,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
84 |
MOM_VI_V_VERTSHEAR |
& +wFld(I,J-1,K,bi,bj)*maskC(i,j-1,Km1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
89 |
MOM_VI_V_VERTSHEAR |
& wFld(I,J,Kp1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
90 |
MOM_VI_V_VERTSHEAR |
& +wFld(I,J-1,Kp1,bi,bj) |
|
./pkg/mom_vecinv/mom_vi_v_vertshear.F |
96 |
MOM_VI_V_VERTSHEAR |
vZm=(vFld(I,J,K,bi,bj)-mask_Km1*vFld(I,J,Km1,bi,bj))*rkSign |
|
./pkg/monitor/mon_advcfl.F |
33 |
MON_ADVCFL |
INTEGER bi,bj,I,J,K |
|
./pkg/monitor/mon_advcfl.F |
42 |
MON_ADVCFL |
DO I=1,sNx |
|
./pkg/monitor/mon_advcfl.F |
43 |
MON_ADVCFL |
tmpVal=abs(U( |