97 |
do L =1,lm |
do L =1,lm |
98 |
do j =1,jm |
do j =1,jm |
99 |
do i =1,im |
do i =1,im |
100 |
cldtot(i,j,L)=min(1.0,max(cldlw(i,j,L),fccave(i,j,L)/imstturb)) |
cldtot(i,j,L)=min(1.0 _d 0,max(cldlw(i,j,L),fccave(i,j,L)/imstturb)) |
101 |
cldmxo(i,j,L) = min( 1.0 , clwmo(i,j,L) ) |
cldmxo(i,j,L) = min( 1.0 _d 0, clwmo(i,j,L) ) |
102 |
lwlz(i,j,L) = lwlz(i,j,L) + qliqave(i,j,L)/imstturb |
lwlz(i,j,L) = lwlz(i,j,L) + qliqave(i,j,L)/imstturb |
103 |
enddo |
enddo |
104 |
enddo |
enddo |
107 |
do L =1,lm |
do L =1,lm |
108 |
do j =1,jm |
do j =1,jm |
109 |
do i =1,im |
do i =1,im |
110 |
cldtot(i,j,L) = min( 1.0,cldlw(i,j,L) ) |
cldtot(i,j,L) = min( 1.0 _d 0,cldlw(i,j,L) ) |
111 |
cldmxo(i,j,L) = min( 1.0,clwmo(i,j,L) ) |
cldmxo(i,j,L) = min( 1.0 _d 0,clwmo(i,j,L) ) |
112 |
enddo |
enddo |
113 |
enddo |
enddo |
114 |
enddo |
enddo |
151 |
IF(NLWCLD.NE.0)THEN |
IF(NLWCLD.NE.0)THEN |
152 |
DO L = 1,lm |
DO L = 1,lm |
153 |
DO I = 1,ISTRIP |
DO I = 1,ISTRIP |
154 |
CLRO(I,L) = min( 1.0,clro(i,L) ) |
CLRO(I,L) = min( 1.0 _d 0,clro(i,L) ) |
155 |
CLMO(I,L) = min( 1.0,clmo(i,L) ) |
CLMO(I,L) = min( 1.0 _d 0,clmo(i,L) ) |
156 |
ENDDO |
ENDDO |
157 |
ENDDO |
ENDDO |
158 |
ENDIF |
ENDIF |
1096 |
|
|
1097 |
if (taux.gt.0.02 .and. fcld(i,j,k).gt.0.01) then |
if (taux.gt.0.02 .and. fcld(i,j,k).gt.0.01) then |
1098 |
|
|
1099 |
reff1=min(reff(i,j,k,1),130.) |
reff1=min(reff(i,j,k,1),130. _d 0) |
1100 |
reff2=min(reff(i,j,k,2),20.0) |
reff2=min(reff(i,j,k,2),20.0 _d 0) |
1101 |
|
|
1102 |
w1=taucl(i,j,k,1)*(aiw(1,ib)+(aiw(2,ib)+(aiw(3,ib) |
w1=taucl(i,j,k,1)*(aiw(1,ib)+(aiw(2,ib)+(aiw(3,ib) |
1103 |
* +aiw(4,ib)*reff1)*reff1)*reff1) |
* +aiw(4,ib)*reff1)*reff1)*reff1) |