Context Navigation

conf_gcm.F90 @ 2526

Last change on this file since 2526 was 2442, checked in by fhourdin, 9 years ago

ngroup passe du statut de paramètre à constante.

Pour pouvoir l'ajuster en fonction de la resolution
ngroup est maintenant lu dans les .def. ngroup=3 par défaut

ngroup changed from parameter to constant

Property copyright set to
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory
Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 32.2 KB

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

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

Download in other formats:

Original Format