Context Navigation

conf_gcm.F @ 1865

Last change on this file since 1865 was 1795, checked in by Ehouarn Millour, 12 years ago

Version testing basee sur la r1794

Testing release based on r1794

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

File size: 29.7 KB

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

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

Download in other formats:

Original Format