snap lowest layer depth to earth radius if close

4 years ago · 46b1b21bbd
parent 6448ca5708
commit 46b1b21bbd
1 changed files with 5 additions and 0 deletions
--- a/src/qpgetinp.f
+++ b/src/qpgetinp.f
@ -619,6 +619,11 @@ c
        endif
      enddo
 c
+      if(abs(dp0(l)-rratmos)/rratmos.lt.1.0d-6)then
+        write(*,*) ' Snapping lowest layer depth to earth model radius.'
+        dp0(l)=rratmos
+      endif
+
      if(dp0(l).gt.rratmos)then
        stop ' Error in qpgetinp: earth radius larger than pre-defined!'
      else if(dp0(l).lt.rratmos)then