 Timestamp:
 Oct 21, 2014, 7:52:35 AM (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/src/advance.f90
r4 r29 463 463 dcwsave=dcwsave+vp*dt 464 464 zt=zt+w*dt*real(ldirect) 465 466 ! HSO/AL: Particle managed to go over highest level > interpolation error in goto 700 467 ! alias interpol_wind (division by zero) 468 if (zt.ge.height(nz)) zt=height(nz)100.*eps 465 469 466 470 if (zt.gt.h) then … … 699 703 endif 700 704 705 ! HSO/AL: Prevent particles from disappearing at the pole 706 !****************************************************************** 707 708 if ( yt.lt.0. ) then 709 xt=mod(xt+180.,360.) 710 yt=yt 711 else if ( yt.gt.real(nymin1) ) then 712 xt=mod(xt+180.,360.) 713 yt=2*real(nymin1)yt 714 endif 701 715 702 716 ! Check position: If trajectory outside model domain, terminate it … … 704 718 705 719 if ((xt.lt.0.).or.(xt.ge.real(nxmin1)).or.(yt.lt.0.).or. & 706 (yt.g e.real(nymin1))) then720 (yt.gt.real(nymin1))) then 707 721 nstop=3 708 722 return … … 859 873 endif 860 874 875 ! HSO/AL: Prevent particles from disappearing at the pole 876 !****************************************************************** 877 878 if ( yt.lt.0. ) then 879 xt=mod(xt+180.,360.) 880 yt=yt 881 else if ( yt.gt.real(nymin1) ) then 882 xt=mod(xt+180.,360.) 883 yt=2*real(nymin1)yt 884 endif 885 861 886 ! Check position: If trajectory outside model domain, terminate it 862 887 !***************************************************************** 863 888 864 889 if ((xt.lt.0.).or.(xt.ge.real(nxmin1)).or.(yt.lt.0.).or. & 865 (yt.g e.real(nymin1))) then890 (yt.gt.real(nymin1))) then 866 891 nstop=3 867 892 return
Note: See TracChangeset
for help on using the changeset viewer.