Context Navigation

conf_gcm.F @ 7

Last change on this file since 7 was 7, checked in by emillour, 14 years ago
Mise a niveau de la dynamique par rapport a la version 1447 de LMDZ5 Voir les details dans chantiers/commit_v7.log
File size: 29.3 KB

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

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

Download in other formats:

Original Format