Context Navigation

conf_gcm.F @ 492

Last change on this file since 492 was 492, checked in by emillour, 14 years ago
Common dynamics: updates to keep up with LMDZ5 Earth (rev 1605) See file "DOC/chantiers/commit_importants.log" for details. EM
File size: 30.1 KB

Rev	Line
[1]	1	!
	2	! $Id: conf_gcm.F 1418 2010-07-19 15:11:24Z jghattas $
	3	!
	4	c
	5	c
[97]	6	SUBROUTINE conf_gcm( tapedef, etatinit )
[1]	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
[492]	15	USE infotrac, ONLY : type_trac
[1]	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 "temps.h"
	37	#include "comconst.h"
	38
	39	! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
	40	! #include "clesphys.h"
	41	#include "iniprint.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
[97]	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
[1]	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
[492]	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
[97]	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
[1]	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
[6]	191	!Config Desc = nombre de subdivisions par pas physique
[1]	192	!Config Def = 1
[6]	193	!Config Help = nombre de subdivisions par pas physique
[1]	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
[270]	241	!Config Key = dissip_period
[1]	242	!Config Desc = periode de la dissipation
[270]	243	!Config Def = 0
[1]	244	!Config Help = periode de la dissipation
[270]	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)
[1]	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.
[108]	315	! Actifs uniquement avec ok_strato=y
[1]	316
[108]	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 )
[1]	324	CALL getin('dissip_deltaz',dissip_deltaz )
[108]	325	CALL getin('dissip_hdelta',dissip_hdelta )
	326	CALL getin('dissip_pupstart',dissip_pupstart )
[1]	327
[108]	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
[1]	334	iflag_top_bound=1
[108]	335	mode_top_bound=3
[1]	336	tau_top_bound=1.e-5
	337	CALL getin('iflag_top_bound',iflag_top_bound)
[108]	338	CALL getin('mode_top_bound',mode_top_bound)
[1]	339	CALL getin('tau_top_bound',tau_top_bound)
	340
	341	!Config Key = coefdis
	342	!Config Desc = coefficient pour gamdissip
	343	!Config Def = 0
	344	!Config Help = coefficient pour gamdissip
	345	coefdis = 0.
	346	CALL getin('coefdis',coefdis)
	347
	348	!Config Key = purmats
	349	!Config Desc = Schema d'integration
	350	!Config Def = n
	351	!Config Help = Choix du schema d'integration temporel.
	352	!Config y = pure Matsuno sinon c'est du Matsuno-leapfrog
	353	purmats = .FALSE.
	354	CALL getin('purmats',purmats)
	355
	356	!Config Key = ok_guide
	357	!Config Desc = Guidage
	358	!Config Def = n
	359	!Config Help = Guidage
	360	ok_guide = .FALSE.
	361	CALL getin('ok_guide',ok_guide)
	362
	363	c ...............................................................
	364
	365	!Config Key = read_start
	366	!Config Desc = Initialize model using a 'start.nc' file
	367	!Config Def = y
	368	!Config Help = y: intialize dynamical fields using a 'start.nc' file
	369	! n: fields are initialized by 'iniacademic' routine
	370	read_start= .true.
	371	CALL getin('read_start',read_start)
	372
	373	!Config Key = iflag_phys
	374	!Config Desc = Avec ls physique
	375	!Config Def = 1
	376	!Config Help = Permet de faire tourner le modele sans
	377	!Config physique.
	378	iflag_phys = 1
	379	CALL getin('iflag_phys',iflag_phys)
	380
	381
	382	!Config Key = iphysiq
	383	!Config Desc = Periode de la physique
	384	!Config Def = 5
	385	!Config Help = Periode de la physique en pas de temps de la dynamique.
	386	iphysiq = 5
	387	CALL getin('iphysiq', iphysiq)
	388
[119]	389	!Config Key = iflag_trac
	390	!Config Desc = traceurs presents ou non
	391	!Config Def = 1
	392	!Config Help = Permet de faire tourner le modele sans traceurs
	393	!Config
	394	iflag_trac = 1
	395	CALL getin('iflag_trac',iflag_trac)
	396
