Context Navigation

conf_gcm.F @ 825

Last change on this file since 825 was 617, checked in by slebonnois, 13 years ago
SL: mise a jour de conf_gcm (cosmetique)
File size: 26.8 KB

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

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

Download in other formats:

Original Format