Context Navigation

conf_gcm.F @ 1384

Last change on this file since 1384 was 1352, checked in by emillour, 11 years ago

Common dynamics:

Added MESU arch files
Bug fix in conf_gcm to enable using FFT filter

File size: 32.7 KB

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

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

Download in other formats:

Original Format