604 |
enddo |
enddo |
605 |
|
|
606 |
if(n0.gt.1) then |
if(n0.gt.1) then |
607 |
IF (prt_level > 9) WRITE(lunout,*) |
IF (prt_level > 9) print *, |
608 |
& 'Nombre de points pour lesquels on advect plus que le' |
& 'Nombre de points pour lesquels on advect plus que le' |
609 |
& ,'contenu de la maille : ',n0 |
& ,'contenu de la maille : ',n0 |
610 |
|
|
648 |
c goto 8888 |
c goto 8888 |
649 |
zz=0.5*(zsig-zsigd(ijq,l))/zsigg(ijq,l) |
zz=0.5*(zsig-zsigd(ijq,l))/zsigg(ijq,l) |
650 |
if(.not.(zz.gt.0..and.zz.le.0.5)) then |
if(.not.(zz.gt.0..and.zz.le.0.5)) then |
651 |
WRITE(lunout,*)'probleme2 au point ij=',ij, |
print *,'probleme2 au point ij=',ij, |
652 |
s ' l=',l |
s ' l=',l |
653 |
WRITE(lunout,*)'zz=',zz |
print *,'zz=',zz |
654 |
stop |
stop |
655 |
endif |
endif |
656 |
u_mq(ij,l)=u_mq(ij,l)+masse(ijq,l)*( |
u_mq(ij,l)=u_mq(ij,l)+masse(ijq,l)*( |
678 |
c goto 9999 |
c goto 9999 |
679 |
zz=0.5*(zsig-zsigg(ijq,l))/zsigd(ijq,l) |
zz=0.5*(zsig-zsigg(ijq,l))/zsigd(ijq,l) |
680 |
if(.not.(zz.gt.0..and.zz.le.0.5)) then |
if(.not.(zz.gt.0..and.zz.le.0.5)) then |
681 |
WRITE(lunout,*)'probleme22 au point ij=',ij |
print *,'probleme22 au point ij=',ij |
682 |
s ,' l=',l |
s ,' l=',l |
683 |
WRITE(lunout,*)'zz=',zz |
print *,'zz=',zz |
684 |
stop |
stop |
685 |
endif |
endif |
686 |
u_mq(ij,l)=u_mq(ij,l)-masse(ijq,l)*( |
u_mq(ij,l)=u_mq(ij,l)-masse(ijq,l)*( |