Context Navigation

conf_gcm.F90 @ 5061

Last change on this file since 5061 was 4669, checked in by Laurent Fairhead, 14 months ago
Merged with trunk revision 4586 corresponding to june 2023 testing
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: 34.2 KB

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

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

Download in other formats:

Original Format