cvsp_check_steady.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\sisyphe\cvsp_check_steady.f
00002 !
00063                       SUBROUTINE CVSP_CHECK_STEADY
00064 !                     ****************************
00065 !
00066      &(J)
00067 !
00068 !***********************************************************************
00069 ! SISYPHE   V6P3                                  12/03/2013
00070 !***********************************************************************
00071 !
00072 !
00073 !
00074 !
00075 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00076 !| J              |<--| INDEX OF A POINT IN MESH
00077 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00078 !
00079       USE DECLARATIONS_SISYPHE
00080 !
00081       IMPLICIT NONE
00082       INTEGER LNG,LU
00083       COMMON/INFO/LNG,LU
00084 !
00085 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00086 !
00087       INTEGER, INTENT(IN) :: J
00088       INTEGER  K, JG
00089       DOUBLE PRECISION AT
00090 !
00091 !-----------------------------------------------------------------------
00092 !
00093       JG = J
00094       IF(NCSIZE > 1) JG = MESH%KNOLG%I(J)
00095 !
00096       AT = DT*LT/PERCOU
00097 !
00098 !-----------------------------------------------------------------------
00099 !
00100       DO K=1,PRO_MAX(J)-1
00101         IF((PRO_D(J,K+1,1) - PRO_D(J,K,1)).LT.0.D0) THEN
00102           WRITE(LU,*) 'ERR: UNSTEADY VSP! ,J,K,AT',
00103      &           JG, K, AT, PRO_D(J,K+1,1), PRO_D(J,K,1)
00104           CALL CVSP_P('./','UNSTEADY_',JG)
00105           CALL PLANTE(1)
00106           STOP
00107         ENDIF
00108       ENDDO
00109 !
00110 !-----------------------------------------------------------------------
00111 !
00112       RETURN
00113       END SUBROUTINE

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0