Context Navigation

conf_gcm.F @ 2111

Last change on this file since 2111 was 2111, checked in by lguez, 10 years ago
Not possible to adjust MPI bands with nudging.
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
File size: 31.9 KB

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

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

Download in other formats:

Original Format