10 |
use comconst |
use comconst |
11 |
use comvert |
use comvert |
12 |
use conf_gcm_m |
use conf_gcm_m |
13 |
use logic |
use conf_gcm_m |
14 |
use comgeom |
use comgeom |
15 |
use temps |
use temps |
16 |
use ener |
use ener |
|
use comdissip |
|
17 |
implicit none |
implicit none |
18 |
|
|
19 |
c--------------------------------------------------- |
c--------------------------------------------------- |
25 |
real massebx(iip1,jjp1,llm),masseby(iip1,jjm,llm) |
real massebx(iip1,jjp1,llm),masseby(iip1,jjm,llm) |
26 |
real w(iip1,jjp1,llm+1) |
real w(iip1,jjp1,llm+1) |
27 |
real fluxwppm(iim,jjp1,llm) |
real fluxwppm(iim,jjp1,llm) |
28 |
real pbaru(iip1,jjp1,llm ) |
real, intent(in):: pbaru(iip1,jjp1,llm ) |
29 |
real pbarv(iip1,jjm,llm) |
real, intent(in):: pbarv(iip1,jjm,llm) |
30 |
real unatppm(iim,jjp1,llm) |
real unatppm(iim,jjp1,llm) |
31 |
real vnatppm(iim,jjp1,llm) |
real vnatppm(iim,jjp1,llm) |
32 |
real psppm(iim,jjp1) |
real psppm(iim,jjp1) |
91 |
do j=1,jjp1 |
do j=1,jjp1 |
92 |
do i=1,iip1 |
do i=1,iip1 |
93 |
fluxw(i,j,l)=w(i,j,l)*g*0.01/aire_2d(i,j) |
fluxw(i,j,l)=w(i,j,l)*g*0.01/aire_2d(i,j) |
|
C print*,i,j,l,'fluxw(i,j,l)=',fluxw(i,j,l), |
|
|
C c 'w(i,j,l)=',w(i,j,l) |
|
94 |
enddo |
enddo |
95 |
enddo |
enddo |
96 |
enddo |
enddo |