Context Navigation

conf_gcm.F90 @ 5118

Last change on this file since 5118 was 5118, checked in by abarral, 8 weeks ago
Replace iniprint.h by lmdz_iniprint.f90 (lint) along the way
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.8 KB

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

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

Download in other formats:

Original Format