The TELEMAC-MASCARET system  trunk
cvsp_check_steady_gaia.f
Go to the documentation of this file.
1 ! ****************************
2  SUBROUTINE cvsp_check_steady_gaia
3 ! ****************************
4 !
5  &(j)
6 !
7 !***********************************************************************
8 ! GAIA V8P1 16/05/2017
9 !***********************************************************************
10 !
12 !
16 !
21 !
25 !
26 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
28 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
29 !
31 
32  IMPLICIT NONE
33 !
34 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
35 !
36  INTEGER, INTENT(IN) :: J
37 !
38 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
39 !
40  INTEGER K, JG
41  DOUBLE PRECISION AT
42 !
43 !-----------------------------------------------------------------------
44 !
45  jg = j
46  IF(ncsize > 1) jg = mesh%KNOLG%I(j)
47 !
48  at = dt*lt/percou
49 !
50 !-----------------------------------------------------------------------
51 !
52  DO k=1,pro_max(j)-1
53  IF((pro_d(j,k+1,1) - pro_d(j,k,1)).LT.0.d0) THEN
54  WRITE(lu,*) 'ERR: UNSTEADY VSP! ,J,K,AT',
55  & jg, k, at, pro_d(j,k+1,1), pro_d(j,k,1),pro_max(j)
56  CALL cvsp_p_gaia('./','UNSTEADY_',jg)
57  CALL layers_p_gaia('./','UNSTEADY_', jg)
58  CALL plante(1)
59  stop
60  ENDIF
61  ENDDO
62 !
63 !-----------------------------------------------------------------------
64 !
65  RETURN
66  END SUBROUTINE
67 
subroutine layers_p_gaia(PATH_PRE, FILE_PRE, JG)
Definition: layers_p_gaia.f:7
double precision, target dt
Time step It may be different from the one in TELEMAC because of the morphological factor...
subroutine cvsp_p_gaia(PATH_PRE, FILE_PRE, JG)
Definition: cvsp_p_gaia.f:7
subroutine cvsp_check_steady_gaia(J)
integer, target lt
Numero du pas de temps.
double precision, dimension(:,:,:), allocatable, target pro_d
Vertical sorting profile: depth for each layer, class, point.
integer percou
COUPLING PERIOD USED IN CVSM TO CALCULATE THE TIME, SHOULD COME FROM TELEMAC ONE DAY.
type(bief_mesh), target mesh
Mesh structure.
integer, dimension(:), allocatable pro_max
Maximum layer number in a vertical sorting profile for each point.