Changeset e52967c in flexpart.git


Ignore:
Timestamp:
Jun 4, 2017, 8:57:12 PM (7 years ago)
Author:
Espen Sollum ATMOS <eso@…>
Branches:
master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
Children:
db375cc
Parents:
72d3a5a
Message:

Quick fix for segmentation fault when particle is at north pole

Location:
src
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • src/advance.f90

    r94735e2 re52967c  
    242242  dt2=real(memtime(2)-itime)
    243243  dtt=1./(dt1+dt2)
     244
     245! eso: Temporary fix for particle exactly at north pole
     246  if (jyp >= nymax) then
     247    ! write(*,*) 'WARNING: advance.f90 jyp >= nymax. xt,yt:',xt,yt
     248    jyp=jyp-1
     249  end if
    244250
    245251  ! Compute maximum mixing height around particle position
  • src/conccalc.f90

    r2bec33e re52967c  
    111111      p4=ddx*ddy
    112112
     113! eso: Temporary fix for particle exactly at north pole
     114      if (jyp >= nymax) then
     115      !  write(*,*) 'WARNING: conccalc.f90 jyp >= nymax'
     116        jyp=jyp-1
     117      end if
     118
    113119      do il=2,nz
    114120        if (height(il).gt.ztra1(i)) then
  • src/conccalc_mpi.f90

    r4c64400 re52967c  
    118118      p4=ddx*ddy
    119119
     120! eso: Temporary fix for particle exactly at north pole
     121      if (jyp >= nymax) then
     122        ! write(*,*) 'WARNING: conccalc.f90 jyp >= nymax'
     123        jyp=jyp-1
     124      end if
     125
    120126      do il=2,nz
    121127        if (height(il).gt.ztra1(i)) then
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG