Context Navigation

conf_gcm.F @ 612

Last change on this file since 612 was 495, checked in by slebonnois, 14 years ago
Mise a jour physique Titan, ajout des forces de marees (dans la dynamique, sous flag titan). SL.
File size: 30.2 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
[495]	341	! FOR TITAN: tidal forces
	342	tidal=.TRUE.
	343	CALL getin('tidal',tidal)
	344
[1]	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
[119]	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
[1]	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	ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ...
	413	c ......... ( modif le 17/04/96 ) .........
	414	c
	415	IF( etatinit ) GO TO 100
	416
	417	!Config Key = clon
	418	!Config Desc = centre du zoom, longitude
	419	!Config Def = 0
	420	!Config Help = longitude en degres du centre
	421	!Config du zoom
	422	clonn = 0.
	423	CALL getin('clon',clonn)
	424
	425	!Config Key = clat
	426	!Config Desc = centre du zoom, latitude
	427	!Config Def = 0
	428	!Config Help = latitude en degres du centre du zoom
	429	!Config
	430	clatt = 0.
	431	CALL getin('clat',clatt)
	432
	433	c
	434	c
	435	IF( ABS(clat - clatt).GE. 0.001 ) THEN
	436	write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
	437	& ' est differente de celle lue sur le fichier start '
	438	STOP
	439	ENDIF
	440
	441	!Config Key = grossismx
	442	!Config Desc = zoom en longitude
	443	!Config Def = 1.0
	444	!Config Help = facteur de grossissement du zoom,
	445	!Config selon la longitude
	446	grossismxx = 1.0
	447	CALL getin('grossismx',grossismxx)
	448
	449
	450	IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN
	451	write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
	452	& 'run.def est differente de celle lue sur le fichier start '
	453	STOP
	454	ENDIF
	455
	456	!Config Key = grossismy
	457	!Config Desc = zoom en latitude
	458	!Config Def = 1.0
	459	!Config Help = facteur de grossissement du zoom,
	460	!Config selon la latitude
	461	grossismyy = 1.0
	462	CALL getin('grossismy',grossismyy)
	463
	464	IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN
	465	write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
	466	& 'run.def est differente de celle lue sur le fichier start '
	467	STOP
	468	ENDIF
	469
	470	IF( grossismx.LT.1. ) THEN
	471	write(lunout,*)
	472	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
	473	STOP
	474	ELSE
	475	alphax = 1. - 1./ grossismx
	476	ENDIF
	477
	478
	479	IF( grossismy.LT.1. ) THEN
	480	write(lunout,*)
	481	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
	482	STOP
	483	ELSE
	484	alphay = 1. - 1./ grossismy
	485	ENDIF
	486
	487	write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
	488	c
	489	c alphax et alphay sont les anciennes formulat. des grossissements
	490	c
	491	c
	492
	493	!Config Key = fxyhypb
	494	!Config Desc = Fonction hyperbolique
	495	!Config Def = y
	496	!Config Help = Fonction f(y) hyperbolique si = .true.
	497	!Config sinon sinusoidale
	498	fxyhypbb = .TRUE.
	499	CALL getin('fxyhypb',fxyhypbb)
	500
	501	IF( .NOT.fxyhypb ) THEN
	502	IF( fxyhypbb ) THEN
	503	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
	504	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
	505	* 'F alors qu il est T sur run.def ***'
	506	STOP
	507	ENDIF
	508	ELSE
	509	IF( .NOT.fxyhypbb ) THEN
	510	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
	511	write(lunout,)' ** fxyhypb lu sur le fichier start est ',
	512	* 'T alors qu il est F sur run.def **** '
	513	STOP
	514	ENDIF
	515	ENDIF
	516	c
	517	!Config Key = dzoomx
	518	!Config Desc = extension en longitude
	519	!Config Def = 0
	520	!Config Help = extension en longitude de la zone du zoom
	521	!Config ( fraction de la zone totale)
	522	dzoomxx = 0.0
	523	CALL getin('dzoomx',dzoomxx)
	524
	525	IF( fxyhypb ) THEN
	526	IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN
	527	write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
	528	* 'run.def est differente de celle lue sur le fichier start '
	529	STOP
	530	ENDIF
	531	ENDIF
	532
	533	!Config Key = dzoomy
	534	!Config Desc = extension en latitude
	535	!Config Def = 0
	536	!Config Help = extension en latitude de la zone du zoom
	537	!Config ( fraction de la zone totale)
	538	dzoomyy = 0.0
	539	CALL getin('dzoomy',dzoomyy)
	540
	541	IF( fxyhypb ) THEN
	542	IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN
	543	write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
	544	* 'run.def est differente de celle lue sur le fichier start '
	545	STOP
	546	ENDIF
	547	ENDIF
	548
	549	!Config Key = taux
	550	!Config Desc = raideur du zoom en X
	551	!Config Def = 3
	552	!Config Help = raideur du zoom en X
	553	tauxx = 3.0
	554	CALL getin('taux',tauxx)
	555
	556	IF( fxyhypb ) THEN
	557	IF( ABS(taux - tauxx).GE. 0.001 ) THEN
	558	write(lunout,*)'conf_gcm: La valeur de taux passee par ',
	559	* 'run.def est differente de celle lue sur le fichier start '
	560	STOP
	561	ENDIF
	562	ENDIF
	563
	564	!Config Key = tauyy
	565	!Config Desc = raideur du zoom en Y
	566	!Config Def = 3
	567	!Config Help = raideur du zoom en Y
	568	tauyy = 3.0
	569	CALL getin('tauy',tauyy)
	570
	571	IF( fxyhypb ) THEN
	572	IF( ABS(tauy - tauyy).GE. 0.001 ) THEN
	573	write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
	574	* 'run.def est differente de celle lue sur le fichier start '
	575	STOP
	576	ENDIF
	577	ENDIF
	578
	579	cc
	580	IF( .NOT.fxyhypb ) THEN
	581
	582	!Config Key = ysinus
	583	!Config IF = !fxyhypb
	584	!Config Desc = Fonction en Sinus
	585	!Config Def = y
	586	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
	587	!Config sinon y = latit.
	588	ysinuss = .TRUE.
	589	CALL getin('ysinus',ysinuss)
	590
	591	IF( .NOT.ysinus ) THEN
	592	IF( ysinuss ) THEN
	593	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
	594	write(lunout,)' ** ysinus lu sur le fichier start est F',
	595	* ' alors qu il est T sur run.def ***'
	596	STOP
	597	ENDIF
	598	ELSE
	599	IF( .NOT.ysinuss ) THEN
	600	write(lunout,)' ***** PBS DANS CONF_GCM ****** '
	601	write(lunout,)' ** ysinus lu sur le fichier start est T',
	602	* ' alors qu il est F sur run.def **** '
	603	STOP
	604	ENDIF
	605	ENDIF
	606	ENDIF ! of IF( .NOT.fxyhypb )
	607	c
	608	!Config Key = offline
	609	!Config Desc = Nouvelle eau liquide
	610	!Config Def = n
	611	!Config Help = Permet de mettre en route la
	612	!Config nouvelle parametrisation de l'eau liquide !
	613	offline = .FALSE.
	614	CALL getin('offline',offline)
