Context Navigation

conf_gcm.F @ 1033

Last change on this file since 1033 was 999, checked in by Laurent Fairhead, 16 years ago

Modifs sur le parallelisme: masquage dans la physique
Inclusion strato
mise en coherence etat0
le mode offline fonctionne maintenant en parallele,
les fichiers de la dynamiques sont correctement sortis et peuvent etre reconstruit avec rebuild
la version parallele de la dynamique peut s'executer sans MPI (sur 1 proc)
L'OPENMP fonctionne maintenant sans la parallelisation MPI.

YM
LF

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 21.6 KB

Line
1	!
2	! $Header$
3	!
4	c
5	c
6	SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
7	c
8	use IOIPSL
9	IMPLICIT NONE
10	c-----------------------------------------------------------------------
11	c Auteurs : L. Fairhead , P. Le Van .
12	c
13	c Arguments :
14	c
15	c tapedef :
16	c etatinit : = TRUE , on ne compare pas les valeurs des para-
17	c -metres du zoom avec celles lues sur le fichier start .
18	c clesphy0 : sortie .
19	c
20	LOGICAL etatinit
21	INTEGER tapedef
22
23	INTEGER longcles
24	PARAMETER( longcles = 20 )
25	REAL clesphy0( longcles )
26	c
27	c Declarations :
28	c --------------
29	#include "dimensions.h"
30	#include "paramet.h"
31	#include "control.h"
32	#include "logic.h"
33	#include "serre.h"
34	#include "comdissnew.h"
35
36	! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
37	! #include "clesphys.h"
38	#include "iniprint.h"
39	c
40	c
41	c local:
42	c ------
43
44	CHARACTER ch172,ch272,ch372,ch412
45	REAL clonn,clatt,grossismxx,grossismyy
46	REAL dzoomxx,dzoomyy, tauxx,tauyy
47	LOGICAL fxyhypbb, ysinuss
48	INTEGER i
49
50	c
51	c -------------------------------------------------------------------
52	c
53	c ......... Version du 29/04/97 ..........
54	c
55	c Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot,
56	c tetatemp ajoutes pour la dissipation .
57	c
58	c Autre parametre ajoute en fin de liste de tapedef : fxyhypb
59	c
60	c Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb.
61	c Sinon , choix de fxynew , a derivee sinusoidale ..
62	c
63	c ...... etatinit = . TRUE. si defrun est appele dans ETAT0_LMD ou
64	c LIMIT_LMD pour l'initialisation de start.dat (dic) et
65	c de limit.dat ( dic) ...........
66	c Sinon etatinit = . FALSE .
67	c
68	c Donc etatinit = .F. si on veut comparer les valeurs de grossismx ,
69	c grossismy,clon,clat, fxyhypb lues sur le fichier start avec
70	c celles passees par run.def , au debut du gcm, apres l'appel a
71	c lectba .
72	c Ces parmetres definissant entre autres la grille et doivent etre
73	c pareils et coherents , sinon il y aura divergence du gcm .
74	c
75	c-----------------------------------------------------------------------
76	c initialisations:
77	c ----------------
78
79	!Config Key = lunout
80	!Config Desc = unite de fichier pour les impressions
81	!Config Def = 6
82	!Config Help = unite de fichier pour les impressions
83	!Config (defaut sortie standard = 6)
84	lunout=6
85	CALL getin('lunout', lunout)
86	IF (lunout /= 5 .and. lunout /= 6) THEN
87	OPEN(lunout,FILE='lmdz.out')
88	ENDIF
89
90	!Config Key = prt_level
91	!Config Desc = niveau d'impressions de débogage
92	!Config Def = 0
93	!Config Help = Niveau d'impression pour le débogage
94	!Config (0 = minimum d'impression)
95	prt_level = 0
96	CALL getin('prt_level',prt_level)
97
98	c-----------------------------------------------------------------------
99	c Parametres de controle du run:
100	c-----------------------------------------------------------------------
101
102	!Config Key = dayref
103	!Config Desc = Jour de l'etat initial
104	!Config Def = 1
105	!Config Help = Jour de l'etat initial ( = 350 si 20 Decembre ,
106	!Config par expl. ,comme ici ) ... A completer
107	dayref=1
108	CALL getin('dayref', dayref)
109
110	!Config Key = anneeref
111	!Config Desc = Annee de l'etat initial
112	!Config Def = 1998
113	!Config Help = Annee de l'etat initial
114	!Config ( avec 4 chiffres ) ... A completer
115	anneeref = 1998
116	CALL getin('anneeref',anneeref)
117
118	!Config Key = raz_date
119	!Config Desc = Remise a zero de la date initiale
120	!Config Def = 0 (pas de remise a zero)
121	!Config Help = Remise a zero de la date initiale
122	!Config 0 pas de remise a zero, on garde la date du fichier restart
123	!Config 1 prise en compte de la date de gcm.def avec remise a zero
124	!Config des compteurs de pas de temps
125	raz_date = 0
126	CALL getin('raz_date', raz_date)
127
128	!Config Key = nday
129	!Config Desc = Nombre de jours d'integration
130	!Config Def = 10
131	!Config Help = Nombre de jours d'integration
132	!Config ... On pourait aussi permettre des mois ou des annees !
133	nday = 10
134	CALL getin('nday',nday)
135
136	!Config Key = day_step
137	!Config Desc = nombre de pas par jour
138	!Config Def = 240
139	!Config Help = nombre de pas par jour (multiple de iperiod) (
140	!Config ici pour dt = 1 min )
141	day_step = 240
142	CALL getin('day_step',day_step)
143
144	!Config Key = iperiod
145	!Config Desc = periode pour le pas Matsuno
146	!Config Def = 5
147	!Config Help = periode pour le pas Matsuno (en pas de temps)
148	iperiod = 5
149	CALL getin('iperiod',iperiod)
150
151	!Config Key = iapp_tracvl
152	!Config Desc = frequence du groupement des flux
153	!Config Def = iperiod
154	!Config Help = frequence du groupement des flux (en pas de temps)
155	iapp_tracvl = iperiod
156	CALL getin('iapp_tracvl',iapp_tracvl)
157
158	!Config Key = iconser
159	!Config Desc = periode de sortie des variables de controle
160	!Config Def = 240
161	!Config Help = periode de sortie des variables de controle
162	!Config (En pas de temps)
163	iconser = 240
164	CALL getin('iconser', iconser)
165
166	!Config Key = iecri
167	!Config Desc = periode d'ecriture du fichier histoire
168	!Config Def = 1
169	!Config Help = periode d'ecriture du fichier histoire (en jour)
170	iecri = 1
171	CALL getin('iecri',iecri)
172
173
174	!Config Key = periodav
175	!Config Desc = periode de stockage fichier histmoy
176	!Config Def = 1
177	!Config Help = periode de stockage fichier histmoy (en jour)
178	periodav = 1.
179	CALL getin('periodav',periodav)
180
181	!Config Key = idissip
182	!Config Desc = periode de la dissipation
183	!Config Def = 10
184	!Config Help = periode de la dissipation
185	!Config (en pas) ... a completer !
186	idissip = 10
187	CALL getin('idissip',idissip)
188
189	ccc .... P. Le Van , modif le 29/04/97 .pour la dissipation ...
190	ccc
191
192	!Config Key = lstardis
193	!Config Desc = choix de l'operateur de dissipation
194	!Config Def = y
195	!Config Help = choix de l'operateur de dissipation
196	!Config 'y' si on veut star et 'n' si on veut non-start !
197	!Config Moi y en a pas comprendre !
198	lstardis = .TRUE.
199	CALL getin('lstardis',lstardis)
200
201
202	!Config Key = nitergdiv
203	!Config Desc = Nombre d'iteration de gradiv
204	!Config Def = 1
205	!Config Help = nombre d'iterations de l'operateur de dissipation
206	!Config gradiv
207	nitergdiv = 1
208	CALL getin('nitergdiv',nitergdiv)
209
210	!Config Key = nitergrot
211	!Config Desc = nombre d'iterations de nxgradrot
212	!Config Def = 2
213	!Config Help = nombre d'iterations de l'operateur de dissipation
214	!Config nxgradrot
215	nitergrot = 2
216	CALL getin('nitergrot',nitergrot)
217
218
219	!Config Key = niterh
220	!Config Desc = nombre d'iterations de divgrad
221	!Config Def = 2
222	!Config Help = nombre d'iterations de l'operateur de dissipation
223	!Config divgrad
224	niterh = 2
225	CALL getin('niterh',niterh)
226
227
228	!Config Key = tetagdiv
229	!Config Desc = temps de dissipation pour div
230	!Config Def = 7200
231	!Config Help = temps de dissipation des plus petites longeur
232	!Config d'ondes pour u,v (gradiv)
233	tetagdiv = 7200.
234	CALL getin('tetagdiv',tetagdiv)
235
236	!Config Key = tetagrot
237	!Config Desc = temps de dissipation pour grad
238	!Config Def = 7200
239	!Config Help = temps de dissipation des plus petites longeur
240	!Config d'ondes pour u,v (nxgradrot)
241	tetagrot = 7200.
242	CALL getin('tetagrot',tetagrot)
243
244	!Config Key = tetatemp
245	!Config Desc = temps de dissipation pour h
246	!Config Def = 7200
247	!Config Help = temps de dissipation des plus petites longeur
248	!Config d'ondes pour h (divgrad)
249	tetatemp = 7200.
250	CALL getin('tetatemp',tetatemp )
251
252	!Config Key = coefdis
253	!Config Desc = coefficient pour gamdissip
254	!Config Def = 0
255	!Config Help = coefficient pour gamdissip
256	coefdis = 0.
257	CALL getin('coefdis',coefdis)
258
259	!Config Key = purmats
260	!Config Desc = Schema d'integration
261	!Config Def = n
262	!Config Help = Choix du schema d'integration temporel.
263	!Config y = pure Matsuno sinon c'est du Matsuno-leapfrog
264	purmats = .FALSE.
265	CALL getin('purmats',purmats)
266
267	!Config Key = ok_guide
268	!Config Desc = Guidage
269	!Config Def = n
270	!Config Help = Guidage
271	ok_guide = .FALSE.
272	CALL getin('ok_guide',ok_guide)
273
274	c ...............................................................
275
276	!Config Key = iflag_phys
277	!Config Desc = Avec ls physique
278	!Config Def = 1
279	!Config Help = Permet de faire tourner le modele sans
280	!Config physique.
281	iflag_phys = 1
282	CALL getin('iflag_phys',iflag_phys)
283
284
285	!Config Key = iphysiq
286	!Config Desc = Periode de la physique
287	!Config Def = 5
288	!Config Help = Periode de la physique en pas de temps de la dynamique.
289	iphysiq = 5
290	CALL getin('iphysiq', iphysiq)
291
292	!Config Key = ip_ebil_dyn
293	!Config Desc = PRINT level for energy conserv. diag.
294	!Config Def = 0
295	!Config Help = PRINT level for energy conservation diag. ;
296	! les options suivantes existent :
297	!Config 0 pas de print
298	!Config 1 pas de print
299	!Config 2 print,
300	ip_ebil_dyn = 0
301	CALL getin('ip_ebil_dyn',ip_ebil_dyn)
302	!
303
304	DO i = 1, longcles
305	clesphy0(i) = 0.
306	ENDDO
307
308	ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ...
309	c ......... ( modif le 17/04/96 ) .........
310	c
311	IF( etatinit ) GO TO 100
312
313	!Config Key = clon
314	!Config Desc = centre du zoom, longitude
315	!Config Def = 0
316	!Config Help = longitude en degres du centre
317	!Config du zoom
318	clonn = 0.
319	CALL getin('clon',clonn)
320
321	!Config Key = clat
322	!Config Desc = centre du zoom, latitude
323	!Config Def = 0
324	!Config Help = latitude en degres du centre du zoom
325	!Config
326	clatt = 0.
327	CALL getin('clat',clatt)
328
329	c
330	c
331	IF( ABS(clat - clatt).GE. 0.001 ) THEN
332	PRINT *,' La valeur de clat passee par run.def est differente de
333	* celle lue sur le fichier start '
334	STOP
335	ENDIF
336
337	!Config Key = grossismx
338	!Config Desc = zoom en longitude
339	!Config Def = 1.0
340	!Config Help = facteur de grossissement du zoom,
341	!Config selon la longitude
342	grossismxx = 1.0
343	CALL getin('grossismx',grossismxx)
344
345
346	IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN
347	PRINT *,' La valeur de grossismx passee par run.def est differente
348	* de celle lue sur le fichier start '
349	STOP
350	ENDIF
351
352	!Config Key = grossismy
353	!Config Desc = zoom en latitude
354	!Config Def = 1.0
355	!Config Help = facteur de grossissement du zoom,
356	!Config selon la latitude
357	grossismyy = 1.0
358	CALL getin('grossismy',grossismyy)
359
360	IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN
361	PRINT *,' La valeur de grossismy passee par run.def est differen
362	* te de celle lue sur le fichier start '
363	STOP
364	ENDIF
365
366	IF( grossismx.LT.1. ) THEN
367	PRINT ,' ATTENTION !! grossismx < 1 . * '
368	STOP
369	ELSE
370	alphax = 1. - 1./ grossismx
371	ENDIF
372
373
374	IF( grossismy.LT.1. ) THEN
375	PRINT ,' ATTENTION !! grossismy < 1 . * '
376	STOP
377	ELSE
378	alphay = 1. - 1./ grossismy
379	ENDIF
380
381	PRINT *,' alphax alphay defrun ',alphax,alphay
382	c
383	c alphax et alphay sont les anciennes formulat. des grossissements
384	c
385	c
386
387	!Config Key = fxyhypb
388	!Config Desc = Fonction hyperbolique
389	!Config Def = y
390	!Config Help = Fonction f(y) hyperbolique si = .true.
391	!Config sinon sinusoidale
392	fxyhypbb = .TRUE.
393	CALL getin('fxyhypb',fxyhypbb)
394
395	IF( .NOT.fxyhypb ) THEN
396	IF( fxyhypbb ) THEN
397	PRINT ,' ***** PBS DANS DEFRUN ****** '
398	PRINT ,' ** fxyhypb lu sur le fichier start est F ',
399	* 'alors qu il est T sur run.def ***'
400	STOP
401	ENDIF
402	ELSE
403	IF( .NOT.fxyhypbb ) THEN
404	PRINT ,' ***** PBS DANS DEFRUN ****** '
405	PRINT ,' ** fxyhypb lu sur le fichier start est T ',
406	* 'alors qu il est F sur run.def **** '
407	STOP
408	ENDIF
409	ENDIF
410	c
411	!Config Key = dzoomx
412	!Config Desc = extension en longitude
413	!Config Def = 0
414	!Config Help = extension en longitude de la zone du zoom
415	!Config ( fraction de la zone totale)
416	dzoomxx = 0.0
417	CALL getin('dzoomx',dzoomxx)
418
419	IF( fxyhypb ) THEN
420	IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN
421	PRINT *,' La valeur de dzoomx passee par run.def est differente
422	* de celle lue sur le fichier start '
423	STOP
424	ENDIF
425	ENDIF
426
427	!Config Key = dzoomy
428	!Config Desc = extension en latitude
429	!Config Def = 0
430	!Config Help = extension en latitude de la zone du zoom
431	!Config ( fraction de la zone totale)
432	dzoomyy = 0.0
433	CALL getin('dzoomy',dzoomyy)
434
435	IF( fxyhypb ) THEN
436	IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN
437	PRINT *,' La valeur de dzoomy passee par run.def est differente
438	* de celle lue sur le fichier start '
439	STOP
440	ENDIF
441	ENDIF
442
443	!Config Key = taux
444	!Config Desc = raideur du zoom en X
445	!Config Def = 3
446	!Config Help = raideur du zoom en X
447	tauxx = 3.0
448	CALL getin('taux',tauxx)
449
450	IF( fxyhypb ) THEN
451	IF( ABS(taux - tauxx).GE. 0.001 ) THEN
452	PRINT *,' La valeur de taux passee par run.def est differente
453	* de celle lue sur le fichier start '
454	STOP
455	ENDIF
456	ENDIF
457
458	!Config Key = tauyy
459	!Config Desc = raideur du zoom en Y
460	!Config Def = 3
461	!Config Help = raideur du zoom en Y
462	tauyy = 3.0
463	CALL getin('tauy',tauyy)
464
465	IF( fxyhypb ) THEN
466	IF( ABS(tauy - tauyy).GE. 0.001 ) THEN
467	PRINT *,' La valeur de tauy passee par run.def est differente
468	* de celle lue sur le fichier start '
469	STOP
470	ENDIF
471	ENDIF
472
473	cc
474	IF( .NOT.fxyhypb ) THEN
475
476	!Config Key = ysinus
477	!Config IF = !fxyhypb
478	!Config Desc = Fonction en Sinus
479	!Config Def = y
480	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
481	!Config sinon y = latit.
482	ysinuss = .TRUE.
483	CALL getin('ysinus',ysinuss)
484
485	IF( .NOT.ysinus ) THEN
486	IF( ysinuss ) THEN
487	PRINT ,' ***** PBS DANS DEFRUN ****** '
488	PRINT ,' ** ysinus lu sur le fichier start est F ',
489	* 'alors qu il est T sur run.def ***'
490	STOP
491	ENDIF
492	ELSE
493	IF( .NOT.ysinuss ) THEN
494	PRINT ,' ***** PBS DANS DEFRUN ****** '
495	PRINT ,' ** ysinus lu sur le fichier start est T ',
496	* 'alors qu il est F sur run.def **** '
497	STOP
498	ENDIF
499	ENDIF
500	ENDIF
501	c
502	!Config Key = offline
503	!Config Desc = Nouvelle eau liquide
504	!Config Def = n
505	!Config Help = Permet de mettre en route la
506	!Config nouvelle parametrisation de l'eau liquide !
507	offline = .FALSE.
508	CALL getin('offline',offline)
509
510	!Config Key = config_inca
511	!Config Desc = Choix de configuration de INCA
512	!Config Def = none
513	!Config Help = Choix de configuration de INCA :
514	!Config 'none' = sans INCA
515	!Config 'chem' = INCA avec calcul de chemie
516	!Config 'aero' = INCA avec calcul des aerosols
517	config_inca = 'none'
518	CALL getin('config_inca',config_inca)
519
520
521	write(lunout,*)' #########################################'
522	write(lunout,*)' Configuration des parametres du gcm: '
523	write(lunout,*)' dayref = ', dayref
524	write(lunout,*)' anneeref = ', anneeref
525	write(lunout,*)' nday = ', nday
526	write(lunout,*)' day_step = ', day_step
527	write(lunout,*)' iperiod = ', iperiod
528	write(lunout,*)' iconser = ', iconser
529	write(lunout,*)' iecri = ', iecri
530	write(lunout,*)' periodav = ', periodav
531	write(lunout,*)' idissip = ', idissip
532	write(lunout,*)' lstardis = ', lstardis
533	write(lunout,*)' nitergdiv = ', nitergdiv
534	write(lunout,*)' nitergrot = ', nitergrot
535	write(lunout,*)' niterh = ', niterh
536	write(lunout,*)' tetagdiv = ', tetagdiv
537	write(lunout,*)' tetagrot = ', tetagrot
538	write(lunout,*)' tetatemp = ', tetatemp
539	write(lunout,*)' coefdis = ', coefdis
540	write(lunout,*)' purmats = ', purmats
541	write(lunout,*)' iflag_phys = ', iflag_phys
542	write(lunout,*)' iphysiq = ', iphysiq
543	write(lunout,*)' clonn = ', clonn
544	write(lunout,*)' clatt = ', clatt
545	write(lunout,*)' grossismx = ', grossismx
546	write(lunout,*)' grossismy = ', grossismy
547	write(lunout,*)' fxyhypbb = ', fxyhypbb
548	write(lunout,*)' dzoomxx = ', dzoomxx
549	write(lunout,*)' dzoomy = ', dzoomyy
550	write(lunout,*)' tauxx = ', tauxx
551	write(lunout,*)' tauyy = ', tauyy
552	write(lunout,*)' offline = ', offline
553	write(lunout,*)' config_inca = ', config_inca
554
555	RETURN
556	c ...............................................
557	c
558	100 CONTINUE
559	!Config Key = clon
560	!Config Desc = centre du zoom, longitude
561	!Config Def = 0
562	!Config Help = longitude en degres du centre
563	!Config du zoom
564	clon = 0.
565	CALL getin('clon',clon)
566
567	!Config Key = clat
568	!Config Desc = centre du zoom, latitude
569	!Config Def = 0
570	!Config Help = latitude en degres du centre du zoom
571	!Config
572	clat = 0.
573	CALL getin('clat',clat)
574
575	!Config Key = grossismx
576	!Config Desc = zoom en longitude
577	!Config Def = 1.0
578	!Config Help = facteur de grossissement du zoom,
579	!Config selon la longitude
580	grossismx = 1.0
581	CALL getin('grossismx',grossismx)
582
583	!Config Key = grossismy
584	!Config Desc = zoom en latitude
585	!Config Def = 1.0
586	!Config Help = facteur de grossissement du zoom,
587	!Config selon la latitude
588	grossismy = 1.0
589	CALL getin('grossismy',grossismy)
590
591	IF( grossismx.LT.1. ) THEN
592	PRINT ,' ATTENTION !! grossismx < 1 . * '
593	STOP
594	ELSE
595	alphax = 1. - 1./ grossismx
596	ENDIF
597
598
599	IF( grossismy.LT.1. ) THEN
600	PRINT ,' ATTENTION !! grossismy < 1 . * '
601	STOP
602	ELSE
603	alphay = 1. - 1./ grossismy
604	ENDIF
605
606	PRINT *,' alphax alphay defrun ',alphax,alphay
607	c
608	c alphax et alphay sont les anciennes formulat. des grossissements
609	c
610	c
611
612	!Config Key = fxyhypb
613	!Config Desc = Fonction hyperbolique
614	!Config Def = y
615	!Config Help = Fonction f(y) hyperbolique si = .true.
616	!Config sinon sinusoidale
617	fxyhypb = .TRUE.
618	CALL getin('fxyhypb',fxyhypb)
619
620	!Config Key = dzoomx
621	!Config Desc = extension en longitude
622	!Config Def = 0
623	!Config Help = extension en longitude de la zone du zoom
624	!Config ( fraction de la zone totale)
625	dzoomx = 0.0
626	CALL getin('dzoomx',dzoomx)
627
628	!Config Key = dzoomy
629	!Config Desc = extension en latitude
630	!Config Def = 0
631	!Config Help = extension en latitude de la zone du zoom
632	!Config ( fraction de la zone totale)
633	dzoomy = 0.0
634	CALL getin('dzoomy',dzoomy)
635
636	!Config Key = taux
637	!Config Desc = raideur du zoom en X
638	!Config Def = 3
639	!Config Help = raideur du zoom en X
640	taux = 3.0
641	CALL getin('taux',taux)
642
643	!Config Key = tauy
644	!Config Desc = raideur du zoom en Y
645	!Config Def = 3
646	!Config Help = raideur du zoom en Y
647	tauy = 3.0
648	CALL getin('tauy',tauy)
649
650	!Config Key = ysinus
651	!Config IF = !fxyhypb
652	!Config Desc = Fonction en Sinus
653	!Config Def = y
654	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
655	!Config sinon y = latit.
656	ysinus = .TRUE.
657	CALL getin('ysinus',ysinus)
658	c
659	!Config Key = offline
660	!Config Desc = Nouvelle eau liquide
661	!Config Def = n
662	!Config Help = Permet de mettre en route la
663	!Config nouvelle parametrisation de l'eau liquide !
664	offline = .FALSE.
665	CALL getin('offline',offline)
666
667	!Config Key = config_inca
668	!Config Desc = Choix de configuration de INCA
669	!Config Def = none
670	!Config Help = Choix de configuration de INCA :
671	!Config 'none' = sans INCA
672	!Config 'chem' = INCA avec calcul de chemie
673	!Config 'aero' = INCA avec calcul des aerosols
674	config_inca = 'none'
675	CALL getin('config_inca',config_inca)
676
677	!Config key = ok_strato
678	!Config Desc = activation de la version strato
679	!Config Def = .FALSE.
680	!Config Help = active la version stratosphérique de LMDZ de F. Lott
681
682	ok_strato=.FALSE.
683	CALL getin('ok_strato',ok_strato)
684
685	write(lunout,*)' #########################################'
686	write(lunout,*)' Configuration des parametres du gcm: '
687	write(lunout,*)' dayref = ', dayref
688	write(lunout,*)' anneeref = ', anneeref
689	write(lunout,*)' nday = ', nday
690	write(lunout,*)' day_step = ', day_step
691	write(lunout,*)' iperiod = ', iperiod
692	write(lunout,*)' iconser = ', iconser
693	write(lunout,*)' iecri = ', iecri
694	write(lunout,*)' periodav = ', periodav
695	write(lunout,*)' idissip = ', idissip
696	write(lunout,*)' lstardis = ', lstardis
697	write(lunout,*)' nitergdiv = ', nitergdiv
698	write(lunout,*)' nitergrot = ', nitergrot
699	write(lunout,*)' niterh = ', niterh
700	write(lunout,*)' tetagdiv = ', tetagdiv
701	write(lunout,*)' tetagrot = ', tetagrot
702	write(lunout,*)' tetatemp = ', tetatemp
703	write(lunout,*)' coefdis = ', coefdis
704	write(lunout,*)' purmats = ', purmats
705	write(lunout,*)' iflag_phys = ', iflag_phys
706	write(lunout,*)' iphysiq = ', iphysiq
707	write(lunout,*)' clonn = ', clonn
708	write(lunout,*)' clatt = ', clatt
709	write(lunout,*)' grossismx = ', grossismx
710	write(lunout,*)' grossismy = ', grossismy
711	write(lunout,*)' fxyhypbb = ', fxyhypbb
712	write(lunout,*)' dzoomx = ', dzoomx
713	write(lunout,*)' dzoomy = ', dzoomy
714	write(lunout,*)' taux = ', taux
715	write(lunout,*)' tauy = ', tauy
716	write(lunout,*)' offline = ', offline
717	write(lunout,*)' config_inca = ', config_inca
718	write(lunout,*)' ok_strato = ', ok_strato
719	c
720	RETURN
721	END

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format