Changeset 467460a in flexpart.git for src/richardson.f90

Timestamp:

Mar 13, 2021, 10:00:37 AM (3 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

GFS_025, dev

Children:

9ca6e38

Parents:

759df5f

Message:

First commit of files from Hachinger

File:

• src/richardson.f90 (modified) (1 diff, 1 prop)

src/richardson.f90

@@ -141 +141 @@
     zold=z
   end do
-  k=k-1 ! ESO: make sure k <= nuvz (ticket #139)
+  ! k=k-1 ! ESO: make sure k <= nuvz (ticket #139)
 20   continue
 
   ! Determine Richardson number between the critical levels
   !********************************************************
+  ! JMA: It may happen that k >= nuvz:
+  ! JMA: In that case: k is set to nuvz 
+
+  if (k .gt. nuvz) k = nuvz ! JMA
 
   zl1=zold