[492]	615
	616	!Config Key = type_trac
	617	!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS
	618	!Config Def = lmdz
	619	!Config Help =
	620	!Config 'lmdz' = pas de couplage, pur LMDZ
	621	!Config 'inca' = model de chime INCA
	622	!Config 'repr' = model de chime REPROBUS
	623	type_trac = 'lmdz'
	624	CALL getin('type_trac',type_trac)
[1]	625
	626	!Config Key = config_inca
	627	!Config Desc = Choix de configuration de INCA
	628	!Config Def = none
	629	!Config Help = Choix de configuration de INCA :
	630	!Config 'none' = sans INCA
	631	!Config 'chem' = INCA avec calcul de chemie
	632	!Config 'aero' = INCA avec calcul des aerosols
	633	config_inca = 'none'
	634	CALL getin('config_inca',config_inca)
	635
	636	!Config Key = ok_dynzon
	637	!Config Desc = calcul et sortie des transports
	638	!Config Def = n
	639	!Config Help = Permet de mettre en route le calcul des transports
	640	!Config
	641	ok_dynzon = .FALSE.
	642	CALL getin('ok_dynzon',ok_dynzon)
	643
	644	!Config Key = ok_dyn_ins
	645	!Config Desc = sorties instantanees dans la dynamique
	646	!Config Def = n
	647	!Config Help =
	648	!Config
	649	ok_dyn_ins = .FALSE.
	650	CALL getin('ok_dyn_ins',ok_dyn_ins)
	651
	652	!Config Key = ok_dyn_ave
	653	!Config Desc = sorties moyennes dans la dynamique
	654	!Config Def = n
	655	!Config Help =
	656	!Config
	657	ok_dyn_ave = .FALSE.
	658	CALL getin('ok_dyn_ave',ok_dyn_ave)
	659
	660
	661	write(lunout,*)' #########################################'
	662	write(lunout,*)' Configuration des parametres du gcm: '
	663	write(lunout,*)' planet_type = ', planet_type
	664	write(lunout,*)' calend = ', calend
	665	write(lunout,*)' dayref = ', dayref
	666	write(lunout,*)' anneeref = ', anneeref
	667	write(lunout,*)' nday = ', nday
