Context Navigation

conf_gcm.F @ 1126

Last change on this file since 1126 was 1056, checked in by slebonnois, 12 years ago
SL: Titan runs ! see DOC/chantiers/commit_importants.log
File size: 30.6 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	use assert_m, only: assert
17	use sponge_mod, only: callsponge,mode_sponge,nsponge,tetasponge
18
19	IMPLICIT NONE
20	c-----------------------------------------------------------------------
21	c Auteurs : L. Fairhead , P. Le Van .
22	c
23	c Arguments :
24	c
25	c tapedef :
26	c etatinit : = TRUE , on ne compare pas les valeurs des para-
27	c -metres du zoom avec celles lues sur le fichier start .
28	c
29	LOGICAL etatinit
30	INTEGER tapedef
31
32	c Declarations :
33	c --------------
34	#include "dimensions.h"
35	#include "paramet.h"
36	#include "logic.h"
37	#include "serre.h"
38	#include "comdissnew.h"
39	#include "iniprint.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	c
46	c
47	c local:
48	c ------
49
50	CHARACTER ch172,ch272,ch372,ch412
51	REAL clonn,clatt,grossismxx,grossismyy
52	REAL dzoomxx,dzoomyy, tauxx,tauyy
53	LOGICAL fxyhypbb, ysinuss
54	INTEGER i
55	LOGICAL use_filtre_fft
56	c
57	c -------------------------------------------------------------------
58	c
59	c ......... Version du 29/04/97 ..........
60	c
61	c Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot,
62	c tetatemp ajoutes pour la dissipation .
63	c
64	c Autre parametre ajoute en fin de liste de tapedef : fxyhypb
65	c
66	c Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb.
67	c Sinon , choix de fxynew , a derivee sinusoidale ..
68	c
69	c ...... etatinit = . TRUE. si defrun est appele dans ETAT0_LMD ou
70	c LIMIT_LMD pour l'initialisation de start.dat (dic) et
71	c de limit.dat ( dic) ...........
72	c Sinon etatinit = . FALSE .
73	c
74	c Donc etatinit = .F. si on veut comparer les valeurs de grossismx ,
75	c grossismy,clon,clat, fxyhypb lues sur le fichier start avec
76	c celles passees par run.def , au debut du gcm, apres l'appel a
77	c lectba .
78	c Ces parmetres definissant entre autres la grille et doivent etre
79	c pareils et coherents , sinon il y aura divergence du gcm .
80	c
81	c-----------------------------------------------------------------------
82	c initialisations:
83	c ----------------
84
85	!Config Key = lunout
86	!Config Desc = unite de fichier pour les impressions
87	!Config Def = 6
88	!Config Help = unite de fichier pour les impressions
89	!Config (defaut sortie standard = 6)
90	lunout=6
91	CALL getin('lunout', lunout)
92	IF (lunout /= 5 .and. lunout /= 6) THEN
93	OPEN(UNIT=lunout,FILE='lmdz.out',ACTION='write',
94	& STATUS='unknown',FORM='formatted')
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 = resetvarc
150	!Config Desc = Reinit des variables de controle
151	!Config Def = n
152	!Config Help = Reinit des variables de controle
153	resetvarc = .false.
154	CALL getin('resetvarc',resetvarc)
155
156	!Config Key = nday
157	!Config Desc = Nombre de jours d'integration
158	!Config Def = 10
159	!Config Help = Nombre de jours d'integration
160	!Config ... On pourait aussi permettre des mois ou des annees !
161	nday = 10
162	CALL getin('nday',nday)
163
164	! alternative to specifying nday (see also 'less1day' and 'fractday'
165	! options below: sopecify numbre of dynamic steps to run:
166	ndynstep = -9999 ! default value ; if ndynstep <0 then option not used.
167	call getin('ndynstep',ndynstep)
168
169	!Config Key = starttime
170	!Config Desc = Heure de depart de la simulation
171	!Config Def = 0
172	!Config Help = Heure de depart de la simulation
173	!Config en jour
174	starttime = 0
175	CALL getin('starttime',starttime)
176
177	!Config Key = less1day
178	!Config Desc = Possibilite d'integrer moins d'un jour
179	!Config Def = n
180	!Config Help = Possibilite d'integrer moins d'un jour
181	less1day = .false.
182	CALL getin('less1day',less1day)
183
184	!Config Key = fractday
185	!Config Desc = integration sur une fraction de jour
186	!Config Def = 0.01
187	!Config Help = integration sur une fraction de jour
188	fractday = 0.01
189	CALL getin('fractday',fractday)
190
191	!Config Key = day_step
192	!Config Desc = nombre de pas par jour
193	!Config Def = 240
194	!Config Help = nombre de pas par jour (multiple de iperiod) (
195	!Config ici pour dt = 1 min )
196	day_step = 240
197	CALL getin('day_step',day_step)
198
199	!Config Key = nsplit_phys
200	!Config Desc = nombre de subdivisions par pas physique
201	!Config Def = 1
202	!Config Help = nombre de subdivisions par pas physique
203	nsplit_phys = 1
204	CALL getin('nsplit_phys',nsplit_phys)
205
206	!Config Key = iperiod
207	!Config Desc = periode pour le pas Matsuno
208	!Config Def = 5
209	!Config Help = periode pour le pas Matsuno (en pas de temps)
210	iperiod = 5
211	CALL getin('iperiod',iperiod)
212
213	!Config Key = iapp_tracvl
214	!Config Desc = frequence du groupement des flux
215	!Config Def = iperiod
216	!Config Help = frequence du groupement des flux (en pas de temps)
217	iapp_tracvl = iperiod
218	CALL getin('iapp_tracvl',iapp_tracvl)
219
220	!Config Key = iconser
221	!Config Desc = periode de sortie des variables de controle
222	!Config Def = 240
223	!Config Help = periode de sortie des variables de controle
224	!Config (En pas de temps)
225	iconser = 240
226	CALL getin('iconser', iconser)
227
228	!Config Key = iecri
229	!Config Desc = periode d'ecriture du fichier histoire
230	!Config Def = 1
231	!Config Help = periode d'ecriture du fichier histoire (en jour)
232	iecri = 1
233	CALL getin('iecri',iecri)
234
235
236	!Config Key = periodav
237	!Config Desc = periode de stockage fichier histmoy
238	!Config Def = 1
239	!Config Help = periode de stockage fichier histmoy (en jour)
240	periodav = 1.
241	CALL getin('periodav',periodav)
242
243	!Config Key = output_grads_dyn
244	!Config Desc = output dynamics diagnostics in 'dyn.dat' file
245	!Config Def = n
246	!Config Help = output dynamics diagnostics in Grads-readable 'dyn.dat' file
247	output_grads_dyn=.false.
248	CALL getin('output_grads_dyn',output_grads_dyn)
249
250	!Config Key = dissip_period
251	!Config Desc = periode de la dissipation
252	!Config Def = 0
253	!Config Help = periode de la dissipation
254	!Config dissip_period=0 => la valeur sera calcule dans inidissip
255	!Config dissip_period>0 => on prend cette valeur
256	dissip_period = 0
257	call getin('idissip',dissip_period) ! old Mars/Genreic model parameter
258	! if there is a "dissip_period" in run.def, it overrides "idissip"
259	CALL getin('dissip_period',dissip_period)
260
261	ccc .... P. Le Van , modif le 29/04/97 .pour la dissipation ...
262	ccc
263
264	!Config Key = lstardis
265	!Config Desc = choix de l'operateur de dissipation
266	!Config Def = y
267	!Config Help = choix de l'operateur de dissipation
268	!Config 'y' si on veut star et 'n' si on veut non-start !
269	!Config Moi y en a pas comprendre !
270	lstardis = .TRUE.
271	CALL getin('lstardis',lstardis)
272
273
274	!Config Key = nitergdiv
275	!Config Desc = Nombre d'iteration de gradiv
276	!Config Def = 1
277	!Config Help = nombre d'iterations de l'operateur de dissipation
278	!Config gradiv
279	nitergdiv = 1
280	CALL getin('nitergdiv',nitergdiv)
281
282	!Config Key = nitergrot
283	!Config Desc = nombre d'iterations de nxgradrot
284	!Config Def = 2
285	!Config Help = nombre d'iterations de l'operateur de dissipation
286	!Config nxgradrot
287	nitergrot = 2
288	CALL getin('nitergrot',nitergrot)
289
290
291	!Config Key = niterh
292	!Config Desc = nombre d'iterations de divgrad
293	!Config Def = 2
294	!Config Help = nombre d'iterations de l'operateur de dissipation
295	!Config divgrad
296	niterh = 2
297	CALL getin('niterh',niterh)
298
299
300	!Config Key = tetagdiv
301	!Config Desc = temps de dissipation pour div
302	!Config Def = 7200
303	!Config Help = temps de dissipation des plus petites longeur
304	!Config d'ondes pour u,v (gradiv)
305	tetagdiv = 7200.
306	CALL getin('tetagdiv',tetagdiv)
307
308	!Config Key = tetagrot
309	!Config Desc = temps de dissipation pour grad
310	!Config Def = 7200
311	!Config Help = temps de dissipation des plus petites longeur
312	!Config d'ondes pour u,v (nxgradrot)
313	tetagrot = 7200.
314	CALL getin('tetagrot',tetagrot)
315
316	!Config Key = tetatemp
317	!Config Desc = temps de dissipation pour h
318	!Config Def = 7200
319	!Config Help = temps de dissipation des plus petites longeur
320	!Config d'ondes pour h (divgrad)
321	tetatemp = 7200.
322	CALL getin('tetatemp',tetatemp )
323
324	! For Earth model only:
325	! Parametres controlant la variation sur la verticale des constantes de
326	! dissipation.
327	! Pour le moment actifs uniquement dans la version a 39 niveaux
328	! avec ok_strato=y
329
330	dissip_factz=4.
331	dissip_deltaz=10.
332	dissip_zref=30.
333	CALL getin('dissip_factz',dissip_factz )
334	CALL getin('dissip_deltaz',dissip_deltaz )
335	CALL getin('dissip_zref',dissip_zref )
336
337	! For other planets:
338	! Parametres controlant la variation sur la verticale des constantes de
339	! dissipation.
340	! Actifs uniquement avec ok_strato=y
341
342	dissip_fac_mid=2.
343	dissip_fac_up=10.
344	dissip_deltaz=10.! Intervalle (km) pour le changement mid / up
345	dissip_hdelta=5. ! scale height (km) dans la zone de la transition(m)
346	dissip_pupstart=1.e3 ! pression (Pa) au bas la transition mid / up
347	CALL getin('dissip_fac_mid',dissip_fac_mid )
348	CALL getin('dissip_fac_up',dissip_fac_up )
349	CALL getin('dissip_deltaz',dissip_deltaz )
350	CALL getin('dissip_hdelta',dissip_hdelta )
351	CALL getin('dissip_pupstart',dissip_pupstart )
352
353	! top_bound sponge: only active if iflag_top_bound!=0
354	! iflag_top_bound=0 for no sponge
355	! iflag_top_bound=1 for sponge over 4 topmost layers
356	! iflag_top_bound=2 for sponge from top to ~1% of top layer pressure
357	iflag_top_bound=0
358	CALL getin('iflag_top_bound',iflag_top_bound)
359
360	! mode_top_bound : fields towards which sponge relaxation will be done:
361	! mode_top_bound=0: no relaxation
362	! mode_top_bound=1: u and v relax towards 0
363	! mode_top_bound=2: u and v relax towards their zonal mean
364	! mode_top_bound=3: u,v and pot. temp. relax towards their zonal mean
365	mode_top_bound=3
366	CALL getin('mode_top_bound',mode_top_bound)
367
368	! top_bound sponge : inverse of charactericstic relaxation time scale for sponge
369	tau_top_bound=1.e-5
370	CALL getin('tau_top_bound',tau_top_bound)
371
372	! the other possible sponge layer (sponge_mod)
373	callsponge=.false. ! default value; don't use the sponge
374	call getin("callsponge",callsponge)
375	! check that user is not trying to use both sponge models
376	if ((iflag_top_bound.ge.1).and.callsponge) then
377	write(lunout,*)'Bad choice of options:'
378	write(lunout,*)' iflag_top_bound=',iflag_top_bound
379	write(lunout,*)' and callsponge=.true.'
380	write(lunout,*)'But both sponge models should not be',
381	& ' used simultaneously!'
382	stop
383	endif
384
385	! nsponge: number of atmospheric layers over which the sponge extends
386	nsponge=3 ! default value
387	call getin("nsponge",nsponge)
388
389	! mode_sponge: (quenching is towards ... over the upper nsponge layers)
390	! 0: (h=hmean,u=v=0)
391	! 1: (h=hmean,u=umean,v=0)
392	! 2: (h=hmean,u=umean,v=vmean)"
393	mode_sponge=2 ! default value
394	call getin("mode_sponge",mode_sponge)
395
396	! tetasponge: characteristic time scale (seconds) at topmost layer
397	! (time scale then doubles with decreasing layer index)."
398	tetasponge=50000.0
399	call getin("tetasponge",tetasponge)
400
401	! FOR TITAN: tidal forces
402	tidal=.TRUE.
403	CALL getin('tidal',tidal)
404
405	!Config Key = coefdis
406	!Config Desc = coefficient pour gamdissip
407	!Config Def = 0
408	!Config Help = coefficient pour gamdissip
409	coefdis = 0.
410	CALL getin('coefdis',coefdis)
411
412	!Config Key = purmats
413	!Config Desc = Schema d'integration
414	!Config Def = n
415	!Config Help = Choix du schema d'integration temporel.
416	!Config y = pure Matsuno sinon c'est du Matsuno-leapfrog
417	purmats = .FALSE.
418	CALL getin('purmats',purmats)
419
420	!Config Key = ok_guide
421	!Config Desc = Guidage
422	!Config Def = n
423	!Config Help = Guidage
424	ok_guide = .FALSE.
425	CALL getin('ok_guide',ok_guide)
426
427	c ...............................................................
428
429	!Config Key = read_start
430	!Config Desc = Initialize model using a 'start.nc' file
431	!Config Def = y
432	!Config Help = y: intialize dynamical fields using a 'start.nc' file
433	! n: fields are initialized by 'iniacademic' routine
434	read_start= .true.
435	CALL getin('read_start',read_start)
436
437	!Config Key = iflag_phys
438	!Config Desc = Avec ls physique
439	!Config Def = 1
440	!Config Help = Permet de faire tourner le modele sans
441	!Config physique.
442	iflag_phys = 1
443	CALL getin('iflag_phys',iflag_phys)
444
445
446	!Config Key = iphysiq
447	!Config Desc = Periode de la physique
448	!Config Def = 5
449	!Config Help = Periode de la physique en pas de temps de la dynamique.
450	iphysiq = 5
451	CALL getin('iphysiq', iphysiq)
452
453	!Config Key = iflag_trac
454	!Config Desc = traceurs presents ou non
455	!Config Def = 1
456	!Config Help = Permet de faire tourner le modele sans traceurs
457	!Config
458	iflag_trac = 1
459	CALL getin('iflag_trac',iflag_trac)
460
461	!Config Key = ip_ebil_dyn
462	!Config Desc = PRINT level for energy conserv. diag.
463	!Config Def = 0
464	!Config Help = PRINT level for energy conservation diag. ;
465	! les options suivantes existent :
466	!Config 0 pas de print
467	!Config 1 pas de print
468	!Config 2 print,
469	ip_ebil_dyn = 0
470	CALL getin('ip_ebil_dyn',ip_ebil_dyn)
471
472	!Config Key = offline
473	!Config Desc = Nouvelle eau liquide
474	!Config Def = n
475	!Config Help = Permet de mettre en route la
476	!Config nouvelle parametrisation de l'eau liquide !
477	offline = .FALSE.
478	CALL getin('offline',offline)
479
480	!Config Key = type_trac
481	!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS
482	!Config Def = lmdz
483	!Config Help =
484	!Config 'lmdz' = pas de couplage, pur LMDZ
485	!Config 'inca' = model de chime INCA
486	!Config 'repr' = model de chime REPROBUS
487	type_trac = 'lmdz'
488	CALL getin('type_trac',type_trac)
489
490	!Config Key = config_inca
491	!Config Desc = Choix de configuration de INCA
492	!Config Def = none
493	!Config Help = Choix de configuration de INCA :
494	!Config 'none' = sans INCA
495	!Config 'chem' = INCA avec calcul de chemie
496	!Config 'aero' = INCA avec calcul des aerosols
497	config_inca = 'none'
498	CALL getin('config_inca',config_inca)
499
500	!Config Key = ok_dynzon
501	!Config Desc = calcul et sortie des transports
502	!Config Def = n
503	!Config Help = Permet de mettre en route le calcul des transports
504	!Config
505	ok_dynzon = .FALSE.
506	CALL getin('ok_dynzon',ok_dynzon)
507
508	!Config Key = ok_dyn_ins
509	!Config Desc = sorties instantanees dans la dynamique
510	!Config Def = n
511	!Config Help =
512	!Config
513	ok_dyn_ins = .FALSE.
514	CALL getin('ok_dyn_ins',ok_dyn_ins)
515
516	!Config Key = ok_dyn_ave
517	!Config Desc = sorties moyennes dans la dynamique
518	!Config Def = n
519	!Config Help =
520	!Config
521	ok_dyn_ave = .FALSE.
522	CALL getin('ok_dyn_ave',ok_dyn_ave)
523
524	!Config Key = use_filtre_fft
525	!Config Desc = flag d'activation des FFT pour le filtre
526	!Config Def = false
527	!Config Help = permet d'activer l'utilisation des FFT pour effectuer
528	!Config le filtrage aux poles.
529	! Le filtre fft n'est pas implemente dans dyn3d
530	use_filtre_fft=.FALSE.
531	CALL getin('use_filtre_fft',use_filtre_fft)
532
533	IF (use_filtre_fft) THEN
534	write(lunout,*)'STOP !!!'
535	write(lunout,*)'use_filtre_fft n est pas implemente dans dyn3d'
536	STOP 1
537	ENDIF
538
539	!Config key = ok_strato
540	!Config Desc = activation de la version strato
541	!Config Def = .FALSE.
542	!Config Help = active la version stratosphérique de LMDZ de F. Lott
543
544	ok_strato=.TRUE.
545	CALL getin('ok_strato',ok_strato)
546
547	! NB: vert_prof_dissip is Earth-specific; should not impact other models
548	if (planet_type=="earth") then
549	vert_prof_dissip = merge(1, 0, ok_strato .and. llm==39)
550	CALL getin('vert_prof_dissip', vert_prof_dissip)
551	call assert(vert_prof_dissip == 0 .or. vert_prof_dissip == 1,
552	$ "bad value for vert_prof_dissip")
553	else
554	vert_prof_dissip=0 ! default for planets !
555	if (planet_type=="mars") then
556	vert_prof_dissip=1 ! use fac_mid & fac_up & startalt & delta
557	endif
558	endif
559
560	!Config Key = ok_gradsfile
561	!Config Desc = activation des sorties grads du guidage
562	!Config Def = n
563	!Config Help = active les sorties grads du guidage
564
565	ok_gradsfile = .FALSE.
566	CALL getin('ok_gradsfile',ok_gradsfile)
567
568	!Config Key = ok_limit
569	!Config Desc = creation des fichiers limit dans create_etat0_limit
570	!Config Def = y
571	!Config Help = production du fichier limit.nc requise
572
573	ok_limit = .TRUE.
574	CALL getin('ok_limit',ok_limit)
575
576	!Config Key = ok_etat0
577	!Config Desc = creation des fichiers etat0 dans create_etat0_limit
578	!Config Def = y
579	!Config Help = production des fichiers start.nc, startphy.nc requise
580
581	ok_etat0 = .TRUE.
582	CALL getin('ok_etat0',ok_etat0)
583
584	!Config Key = grilles_gcm_netcdf
585	!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit
586	!Config Def = n
587	grilles_gcm_netcdf = .FALSE.
588	CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf)
589
590	c----------------------------------------
591	c Parameters for zonal averages in the case of Titan
592	moyzon_mu = .false.
593	moyzon_ch = .false.
594	if (planet_type=="titan") then
595	CALL getin('moyzon_mu', moyzon_mu)
596	CALL getin('moyzon_ch', moyzon_ch)
597	endif
598	c----------------------------------------
599
600	c----------------------------------------
601	ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ...
602	c ......... ( modif le 17/04/96 ) .........
603	c
604	C ZOOM PARAMETERS ... the ones read in start.nc prevail anyway ! (SL, 2012)
605	c
606	c----------------------------------------
607	IF( etatinit ) then
608
609	!Config Key = clon
610	!Config Desc = centre du zoom, longitude
611	!Config Def = 0
612	!Config Help = longitude en degres du centre
613	!Config du zoom
614	clon = 0.
615	CALL getin('clon',clon)
616
617	!Config Key = clat
618	!Config Desc = centre du zoom, latitude
619	!Config Def = 0
620	!Config Help = latitude en degres du centre du zoom
621	!Config
622	clat = 0.
623	CALL getin('clat',clat)
624
625	!Config Key = grossismx
626	!Config Desc = zoom en longitude
627	!Config Def = 1.0
628	!Config Help = facteur de grossissement du zoom,
629	!Config selon la longitude
630	grossismx = 1.0
631	CALL getin('grossismx',grossismx)
632
633	!Config Key = grossismy
634	!Config Desc = zoom en latitude
635	!Config Def = 1.0
636	!Config Help = facteur de grossissement du zoom,
637	!Config selon la latitude
638	grossismy = 1.0
639	CALL getin('grossismy',grossismy)
640
641	IF( grossismx.LT.1. ) THEN
642	write(lunout,*)
643	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
644	STOP
645	ELSE
646	alphax = 1. - 1./ grossismx
647	ENDIF
648
649
650	IF( grossismy.LT.1. ) THEN
651	write(lunout,*)
652	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
653	STOP
654	ELSE
655	alphay = 1. - 1./ grossismy
656	ENDIF
657
658	write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
659	c
660	c alphax et alphay sont les anciennes formulat. des grossissements
661	c
662	c
663
664	!Config Key = fxyhypb
665	!Config Desc = Fonction hyperbolique
666	!Config Def = y
667	!Config Help = Fonction f(y) hyperbolique si = .true.
668	!Config sinon sinusoidale
669	fxyhypb = .TRUE.
670	CALL getin('fxyhypb',fxyhypb)
671
672	!Config Key = dzoomx
673	!Config Desc = extension en longitude
674	!Config Def = 0
675	!Config Help = extension en longitude de la zone du zoom
676	!Config ( fraction de la zone totale)
677	dzoomx = 0.0
678	CALL getin('dzoomx',dzoomx)
679
680	!Config Key = dzoomy
681	!Config Desc = extension en latitude
682	!Config Def = 0
683	!Config Help = extension en latitude de la zone du zoom
684	!Config ( fraction de la zone totale)
685	dzoomy = 0.0
686	CALL getin('dzoomy',dzoomy)
687
688	!Config Key = taux
689	!Config Desc = raideur du zoom en X
690	!Config Def = 3
691	!Config Help = raideur du zoom en X
692	taux = 3.0
693	CALL getin('taux',taux)
694
695	!Config Key = tauy
696	!Config Desc = raideur du zoom en Y
697	!Config Def = 3
698	!Config Help = raideur du zoom en Y
699	tauy = 3.0
700	CALL getin('tauy',tauy)
701
702	!Config Key = ysinus
703	!Config IF = !fxyhypb
704	!Config Desc = Fonction en Sinus
705	!Config Def = y
706	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
707	!Config sinon y = latit.
708	ysinus = .TRUE.
709	CALL getin('ysinus',ysinus)
710	c
711	c----------------------------------------
712	else ! etatinit=false
713	c----------------------------------------
714
715	!Config Key = clon
716	!Config Desc = centre du zoom, longitude
717	!Config Def = 0
718	!Config Help = longitude en degres du centre
719	!Config du zoom
720	clonn = 0.
721	CALL getin('clon',clonn)
722
723	!Config Key = clat
724	!Config Desc = centre du zoom, latitude
725	!Config Def = 0
726	!Config Help = latitude en degres du centre du zoom
727	!Config
728	clatt = 0.
729	CALL getin('clat',clatt)
730
731	c
732	c
733	IF( ABS(clat - clatt).GE. 0.001 ) THEN
734	write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
735	& ' est differente de celle lue sur le fichier start '
736	STOP
737	ENDIF
738
739	!Config Key = grossismx
740	!Config Desc = zoom en longitude
741	!Config Def = 1.0
742	!Config Help = facteur de grossissement du zoom,
743	!Config selon la longitude
744	grossismxx = 1.0
745	CALL getin('grossismx',grossismxx)
746
747
748	IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN
749	write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
750	& 'run.def est differente de celle lue sur le fichier start '
751	STOP
752	ENDIF
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	grossismyy = 1.0
760	CALL getin('grossismy',grossismyy)
761
762	IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN
763	write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
764	& 'run.def est differente de celle lue sur le fichier start '
765	STOP
766	ENDIF
767
768	IF( grossismx.LT.1. ) THEN
769	write(lunout,*)
770	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
771	STOP
772	ELSE
773	alphax = 1. - 1./ grossismx
774	ENDIF
775
776
777	IF( grossismy.LT.1. ) THEN
778	write(lunout,*)
779	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
780	STOP
781	ELSE
782	alphay = 1. - 1./ grossismy
783	ENDIF
784
785	write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
786	c
787	c alphax et alphay sont les anciennes formulat. des grossissements
788	c
789	c
790
791	!Config Key = fxyhypb
792	!Config Desc = Fonction hyperbolique
793	!Config Def = y
794	!Config Help = Fonction f(y) hyperbolique si = .true.
795	!Config sinon sinusoidale
796	fxyhypbb = .TRUE.
797	CALL getin('fxyhypb',fxyhypbb)
798
799	IF( .NOT.fxyhypb ) THEN
800	IF( fxyhypbb ) THEN
801	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
802	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
803	* 'F alors qu il est T sur run.def ***'
804	STOP
805	ENDIF
806	ELSE
807	IF( .NOT.fxyhypbb ) THEN
808	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
809	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
810	* 'T alors qu il est F sur run.def **** '
811	STOP
812	ENDIF
813	ENDIF
814	c
815	!Config Key = dzoomx
816	!Config Desc = extension en longitude
817	!Config Def = 0
818	!Config Help = extension en longitude de la zone du zoom
819	!Config ( fraction de la zone totale)
820	dzoomxx = 0.0
821	CALL getin('dzoomx',dzoomxx)
822
823	IF( fxyhypb ) THEN
824	IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN
825	write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
826	* 'run.def est differente de celle lue sur le fichier start '
827	STOP
828	ENDIF
829	ENDIF
830
831	!Config Key = dzoomy
832	!Config Desc = extension en latitude
833	!Config Def = 0
834	!Config Help = extension en latitude de la zone du zoom
835	!Config ( fraction de la zone totale)
836	dzoomyy = 0.0
837	CALL getin('dzoomy',dzoomyy)
838
839	IF( fxyhypb ) THEN
840	IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN
841	write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
842	* 'run.def est differente de celle lue sur le fichier start '
843	STOP
844	ENDIF
845	ENDIF
846
847	!Config Key = taux
848	!Config Desc = raideur du zoom en X
849	!Config Def = 3
850	!Config Help = raideur du zoom en X
851	tauxx = 3.0
852	CALL getin('taux',tauxx)
853
854	IF( fxyhypb ) THEN
855	IF( ABS(taux - tauxx).GE. 0.001 ) THEN
856	write(lunout,*)'conf_gcm: La valeur de taux passee par ',
857	* 'run.def est differente de celle lue sur le fichier start '
858	STOP
859	ENDIF
860	ENDIF
861
862	!Config Key = tauyy
863	!Config Desc = raideur du zoom en Y
864	!Config Def = 3
865	!Config Help = raideur du zoom en Y
866	tauyy = 3.0
867	CALL getin('tauy',tauyy)
868
869	IF( fxyhypb ) THEN
870	IF( ABS(tauy - tauyy).GE. 0.001 ) THEN
871	write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
872	* 'run.def est differente de celle lue sur le fichier start '
873	STOP
874	ENDIF
875	ENDIF
876
877	cc
878	IF( .NOT.fxyhypb ) THEN
879
880	!Config Key = ysinus
881	!Config IF = !fxyhypb
882	!Config Desc = Fonction en Sinus
883	!Config Def = y
884	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
885	!Config sinon y = latit.
886	ysinuss = .TRUE.
887	CALL getin('ysinus',ysinuss)
888
889	IF( .NOT.ysinus ) THEN
890	IF( ysinuss ) THEN
891	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
892	write(lunout,)' ** ysinus lu sur le fichier start est F',
893	* ' alors qu il est T sur run.def ***'
894	STOP
895	ENDIF
896	ELSE
897	IF( .NOT.ysinuss ) THEN
898	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
899	write(lunout,)' ** ysinus lu sur le fichier start est T',
900	* ' alors qu il est F sur run.def **** '
901	STOP
902	ENDIF
903	ENDIF
904	ENDIF ! of IF( .NOT.fxyhypb )
905
906	endif ! etatinit
907	c----------------------------------------
908
909
910	write(lunout,*)' #########################################'
911	write(lunout,*)' Configuration des parametres lus via run.def '
912	write(lunout,*)' planet_type = ', planet_type
913	write(lunout,*)' calend = ', calend
914	write(lunout,*)' dayref = ', dayref
915	write(lunout,*)' anneeref = ', anneeref
916	write(lunout,*)' nday = ', nday
917	if (ndynstep.ne.-9999) write(lunout,*)' ndynstep = ', ndynstep
918	if (less1day) write(lunout,*)' fractday = ', fractday
919	write(lunout,*)' day_step = ', day_step
920	write(lunout,*)' iperiod = ', iperiod
921	write(lunout,*)' nsplit_phys = ', nsplit_phys
922	write(lunout,*)' iconser = ', iconser
923	write(lunout,*)' iecri = ', iecri
924	write(lunout,*)' periodav = ', periodav
925	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
926	write(lunout,*)' dissip_period = ', dissip_period
927	write(lunout,*)' lstardis = ', lstardis
928	write(lunout,*)' nitergdiv = ', nitergdiv
929	write(lunout,*)' nitergrot = ', nitergrot
930	write(lunout,*)' niterh = ', niterh
931	write(lunout,*)' tetagdiv = ', tetagdiv
932	write(lunout,*)' tetagrot = ', tetagrot
933	write(lunout,*)' tetatemp = ', tetatemp
934	write(lunout,*)' coefdis = ', coefdis
935	write(lunout,*)' purmats = ', purmats
936	write(lunout,*)' read_start = ', read_start
937	write(lunout,*)' iflag_phys = ', iflag_phys
938	write(lunout,*)' iphysiq = ', iphysiq
939	write(lunout,*)' iflag_trac = ', iflag_trac
940	write(lunout,*)' clon = ', clon
941	write(lunout,*)' clat = ', clat
942	write(lunout,*)' grossismx = ', grossismx
943	write(lunout,*)' grossismy = ', grossismy
944	write(lunout,*)' fxyhypb = ', fxyhypb
945	write(lunout,*)' dzoomx = ', dzoomx
946	write(lunout,*)' dzoomy = ', dzoomy
947	write(lunout,*)' taux = ', taux
948	write(lunout,*)' tauy = ', tauy
949	write(lunout,*)' offline = ', offline
950	write(lunout,*)' type_trac = ', type_trac
951	write(lunout,*)' config_inca = ', config_inca
952	write(lunout,*)' ok_dynzon = ', ok_dynzon
953	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
954	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
955	write(lunout,*)' ok_strato = ', ok_strato
956	write(lunout,*)' ok_gradsfile = ', ok_gradsfile
957	write(lunout,*)' ok_limit = ', ok_limit
958	write(lunout,*)' ok_etat0 = ', ok_etat0
959	write(lunout,*)' grilles_gcm_netcdf = ', grilles_gcm_netcdf
960	if (planet_type=="titan") then
961	write(lunout,*)' moyzon_mu = ', moyzon_mu
962	write(lunout,*)' moyzon_ch = ', moyzon_ch
963	endif
964
965	RETURN
966	END

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

Download in other formats:

Original Format