Context Navigation

conf_gcm.F @ 1664

Last change on this file since 1664 was 1664, checked in by Laurent Fairhead, 12 years ago

Version testing basée sur la r1575

Testing release based on r1575

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

File size: 28.5 KB

Line
1	!
2	! $Id: conf_gcm.F 1664 2012-10-09 13:29:15Z fairhead $
3	!
4	c
5	c
6	SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
7	c
8	USE control_mod
9	#ifdef CPP_IOIPSL
10	use IOIPSL
11	#else
12	! if not using IOIPSL, we still need to use (a local version of) getin
13	use ioipsl_getincom
14	#endif
15	USE infotrac, ONLY : type_trac
16	IMPLICIT NONE
17	c-----------------------------------------------------------------------
18	c Auteurs : L. Fairhead , P. Le Van .
19	c
20	c Arguments :
21	c
22	c tapedef :
23	c etatinit : = TRUE , on ne compare pas les valeurs des para-
24	c -metres du zoom avec celles lues sur le fichier start .
25	c clesphy0 : sortie .
26	c
27	LOGICAL etatinit
28	INTEGER tapedef
29
30	INTEGER longcles
31	PARAMETER( longcles = 20 )
32	REAL clesphy0( longcles )
33	c
34	c Declarations :
35	c --------------
36	#include "dimensions.h"
37	#include "paramet.h"
38	#include "logic.h"
39	#include "serre.h"
40	#include "comdissnew.h"
41	#include "temps.h"
42	#include "comconst.h"
43
44	! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
45	! #include "clesphys.h"
46	#include "iniprint.h"
47	c
48	c
49	c local:
50	c ------
51
52	CHARACTER ch172,ch272,ch372,ch412
53	REAL clonn,clatt,grossismxx,grossismyy
54	REAL dzoomxx,dzoomyy, tauxx,tauyy
55	LOGICAL fxyhypbb, ysinuss
56	INTEGER i
57	LOGICAL use_filtre_fft
58	c
59	c -------------------------------------------------------------------
60	c
61	c ......... Version du 29/04/97 ..........
62	c
63	c Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot,
64	c tetatemp ajoutes pour la dissipation .
65	c
66	c Autre parametre ajoute en fin de liste de tapedef : fxyhypb
67	c
68	c Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb.
69	c Sinon , choix de fxynew , a derivee sinusoidale ..
70	c
71	c ...... etatinit = . TRUE. si defrun est appele dans ETAT0_LMD ou
72	c LIMIT_LMD pour l'initialisation de start.dat (dic) et
73	c de limit.dat ( dic) ...........
74	c Sinon etatinit = . FALSE .
75	c
76	c Donc etatinit = .F. si on veut comparer les valeurs de grossismx ,
77	c grossismy,clon,clat, fxyhypb lues sur le fichier start avec
78	c celles passees par run.def , au debut du gcm, apres l'appel a
79	c lectba .
80	c Ces parmetres definissant entre autres la grille et doivent etre
81	c pareils et coherents , sinon il y aura divergence du gcm .
82	c
83	c-----------------------------------------------------------------------
84	c initialisations:
85	c ----------------
86
87	!Config Key = lunout
88	!Config Desc = unite de fichier pour les impressions
89	!Config Def = 6
90	!Config Help = unite de fichier pour les impressions
91	!Config (defaut sortie standard = 6)
92	lunout=6
93	CALL getin('lunout', lunout)
94	IF (lunout /= 5 .and. lunout /= 6) THEN
95	OPEN(lunout,FILE='lmdz.out')
96	ENDIF
97
98	!Config Key = prt_level
99	!Config Desc = niveau d'impressions de débogage
100	!Config Def = 0
101	!Config Help = Niveau d'impression pour le débogage
102	!Config (0 = minimum d'impression)
103	prt_level = 0
104	CALL getin('prt_level',prt_level)
105
106	c-----------------------------------------------------------------------
107	c Parametres de controle du run:
108	c-----------------------------------------------------------------------
109	!Config Key = planet_type
110	!Config Desc = planet type ("earth", "mars", "venus", ...)
111	!Config Def = earth
112	!Config Help = this flag sets the type of atymosphere that is considered
113	planet_type="earth"
114	CALL getin('planet_type',planet_type)
115
116	!Config Key = calend
117	!Config Desc = type de calendrier utilise
118	!Config Def = earth_360d
119	!Config Help = valeur possible: earth_360d, earth_365d, earth_366d
120	!Config
121	calend = 'earth_360d'
122	CALL getin('calend', calend)
123
124	!Config Key = dayref
125	!Config Desc = Jour de l'etat initial
126	!Config Def = 1
127	!Config Help = Jour de l'etat initial ( = 350 si 20 Decembre ,
128	!Config par expl. ,comme ici ) ... A completer
129	dayref=1
130	CALL getin('dayref', dayref)
131
132	!Config Key = anneeref
133	!Config Desc = Annee de l'etat initial
134	!Config Def = 1998
135	!Config Help = Annee de l'etat initial
136	!Config ( avec 4 chiffres ) ... A completer
137	anneeref = 1998
138	CALL getin('anneeref',anneeref)
139
140	!Config Key = raz_date
141	!Config Desc = Remise a zero de la date initiale
142	!Config Def = 0 (pas de remise a zero)
143	!Config Help = Remise a zero de la date initiale
144	!Config 0 pas de remise a zero, on garde la date du fichier restart
145	!Config 1 prise en compte de la date de gcm.def avec remise a zero
146	!Config des compteurs de pas de temps
147	raz_date = 0
148	CALL getin('raz_date', raz_date)
149
150	!Config Key = nday
151	!Config Desc = Nombre de jours d'integration
152	!Config Def = 10
153	!Config Help = Nombre de jours d'integration
154	!Config ... On pourait aussi permettre des mois ou des annees !
155	nday = 10
156	CALL getin('nday',nday)
157
158	!Config Key = day_step
159	!Config Desc = nombre de pas par jour
160	!Config Def = 240
161	!Config Help = nombre de pas par jour (multiple de iperiod) (
162	!Config ici pour dt = 1 min )
163	day_step = 240
164	CALL getin('day_step',day_step)
165
166	!Config Key = nsplit_phys
167	!Config Desc = nombre de pas par jour
168	!Config Def = 1
169	!Config Help = nombre de pas par jour (multiple de iperiod) (
170	!Config ici pour dt = 1 min )
171	nsplit_phys = 1
172	CALL getin('nsplit_phys',nsplit_phys)
173
174	!Config Key = iperiod
175	!Config Desc = periode pour le pas Matsuno
176	!Config Def = 5
177	!Config Help = periode pour le pas Matsuno (en pas de temps)
178	iperiod = 5
179	CALL getin('iperiod',iperiod)
180
181	!Config Key = iapp_tracvl
182	!Config Desc = frequence du groupement des flux
183	!Config Def = iperiod
184	!Config Help = frequence du groupement des flux (en pas de temps)
185	iapp_tracvl = iperiod
186	CALL getin('iapp_tracvl',iapp_tracvl)
187
188	!Config Key = iconser
189	!Config Desc = periode de sortie des variables de controle
190	!Config Def = 240
191	!Config Help = periode de sortie des variables de controle
192	!Config (En pas de temps)
193	iconser = 240
194	CALL getin('iconser', iconser)
195
196	!Config Key = iecri
197	!Config Desc = periode d'ecriture du fichier histoire
198	!Config Def = 1
199	!Config Help = periode d'ecriture du fichier histoire (en jour)
200	iecri = 1
201	CALL getin('iecri',iecri)
202
203
204	!Config Key = periodav
205	!Config Desc = periode de stockage fichier histmoy
206	!Config Def = 1
207	!Config Help = periode de stockage fichier histmoy (en jour)
208	periodav = 1.
209	CALL getin('periodav',periodav)
210
211	!Config Key = output_grads_dyn
212	!Config Desc = output dynamics diagnostics in 'dyn.dat' file
213	!Config Def = n
214	!Config Help = output dynamics diagnostics in Grads-readable 'dyn.dat' file
215	output_grads_dyn=.false.
216	CALL getin('output_grads_dyn',output_grads_dyn)
217
218	!Config Key = dissip_period
219	!Config Desc = periode de la dissipation
220	!Config Def = 0
221	!Config Help = periode de la dissipation
222	!Config dissip_period=0 => la valeur sera calcule dans inidissip
223	!Config dissip_period>0 => on prend cette valeur
224	dissip_period = 0
225	CALL getin('dissip_period',dissip_period)
226
227	ccc .... P. Le Van , modif le 29/04/97 .pour la dissipation ...
228	ccc
229
230	!Config Key = lstardis
231	!Config Desc = choix de l'operateur de dissipation
232	!Config Def = y
233	!Config Help = choix de l'operateur de dissipation
234	!Config 'y' si on veut star et 'n' si on veut non-start !
235	!Config Moi y en a pas comprendre !
236	lstardis = .TRUE.
237	CALL getin('lstardis',lstardis)
238
239
240	!Config Key = nitergdiv
241	!Config Desc = Nombre d'iteration de gradiv
242	!Config Def = 1
243	!Config Help = nombre d'iterations de l'operateur de dissipation
244	!Config gradiv
245	nitergdiv = 1
246	CALL getin('nitergdiv',nitergdiv)
247
248	!Config Key = nitergrot
249	!Config Desc = nombre d'iterations de nxgradrot
250	!Config Def = 2
251	!Config Help = nombre d'iterations de l'operateur de dissipation
252	!Config nxgradrot
253	nitergrot = 2
254	CALL getin('nitergrot',nitergrot)
255
256
257	!Config Key = niterh
258	!Config Desc = nombre d'iterations de divgrad
259	!Config Def = 2
260	!Config Help = nombre d'iterations de l'operateur de dissipation
261	!Config divgrad
262	niterh = 2
263	CALL getin('niterh',niterh)
264
265
266	!Config Key = tetagdiv
267	!Config Desc = temps de dissipation pour div
268	!Config Def = 7200
269	!Config Help = temps de dissipation des plus petites longeur
270	!Config d'ondes pour u,v (gradiv)
271	tetagdiv = 7200.
272	CALL getin('tetagdiv',tetagdiv)
273
274	!Config Key = tetagrot
275	!Config Desc = temps de dissipation pour grad
276	!Config Def = 7200
277	!Config Help = temps de dissipation des plus petites longeur
278	!Config d'ondes pour u,v (nxgradrot)
279	tetagrot = 7200.
280	CALL getin('tetagrot',tetagrot)
281
282	!Config Key = tetatemp
283	!Config Desc = temps de dissipation pour h
284	!Config Def = 7200
285	!Config Help = temps de dissipation des plus petites longeur
286	!Config d'ondes pour h (divgrad)
287	tetatemp = 7200.
288	CALL getin('tetatemp',tetatemp )
289
290	! Parametres controlant la variation sur la verticale des constantes de
291	! dissipation.
292	! Pour le moment actifs uniquement dans la version a 39 niveaux
293	! avec ok_strato=y
294
295	dissip_factz=4.
296	dissip_deltaz=10.
297	dissip_zref=30.
298	CALL getin('dissip_factz',dissip_factz )
299	CALL getin('dissip_deltaz',dissip_deltaz )
300	CALL getin('dissip_zref',dissip_zref )
301
302	iflag_top_bound=1
303	tau_top_bound=1.e-5
304	CALL getin('iflag_top_bound',iflag_top_bound)
305	CALL getin('tau_top_bound',tau_top_bound)
306
307	!Config Key = coefdis
308	!Config Desc = coefficient pour gamdissip
309	!Config Def = 0
310	!Config Help = coefficient pour gamdissip
311	coefdis = 0.
312	CALL getin('coefdis',coefdis)
313
314	!Config Key = purmats
315	!Config Desc = Schema d'integration
316	!Config Def = n
317	!Config Help = Choix du schema d'integration temporel.
318	!Config y = pure Matsuno sinon c'est du Matsuno-leapfrog
319	purmats = .FALSE.
320	CALL getin('purmats',purmats)
321
322	!Config Key = ok_guide
323	!Config Desc = Guidage
324	!Config Def = n
325	!Config Help = Guidage
326	ok_guide = .FALSE.
327	CALL getin('ok_guide',ok_guide)
328
329	c ...............................................................
330
331	!Config Key = read_start
332	!Config Desc = Initialize model using a 'start.nc' file
333	!Config Def = y
334	!Config Help = y: intialize dynamical fields using a 'start.nc' file
335	! n: fields are initialized by 'iniacademic' routine
336	read_start= .true.
337	CALL getin('read_start',read_start)
338
339	!Config Key = iflag_phys
340	!Config Desc = Avec ls physique
341	!Config Def = 1
342	!Config Help = Permet de faire tourner le modele sans
343	!Config physique.
344	iflag_phys = 1
345	CALL getin('iflag_phys',iflag_phys)
346
347
348	!Config Key = iphysiq
349	!Config Desc = Periode de la physique
350	!Config Def = 5
351	!Config Help = Periode de la physique en pas de temps de la dynamique.
352	iphysiq = 5
353	CALL getin('iphysiq', iphysiq)
354
355	!Config Key = ip_ebil_dyn
356	!Config Desc = PRINT level for energy conserv. diag.
357	!Config Def = 0
358	!Config Help = PRINT level for energy conservation diag. ;
359	! les options suivantes existent :
360	!Config 0 pas de print
361	!Config 1 pas de print
362	!Config 2 print,
363	ip_ebil_dyn = 0
364	CALL getin('ip_ebil_dyn',ip_ebil_dyn)
365	!
366
367	DO i = 1, longcles
368	clesphy0(i) = 0.
369	ENDDO
370
371	ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ...
372	c ......... ( modif le 17/04/96 ) .........
373	c
374	IF( etatinit ) GO TO 100
375
376	!Config Key = clon
377	!Config Desc = centre du zoom, longitude
378	!Config Def = 0
379	!Config Help = longitude en degres du centre
380	!Config du zoom
381	clonn = 0.
382	CALL getin('clon',clonn)
383
384	!Config Key = clat
385	!Config Desc = centre du zoom, latitude
386	!Config Def = 0
387	!Config Help = latitude en degres du centre du zoom
388	!Config
389	clatt = 0.
390	CALL getin('clat',clatt)
391
392	c
393	c
394	IF( ABS(clat - clatt).GE. 0.001 ) THEN
395	write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
396	& ' est differente de celle lue sur le fichier start '
397	STOP
398	ENDIF
399
400	!Config Key = grossismx
401	!Config Desc = zoom en longitude
402	!Config Def = 1.0
403	!Config Help = facteur de grossissement du zoom,
404	!Config selon la longitude
405	grossismxx = 1.0
406	CALL getin('grossismx',grossismxx)
407
408
409	IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN
410	write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
411	& 'run.def est differente de celle lue sur le fichier start '
412	STOP
413	ENDIF
414
415	!Config Key = grossismy
416	!Config Desc = zoom en latitude
417	!Config Def = 1.0
418	!Config Help = facteur de grossissement du zoom,
419	!Config selon la latitude
420	grossismyy = 1.0
421	CALL getin('grossismy',grossismyy)
422
423	IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN
424	write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
425	& 'run.def est differente de celle lue sur le fichier start '
426	STOP
427	ENDIF
428
429	IF( grossismx.LT.1. ) THEN
430	write(lunout,*)
431	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
432	STOP
433	ELSE
434	alphax = 1. - 1./ grossismx
435	ENDIF
436
437
438	IF( grossismy.LT.1. ) THEN
439	write(lunout,*)
440	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
441	STOP
442	ELSE
443	alphay = 1. - 1./ grossismy
444	ENDIF
445
446	write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
447	c
448	c alphax et alphay sont les anciennes formulat. des grossissements
449	c
450	c
451
452	!Config Key = fxyhypb
453	!Config Desc = Fonction hyperbolique
454	!Config Def = y
455	!Config Help = Fonction f(y) hyperbolique si = .true.
456	!Config sinon sinusoidale
457	fxyhypbb = .TRUE.
458	CALL getin('fxyhypb',fxyhypbb)
459
460	IF( .NOT.fxyhypb ) THEN
461	IF( fxyhypbb ) THEN
462	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
463	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
464	* 'F alors qu il est T sur run.def ***'
465	STOP
466	ENDIF
467	ELSE
468	IF( .NOT.fxyhypbb ) THEN
469	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
470	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
471	* 'T alors qu il est F sur run.def **** '
472	STOP
473	ENDIF
474	ENDIF
475	c
476	!Config Key = dzoomx
477	!Config Desc = extension en longitude
478	!Config Def = 0
479	!Config Help = extension en longitude de la zone du zoom
480	!Config ( fraction de la zone totale)
481	dzoomxx = 0.0
482	CALL getin('dzoomx',dzoomxx)
483
484	IF( fxyhypb ) THEN
485	IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN
486	write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
487	* 'run.def est differente de celle lue sur le fichier start '
488	STOP
489	ENDIF
490	ENDIF
491
492	!Config Key = dzoomy
493	!Config Desc = extension en latitude
494	!Config Def = 0
495	!Config Help = extension en latitude de la zone du zoom
496	!Config ( fraction de la zone totale)
497	dzoomyy = 0.0
498	CALL getin('dzoomy',dzoomyy)
499
500	IF( fxyhypb ) THEN
501	IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN
502	write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
503	* 'run.def est differente de celle lue sur le fichier start '
504	STOP
505	ENDIF
506	ENDIF
507
508	!Config Key = taux
509	!Config Desc = raideur du zoom en X
510	!Config Def = 3
511	!Config Help = raideur du zoom en X
512	tauxx = 3.0
513	CALL getin('taux',tauxx)
514
515	IF( fxyhypb ) THEN
516	IF( ABS(taux - tauxx).GE. 0.001 ) THEN
517	write(lunout,*)'conf_gcm: La valeur de taux passee par ',
518	* 'run.def est differente de celle lue sur le fichier start '
519	STOP
520	ENDIF
521	ENDIF
522
523	!Config Key = tauyy
524	!Config Desc = raideur du zoom en Y
525	!Config Def = 3
526	!Config Help = raideur du zoom en Y
527	tauyy = 3.0
528	CALL getin('tauy',tauyy)
529
530	IF( fxyhypb ) THEN
531	IF( ABS(tauy - tauyy).GE. 0.001 ) THEN
532	write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
533	* 'run.def est differente de celle lue sur le fichier start '
534	STOP
535	ENDIF
536	ENDIF
537
538	cc
539	IF( .NOT.fxyhypb ) THEN
540
541	!Config Key = ysinus
542	!Config IF = !fxyhypb
543	!Config Desc = Fonction en Sinus
544	!Config Def = y
545	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
546	!Config sinon y = latit.
547	ysinuss = .TRUE.
548	CALL getin('ysinus',ysinuss)
549
550	IF( .NOT.ysinus ) THEN
551	IF( ysinuss ) THEN
552	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
553	write(lunout,)' ** ysinus lu sur le fichier start est F',
554	* ' alors qu il est T sur run.def ***'
555	STOP
556	ENDIF
557	ELSE
558	IF( .NOT.ysinuss ) THEN
559	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
560	write(lunout,)' ** ysinus lu sur le fichier start est T',
561	* ' alors qu il est F sur run.def **** '
562	STOP
563	ENDIF
564	ENDIF
565	ENDIF ! of IF( .NOT.fxyhypb )
566	c
567	!Config Key = offline
568	!Config Desc = Nouvelle eau liquide
569	!Config Def = n
570	!Config Help = Permet de mettre en route la
571	!Config nouvelle parametrisation de l'eau liquide !
572	offline = .FALSE.
573	CALL getin('offline',offline)
574
575	!Config Key = type_trac
576	!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS
577	!Config Def = lmdz
578	!Config Help =
579	!Config 'lmdz' = pas de couplage, pur LMDZ
580	!Config 'inca' = model de chime INCA
581	!Config 'repr' = model de chime REPROBUS
582	type_trac = 'lmdz'
583	CALL getin('type_trac',type_trac)
584
585	!Config Key = config_inca
586	!Config Desc = Choix de configuration de INCA
587	!Config Def = none
588	!Config Help = Choix de configuration de INCA :
589	!Config 'none' = sans INCA
590	!Config 'chem' = INCA avec calcul de chemie
591	!Config 'aero' = INCA avec calcul des aerosols
592	config_inca = 'none'
593	CALL getin('config_inca',config_inca)
594
595	!Config Key = ok_dynzon
596	!Config Desc = calcul et sortie des transports
597	!Config Def = n
598	!Config Help = Permet de mettre en route le calcul des transports
599	!Config
600	ok_dynzon = .FALSE.
601	CALL getin('ok_dynzon',ok_dynzon)
602
603	!Config Key = ok_dyn_ins
604	!Config Desc = sorties instantanees dans la dynamique
605	!Config Def = n
606	!Config Help =
607	!Config
608	ok_dyn_ins = .FALSE.
609	CALL getin('ok_dyn_ins',ok_dyn_ins)
610
611	!Config Key = ok_dyn_ave
612	!Config Desc = sorties moyennes dans la dynamique
613	!Config Def = n
614	!Config Help =
615	!Config
616	ok_dyn_ave = .FALSE.
617	CALL getin('ok_dyn_ave',ok_dyn_ave)
618
619
620	write(lunout,*)' #########################################'
621	write(lunout,*)' Configuration des parametres du gcm: '
622	write(lunout,*)' planet_type = ', planet_type
623	write(lunout,*)' calend = ', calend
624	write(lunout,*)' dayref = ', dayref
625	write(lunout,*)' anneeref = ', anneeref
626	write(lunout,*)' nday = ', nday
627	write(lunout,*)' day_step = ', day_step
628	write(lunout,*)' iperiod = ', iperiod
629	write(lunout,*)' iconser = ', iconser
630	write(lunout,*)' iecri = ', iecri
631	write(lunout,*)' periodav = ', periodav
632	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
633	write(lunout,*)' dissip_period = ', dissip_period
634	write(lunout,*)' lstardis = ', lstardis
635	write(lunout,*)' nitergdiv = ', nitergdiv
636	write(lunout,*)' nitergrot = ', nitergrot
637	write(lunout,*)' niterh = ', niterh
638	write(lunout,*)' tetagdiv = ', tetagdiv
639	write(lunout,*)' tetagrot = ', tetagrot
640	write(lunout,*)' tetatemp = ', tetatemp
641	write(lunout,*)' coefdis = ', coefdis
642	write(lunout,*)' purmats = ', purmats
643	write(lunout,*)' read_start = ', read_start
644	write(lunout,*)' iflag_phys = ', iflag_phys
645	write(lunout,*)' iphysiq = ', iphysiq
646	write(lunout,*)' clonn = ', clonn
647	write(lunout,*)' clatt = ', clatt
648	write(lunout,*)' grossismx = ', grossismx
649	write(lunout,*)' grossismy = ', grossismy
650	write(lunout,*)' fxyhypbb = ', fxyhypbb
651	write(lunout,*)' dzoomxx = ', dzoomxx
652	write(lunout,*)' dzoomy = ', dzoomyy
653	write(lunout,*)' tauxx = ', tauxx
654	write(lunout,*)' tauyy = ', tauyy
655	write(lunout,*)' offline = ', offline
656	write(lunout,*)' type_trac = ', type_trac
657	write(lunout,*)' config_inca = ', config_inca
658	write(lunout,*)' ok_dynzon = ', ok_dynzon
659	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
660	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
661
662	RETURN
663	c ...............................................
664	c
665	100 CONTINUE
666	!Config Key = clon
667	!Config Desc = centre du zoom, longitude
668	!Config Def = 0
669	!Config Help = longitude en degres du centre
670	!Config du zoom
671	clon = 0.
672	CALL getin('clon',clon)
673
674	!Config Key = clat
675	!Config Desc = centre du zoom, latitude
676	!Config Def = 0
677	!Config Help = latitude en degres du centre du zoom
678	!Config
679	clat = 0.
680	CALL getin('clat',clat)
681
682	!Config Key = grossismx
683	!Config Desc = zoom en longitude
684	!Config Def = 1.0
685	!Config Help = facteur de grossissement du zoom,
686	!Config selon la longitude
687	grossismx = 1.0
688	CALL getin('grossismx',grossismx)
689
690	!Config Key = grossismy
691	!Config Desc = zoom en latitude
692	!Config Def = 1.0
693	!Config Help = facteur de grossissement du zoom,
694	!Config selon la latitude
695	grossismy = 1.0
696	CALL getin('grossismy',grossismy)
697
698	IF( grossismx.LT.1. ) THEN
699	write(lunout,*)
700	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
701	STOP
702	ELSE
703	alphax = 1. - 1./ grossismx
704	ENDIF
705
706
707	IF( grossismy.LT.1. ) THEN
708	write(lunout,*)
709	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
710	STOP
711	ELSE
712	alphay = 1. - 1./ grossismy
713	ENDIF
714
715	write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
716	c
717	c alphax et alphay sont les anciennes formulat. des grossissements
718	c
719	c
720
721	!Config Key = fxyhypb
722	!Config Desc = Fonction hyperbolique
723	!Config Def = y
724	!Config Help = Fonction f(y) hyperbolique si = .true.
725	!Config sinon sinusoidale
726	fxyhypb = .TRUE.
727	CALL getin('fxyhypb',fxyhypb)
728
729	!Config Key = dzoomx
730	!Config Desc = extension en longitude
731	!Config Def = 0
732	!Config Help = extension en longitude de la zone du zoom
733	!Config ( fraction de la zone totale)
734	dzoomx = 0.0
735	CALL getin('dzoomx',dzoomx)
736
737	!Config Key = dzoomy
738	!Config Desc = extension en latitude
739	!Config Def = 0
740	!Config Help = extension en latitude de la zone du zoom
741	!Config ( fraction de la zone totale)
742	dzoomy = 0.0
743	CALL getin('dzoomy',dzoomy)
744
745	!Config Key = taux
746	!Config Desc = raideur du zoom en X
747	!Config Def = 3
748	!Config Help = raideur du zoom en X
749	taux = 3.0
750	CALL getin('taux',taux)
751
752	!Config Key = tauy
753	!Config Desc = raideur du zoom en Y
754	!Config Def = 3
755	!Config Help = raideur du zoom en Y
756	tauy = 3.0
757	CALL getin('tauy',tauy)
758
759	!Config Key = ysinus
760	!Config IF = !fxyhypb
761	!Config Desc = Fonction en Sinus
762	!Config Def = y
763	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
764	!Config sinon y = latit.
765	ysinus = .TRUE.
766	CALL getin('ysinus',ysinus)
767	c
768	!Config Key = offline
769	!Config Desc = Nouvelle eau liquide
770	!Config Def = n
771	!Config Help = Permet de mettre en route la
772	!Config nouvelle parametrisation de l'eau liquide !
773	offline = .FALSE.
774	CALL getin('offline',offline)
775
776	!Config Key = type_trac
777	!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS
778	!Config Def = lmdz
779	!Config Help =
780	!Config 'lmdz' = pas de couplage, pur LMDZ
781	!Config 'inca' = model de chime INCA
782	!Config 'repr' = model de chime REPROBUS
783	type_trac = 'lmdz'
784	CALL getin('type_trac',type_trac)
785
786	!Config Key = config_inca
787	!Config Desc = Choix de configuration de INCA
788	!Config Def = none
789	!Config Help = Choix de configuration de INCA :
790	!Config 'none' = sans INCA
791	!Config 'chem' = INCA avec calcul de chemie
792	!Config 'aero' = INCA avec calcul des aerosols
793	config_inca = 'none'
794	CALL getin('config_inca',config_inca)
795
796	!Config Key = ok_dynzon
797	!Config Desc = sortie des transports zonaux dans la dynamique
798	!Config Def = n
799	!Config Help =
800	!Config
801	ok_dynzon = .FALSE.
802	CALL getin('ok_dynzon',ok_dynzon)
803
804	!Config Key = ok_dyn_ins
805	!Config Desc = sorties instantanees dans la dynamique
806	!Config Def = n
807	!Config Help =
808	!Config
809	ok_dyn_ins = .FALSE.
810	CALL getin('ok_dyn_ins',ok_dyn_ins)
811
812	!Config Key = ok_dyn_ave
813	!Config Desc = sorties moyennes dans la dynamique
814	!Config Def = n
815	!Config Help =
816	!Config
817	ok_dyn_ave = .FALSE.
818	CALL getin('ok_dyn_ave',ok_dyn_ave)
819
820	!Config Key = use_filtre_fft
821	!Config Desc = flag d'activation des FFT pour le filtre
822	!Config Def = false
823	!Config Help = permet d'activer l'utilisation des FFT pour effectuer
824	!Config le filtrage aux poles.
825	! Le filtre fft n'est pas implemente dans dyn3d
826	use_filtre_fft=.FALSE.
827	CALL getin('use_filtre_fft',use_filtre_fft)
828
829	IF (use_filtre_fft) THEN
830	write(lunout,*)'STOP !!!'
831	write(lunout,*)'use_filtre_fft n est pas implemente dans dyn3d'
832	STOP
833	ENDIF
834
835	!Config key = ok_strato
836	!Config Desc = activation de la version strato
837	!Config Def = .FALSE.
838	!Config Help = active la version stratosphérique de LMDZ de F. Lott
839
840	ok_strato=.FALSE.
841	CALL getin('ok_strato',ok_strato)
842
843	!Config Key = ok_gradsfile
844	!Config Desc = activation des sorties grads du guidage
845	!Config Def = n
846	!Config Help = active les sorties grads du guidage
847
848	ok_gradsfile = .FALSE.
849	CALL getin('ok_gradsfile',ok_gradsfile)
850
851	!Config Key = ok_limit
852	!Config Desc = creation des fichiers limit dans create_etat0_limit
853	!Config Def = y
854	!Config Help = production du fichier limit.nc requise
855
856	ok_limit = .TRUE.
857	CALL getin('ok_limit',ok_limit)
858
859	!Config Key = ok_etat0
860	!Config Desc = creation des fichiers etat0 dans create_etat0_limit
861	!Config Def = y
862	!Config Help = production des fichiers start.nc, startphy.nc requise
863
864	ok_etat0 = .TRUE.
865	CALL getin('ok_etat0',ok_etat0)
866
867	!Config Key = grilles_gcm_netcdf
868	!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit
869	!Config Def = n
870	grilles_gcm_netcdf = .FALSE.
871	CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf)
872
873	write(lunout,*)' #########################################'
874	write(lunout,*)' Configuration des parametres de cel0'
875	& //'_limit: '
876	write(lunout,*)' planet_type = ', planet_type
877	write(lunout,*)' calend = ', calend
878	write(lunout,*)' dayref = ', dayref
879	write(lunout,*)' anneeref = ', anneeref
880	write(lunout,*)' nday = ', nday
881	write(lunout,*)' day_step = ', day_step
882	write(lunout,*)' iperiod = ', iperiod
883	write(lunout,*)' iconser = ', iconser
884	write(lunout,*)' iecri = ', iecri
885	write(lunout,*)' periodav = ', periodav
886	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
887	write(lunout,*)' dissip_period = ', dissip_period
888	write(lunout,*)' lstardis = ', lstardis
889	write(lunout,*)' nitergdiv = ', nitergdiv
890	write(lunout,*)' nitergrot = ', nitergrot
891	write(lunout,*)' niterh = ', niterh
892	write(lunout,*)' tetagdiv = ', tetagdiv
893	write(lunout,*)' tetagrot = ', tetagrot
894	write(lunout,*)' tetatemp = ', tetatemp
895	write(lunout,*)' coefdis = ', coefdis
896	write(lunout,*)' purmats = ', purmats
897	write(lunout,*)' read_start = ', read_start
898	write(lunout,*)' iflag_phys = ', iflag_phys
899	write(lunout,*)' iphysiq = ', iphysiq
900	write(lunout,*)' clon = ', clon
901	write(lunout,*)' clat = ', clat
902	write(lunout,*)' grossismx = ', grossismx
903	write(lunout,*)' grossismy = ', grossismy
904	write(lunout,*)' fxyhypb = ', fxyhypb
905	write(lunout,*)' dzoomx = ', dzoomx
906	write(lunout,*)' dzoomy = ', dzoomy
907	write(lunout,*)' taux = ', taux
908	write(lunout,*)' tauy = ', tauy
909	write(lunout,*)' offline = ', offline
910	write(lunout,*)' type_trac = ', type_trac
911	write(lunout,*)' config_inca = ', config_inca
912	write(lunout,*)' ok_dynzon = ', ok_dynzon
913	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
914	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
915	write(lunout,*)' ok_strato = ', ok_strato
916	write(lunout,*)' ok_gradsfile = ', ok_gradsfile
917	write(lunout,*)' ok_limit = ', ok_limit
918	write(lunout,*)' ok_etat0 = ', ok_etat0
919	write(lunout,*)' grilles_gcm_netcdf = ', grilles_gcm_netcdf
920	c
921	RETURN
922	END

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

Download in other formats:

Original Format