Context Navigation

conf_gcm.F @ 119

Last change on this file since 119 was 119, checked in by slebonnois, 14 years ago
Sebastien Lebonnois: apres validation des versions Venus et Titan, correction d'un certain nombre de bugs.
File size: 30.8 KB

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

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

Download in other formats:

Original Format