Context Navigation

conf_gcm.F @ 1243

Last change on this file since 1243 was 1189, checked in by emillour, 11 years ago

Common dynamics: a couple of bug fixes:

Correctly account for the change in pressure, mass, etc. after modifying surface pressure following a call to the physics.
Corrected tracer advection, which is computed using values at the beginning of the time step, so it is done at Matsuno forward step.

File size: 33.0 KB

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

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

Download in other formats:

Original Format