Context Navigation

conf_gcm.F @ 1200

Last change on this file since 1200 was 1190, checked in by yann meurdesoif, 15 years ago

Correction de FH sur la sponge layer de top_bound + parallelisation des corrections.

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

File size: 26.4 KB

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

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

Download in other formats:

Original Format