Context Navigation

conf_gcm.F @ 1893

Last change on this file since 1893 was 1864, checked in by Laurent Fairhead, 11 years ago

Création d'une nouvelle testing:

merge des modifications du trunk entre r1796 et r1860

New testing version

merged modifications between r1796 and r1860 from the trunk

i.e.
svn merge -r1796:1860 http://svn.lmd.jussieu.fr/LMDZ/LMDZ5/trunk

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

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

Download in other formats:

Original Format