cvsp_check_steady.f

Go to the documentation of this file.

!                     ****************************
                      SUBROUTINE cvsp_check_steady
!                     ****************************
!
     &(j)
!
!***********************************************************************
! SISYPHE   V7P2                                  16/05/2017
!***********************************************************************
!
!brief   CHECKS VERTICAL SORTING PROFILE TO BE STEADY IN PRO_D
!
!history UWE MERKEL
!+        2012
!+       V6P2
!+
!
!history  P. A. TASSI (EDF R&D, LNHE)
!+        12/03/2013
!+        V6P3
!+   Cleaning, cosmetic
!
!history UWE MERKEL, R KOPMANN (BAW)
!+       18/06/2016, 2017
!+       V6P3, V7P2
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| J              |<--| INDEX OF A POINT IN MESH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_sisyphe

      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: J
      INTEGER  K, JG
      DOUBLE PRECISION AT
!
!-----------------------------------------------------------------------
!
      jg = j
      IF(ncsize > 1) jg = mesh%KNOLG%I(j)
!
      at = dt*lt/percou
!
!-----------------------------------------------------------------------
!
      DO k=1,pro_max(j)-1
        IF((pro_d(j,k+1,1) - pro_d(j,k,1)).LT.0.d0) THEN
          WRITE(lu,*) 'ERR: UNSTEADY VSP! ,J,K,AT',
     &           jg, k, at, pro_d(j,k+1,1), pro_d(j,k,1)
          CALL cvsp_p('./','UNSTEADY_',jg)
          CALL layers_p('./','UNSTEADY_', jg)
          CALL plante(1)
          stop
        ENDIF
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END SUBROUTINE