Context Navigation

conf_gcm.F @ 1664

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

Version testing basée sur la r1575

Testing release based on r1575

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

File size: 30.5 KB

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

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

Download in other formats:

Original Format