[97]	668	if (less1day) then
	669	write(lunout,*)' Run only for a fraction of day ! '
	670	write(lunout,*)' fractday = ', fractday
	671	endif
[1]	672	write(lunout,*)' day_step = ', day_step
	673	write(lunout,*)' iperiod = ', iperiod
	674	write(lunout,*)' iconser = ', iconser
	675	write(lunout,*)' iecri = ', iecri
	676	write(lunout,*)' periodav = ', periodav
	677	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
[270]	678	write(lunout,*)' dissip_period = ', dissip_period
[1]	679	write(lunout,*)' lstardis = ', lstardis
	680	write(lunout,*)' nitergdiv = ', nitergdiv
	681	write(lunout,*)' nitergrot = ', nitergrot
	682	write(lunout,*)' niterh = ', niterh
	683	write(lunout,*)' tetagdiv = ', tetagdiv
	684	write(lunout,*)' tetagrot = ', tetagrot
	685	write(lunout,*)' tetatemp = ', tetatemp
	686	write(lunout,*)' coefdis = ', coefdis
	687	write(lunout,*)' purmats = ', purmats
	688	write(lunout,*)' read_start = ', read_start
	689	write(lunout,*)' iflag_phys = ', iflag_phys
	690	write(lunout,*)' iphysiq = ', iphysiq
[119]	691	write(lunout,*)' iflag_trac = ', iflag_trac
[1]	692	write(lunout,*)' clonn = ', clonn
	693	write(lunout,*)' clatt = ', clatt
	694	write(lunout,*)' grossismx = ', grossismx
	695	write(lunout,*)' grossismy = ', grossismy
	696	write(lunout,*)' fxyhypbb = ', fxyhypbb
	697	write(lunout,*)' dzoomxx = ', dzoomxx
	698	write(lunout,*)' dzoomy = ', dzoomyy
	699	write(lunout,*)' tauxx = ', tauxx
	700	write(lunout,*)' tauyy = ', tauyy
	701	write(lunout,*)' offline = ', offline
[492]	702	write(lunout,*)' type_trac = ', type_trac
[1]	703	write(lunout,*)' config_inca = ', config_inca
	704	write(lunout,*)' ok_dynzon = ', ok_dynzon
	705	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
	706	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
	707
	708	RETURN
	709	c ...............................................
	710	c
	711	100 CONTINUE
	712	!Config Key = clon
	713	!Config Desc = centre du zoom, longitude
	714	!Config Def = 0
	715	!Config Help = longitude en degres du centre
	716	!Config du zoom
	717	clon = 0.
	718	CALL getin('clon',clon)
	719
	720	!Config Key = clat
	721	!Config Desc = centre du zoom, latitude
	722	!Config Def = 0
	723	!Config Help = latitude en degres du centre du zoom
	724	!Config
	725	clat = 0.
	726	CALL getin('clat',clat)
	727
	728	!Config Key = grossismx
	729	!Config Desc = zoom en longitude
	730	!Config Def = 1.0
	731	!Config Help = facteur de grossissement du zoom,
	732	!Config selon la longitude
	733	grossismx = 1.0
	734	CALL getin('grossismx',grossismx)
	735
	736	!Config Key = grossismy
	737	!Config Desc = zoom en latitude
	738	!Config Def = 1.0
	739	!Config Help = facteur de grossissement du zoom,
	740	!Config selon la latitude
	741	grossismy = 1.0
	742	CALL getin('grossismy',grossismy)
	743
	744	IF( grossismx.LT.1. ) THEN
	745	write(lunout,*)
	746	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
	747	STOP
	748	ELSE
	749	alphax = 1. - 1./ grossismx
	750	ENDIF
	751
	752
	753	IF( grossismy.LT.1. ) THEN
	754	write(lunout,*)
	755	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
	756	STOP
	757	ELSE
	758	alphay = 1. - 1./ grossismy
	759	ENDIF
	760
	761	write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
	762	c
	763	c alphax et alphay sont les anciennes formulat. des grossissements
	764	c
	765	c
	766
	767	!Config Key = fxyhypb
	768	!Config Desc = Fonction hyperbolique
	769	!Config Def = y
	770	!Config Help = Fonction f(y) hyperbolique si = .true.
	771	!Config sinon sinusoidale
	772	fxyhypb = .TRUE.
	773	CALL getin('fxyhypb',fxyhypb)
	774
	775	!Config Key = dzoomx
	776	!Config Desc = extension en longitude
	777	!Config Def = 0
	778	!Config Help = extension en longitude de la zone du zoom
	779	!Config ( fraction de la zone totale)
	780	dzoomx = 0.0
	781	CALL getin('dzoomx',dzoomx)
	782
	783	!Config Key = dzoomy
	784	!Config Desc = extension en latitude
	785	!Config Def = 0
	786	!Config Help = extension en latitude de la zone du zoom
	787	!Config ( fraction de la zone totale)
	788	dzoomy = 0.0
	789	CALL getin('dzoomy',dzoomy)
	790
	791	!Config Key = taux
	792	!Config Desc = raideur du zoom en X
	793	!Config Def = 3
	794	!Config Help = raideur du zoom en X
	795	taux = 3.0
	796	CALL getin('taux',taux)
	797
	798	!Config Key = tauy
	799	!Config Desc = raideur du zoom en Y
	800	!Config Def = 3
	801	!Config Help = raideur du zoom en Y
	802	tauy = 3.0
	803	CALL getin('tauy',tauy)
	804
	805	!Config Key = ysinus
	806	!Config IF = !fxyhypb
	807	!Config Desc = Fonction en Sinus
	808	!Config Def = y
	809	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
	810	!Config sinon y = latit.
	811	ysinus = .TRUE.
	812	CALL getin('ysinus',ysinus)
	813	c
	814	!Config Key = offline
	815	!Config Desc = Nouvelle eau liquide
	816	!Config Def = n
	817	!Config Help = Permet de mettre en route la
	818	!Config nouvelle parametrisation de l'eau liquide !
	819	offline = .FALSE.
	820	CALL getin('offline',offline)
	821
