Changeset 2296 for trunk/LMDZ.COMMON/libf/dyn3d
- Timestamp:
- Apr 24, 2020, 6:55:59 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/dyn3d/vlsplt.F
r1508 r2296 456 456 DO ij=iip2,ip1jm 457 457 ! On a besoin de q et masse seulement entre iip2 et ip1jm 458 masseq(ij,l,iq2)=masse(ij,l,iq)*q(ij,l,iq) 459 Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 458 !masseq(ij,l,iq2)=masse(ij,l,iq)*q(ij,l,iq) 459 !Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 460 !MVals: veiller a ce qu'on n'ait pas de denominateur nul 461 masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),1e-16) 462 if (q(ij,l,iq).gt.1e-16) then 463 Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 464 else 465 Ratio(ij,l,iq2)=0. 466 endif 460 467 enddo 461 468 enddo … … 473 480 DO l=1,llm 474 481 DO ij=iip2+1,ip1jm 475 new_m=masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l) 482 !MVals: veiller a ce qu'on n'ait pas de denominateur nul 483 new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),1e-16) 476 484 q(ij,l,iq)=(q(ij,l,iq)*masse(ij,l,iq)+ 477 485 & u_mq(ij-1,l)-u_mq(ij,l)) … … 777 785 ! attention, chaque fils doit avoir son masseq, sinon, le 1er 778 786 ! fils ecrase le masseq de ses freres. 779 masseq(ij,l,iq2)=masse(ij,l,iq)*q(ij,l,iq) 780 Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 787 !masseq(ij,l,iq2)=masse(ij,l,iq)*q(ij,l,iq) 788 !Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 789 !MVals: veiller a ce qu'on n'ait pas de denominateur nul 790 masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),1e-16) 791 if (q(ij,l,iq).gt.1e-16) then 792 Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 793 else 794 Ratio(ij,l,iq2)=0. 795 endif 781 796 enddo 782 797 enddo … … 997 1012 DO l=1,llm 998 1013 DO ij=1,ip1jmp1 999 masseq(ij,l,iq2)=masse(ij,l,iq)*q(ij,l,iq) 1000 Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 1014 !masseq(ij,l,iq2)=masse(ij,l,iq)*q(ij,l,iq) 1015 !Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 1016 !MVals: veiller a ce qu'on n'ait pas de denominateur nul 1017 masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),1e-16) 1018 if (q(ij,l,iq).gt.1e-16) then 1019 Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq) 1020 else 1021 Ratio(ij,l,iq2)=0. 1022 endif 1001 1023 enddo 1002 1024 enddo
Note: See TracChangeset
for help on using the changeset viewer.