[1]	397	!Config Key = ip_ebil_dyn
	398	!Config Desc = PRINT level for energy conserv. diag.
	399	!Config Def = 0
	400	!Config Help = PRINT level for energy conservation diag. ;
	401	! les options suivantes existent :
	402	!Config 0 pas de print
	403	!Config 1 pas de print
	404	!Config 2 print,
	405	ip_ebil_dyn = 0
	406	CALL getin('ip_ebil_dyn',ip_ebil_dyn)
	407
	408	ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ...
	409	c ......... ( modif le 17/04/96 ) .........
	410	c
	411	IF( etatinit ) GO TO 100
	412
	413	!Config Key = clon
	414	!Config Desc = centre du zoom, longitude
	415	!Config Def = 0
	416	!Config Help = longitude en degres du centre
	417	!Config du zoom
	418	clonn = 0.
	419	CALL getin('clon',clonn)
	420
	421	!Config Key = clat
	422	!Config Desc = centre du zoom, latitude
	423	!Config Def = 0
	424	!Config Help = latitude en degres du centre du zoom
	425	!Config
	426	clatt = 0.
	427	CALL getin('clat',clatt)
	428
	429	c
	430	c
	431	IF( ABS(clat - clatt).GE. 0.001 ) THEN
	432	write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
	433	& ' est differente de celle lue sur le fichier start '
	434	STOP
	435	ENDIF
	436
	437	!Config Key = grossismx
	438	!Config Desc = zoom en longitude
	439	!Config Def = 1.0
	440	!Config Help = facteur de grossissement du zoom,
	441	!Config selon la longitude
	442	grossismxx = 1.0
	443	CALL getin('grossismx',grossismxx)
	444
	445
	446	IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN
	447	write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
	448	& 'run.def est differente de celle lue sur le fichier start '
	449	STOP
	450	ENDIF
	451
	452	!Config Key = grossismy
	453	!Config Desc = zoom en latitude
	454	!Config Def = 1.0
	455	!Config Help = facteur de grossissement du zoom,
	456	!Config selon la latitude
	457	grossismyy = 1.0
	458	CALL getin('grossismy',grossismyy)
	459
	460	IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN
	461	write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
	462	& 'run.def est differente de celle lue sur le fichier start '
	463	STOP
	464	ENDIF
	465
	466	IF( grossismx.LT.1. ) THEN
	467	write(lunout,*)
	468	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
	469	STOP
	470	ELSE
	471	alphax = 1. - 1./ grossismx
	472	ENDIF
	473
	474
	475	IF( grossismy.LT.1. ) THEN
	476	write(lunout,*)
	477	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
	478	STOP
	479	ELSE
	480	alphay = 1. - 1./ grossismy
	481	ENDIF
	482
	483	write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
	484	c
	485	c alphax et alphay sont les anciennes formulat. des grossissements
	486	c
	487	c
	488
	489	!Config Key = fxyhypb
	490	!Config Desc = Fonction hyperbolique
	491	!Config Def = y
	492	!Config Help = Fonction f(y) hyperbolique si = .true.
	493	!Config sinon sinusoidale
	494	fxyhypbb = .TRUE.
	495	CALL getin('fxyhypb',fxyhypbb)
	496
	497	IF( .NOT.fxyhypb ) THEN
	498	IF( fxyhypbb ) THEN
	499	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
	500	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
	501	* 'F alors qu il est T sur run.def ***'
	502	STOP
	503	ENDIF
	504	ELSE
	505	IF( .NOT.fxyhypbb ) THEN
	506	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
	507	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
	508	* 'T alors qu il est F sur run.def **** '
	509	STOP
	510	ENDIF
	511	ENDIF
	512	c
	513	!Config Key = dzoomx
	514	!Config Desc = extension en longitude
	515	!Config Def = 0
	516	!Config Help = extension en longitude de la zone du zoom
	517	!Config ( fraction de la zone totale)
	518	dzoomxx = 0.0
	519	CALL getin('dzoomx',dzoomxx)
	520
	521	IF( fxyhypb ) THEN
	522	IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN
	523	write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
	524	* 'run.def est differente de celle lue sur le fichier start '
	525	STOP
	526	ENDIF
	527	ENDIF
	528
	529	!Config Key = dzoomy
	530	!Config Desc = extension en latitude
	531	!Config Def = 0
	532	!Config Help = extension en latitude de la zone du zoom
	533	!Config ( fraction de la zone totale)
	534	dzoomyy = 0.0
	535	CALL getin('dzoomy',dzoomyy)
	536
	537	IF( fxyhypb ) THEN
	538	IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN
	539	write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
	540	* 'run.def est differente de celle lue sur le fichier start '
	541	STOP
	542	ENDIF
	543	ENDIF
	544
	545	!Config Key = taux
	546	!Config Desc = raideur du zoom en X
	547	!Config Def = 3
	548	!Config Help = raideur du zoom en X
	549	tauxx = 3.0
	550	CALL getin('taux',tauxx)
	551
	552	IF( fxyhypb ) THEN
	553	IF( ABS(taux - tauxx).GE. 0.001 ) THEN
	554	write(lunout,*)'conf_gcm: La valeur de taux passee par ',
	555	* 'run.def est differente de celle lue sur le fichier start '
	556	STOP
	557	ENDIF
	558	ENDIF
	559
	560	!Config Key = tauyy
	561	!Config Desc = raideur du zoom en Y
	562	!Config Def = 3
	563	!Config Help = raideur du zoom en Y
	564	tauyy = 3.0
	565	CALL getin('tauy',tauyy)
	566
	567	IF( fxyhypb ) THEN
	568	IF( ABS(tauy - tauyy).GE. 0.001 ) THEN
	569	write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
	570	* 'run.def est differente de celle lue sur le fichier start '
	571	STOP
	572	ENDIF
	573	ENDIF
	574
	575	cc
	576	IF( .NOT.fxyhypb ) THEN
	577
	578	!Config Key = ysinus
	579	!Config IF = !fxyhypb
	580	!Config Desc = Fonction en Sinus
	581	!Config Def = y
	582	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
	583	!Config sinon y = latit.
	584	ysinuss = .TRUE.
	585	CALL getin('ysinus',ysinuss)
	586
	587	IF( .NOT.ysinus ) THEN
	588	IF( ysinuss ) THEN
	589	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
	590	write(lunout,)' ** ysinus lu sur le fichier start est F',
	591	* ' alors qu il est T sur run.def ***'
	592	STOP
	593	ENDIF
	594	ELSE
	595	IF( .NOT.ysinuss ) THEN
	596	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
	597	write(lunout,)' ** ysinus lu sur le fichier start est T',
	598	* ' alors qu il est F sur run.def **** '
	599	STOP
	600	ENDIF
	601	ENDIF
	602	ENDIF ! of IF( .NOT.fxyhypb )
	603	c
	604	!Config Key = offline
	605	!Config Desc = Nouvelle eau liquide
	606	!Config Def = n
	607	!Config Help = Permet de mettre en route la
	608	!Config nouvelle parametrisation de l'eau liquide !
	609	offline = .FALSE.
	610	CALL getin('offline',offline)