[492]	822	!Config Key = type_trac
	823	!Config Desc = Choix de couplage avec model de chimie INCA ou REPROBUS
	824	!Config Def = lmdz
	825	!Config Help =
	826	!Config 'lmdz' = pas de couplage, pur LMDZ
	827	!Config 'inca' = model de chime INCA
	828	!Config 'repr' = model de chime REPROBUS
	829	type_trac = 'lmdz'
	830	CALL getin('type_trac',type_trac)
	831
[1]	832	!Config Key = config_inca
	833	!Config Desc = Choix de configuration de INCA
	834	!Config Def = none
	835	!Config Help = Choix de configuration de INCA :
	836	!Config 'none' = sans INCA
	837	!Config 'chem' = INCA avec calcul de chemie
	838	!Config 'aero' = INCA avec calcul des aerosols
	839	config_inca = 'none'
	840	CALL getin('config_inca',config_inca)
	841
	842	!Config Key = ok_dynzon
	843	!Config Desc = sortie des transports zonaux dans la dynamique
	844	!Config Def = n
	845	!Config Help =
	846	!Config
	847	ok_dynzon = .FALSE.
	848	CALL getin('ok_dynzon',ok_dynzon)
	849
	850	!Config Key = ok_dyn_ins
	851	!Config Desc = sorties instantanees dans la dynamique
	852	!Config Def = n
	853	!Config Help =
	854	!Config
	855	ok_dyn_ins = .FALSE.
	856	CALL getin('ok_dyn_ins',ok_dyn_ins)
	857
	858	!Config Key = ok_dyn_ave
	859	!Config Desc = sorties moyennes dans la dynamique
	860	!Config Def = n
	861	!Config Help =
	862	!Config
	863	ok_dyn_ave = .FALSE.
	864	CALL getin('ok_dyn_ave',ok_dyn_ave)
	865
	866	!Config Key = use_filtre_fft
	867	!Config Desc = flag d'activation des FFT pour le filtre
	868	!Config Def = false
	869	!Config Help = permet d'activer l'utilisation des FFT pour effectuer
	870	!Config le filtrage aux poles.
	871	! Le filtre fft n'est pas implemente dans dyn3d
	872	use_filtre_fft=.FALSE.
	873	CALL getin('use_filtre_fft',use_filtre_fft)
	874
	875	IF (use_filtre_fft) THEN
	876	write(lunout,*)'STOP !!!'
	877	write(lunout,*)'use_filtre_fft n est pas implemente dans dyn3d'
	878	STOP
	879	ENDIF
	880
	881	!Config key = ok_strato
	882	!Config Desc = activation de la version strato
	883	!Config Def = .FALSE.
	884	!Config Help = active la version stratosphérique de LMDZ de F. Lott
	885
