/[lmdze]/trunk/dyn3d/advn.f
ViewVC logotype

Diff of /trunk/dyn3d/advn.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 12 by guez, Mon Jul 21 16:05:07 2008 UTC
# Line 604  c   tris des regions a traiter Line 604  c   tris des regions a traiter
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    
# Line 648  c               u_mq(ij,l)=u_mq(ij,l)+zu Line 648  c               u_mq(ij,l)=u_mq(ij,l)+zu
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)*(
# Line 678  c               u_mq(ij,l)=u_mq(ij,l)+zu Line 678  c               u_mq(ij,l)=u_mq(ij,l)+zu
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)*(
Legend:
Removed from v.3  
changed lines
  Added in v.12
  ViewVC Help
Powered by ViewVC 1.1.21