Context Navigation

conf_gcm.F @ 1664

Last change on this file since 1664 was 1664, checked in by Laurent Fairhead, 12 years ago

Version testing basée sur la r1575

Testing release based on r1575

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 30.5 KB

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

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

Download in other formats:

Original Format