Context Navigation

conf_gcm.F @ 2124

Last change on this file since 2124 was 2111, checked in by lguez, 10 years ago
Not possible to adjust MPI bands with nudging.
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 31.8 KB

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

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

Download in other formats:

Original Format