[492]	611
	612	!Config Key = type_trac
	613	!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS
	614	!Config Def = lmdz
	615	!Config Help =
	616	!Config 'lmdz' = pas de couplage, pur LMDZ
	617	!Config 'inca' = model de chime INCA
	618	!Config 'repr' = model de chime REPROBUS
	619	type_trac = 'lmdz'
	620	CALL getin('type_trac',type_trac)
[1]	621
	622	!Config Key = config_inca
	623	!Config Desc = Choix de configuration de INCA
	624	!Config Def = none
	625	!Config Help = Choix de configuration de INCA :
	626	!Config 'none' = sans INCA
	627	!Config 'chem' = INCA avec calcul de chemie
	628	!Config 'aero' = INCA avec calcul des aerosols
	629	config_inca = 'none'
	630	CALL getin('config_inca',config_inca)
	631
	632	!Config Key = ok_dynzon
	633	!Config Desc = calcul et sortie des transports
	634	!Config Def = n
	635	!Config Help = Permet de mettre en route le calcul des transports
	636	!Config
	637	ok_dynzon = .FALSE.
	638	CALL getin('ok_dynzon',ok_dynzon)
	639
	640	!Config Key = ok_dyn_ins
	641	!Config Desc = sorties instantanees dans la dynamique
	642	!Config Def = n
	643	!Config Help =
	644	!Config
	645	ok_dyn_ins = .FALSE.
	646	CALL getin('ok_dyn_ins',ok_dyn_ins)
	647
	648	!Config Key = ok_dyn_ave
	649	!Config Desc = sorties moyennes dans la dynamique
	650	!Config Def = n
	651	!Config Help =
	652	!Config
	653	ok_dyn_ave = .FALSE.
	654	CALL getin('ok_dyn_ave',ok_dyn_ave)
	655
	656
	657	write(lunout,*)' #########################################'
	658	write(lunout,*)' Configuration des parametres du gcm: '
	659	write(lunout,*)' planet_type = ', planet_type
	660	write(lunout,*)' calend = ', calend
	661	write(lunout,*)' dayref = ', dayref
	662	write(lunout,*)' anneeref = ', anneeref
	663	write(lunout,*)' nday = ', nday
