Context Navigation

conf_gcm.F @ 271

Last change on this file since 271 was 271, checked in by emillour, 13 years ago
Petit oubli lors de la mise a jour de la dynamique de la ver 270. Ehouarn
File size: 31.2 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 = dissip_period
244	!Config Desc = periode de la dissipation
245	!Config Def = 0
246	!Config Help = periode de la dissipation
247	!Config dissip_period=0 => la valeur sera calcule dans inidissip
248	!Config dissip_period>0 => on prend cette valeur
249	dissip_period = 0
250	CALL getin('dissip_period',dissip_period)
251
252	ccc .... P. Le Van , modif le 29/04/97 .pour la dissipation ...
253	ccc
254
255	!Config Key = lstardis
256	!Config Desc = choix de l'operateur de dissipation
257	!Config Def = y
258	!Config Help = choix de l'operateur de dissipation
259	!Config 'y' si on veut star et 'n' si on veut non-start !
260	!Config Moi y en a pas comprendre !
261	lstardis = .TRUE.
262	CALL getin('lstardis',lstardis)
263
264
265	!Config Key = nitergdiv
266	!Config Desc = Nombre d'iteration de gradiv
267	!Config Def = 1
268	!Config Help = nombre d'iterations de l'operateur de dissipation
269	!Config gradiv
270	nitergdiv = 1
271	CALL getin('nitergdiv',nitergdiv)
272
273	!Config Key = nitergrot
274	!Config Desc = nombre d'iterations de nxgradrot
275	!Config Def = 2
276	!Config Help = nombre d'iterations de l'operateur de dissipation
277	!Config nxgradrot
278	nitergrot = 2
279	CALL getin('nitergrot',nitergrot)
280
281
282	!Config Key = niterh
283	!Config Desc = nombre d'iterations de divgrad
284	!Config Def = 2
285	!Config Help = nombre d'iterations de l'operateur de dissipation
286	!Config divgrad
287	niterh = 2
288	CALL getin('niterh',niterh)
289
290
291	!Config Key = tetagdiv
292	!Config Desc = temps de dissipation pour div
293	!Config Def = 7200
294	!Config Help = temps de dissipation des plus petites longeur
295	!Config d'ondes pour u,v (gradiv)
296	tetagdiv = 7200.
297	CALL getin('tetagdiv',tetagdiv)
298
299	!Config Key = tetagrot
300	!Config Desc = temps de dissipation pour grad
301	!Config Def = 7200
302	!Config Help = temps de dissipation des plus petites longeur
303	!Config d'ondes pour u,v (nxgradrot)
304	tetagrot = 7200.
305	CALL getin('tetagrot',tetagrot)
306
307	!Config Key = tetatemp
308	!Config Desc = temps de dissipation pour h
309	!Config Def = 7200
310	!Config Help = temps de dissipation des plus petites longeur
311	!Config d'ondes pour h (divgrad)
312	tetatemp = 7200.
313	CALL getin('tetatemp',tetatemp )
314
315	! Parametres controlant la variation sur la verticale des constantes de
316	! dissipation.
317	! Actifs uniquement avec ok_strato=y
318
319	dissip_fac_mid=2.
320	dissip_fac_up=10.
321	dissip_deltaz=10.! Intervalle (km) pour le changement mid / up
322	dissip_hdelta=5. ! scale height (km) dans la zone de la transition(m)
323	dissip_pupstart=1.e3 ! pression (Pa) au bas la transition mid / up
324	CALL getin('dissip_fac_mid',dissip_fac_mid )
325	CALL getin('dissip_fac_up',dissip_fac_up )
326	CALL getin('dissip_deltaz',dissip_deltaz )
327	CALL getin('dissip_hdelta',dissip_hdelta )
328	CALL getin('dissip_pupstart',dissip_pupstart )
329
330	! Parametres controlant la couche eponge
331	! Actifs uniquement avec ok_strato=y
332	! mode = 0 : pas de sponge
333	! mode = 1 : u et v -> 0
334	! mode = 2 : u et v -> moyenne zonale
335	! mode = 3 : u, v et h -> moyenne zonale
336	iflag_top_bound=1
337	mode_top_bound=3
338	tau_top_bound=1.e-5
339	CALL getin('iflag_top_bound',iflag_top_bound)
340	CALL getin('mode_top_bound',mode_top_bound)
341	CALL getin('tau_top_bound',tau_top_bound)
342
343	!
344	!Config Key = coefdis
345	!Config Desc = coefficient pour gamdissip
346	!Config Def = 0
347	!Config Help = coefficient pour gamdissip
348	coefdis = 0.
349	CALL getin('coefdis',coefdis)
350
351	!Config Key = purmats
352	!Config Desc = Schema d'integration
353	!Config Def = n
354	!Config Help = Choix du schema d'integration temporel.
355	!Config y = pure Matsuno sinon c'est du Matsuno-leapfrog
356	purmats = .FALSE.
357	CALL getin('purmats',purmats)
358
359	!Config Key = ok_guide
360	!Config Desc = Guidage
361	!Config Def = n
362	!Config Help = Guidage
363	ok_guide = .FALSE.
364	CALL getin('ok_guide',ok_guide)
365
366	c ...............................................................
367
368	!Config Key = read_start
369	!Config Desc = Initialize model using a 'start.nc' file
370	!Config Def = y
371	!Config Help = y: intialize dynamical fields using a 'start.nc' file
372	! n: fields are initialized by 'iniacademic' routine
373	read_start= .true.
374	CALL getin('read_start',read_start)
375
376	!Config Key = iflag_phys
377	!Config Desc = Avec ls physique
378	!Config Def = 1
379	!Config Help = Permet de faire tourner le modele sans
380	!Config physique.
381	iflag_phys = 1
382	CALL getin('iflag_phys',iflag_phys)
383
384
385	!Config Key = iphysiq
386	!Config Desc = Periode de la physique
387	!Config Def = 5
388	!Config Help = Periode de la physique en pas de temps de la dynamique.
389	iphysiq = 5
390	CALL getin('iphysiq', iphysiq)
391
392	!Config Key = iflag_trac
393	!Config Desc = traceurs presents ou non
394	!Config Def = 1
395	!Config Help = Permet de faire tourner le modele sans traceurs
396	!Config
397	iflag_trac = 1
398	CALL getin('iflag_trac',iflag_trac)
399
400	!Config Key = ip_ebil_dyn
401	!Config Desc = PRINT level for energy conserv. diag.
402	!Config Def = 0
403	!Config Help = PRINT level for energy conservation diag. ;
404	! les options suivantes existent :
405	!Config 0 pas de print
406	!Config 1 pas de print
407	!Config 2 print,
408	ip_ebil_dyn = 0
409	CALL getin('ip_ebil_dyn',ip_ebil_dyn)
410	!
411
412
413	ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ...
414	c ......... ( modif le 17/04/96 ) .........
415	c
416	IF( etatinit ) GO TO 100
417
418	!Config Key = clon
419	!Config Desc = centre du zoom, longitude
420	!Config Def = 0
421	!Config Help = longitude en degres du centre
422	!Config du zoom
423	clonn = 0.
424	CALL getin('clon',clonn)
425
426	!Config Key = clat
427	!Config Desc = centre du zoom, latitude
428	!Config Def = 0
429	!Config Help = latitude en degres du centre du zoom
430	!Config
431	clatt = 0.
432	CALL getin('clat',clatt)
433
434	c
435	c
436	IF( ABS(clat - clatt).GE. 0.001 ) THEN
437	write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
438	& ' est differente de celle lue sur le fichier start '
439	STOP
440	ENDIF
441
442	!Config Key = grossismx
443	!Config Desc = zoom en longitude
444	!Config Def = 1.0
445	!Config Help = facteur de grossissement du zoom,
446	!Config selon la longitude
447	grossismxx = 1.0
448	CALL getin('grossismx',grossismxx)
449
450
451	IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN
452	write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
453	& 'run.def est differente de celle lue sur le fichier start '
454	STOP
455	ENDIF
456
457	!Config Key = grossismy
458	!Config Desc = zoom en latitude
459	!Config Def = 1.0
460	!Config Help = facteur de grossissement du zoom,
461	!Config selon la latitude
462	grossismyy = 1.0
463	CALL getin('grossismy',grossismyy)
464
465	IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN
466	write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
467	& 'run.def est differente de celle lue sur le fichier start '
468	STOP
469	ENDIF
470
471	IF( grossismx.LT.1. ) THEN
472	write(lunout,*)
473	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
474	STOP
475	ELSE
476	alphax = 1. - 1./ grossismx
477	ENDIF
478
479
480	IF( grossismy.LT.1. ) THEN
481	write(lunout,*)
482	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
483	STOP
484	ELSE
485	alphay = 1. - 1./ grossismy
486	ENDIF
487
488	write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
489	c
490	c alphax et alphay sont les anciennes formulat. des grossissements
491	c
492	c
493
494	!Config Key = fxyhypb
495	!Config Desc = Fonction hyperbolique
496	!Config Def = y
497	!Config Help = Fonction f(y) hyperbolique si = .true.
498	!Config sinon sinusoidale
499	fxyhypbb = .TRUE.
500	CALL getin('fxyhypb',fxyhypbb)
501
502	IF( .NOT.fxyhypb ) THEN
503	IF( fxyhypbb ) THEN
504	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
505	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
506	* 'F alors qu il est T sur run.def ***'
507	STOP
508	ENDIF
509	ELSE
510	IF( .NOT.fxyhypbb ) THEN
511	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
512	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
513	* 'T alors qu il est F sur run.def **** '
514	STOP
515	ENDIF
516	ENDIF
517	c
518	!Config Key = dzoomx
519	!Config Desc = extension en longitude
520	!Config Def = 0
521	!Config Help = extension en longitude de la zone du zoom
522	!Config ( fraction de la zone totale)
523	dzoomxx = 0.0
524	CALL getin('dzoomx',dzoomxx)
525
526	IF( fxyhypb ) THEN
527	IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN
528	write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
529	* 'run.def est differente de celle lue sur le fichier start '
530	STOP
531	ENDIF
532	ENDIF
533
534	!Config Key = dzoomy
535	!Config Desc = extension en latitude
536	!Config Def = 0
537	!Config Help = extension en latitude de la zone du zoom
538	!Config ( fraction de la zone totale)
539	dzoomyy = 0.0
540	CALL getin('dzoomy',dzoomyy)
541
542	IF( fxyhypb ) THEN
543	IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN
544	write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
545	* 'run.def est differente de celle lue sur le fichier start '
546	STOP
547	ENDIF
548	ENDIF
549
550	!Config Key = taux
551	!Config Desc = raideur du zoom en X
552	!Config Def = 3
553	!Config Help = raideur du zoom en X
554	tauxx = 3.0
555	CALL getin('taux',tauxx)
556
557	IF( fxyhypb ) THEN
558	IF( ABS(taux - tauxx).GE. 0.001 ) THEN
559	write(lunout,*)'conf_gcm: La valeur de taux passee par ',
560	* 'run.def est differente de celle lue sur le fichier start '
561	STOP
562	ENDIF
563	ENDIF
564
565	!Config Key = tauyy
566	!Config Desc = raideur du zoom en Y
567	!Config Def = 3
568	!Config Help = raideur du zoom en Y
569	tauyy = 3.0
570	CALL getin('tauy',tauyy)
571
572	IF( fxyhypb ) THEN
573	IF( ABS(tauy - tauyy).GE. 0.001 ) THEN
574	write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
575	* 'run.def est differente de celle lue sur le fichier start '
576	STOP
577	ENDIF
578	ENDIF
579
580	cc
581	IF( .NOT.fxyhypb ) THEN
582
583	!Config Key = ysinus
584	!Config IF = !fxyhypb
585	!Config Desc = Fonction en Sinus
586	!Config Def = y
587	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
588	!Config sinon y = latit.
589	ysinuss = .TRUE.
590	CALL getin('ysinus',ysinuss)
591
592	IF( .NOT.ysinus ) THEN
593	IF( ysinuss ) THEN
594	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
595	write(lunout,)' ** ysinus lu sur le fichier start est F',
596	* ' alors qu il est T sur run.def ***'
597	STOP
598	ENDIF
599	ELSE
600	IF( .NOT.ysinuss ) THEN
601	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
602	write(lunout,)' ** ysinus lu sur le fichier start est T',
603	* ' alors qu il est F sur run.def **** '
604	STOP
605	ENDIF
606	ENDIF
607	ENDIF ! of IF( .NOT.fxyhypb )
608	c
609	!Config Key = offline
610	!Config Desc = Nouvelle eau liquide
611	!Config Def = n
612	!Config Help = Permet de mettre en route la
613	!Config nouvelle parametrisation de l'eau liquide !
614	offline = .FALSE.
615	CALL getin('offline',offline)
616	IF (offline .AND. adjust) THEN
617	WRITE(lunout,*)
618	& 'WARNING : option offline does not work with adjust=y :'
619	WRITE(lunout,*) 'the files defstoke.nc, fluxstoke.nc ',
620	& 'and fluxstokev.nc will not be created'
621	WRITE(lunout,*)
622	& 'only the file phystoke.nc will still be created '
623	END IF
624
625	!Config Key = config_inca
626	!Config Desc = Choix de configuration de INCA
627	!Config Def = none
628	!Config Help = Choix de configuration de INCA :
629	!Config 'none' = sans INCA
630	!Config 'chem' = INCA avec calcul de chemie
631	!Config 'aero' = INCA avec calcul des aerosols
632	config_inca = 'none'
633	CALL getin('config_inca',config_inca)
634
635	!Config Key = ok_dynzon
636	!Config Desc = calcul et sortie des transports
637	!Config Def = n
638	!Config Help = Permet de mettre en route le calcul des transports
639	!Config
640	ok_dynzon = .FALSE.
641	CALL getin('ok_dynzon',ok_dynzon)
642
643	!Config Key = ok_dyn_ins
644	!Config Desc = sorties instantanees dans la dynamique
645	!Config Def = n
646	!Config Help =
647	!Config
648	ok_dyn_ins = .FALSE.
649	CALL getin('ok_dyn_ins',ok_dyn_ins)
650
651	!Config Key = ok_dyn_ave
652	!Config Desc = sorties moyennes dans la dynamique
653	!Config Def = n
654	!Config Help =
655	!Config
656	ok_dyn_ave = .FALSE.
657	CALL getin('ok_dyn_ave',ok_dyn_ave)
658
659	write(lunout,*)' #########################################'
660	write(lunout,*)' Configuration des parametres du gcm: '
661	write(lunout,*)' planet_type = ', planet_type
662	write(lunout,*)' calend = ', calend
663	write(lunout,*)' dayref = ', dayref
664	write(lunout,*)' anneeref = ', anneeref
665	write(lunout,*)' nday = ', nday
666	if (less1day) then
667	write(lunout,*)' Run only for a fraction of day ! '
668	write(lunout,*)' fractday = ', fractday
669	endif
670	write(lunout,*)' day_step = ', day_step
671	write(lunout,*)' iperiod = ', iperiod
672	write(lunout,*)' nsplit_phys = ', nsplit_phys
673	write(lunout,*)' iconser = ', iconser
674	write(lunout,*)' iecri = ', iecri
675	write(lunout,*)' periodav = ', periodav
676	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
677	write(lunout,*)' dissip_period = ', dissip_period
678	write(lunout,*)' lstardis = ', lstardis
679	write(lunout,*)' nitergdiv = ', nitergdiv
680	write(lunout,*)' nitergrot = ', nitergrot
681	write(lunout,*)' niterh = ', niterh
682	write(lunout,*)' tetagdiv = ', tetagdiv
683	write(lunout,*)' tetagrot = ', tetagrot
684	write(lunout,*)' tetatemp = ', tetatemp
685	write(lunout,*)' coefdis = ', coefdis
686	write(lunout,*)' purmats = ', purmats
687	write(lunout,*)' read_start = ', read_start
688	write(lunout,*)' iflag_phys = ', iflag_phys
689	write(lunout,*)' iphysiq = ', iphysiq
690	write(lunout,*)' iflag_trac = ', iflag_trac
691	write(lunout,*)' clonn = ', clonn
692	write(lunout,*)' clatt = ', clatt
693	write(lunout,*)' grossismx = ', grossismx
694	write(lunout,*)' grossismy = ', grossismy
695	write(lunout,*)' fxyhypbb = ', fxyhypbb
696	write(lunout,*)' dzoomxx = ', dzoomxx
697	write(lunout,*)' dzoomy = ', dzoomyy
698	write(lunout,*)' tauxx = ', tauxx
699	write(lunout,*)' tauyy = ', tauyy
700	write(lunout,*)' offline = ', offline
701	write(lunout,*)' config_inca = ', config_inca
702	write(lunout,*)' ok_dynzon = ', ok_dynzon
703	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
704	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
705
706	RETURN
707	c ...............................................
708	c
709	100 CONTINUE
710	!Config Key = clon
711	!Config Desc = centre du zoom, longitude
712	!Config Def = 0
713	!Config Help = longitude en degres du centre
714	!Config du zoom
715	clon = 0.
716	CALL getin('clon',clon)
717
718	!Config Key = clat
719	!Config Desc = centre du zoom, latitude
720	!Config Def = 0
721	!Config Help = latitude en degres du centre du zoom
722	!Config
723	clat = 0.
724	CALL getin('clat',clat)
725
726	!Config Key = grossismx
727	!Config Desc = zoom en longitude
728	!Config Def = 1.0
729	!Config Help = facteur de grossissement du zoom,
730	!Config selon la longitude
731	grossismx = 1.0
732	CALL getin('grossismx',grossismx)
733
734	!Config Key = grossismy
735	!Config Desc = zoom en latitude
736	!Config Def = 1.0
737	!Config Help = facteur de grossissement du zoom,
738	!Config selon la latitude
739	grossismy = 1.0
740	CALL getin('grossismy',grossismy)
741
742	IF( grossismx.LT.1. ) THEN
743	write(lunout,*)
744	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
745	STOP
746	ELSE
747	alphax = 1. - 1./ grossismx
748	ENDIF
749
750
751	IF( grossismy.LT.1. ) THEN
752	write(lunout,*)
753	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
754	STOP
755	ELSE
756	alphay = 1. - 1./ grossismy
757	ENDIF
758
759	write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
760	c
761	c alphax et alphay sont les anciennes formulat. des grossissements
762	c
763	c
764
765	!Config Key = fxyhypb
766	!Config Desc = Fonction hyperbolique
767	!Config Def = y
768	!Config Help = Fonction f(y) hyperbolique si = .true.
769	!Config sinon sinusoidale
770	fxyhypb = .TRUE.
771	CALL getin('fxyhypb',fxyhypb)
772
773	!Config Key = dzoomx
774	!Config Desc = extension en longitude
775	!Config Def = 0
776	!Config Help = extension en longitude de la zone du zoom
777	!Config ( fraction de la zone totale)
778	dzoomx = 0.0
779	CALL getin('dzoomx',dzoomx)
780
781	!Config Key = dzoomy
782	!Config Desc = extension en latitude
783	!Config Def = 0
784	!Config Help = extension en latitude de la zone du zoom
785	!Config ( fraction de la zone totale)
786	dzoomy = 0.0
787	CALL getin('dzoomy',dzoomy)
788
789	!Config Key = taux
790	!Config Desc = raideur du zoom en X
791	!Config Def = 3
792	!Config Help = raideur du zoom en X
793	taux = 3.0
794	CALL getin('taux',taux)
795
796	!Config Key = tauy
797	!Config Desc = raideur du zoom en Y
798	!Config Def = 3
799	!Config Help = raideur du zoom en Y
800	tauy = 3.0
801	CALL getin('tauy',tauy)
802
803	!Config Key = ysinus
804	!Config IF = !fxyhypb
805	!Config Desc = Fonction en Sinus
806	!Config Def = y
807	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
808	!Config sinon y = latit.
809	ysinus = .TRUE.
810	CALL getin('ysinus',ysinus)
811	c
812	!Config Key = offline
813	!Config Desc = Nouvelle eau liquide
814	!Config Def = n
815	!Config Help = Permet de mettre en route la
816	!Config nouvelle parametrisation de l'eau liquide !
817	offline = .FALSE.
818	CALL getin('offline',offline)
819	IF (offline .AND. adjust) THEN
820	WRITE(lunout,*)
821	& 'WARNING : option offline does not work with adjust=y :'
822	WRITE(lunout,*) 'the files defstoke.nc, fluxstoke.nc ',
823	& 'and fluxstokev.nc will not be created'
824	WRITE(lunout,*)
825	& 'only the file phystoke.nc will still be created '
826	END IF
827
828	!Config Key = config_inca
829	!Config Desc = Choix de configuration de INCA
830	!Config Def = none
831	!Config Help = Choix de configuration de INCA :
832	!Config 'none' = sans INCA
833	!Config 'chem' = INCA avec calcul de chemie
834	!Config 'aero' = INCA avec calcul des aerosols
835	config_inca = 'none'
836	CALL getin('config_inca',config_inca)
837
838	!Config Key = ok_dynzon
839	!Config Desc = calcul et sortie des transports
840	!Config Def = n
841	!Config Help = Permet de mettre en route le calcul des transports
842	!Config
843	ok_dynzon = .FALSE.
844	CALL getin('ok_dynzon',ok_dynzon)
845
846	!Config Key = ok_dyn_ins
847	!Config Desc = sorties instantanees dans la dynamique
848	!Config Def = n
849	!Config Help =
850	!Config
851	ok_dyn_ins = .FALSE.
852	CALL getin('ok_dyn_ins',ok_dyn_ins)
853
854	!Config Key = ok_dyn_ave
855	!Config Desc = sorties moyennes dans la dynamique
856	!Config Def = n
857	!Config Help =
858	!Config
859	ok_dyn_ave = .FALSE.
860	CALL getin('ok_dyn_ave',ok_dyn_ave)
861
862	!Config Key = use_filtre_fft
863	!Config Desc = flag d'activation des FFT pour le filtre
864	!Config Def = false
865	!Config Help = permet d'activer l'utilisation des FFT pour effectuer
866	!Config le filtrage aux poles.
867	use_filtre_fft=.FALSE.
868	CALL getin('use_filtre_fft',use_filtre_fft)
869
870	IF (use_filtre_fft .AND. grossismx /= 1.0) THEN
871	write(lunout,*)'WARNING !!! '
872	write(lunout,*)"Le zoom en longitude est incompatible",
873	& " avec l'utilisation du filtre FFT ",
874	& "---> filtre FFT dÃ©sactivÃ© "
875	use_filtre_fft=.FALSE.
876	ENDIF
877
878
879
880	!Config Key = use_mpi_alloc
881	!Config Desc = Utilise un buffer MPI en mï¿½moire globale
882	!Config Def = false
883	!Config Help = permet d'activer l'utilisation d'un buffer MPI
884	!Config en mï¿½moire globale a l'aide de la fonction MPI_ALLOC.
885	!Config Cela peut amï¿½liorer la bande passante des transferts MPI
886	!Config d'un facteur 2
887	use_mpi_alloc=.FALSE.
888	CALL getin('use_mpi_alloc',use_mpi_alloc)
889
890	!Config Key = omp_chunk
891	!Config Desc = taille des blocs openmp
892	!Config Def = 1
893	!Config Help = defini la taille des packets d'itï¿½ration openmp
894	!Config distribuï¿½e ï¿½ chaque tï¿½che lors de l'entrï¿½e dans une
895	!Config boucle parallï¿½lisï¿½e
896
897	omp_chunk=1
898	CALL getin('omp_chunk',omp_chunk)
899
900	!Config key = ok_strato
901	!Config Desc = activation de la version strato
902	!Config Def = .FALSE.
903	!Config Help = active la version stratosphérique de LMDZ de F. Lott
904
905	ok_strato=.TRUE.
906	CALL getin('ok_strato',ok_strato)
907
908	!Config Key = ok_gradsfile
909	!Config Desc = activation des sorties grads du guidage
910	!Config Def = n
911	!Config Help = active les sorties grads du guidage
912
913	ok_gradsfile = .FALSE.
914	CALL getin('ok_gradsfile',ok_gradsfile)
915
916	!Config Key = ok_limit
917	!Config Desc = creation des fichiers limit dans create_etat0_limit
918	!Config Def = y
919	!Config Help = production du fichier limit.nc requise
920
921	ok_limit = .TRUE.
922	CALL getin('ok_limit',ok_limit)
923
924	!Config Key = ok_etat0
925	!Config Desc = creation des fichiers etat0 dans create_etat0_limit
926	!Config Def = y
927	!Config Help = production des fichiers start.nc, startphy.nc requise
928
929	ok_etat0 = .TRUE.
930	CALL getin('ok_etat0',ok_etat0)
931
932	!Config Key = grilles_gcm_netcdf
933	!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit
934	!Config Def = n
935	grilles_gcm_netcdf = .FALSE.
936	CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf)
937
938	write(lunout,*)' #########################################'
939	write(lunout,*)' Configuration des parametres de cel0'
940	& //'_limit: '
941	write(lunout,*)' planet_type = ', planet_type
942	write(lunout,*)' calend = ', calend
943	write(lunout,*)' dayref = ', dayref
944	write(lunout,*)' anneeref = ', anneeref
945	write(lunout,*)' nday = ', nday
946	write(lunout,*)' day_step = ', day_step
947	write(lunout,*)' iperiod = ', iperiod
948	write(lunout,*)' iconser = ', iconser
949	write(lunout,*)' iecri = ', iecri
950	write(lunout,*)' periodav = ', periodav
951	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
952	write(lunout,*)' dissip_period = ', dissip_period
953	write(lunout,*)' lstardis = ', lstardis
954	write(lunout,*)' nitergdiv = ', nitergdiv
955	write(lunout,*)' nitergrot = ', nitergrot
956	write(lunout,*)' niterh = ', niterh
957	write(lunout,*)' tetagdiv = ', tetagdiv
958	write(lunout,*)' tetagrot = ', tetagrot
959	write(lunout,*)' tetatemp = ', tetatemp
960	write(lunout,*)' coefdis = ', coefdis
961	write(lunout,*)' purmats = ', purmats
962	write(lunout,*)' read_start = ', read_start
963	write(lunout,*)' iflag_phys = ', iflag_phys
964	write(lunout,*)' iphysiq = ', iphysiq
965	write(lunout,*)' iflag_trac = ', iflag_trac
966	write(lunout,*)' clon = ', clon
967	write(lunout,*)' clat = ', clat
968	write(lunout,*)' grossismx = ', grossismx
969	write(lunout,*)' grossismy = ', grossismy
970	write(lunout,*)' fxyhypb = ', fxyhypb
971	write(lunout,*)' dzoomx = ', dzoomx
972	write(lunout,*)' dzoomy = ', dzoomy
973	write(lunout,*)' taux = ', taux
974	write(lunout,*)' tauy = ', tauy
975	write(lunout,*)' offline = ', offline
976	write(lunout,*)' config_inca = ', config_inca
977	write(lunout,*)' ok_dynzon = ', ok_dynzon
978	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
979	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
980	write(lunout,*)' use_filtre_fft = ', use_filtre_fft
981	write(lunout,*)' use_mpi_alloc = ', use_mpi_alloc
982	write(lunout,*)' omp_chunk = ', omp_chunk
983	write(lunout,*)' ok_strato = ', ok_strato
984	write(lunout,*)' ok_gradsfile = ', ok_gradsfile
985	write(lunout,*)' ok_limit = ', ok_limit
986	write(lunout,*)' ok_etat0 = ', ok_etat0
987	write(lunout,*)' grilles_gcm_netcdf = ', grilles_gcm_netcdf
988	c
989	RETURN
990	END

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

Download in other formats:

Original Format