Context Navigation

conf_gcm.F @ 5451

Last change on this file since 5451 was 1438, checked in by jghattas, 14 years ago
Bug corrections for option offline. This option does not work at the same time as adjust=y.
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 29.3 KB

Rev	Line
[630]	1	!
[1279]	2	! $Id: conf_gcm.F 1438 2010-10-08 10:19:34Z fhourdin $
[630]	3	!
	4	c
	5	c
	6	SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
	7	c
[1146]	8	#ifdef CPP_IOIPSL
[630]	9	use IOIPSL
[1146]	10	#else
	11	! if not using IOIPSL, we still need to use (a local version of) getin
	12	use ioipsl_getincom
	13	#endif
[630]	14	use misc_mod
[985]	15	use mod_filtre_fft, ONLY : use_filtre_fft
	16	use mod_hallo, ONLY : use_mpi_alloc
	17	use parallel, ONLY : omp_chunk
[1403]	18	USE control_mod
[630]	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"
[985]	44	!#include "clesphys.h"
[630]	45	#include "iniprint.h"
[1279]	46	#include "temps.h"
	47	#include "comconst.h"
[985]	48
	49	! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
[630]	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-----------------------------------------------------------------------
[1146]	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)
[630]	124
[1279]	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
[630]	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
[1403]	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
[630]	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
[1146]	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
[630]	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
[1279]	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	!
[630]	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
[1146]	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
[630]	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
[764]	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	!
[630]	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
[1146]	401	write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
	402	& ' est differente de celle lue sur le fichier start '
[630]	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
[1146]	416	write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
	417	& 'run.def est differente de celle lue sur le fichier start '
[630]	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
[1146]	430	write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
	431	& 'run.def est differente de celle lue sur le fichier start '
[630]	432	STOP
	433	ENDIF
	434
	435	IF( grossismx.LT.1. ) THEN
[1146]	436	write(lunout,*)
	437	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
[630]	438	STOP
	439	ELSE
	440	alphax = 1. - 1./ grossismx
	441	ENDIF
	442
	443
	444	IF( grossismy.LT.1. ) THEN
[1146]	445	write(lunout,*)
	446	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
[630]	447	STOP
	448	ELSE
	449	alphay = 1. - 1./ grossismy
	450	ENDIF
	451
[1146]	452	write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
[630]	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
[1146]	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 ***'
[630]	471	STOP
[1146]	472	ENDIF
[630]	473	ELSE
[1146]	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 **** '
[630]	478	STOP
[1146]	479	ENDIF
[630]	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
[1146]	492	write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
	493	* 'run.def est differente de celle lue sur le fichier start '
[630]	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
[1146]	508	write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
	509	* 'run.def est differente de celle lue sur le fichier start '
[630]	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
[1146]	523	write(lunout,*)'conf_gcm: La valeur de taux passee par ',
	524	* 'run.def est differente de celle lue sur le fichier start '
[630]	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
[1146]	538	write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
	539	* 'run.def est differente de celle lue sur le fichier start '
[630]	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
[1146]	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
[630]	563	ELSE
[1146]	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 **** '
[630]	568	STOP
[1146]	569	ENDIF
[630]	570	ENDIF
[1146]	571	ENDIF ! of IF( .NOT.fxyhypb )
[630]	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)
[1438]	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
[960]	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)
[630]	598
[1146]	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)
[960]	606
[1403]	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)
[1146]	614
[1403]	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
[630]	623	write(lunout,*)' #########################################'
	624	write(lunout,*)' Configuration des parametres du gcm: '
[1146]	625	write(lunout,*)' planet_type = ', planet_type
[1279]	626	write(lunout,*)' calend = ', calend
[630]	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
[1403]	632	write(lunout,*)' nsplit_phys = ', nsplit_phys
[630]	633	write(lunout,*)' iconser = ', iconser
	634	write(lunout,*)' iecri = ', iecri
	635	write(lunout,*)' periodav = ', periodav
