Quick fix for segmentation fault when particle is at north pole

src/advance.f90

@@ -242,6 +242,12 @@ subroutine advance(itime,nrelpoint,ldt,up,vp,wp, &
   dt2=real(memtime(2)-itime)
   dtt=1./(dt1+dt2)
 
+! eso: Temporary fix for particle exactly at north pole
+  if (jyp >= nymax) then
+    ! write(*,*) 'WARNING: advance.f90 jyp >= nymax. xt,yt:',xt,yt
+    jyp=jyp-1
+  end if
+
   ! Compute maximum mixing height around particle position
   !*******************************************************
 

src/conccalc.f90

@@ -110,6 +110,12 @@ subroutine conccalc(itime,weight)
       p3=rddx*ddy
       p4=ddx*ddy
 
+! eso: Temporary fix for particle exactly at north pole
+      if (jyp >= nymax) then
+      !  write(*,*) 'WARNING: conccalc.f90 jyp >= nymax'
+        jyp=jyp-1
+      end if
+
       do il=2,nz
         if (height(il).gt.ztra1(i)) then
           indz=il-1

src/conccalc_mpi.f90

@@ -117,6 +117,12 @@ subroutine conccalc(itime,weight)
       p3=rddx*ddy
       p4=ddx*ddy
 
+! eso: Temporary fix for particle exactly at north pole
+      if (jyp >= nymax) then
+        ! write(*,*) 'WARNING: conccalc.f90 jyp >= nymax'
+        jyp=jyp-1
+      end if
+
       do il=2,nz
         if (height(il).gt.ztra1(i)) then
           indz=il-1