[6]	886	ok_strato=.TRUE.
[1]	887	CALL getin('ok_strato',ok_strato)
	888
	889	!Config Key = ok_gradsfile
	890	!Config Desc = activation des sorties grads du guidage
	891	!Config Def = n
	892	!Config Help = active les sorties grads du guidage
	893
	894	ok_gradsfile = .FALSE.
	895	CALL getin('ok_gradsfile',ok_gradsfile)
	896
	897	!Config Key = ok_limit
	898	!Config Desc = creation des fichiers limit dans create_etat0_limit
	899	!Config Def = y
	900	!Config Help = production du fichier limit.nc requise
	901
	902	ok_limit = .TRUE.
	903	CALL getin('ok_limit',ok_limit)
	904
	905	!Config Key = ok_etat0
	906	!Config Desc = creation des fichiers etat0 dans create_etat0_limit
	907	!Config Def = y
	908	!Config Help = production des fichiers start.nc, startphy.nc requise
	909
	910	ok_etat0 = .TRUE.
	911	CALL getin('ok_etat0',ok_etat0)
	912
[270]	913	!Config Key = grilles_gcm_netcdf
	914	!Config Desc = creation de fichier grilles_gcm.nc dans create_etat0_limit
	915	!Config Def = n
	916	grilles_gcm_netcdf = .FALSE.
	917	CALL getin('grilles_gcm_netcdf',grilles_gcm_netcdf)
	918
[1]	919	write(lunout,*)' #########################################'
	920	write(lunout,*)' Configuration des parametres de cel0'
	921	& //'_limit: '
	922	write(lunout,*)' planet_type = ', planet_type
	923	write(lunout,*)' calend = ', calend
	924	write(lunout,*)' dayref = ', dayref
	925	write(lunout,*)' anneeref = ', anneeref
	926	write(lunout,*)' nday = ', nday
	927	write(lunout,*)' day_step = ', day_step
	928	write(lunout,*)' iperiod = ', iperiod
	929	write(lunout,*)' iconser = ', iconser
	930	write(lunout,*)' iecri = ', iecri
	931	write(lunout,*)' periodav = ', periodav
	932	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
[270]	933	write(lunout,*)' dissip_period = ', dissip_period
[1]	934	write(lunout,*)' lstardis = ', lstardis
	935	write(lunout,*)' nitergdiv = ', nitergdiv
	936	write(lunout,*)' nitergrot = ', nitergrot
	937	write(lunout,*)' niterh = ', niterh
	938	write(lunout,*)' tetagdiv = ', tetagdiv
	939	write(lunout,*)' tetagrot = ', tetagrot
	940	write(lunout,*)' tetatemp = ', tetatemp
	941	write(lunout,*)' coefdis = ', coefdis
	942	write(lunout,*)' purmats = ', purmats
	943	write(lunout,*)' read_start = ', read_start
	944	write(lunout,*)' iflag_phys = ', iflag_phys
	945	write(lunout,*)' iphysiq = ', iphysiq
[119]	946	write(lunout,*)' iflag_trac = ', iflag_trac
[1]	947	write(lunout,*)' clon = ', clon
	948	write(lunout,*)' clat = ', clat
	949	write(lunout,*)' grossismx = ', grossismx
	950	write(lunout,*)' grossismy = ', grossismy
	951	write(lunout,*)' fxyhypb = ', fxyhypb
	952	write(lunout,*)' dzoomx = ', dzoomx
	953	write(lunout,*)' dzoomy = ', dzoomy
	954	write(lunout,*)' taux = ', taux
	955	write(lunout,*)' tauy = ', tauy
	956	write(lunout,*)' offline = ', offline
[492]	957	write(lunout,*)' type_trac = ', type_trac
[1]	958	write(lunout,*)' config_inca = ', config_inca
	959	write(lunout,*)' ok_dynzon = ', ok_dynzon
	960	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
	961	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
	962	write(lunout,*)' ok_strato = ', ok_strato
	963	write(lunout,*)' ok_gradsfile = ', ok_gradsfile
	964	write(lunout,*)' ok_limit = ', ok_limit
	965	write(lunout,*)' ok_etat0 = ', ok_etat0
[270]	966	write(lunout,*)' grilles_gcm_netcdf = ', grilles_gcm_netcdf
[1]	967	c
	968	RETURN
	969	END

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

Download in other formats:

Original Format