source: LMDZ.3.3/branches/rel-LF/libf/dyn3d/defrun.F @ 4122

Last change on this file since 4122 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 5.3 KB
Line 
1      SUBROUTINE defrun(tapedef,creinit)
2      IMPLICIT NONE
3c-----------------------------------------------------------------------
4c   declarations:
5c   -------------
6
7#include "dimensions.h"
8#include "paramet.h"
9#include "control.h"
10#include "logic.h"
11#include "serre.h"
12#include "comdissip.h"
13      LOGICAL creinit
14
15c   local:
16c   ------
17
18      CHARACTER ch1*72,ch2*72,ch3*72,ch4*8
19      INTEGER tapedef,tapeout
20      REAL clonn,clatt,alphaxx,alphayy
21c
22c-----------------------------------------------------------------------
23c   initialisations:
24c   ----------------
25
26      tapeout=6
27
28c-----------------------------------------------------------------------
29c  Parametres de controle du run:
30c-----------------------------------------------------------------------
31
32      READ (tapedef,9000) ch1,ch2,ch3
33      WRITE(tapeout,9000) ch1,ch2,ch3
34
35      READ (tapedef,9001) ch1,ch4
36      READ (tapedef,*)    nday
37      WRITE(tapeout,9001) ch1,'nday'
38      WRITE(tapeout,*)    nday
39
40      READ (tapedef,9001) ch1,ch4
41      READ (tapedef,*)    day_step
42      WRITE(tapeout,9001) ch1,'day_step'
43      WRITE(tapeout,*)    day_step
44
45      READ (tapedef,9001) ch1,ch4
46      READ (tapedef,*)    iperiod
47      WRITE(tapeout,9001) ch1,'iperiod'
48      WRITE(tapeout,*)    iperiod
49
50      READ (tapedef,9001) ch1,ch4
51      READ (tapedef,*)    iconser
52      WRITE(tapeout,9001) ch1,'iconser'
53      WRITE(tapeout,*)    iconser
54
55      READ (tapedef,9001) ch1,ch4
56      READ (tapedef,*)    iecri
57      WRITE(tapeout,9001) ch1,'iecri'
58      WRITE(tapeout,*)    iecri
59
60      READ (tapedef,9001) ch1,ch4
61      READ (tapedef,*)    periodav
62      WRITE(tapeout,9001) ch1,'periodav'
63      WRITE(tapeout,*)    periodav
64
65      READ (tapedef,9001) ch1,ch4
66      READ (tapedef,*)    idissip
67      WRITE(tapeout,9001) ch1,'idissip'
68      WRITE(tapeout,*)    idissip
69
70ccc  ....   P. Le Van , ajout  le 5/09/95 .pour la dissipation  ...
71ccc
72      READ (tapedef,9001) ch1,ch4
73      READ (tapedef,*)    lstardis
74      WRITE(tapeout,9001) ch1,'lstardis'
75      WRITE(tapeout,*)    lstardis
76
77      READ (tapedef,9001) ch1,ch4
78      READ (tapedef,*)    niterdis
79      WRITE(tapeout,9001) ch1,'niterdis'
80      WRITE(tapeout,*)    niterdis
81
82      READ (tapedef,9001) ch1,ch4
83      READ (tapedef,*)    coefdis
84      WRITE(tapeout,9001) ch1,'coefdis'
85      WRITE(tapeout,*)    coefdis
86
87      READ (tapedef,9001) ch1,ch4
88      READ (tapedef,*)    tetavel
89      WRITE(tapeout,9001) ch1,'tetavel'
90      WRITE(tapeout,*)    tetavel
91
92      READ (tapedef,9001) ch1,ch4
93      READ (tapedef,*)    tetatemp
94      WRITE(tapeout,9001) ch1,'tetatemp'
95      WRITE(tapeout,*)    tetatemp
96
97      READ (tapedef,9001) ch1,ch4
98      READ (tapedef,*)    iphysiq
99      WRITE(tapeout,9001) ch1,'iphysiq'
100      WRITE(tapeout,*)    iphysiq
101
102      READ (tapedef,9001) ch1,ch4
103      READ (tapedef,*)    purmats
104      WRITE(tapeout,9001) ch1,'purmats'
105      WRITE(tapeout,*)    purmats
106
107      READ (tapedef,9001) ch1,ch4
108      READ (tapedef,*)    physic
109      WRITE(tapeout,9001) ch1,'physic'
110      WRITE(tapeout,*)    physic
111
112ccc   .... P.Le Van, ajout le 03/01/96 pour l'ecriture phys ...
113c
114      READ (tapedef,9001) ch1,ch4
115      READ (tapedef,*)    ecritphy
116      WRITE(tapeout,9001) ch1,'ecritphy'
117      WRITE(tapeout,*)    ecritphy
118
119ccc  ....   P. Le Van , ajout  le 7/03/95 .pour le zoom ...
120
121      IF( creinit ) go to 100
122
123      READ (tapedef,9001) ch1,ch4
124      READ (tapedef,*)    clonn
125      WRITE(tapeout,9001) ch1,'clon'
126      WRITE(tapeout,*)    clonn
127      IF( clon - clonn .gt.0.001)  THEN
128        PRINT *,' La valeur de clon passee par run.def est differente de
129     *  celle lue sur le fichier  start '
130        STOP
131      ENDIF
132
133      READ (tapedef,9001) ch1,ch4
134      READ (tapedef,*)    clatt
135      WRITE(tapeout,9001) ch1,'clat'
136      WRITE(tapeout,*)    clatt
137
138      IF( clat -clatt.gt.0.001 )  THEN
139        PRINT *,' La valeur de clat passee par run.def est differente de
140     *  celle lue sur le fichier  start '
141        STOP
142      ENDIF
143
144      READ (tapedef,9001) ch1,ch4
145      READ (tapedef,*)    alphaxx
146      WRITE(tapeout,9001) ch1,'alphax'
147      WRITE(tapeout,*)    alphaxx
148
149      IF( alphax-alphaxx.gt.0.001 )  THEN
150        PRINT *,' La valeur de alphax passee par run.def est differente 
151     *  de celle lue sur le fichier  start '
152        STOP
153      ENDIF
154
155      READ (tapedef,9001) ch1,ch4
156      READ (tapedef,*)    alphayy
157      WRITE(tapeout,9001) ch1,'alphay'
158      WRITE(tapeout,*)    alphayy
159
160      IF( alphay-alphayy.gt.0.001 )  THEN
161        PRINT *,' La valeur de alphay passee par run.def est differente 
162     * de celle lue sur le fichier  start '
163        STOP
164      ENDIF
165
166      RETURN
167c
168c   
169100   CONTINUE
170      READ (tapedef,9001) ch1,ch4
171      READ (tapedef,*)    clon
172      WRITE(tapeout,9001) ch1,'clon'
173      WRITE(tapeout,*)    clon
174
175      READ (tapedef,9001) ch1,ch4
176      READ (tapedef,*)    clat
177      WRITE(tapeout,9001) ch1,'clat'
178      WRITE(tapeout,*)    clat
179
180      READ (tapedef,9001) ch1,ch4
181      READ (tapedef,*)    alphax
182      WRITE(tapeout,9001) ch1,'alphax'
183      WRITE(tapeout,*)    alphax
184
185      READ (tapedef,9001) ch1,ch4
186      READ (tapedef,*)    alphay
187      WRITE(tapeout,9001) ch1,'alphay'
188      WRITE(tapeout,*)    alphay
189
190c-----------------------------------------------------------------------
191
192      RETURN
1939000  FORMAT(3(/,a72))
1949001  FORMAT(/,a72,/,a8)
195      END
Note: See TracBrowser for help on using the repository browser.