Context Navigation

conf_gcm.F @ 5427

Last change on this file since 5427 was 1323, checked in by Laurent Fairhead, 15 years ago

Changes made in r1293 are integrated into the trunk
Start files are identical between r1293 and this version

Les modifications de la r1293 sont intégrées à la trunk
Les fichiers start et startphy sont identiques entre la version 1293 et celle-ci

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 25.5 KB

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

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

Download in other formats:

Original Format