[97]	664	if (less1day) then
	665	write(lunout,*)' Run only for a fraction of day ! '
	666	write(lunout,*)' fractday = ', fractday
	667	endif
[1]	668	write(lunout,*)' day_step = ', day_step
	669	write(lunout,*)' iperiod = ', iperiod
	670	write(lunout,*)' iconser = ', iconser
	671	write(lunout,*)' iecri = ', iecri
	672	write(lunout,*)' periodav = ', periodav
	673	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
[270]	674	write(lunout,*)' dissip_period = ', dissip_period
[1]	675	write(lunout,*)' lstardis = ', lstardis
	676	write(lunout,*)' nitergdiv = ', nitergdiv
	677	write(lunout,*)' nitergrot = ', nitergrot
	678	write(lunout,*)' niterh = ', niterh
	679	write(lunout,*)' tetagdiv = ', tetagdiv
	680	write(lunout,*)' tetagrot = ', tetagrot
	681	write(lunout,*)' tetatemp = ', tetatemp
	682	write(lunout,*)' coefdis = ', coefdis
	683	write(lunout,*)' purmats = ', purmats
	684	write(lunout,*)' read_start = ', read_start
	685	write(lunout,*)' iflag_phys = ', iflag_phys
	686	write(lunout,*)' iphysiq = ', iphysiq
[119]	687	write(lunout,*)' iflag_trac = ', iflag_trac
[1]	688	write(lunout,*)' clonn = ', clonn
	689	write(lunout,*)' clatt = ', clatt
	690	write(lunout,*)' grossismx = ', grossismx
	691	write(lunout,*)' grossismy = ', grossismy
	692	write(lunout,*)' fxyhypbb = ', fxyhypbb
	693	write(lunout,*)' dzoomxx = ', dzoomxx
	694	write(lunout,*)' dzoomy = ', dzoomyy
	695	write(lunout,*)' tauxx = ', tauxx
	696	write(lunout,*)' tauyy = ', tauyy
	697	write(lunout,*)' offline = ', offline
[492]	698	write(lunout,*)' type_trac = ', type_trac
[1]	699	write(lunout,*)' config_inca = ', config_inca
	700	write(lunout,*)' ok_dynzon = ', ok_dynzon
	701	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
	702	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
	703
	704	RETURN
	705	c ...............................................
	706	c
	707	100 CONTINUE
	708	!Config Key = clon
	709	!Config Desc = centre du zoom, longitude
	710	!Config Def = 0
	711	!Config Help = longitude en degres du centre
	712	!Config du zoom
	713	clon = 0.
	714	CALL getin('clon',clon)
	715
	716	!Config Key = clat
	717	!Config Desc = centre du zoom, latitude
	718	!Config Def = 0
	719	!Config Help = latitude en degres du centre du zoom
	720	!Config
	721	clat = 0.
	722	CALL getin('clat',clat)
	723
	724	!Config Key = grossismx
	725	!Config Desc = zoom en longitude
	726	!Config Def = 1.0
	727	!Config Help = facteur de grossissement du zoom,
	728	!Config selon la longitude
	729	grossismx = 1.0
	730	CALL getin('grossismx',grossismx)
	731
	732	!Config Key = grossismy
	733	!Config Desc = zoom en latitude
	734	!Config Def = 1.0
	735	!Config Help = facteur de grossissement du zoom,
	736	!Config selon la latitude
	737	grossismy = 1.0
	738	CALL getin('grossismy',grossismy)
	739
	740	IF( grossismx.LT.1. ) THEN
	741	write(lunout,*)
	742	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
	743	STOP
	744	ELSE
	745	alphax = 1. - 1./ grossismx
	746	ENDIF
	747
	748
	749	IF( grossismy.LT.1. ) THEN
	750	write(lunout,*)
	751	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
	752	STOP
	753	ELSE
	754	alphay = 1. - 1./ grossismy
	755	ENDIF
	756
	757	write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
	758	c
	759	c alphax et alphay sont les anciennes formulat. des grossissements
	760	c
	761	c
	762
	763	!Config Key = fxyhypb
	764	!Config Desc = Fonction hyperbolique
	765	!Config Def = y
	766	!Config Help = Fonction f(y) hyperbolique si = .true.
	767	!Config sinon sinusoidale
	768	fxyhypb = .TRUE.
	769	CALL getin('fxyhypb',fxyhypb)
	770
	771	!Config Key = dzoomx
	772	!Config Desc = extension en longitude
	773	!Config Def = 0
	774	!Config Help = extension en longitude de la zone du zoom
	775	!Config ( fraction de la zone totale)
	776	dzoomx = 0.0
	777	CALL getin('dzoomx',dzoomx)
	778
	779	!Config Key = dzoomy
	780	!Config Desc = extension en latitude
	781	!Config Def = 0
	782	!Config Help = extension en latitude de la zone du zoom
	783	!Config ( fraction de la zone totale)
	784	dzoomy = 0.0
	785	CALL getin('dzoomy',dzoomy)
	786
	787	!Config Key = taux
	788	!Config Desc = raideur du zoom en X
	789	!Config Def = 3
	790	!Config Help = raideur du zoom en X
	791	taux = 3.0
	792	CALL getin('taux',taux)
	793
	794	!Config Key = tauy
	795	!Config Desc = raideur du zoom en Y
	796	!Config Def = 3
	797	!Config Help = raideur du zoom en Y
	798	tauy = 3.0
	799	CALL getin('tauy',tauy)
	800
	801	!Config Key = ysinus
	802	!Config IF = !fxyhypb
	803	!Config Desc = Fonction en Sinus
	804	!Config Def = y
	805	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
	806	!Config sinon y = latit.
	807	ysinus = .TRUE.
	808	CALL getin('ysinus',ysinus)
	809	c
	810	!Config Key = offline
	811	!Config Desc = Nouvelle eau liquide
	812	!Config Def = n
	813	!Config Help = Permet de mettre en route la
	814	!Config nouvelle parametrisation de l'eau liquide !
	815	offline = .FALSE.
	816	CALL getin('offline',offline)
	817
