Context Navigation

conf_gcm.F @ 5451

Last change on this file since 5451 was 1438, checked in by jghattas, 14 years ago
Bug corrections for option offline. This option does not work at the same time as adjust=y.
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 29.3 KB

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

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

Download in other formats:

Original Format