Philippe 23/02/2016: lfi2cdf: modif for JPHEXT/=1
[MNH-git_open_source-lfs.git] / LIBTOOLS / tools / diachro / src / DIAPRO / axelogpres.f90
1 !     ######spl
2       SUBROUTINE AXELOGPRES(PHMIN,PHMAX)
3 !     ##################################
4 !
5 !!****  *AXELOGPRES* - 
6 !!****    
7 !!
8 !!    PURPOSE
9 !!    -------
10 !
11 !!**  METHOD
12 !!    ------
13 !!
14 !!    EXTERNAL
15 !!    --------
16 !!
17 !!    IMPLICIT ARGUMENTS
18 !!    ------------------
19 !!
20 !!    REFERENCE
21 !!    ---------
22 !!
23 !!    AUTHOR
24 !!    ------
25 !!      J. Duron    * Laboratoire d'Aerologie *
26 !!
27 !!    MODIFICATIONS
28 !!    -------------
29 !!      Original       20/10/2000
30 !!      Updated   PM   
31 !-------------------------------------------------------------------------------
32 !
33 !*       0.    DECLARATIONS
34 !              ------------
35 !
36 !USE MODD_RESOLVCAR
37 USE MODD_PVT
38 !
39 IMPLICIT NONE
40 !
41 !*       0.1  Dummy arguments and results
42 !
43 REAL :: PHMIN,PHMAX
44 !
45 !*       0.2  Local variables
46 !
47 INTEGER             :: J, JA, ID
48 !
49 REAL :: ZVL,ZVR,ZVB,ZVT,ZWL,ZWR,ZWB,ZWT
50 CHARACTER(LEN=5) :: YCAR
51 !
52 !-------------------------------------------------------------------------------
53 !
54 CALL GETSET(ZVL,ZVR,ZVB,ZVT,ZWL,ZWR,ZWB,ZWT,ID)
55 IF(LPRESY)THEN
56 IF(XPMAX /= 0. .AND. XPMIN /= 0. .AND. XPINT /= 0.)THEN
57   IF(XPMIN < 1300.)THEN
58     XPMAX=XPMAX*100.
59     XPMIN=XPMIN*100.
60     XPINT=XPINT*100.
61   ENDIF
62 DO J=INT(XPMIN),INT(XPMAX),-INT((ABS(XPINT)))
63   IF(FLOAT(J) >= ANINT(ZWT) .AND. FLOAT(J) <= ANINT(ZWB))THEN
64   YCAR=' '
65   IF(XPINT > 1000.)THEN
66     WRITE(YCAR,'(F5.0)')FLOAT(J)/100.
67   ELSE
68     WRITE(YCAR,'(F5.0)')FLOAT(J)
69   ENDIF
70   YCAR=ADJUSTR(YCAR)
71   CALL PLCHHQ(ZWL-ZWL/110.,FLOAT(J),YCAR,13.,0.,1.)
72   CALL FRSTPT(ZWL,FLOAT(J))
73   CALL VECTOR(ZWL+(ZWR-ZWL)/(ZVR-ZVL)*.015,FLOAT(J))
74   ENDIF
75 ENDDO
76 ELSE
77   IF(PHMIN < 1300)THEN
78     PHMIN=PHMIN*100
79     PHMAX=PHMAX*100
80   ENDIF
81 DO J=INT(PHMIN),INT(PHMAX),-10000
82   IF(FLOAT(J) >= ANINT(ZWT) .AND. FLOAT(J) <= ANINT(ZWB))THEN
83   YCAR=' '
84   IF(PHMAX > 1300.)THEN
85     WRITE(YCAR,'(F5.0)')FLOAT(J)/100.
86   ELSE
87     WRITE(YCAR,'(F5.0)')FLOAT(J)
88   ENDIF
89   YCAR=ADJUSTR(YCAR)
90   print *,' **axelogpres PHMIN,PHMAX ',PHMIN,PHMAX
91   print *,' **axelogpres ZWL-ZWL/20.,FLOAT(J),YCAR ',ZWL-ZWL/20.,FLOAT(J),YCAR 
92   CALL PLCHHQ(ZWL-ZWL/100.,FLOAT(J),YCAR,13.,0.,1.)
93   CALL FRSTPT(ZWL,FLOAT(J))
94   CALL VECTOR(ZWL+(ZWR-ZWL)/(ZVR-ZVL)*.015,FLOAT(J))
95   ENDIF
96 ENDDO
97 ENDIF
98 ELSE
99 ENDIF
100 !*        2.     EXIT
101 !                ----
102 !
103 RETURN
104 END SUBROUTINE AXELOGPRES