[492]	818	!Config Key = type_trac
	819	!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS
	820	!Config Def = lmdz
	821	!Config Help =
	822	!Config 'lmdz' = pas de couplage, pur LMDZ
	823	!Config 'inca' = model de chime INCA
	824	!Config 'repr' = model de chime REPROBUS
	825	type_trac = 'lmdz'
	826	CALL getin('type_trac',type_trac)
	827
[1]	828	!Config Key = config_inca
	829	!Config Desc = Choix de configuration de INCA
	830	!Config Def = none
	831	!Config Help = Choix de configuration de INCA :
	832	!Config 'none' = sans INCA
	833	!Config 'chem' = INCA avec calcul de chemie
	834	!Config 'aero' = INCA avec calcul des aerosols
	835	config_inca = 'none'
	836	CALL getin('config_inca',config_inca)
	837
	838	!Config Key = ok_dynzon
	839	!Config Desc = sortie des transports zonaux dans la dynamique
	840	!Config Def = n
	841	!Config Help =
	842	!Config
	843	ok_dynzon = .FALSE.
	844	CALL getin('ok_dynzon',ok_dynzon)
	845
	846	!Config Key = ok_dyn_ins
	847	!Config Desc = sorties instantanees dans la dynamique
	848	!Config Def = n
	849	!Config Help =
	850	!Config
	851	ok_dyn_ins = .FALSE.
	852	CALL getin('ok_dyn_ins',ok_dyn_ins)
	853
	854	!Config Key = ok_dyn_ave
	855	!Config Desc = sorties moyennes dans la dynamique
	856	!Config Def = n
	857	!Config Help =
	858	!Config
	859	ok_dyn_ave = .FALSE.
	860	CALL getin('ok_dyn_ave',ok_dyn_ave)
	861
	862	!Config Key = use_filtre_fft
	863	!Config Desc = flag d'activation des FFT pour le filtre
	864	!Config Def = false
	865	!Config Help = permet d'activer l'utilisation des FFT pour effectuer
	866	!Config le filtrage aux poles.
	867	! Le filtre fft n'est pas implemente dans dyn3d
	868	use_filtre_fft=.FALSE.
	869	CALL getin('use_filtre_fft',use_filtre_fft)
	870
	871	IF (use_filtre_fft) THEN
	872	write(lunout,*)'STOP !!!'
	873	write(lunout,*)'use_filtre_fft n est pas implemente dans dyn3d'
	874	STOP
	875	ENDIF
	876
	877	!Config key = ok_strato
	878	!Config Desc = activation de la version strato
	879	!Config Def = .FALSE.
	880	!Config Help = active la version stratosphérique de LMDZ de F. Lott
	881