[1146]	636	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
[630]	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
[1146]	647	write(lunout,*)' read_start = ', read_start
[630]	648	write(lunout,*)' iflag_phys = ', iflag_phys
[1279]	649	write(lunout,*)' iphysiq = ', iphysiq
[630]	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
[960]	660	write(lunout,*)' config_inca = ', config_inca
[1146]	661	write(lunout,*)' ok_dynzon = ', ok_dynzon
[1403]	662	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
	663	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
[630]	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
[1146]	702	write(lunout,*)
	703	& 'conf_gcm: * ATTENTION !! grossismx < 1 . * '
[630]	704	STOP
	705	ELSE
	706	alphax = 1. - 1./ grossismx
	707	ENDIF
	708
	709
	710	IF( grossismy.LT.1. ) THEN
[1146]	711	write(lunout,*)
	712	& 'conf_gcm: * ATTENTION !! grossismy < 1 . * '
[630]	713	STOP
	714	ELSE
	715	alphay = 1. - 1./ grossismy
	716	ENDIF
	717
[1146]	718	write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
[630]	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)
[1438]	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
[630]	786
[960]	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)
[630]	796
[1146]	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
[1403]	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
[985]	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)
[1146]	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
[985]	836
[1146]	837
	838
[985]	839	!Config Key = use_mpi_alloc
[1146]	840	!Config Desc = Utilise un buffer MPI en mï¿½moire globale
[985]	841	!Config Def = false
	842	!Config Help = permet d'activer l'utilisation d'un buffer MPI
[1146]	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
[985]	845	!Config d'un facteur 2
	846	use_mpi_alloc=.FALSE.
	847	CALL getin('use_mpi_alloc',use_mpi_alloc)
[960]	848
[985]	849	!Config Key = omp_chunk
	850	!Config Desc = taille des blocs openmp
	851	!Config Def = 1
[1146]	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
[985]	855
	856	omp_chunk=1
	857	CALL getin('omp_chunk',omp_chunk)
	858
[1000]	859	!Config key = ok_strato
	860	!Config Desc = activation de la version strato
	861	!Config Def = .FALSE.
[1279]	862	!Config Help = active la version stratosphérique de LMDZ de F. Lott
[1000]	863
	864	ok_strato=.FALSE.
	865	CALL getin('ok_strato',ok_strato)
	866
[1046]	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
[1319]	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
[630]	891	write(lunout,*)' #########################################'
[1323]	892	write(lunout,*)' Configuration des parametres de cel0'
[1319]	893	& //'_limit: '
[1146]	894	write(lunout,*)' planet_type = ', planet_type
[1279]	895	write(lunout,*)' calend = ', calend
[630]	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
[1146]	904	write(lunout,*)' output_grads_dyn = ', output_grads_dyn
[630]	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
[1146]	915	write(lunout,*)' read_start = ', read_start
[630]	916	write(lunout,*)' iflag_phys = ', iflag_phys
[1279]	917	write(lunout,*)' iphysiq = ', iphysiq
	918	write(lunout,*)' clon = ', clon
[985]	919	write(lunout,*)' clat = ', clat
[630]	920	write(lunout,*)' grossismx = ', grossismx
	921	write(lunout,*)' grossismy = ', grossismy
[1146]	922	write(lunout,*)' fxyhypb = ', fxyhypb
[630]	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
[960]	928	write(lunout,*)' config_inca = ', config_inca
[1146]	929	write(lunout,*)' ok_dynzon = ', ok_dynzon
[1403]	930	write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
	931	write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
[985]	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
[1000]	935	write(lunout,*)' ok_strato = ', ok_strato
[1046]	936	write(lunout,*)' ok_gradsfile = ', ok_gradsfile
[1319]	937	write(lunout,*)' ok_limit = ', ok_limit
	938	write(lunout,*)' ok_etat0 = ', ok_etat0
[630]	939	c
	940	RETURN
	941	END

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

Download in other formats:

Original Format