Context Navigation

conf_gcm.F90 @ 2598

Last change on this file since 2598 was 2598, checked in by Ehouarn Millour, 8 years ago
Cleanup in the dynamics: turn serre.h into module serre_mod.F90 EM
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.5 KB

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

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

Download in other formats:

Original Format