[6]	882	ok_strato=.TRUE.
[1]	883	CALL getin('ok_strato',ok_strato)
	884
	885	!Config Key = ok_gradsfile
	886	!Config Desc = activation des sorties grads du guidage
	887	!Config Def = n
	888	!Config Help = active les sorties grads du guidage
	889
	890	ok_gradsfile = .FALSE.
	891	CALL getin('ok_gradsfile',ok_gradsfile)
	892
	893	!Config Key = ok_limit
	894	!Config Desc = creation des fichiers limit dans create_etat0_limit
	895	!Config Def = y
	896	!Config Help = production du fichier limit.nc requise
	897
	898	ok_limit = .TRUE.
	899	CALL getin('ok_limit',ok_limit)
	900
	901	!Config Key = ok_etat0
	902	!Config Desc = creation des fichiers etat0 dans create_etat0_limit
	903	!Config Def = y
	904	!Config Help = production des fichiers start.nc, startphy.nc requise
	905
	906	ok_etat0 = .TRUE.
	907	CALL getin('ok_etat0',ok_etat0)
	908
[270]	909	!Config Key = grilles_gcm_netcdf
	910	!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit
	911	!Config Def = n
	912	grilles_gcm_netcdf = .FALSE.
	913	CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf)
	914
[1]	915	write(lunout,*)' #########################################'
	916	write(lunout,*)' Configuration des parametres de cel0'
	917	& //'_limit: '
	918	write(lunout,*)' planet_type = ', planet_type
	919	write(lunout,*)' calend = ', calend
	920	write(lunout,*)' dayref = ', dayref
	921	write(lunout,*)' anneeref = ', anneeref
	922	write(lunout,*)' nday = ', nday
	923	write(lunout,*)' day_step = ', day_step
	924	write(lunout,*)' iperiod = ', iperiod
	925	write(lunout,*)' iconser = ', iconser
	926	write(lunout,*)' iecri = ', iecri
	927	write(lunout,*)' periodav = ', periodav
	928	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
[270]	929	write(lunout,*)' dissip_period = ', dissip_period
[1]	930	write(lunout,*)' lstardis = ', lstardis
	931	write(lunout,*)' nitergdiv = ', nitergdiv
	932	write(lunout,*)' nitergrot = ', nitergrot
	933	write(lunout,*)' niterh = ', niterh
	934	write(lunout,*)' tetagdiv = ', tetagdiv
	935	write(lunout,*)' tetagrot = ', tetagrot
	936	write(lunout,*)' tetatemp = ', tetatemp
	937	write(lunout,*)' coefdis = ', coefdis
	938	write(lunout,*)' purmats = ', purmats
	939	write(lunout,*)' read_start = ', read_start
	940	write(lunout,*)' iflag_phys = ', iflag_phys
	941	write(lunout,*)' iphysiq = ', iphysiq
[119]	942	write(lunout,*)' iflag_trac = ', iflag_trac
[1]	943	write(lunout,*)' clon = ', clon
	944	write(lunout,*)' clat = ', clat
	945	write(lunout,*)' grossismx = ', grossismx
	946	write(lunout,*)' grossismy = ', grossismy
	947	write(lunout,*)' fxyhypb = ', fxyhypb
	948	write(lunout,*)' dzoomx = ', dzoomx
	949	write(lunout,*)' dzoomy = ', dzoomy
	950	write(lunout,*)' taux = ', taux
	951	write(lunout,*)' tauy = ', tauy
	952	write(lunout,*)' offline = ', offline
[492]	953	write(lunout,*)' type_trac = ', type_trac
[1]	954	write(lunout,*)' config_inca = ', config_inca
	955	write(lunout,*)' ok_dynzon = ', ok_dynzon
	956	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
	957	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
	958	write(lunout,*)' ok_strato = ', ok_strato
	959	write(lunout,*)' ok_gradsfile = ', ok_gradsfile
	960	write(lunout,*)' ok_limit = ', ok_limit
	961	write(lunout,*)' ok_etat0 = ', ok_etat0
[270]	962	write(lunout,*)' grilles_gcm_netcdf = ', grilles_gcm_netcdf
[1]	963	c
	964	RETURN
	965	END

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

Download in other formats:

Original Format