125 |
c ---------------------------------------------- |
c ---------------------------------------------- |
126 |
do j=1,jm |
do j=1,jm |
127 |
do i=1,im |
do i=1,im |
128 |
phis_std(i,j) = min( 400.0, sqrt( max(0.0,phis_var(i,j)) )/grav ) |
phis_std(i,j) = min( 400.0 _d 0, sqrt( max(0.0 _d 0,phis_var(i,j)) )/grav ) |
129 |
enddo |
enddo |
130 |
enddo |
enddo |
131 |
|
|
183 |
do L = 1,Lm |
do L = 1,Lm |
184 |
do j = 1,jm |
do j = 1,jm |
185 |
do i = 1,im |
do i = 1,im |
186 |
dragu(i,j,L) = sign( min(0.006,abs(dragu(i,j,L))),dragu(i,j,L) ) |
dragu(i,j,L) = sign( min(0.006 _d 0,abs(dragu(i,j,L))),dragu(i,j,L) ) |
187 |
dragv(i,j,L) = sign( min(0.006,abs(dragv(i,j,L))),dragv(i,j,L) ) |
dragv(i,j,L) = sign( min(0.006 _d 0,abs(dragv(i,j,L))),dragv(i,j,L) ) |
188 |
dragt(i,j,L) = -( uz(i,j,L)*dragu(i,j,L)+vz(i,j,L)*dragv(i,j,L) ) |
dragt(i,j,L) = -( uz(i,j,L)*dragu(i,j,L)+vz(i,j,L)*dragv(i,j,L) ) |
189 |
. *cpinv |
. *cpinv |
190 |
dudt(i,j,L) = dudt(i,j,L) + dragu(i,j,L) |
dudt(i,j,L) = dudt(i,j,L) + dragu(i,j,L) |
524 |
c ---------------------------------- |
c ---------------------------------- |
525 |
crifro = 1.0 - 0.25/richsn |
crifro = 1.0 - 0.25/richsn |
526 |
crif2 = crifro*crifro |
crif2 = crifro*crifro |
527 |
if( l.eq.2 ) crif2 = min(0.7,crif2) |
if( l.eq.2 ) crif2 = min(0.7 _d 0,crif2) |
528 |
|
|
529 |
if( fro2.gt.crif2 ) then |
if( fro2.gt.crif2 ) then |
530 |
tensio(i,L) = crif2/fro2*tensio(i,L-1) |
tensio(i,L) = crif2/fro2*tensio(i,L-1) |