Context Navigation

conf_gcm.F @ 113

Last change on this file since 113 was 108, checked in by slebonnois, 14 years ago
Sebastien Lebonnois: sponge layer et dissip horizontale.
File size: 30.5 KB

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

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

Download in other formats:

Original Format