Context Navigation

conf_gcm.F @ 1707

Last change on this file since 1707 was 1707, checked in by Laurent Fairhead, 11 years ago

Version testing basée sur la r1706

Testing release based on r1706

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

File size: 30.9 KB

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

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

Download in other formats:

Original Format