5 |
use dimens_m |
use dimens_m |
6 |
use paramet_m |
use paramet_m |
7 |
use comconst |
use comconst |
8 |
use comvert |
use disvert_m |
9 |
use comgeom |
use comgeom |
10 |
|
USE nr_util, ONLY : pi |
11 |
IMPLICIT NONE |
IMPLICIT NONE |
12 |
|
|
13 |
c======================================================================= |
c======================================================================= |
28 |
c Arguments: |
c Arguments: |
29 |
c ---------- |
c ---------- |
30 |
INTEGER iq,nt |
INTEGER iq,nt |
31 |
REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ) |
REAL, intent(in):: pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ) |
32 |
REAL masse(iip1,jjp1,llm) |
REAL masse(iip1,jjp1,llm) |
33 |
REAL q( iip1,jjp1,llm,0:9) |
REAL q( iip1,jjp1,llm,0:9) |
34 |
REAL w( ip1jmp1,llm ) |
REAL w( ip1jmp1,llm ) |
249 |
q( i,1,llm+1-l,0)=qpn |
q( i,1,llm+1-l,0)=qpn |
250 |
q( i,jjp1,llm+1-l,0)=qps |
q( i,jjp1,llm+1-l,0)=qps |
251 |
enddo |
enddo |
|
c enddo |
|
|
c print*,'qpn',qpn,'qps',qps |
|
|
c print*,'dqzpn',dqzpn,'dqzps',dqzps |
|
|
c enddo |
|
252 |
dyn1=0. |
dyn1=0. |
253 |
dys1=0. |
dys1=0. |
254 |
dyn2=0. |
dyn2=0. |
333 |
PRINT*,'SY(',i,j,l,')=',q(i,j,l,2), |
PRINT*,'SY(',i,j,l,')=',q(i,j,l,2), |
334 |
$ q(i,j-1,l,2) |
$ q(i,j-1,l,2) |
335 |
PRINT*,'SZ(',i,j,l,')=',q(i,j,l,3) |
PRINT*,'SZ(',i,j,l,')=',q(i,j,l,3) |
|
c PRINT*,' PBL EN SORTIE D'' ADVZP' |
|
336 |
q(i,j,l,0)=0. |
q(i,j,l,0)=0. |
|
c STOP |
|
337 |
ENDIF |
ENDIF |
338 |
ENDDO |
ENDDO |
339 |
ENDDO |
ENDDO |