Context Navigation

conf_gcm.F @ 377

Last change on this file since 377 was 351, checked in by lmdzadmin, 23 years ago
Changement de la date de depart par defaut (1/01 plutot que 1/07) LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 21.8 KB

Rev	Line
[237]	1	c
	2	c $Header$
	3	c
	4	SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
	5	c
	6	use IOIPSL
	7	IMPLICIT NONE
	8	c-----------------------------------------------------------------------
	9	c Auteurs : L. Fairhead , P. Le Van .
	10	c
	11	c Arguments :
	12	c
	13	c tapedef :
	14	c etatinit : = TRUE , on ne compare pas les valeurs des para-
	15	c -metres du zoom avec celles lues sur le fichier start .
	16	c clesphy0 : sortie .
	17	c
	18	LOGICAL etatinit
	19	INTEGER tapedef
	20
	21	INTEGER longcles
	22	PARAMETER( longcles = 20 )
	23	REAL clesphy0( longcles )
	24	c
	25	c Declarations :
	26	c --------------
	27	#include "dimensions.h"
	28	#include "paramet.h"
	29	#include "control.h"
	30	#include "logic.h"
	31	#include "serre.h"
	32	#include "comdissnew.h"
[351]	33	#include "clesphys.h"
[237]	34	c
	35	c
	36	c local:
	37	c ------
	38
	39	CHARACTER ch172,ch272,ch372,ch412
	40	INTEGER tapeout
	41	REAL clonn,clatt,grossismxx,grossismyy
[240]	42	REAL dzoomxx,dzoomyy, tauxx,tauyy
[237]	43	LOGICAL fxyhypbb, ysinuss
	44	INTEGER i
	45
	46	c
	47	c -------------------------------------------------------------------
	48	c
	49	c ......... Version du 29/04/97 ..........
	50	c
	51	c Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot,
	52	c tetatemp ajoutes pour la dissipation .
	53	c
	54	c Autre parametre ajoute en fin de liste de tapedef : fxyhypb
	55	c
	56	c Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb.
	57	c Sinon , choix de fxynew , a derivee sinusoidale ..
	58	c
	59	c ...... etatinit = . TRUE. si defrun est appele dans ETAT0_LMD ou
	60	c LIMIT_LMD pour l'initialisation de start.dat (dic) et
	61	c de limit.dat ( dic) ...........
	62	c Sinon etatinit = . FALSE .
	63	c
	64	c Donc etatinit = .F. si on veut comparer les valeurs de grossismx ,
	65	c grossismy,clon,clat, fxyhypb lues sur le fichier start avec
	66	c celles passees par run.def , au debut du gcm, apres l'appel a
	67	c lectba .
	68	c Ces parmetres definissant entre autres la grille et doivent etre
	69	c pareils et coherents , sinon il y aura divergence du gcm .
	70	c
	71	c-----------------------------------------------------------------------
	72	c initialisations:
	73	c ----------------
	74
	75	tapeout = 6
	76
	77	c-----------------------------------------------------------------------
	78	c Parametres de controle du run:
	79	c-----------------------------------------------------------------------
	80
	81	!Config Key = dayref
	82	!Config Desc = Jour de l'etat initial
[351]	83	!Config Def = 1
[237]	84	!Config Help = Jour de l'etat initial ( = 350 si 20 Decembre ,
	85	!Config par expl. ,comme ici ) ... A completer
[351]	86	dayref=1
[237]	87	CALL getin('dayref', dayref)
	88
	89	!Config Key = anneeref
	90	!Config Desc = Annee de l'etat initial
	91	!Config Def = 1998
	92	!Config Help = Annee de l'etat initial
	93	!Config ( avec 4 chiffres ) ... A completer
	94	anneeref = 1998
	95	CALL getin('anneeref',anneeref)
	96
	97	!Config Key = nday
	98	!Config Desc = Nombre de jours d'integration
	99	!Config Def = 10
	100	!Config Help = Nombre de jours d'integration
	101	!Config ... On pourait aussi permettre des mois ou des annees !
	102	nday = 10
	103	CALL getin('nday',nday)
	104
	105	!Config Key = day_step
	106	!Config Desc = nombre de pas par jour
	107	!Config Def = 240
	108	!Config Help = nombre de pas par jour (multiple de iperiod) (
	109	!Config ici pour dt = 1 min )
	110	day_step = 240
	111	CALL getin('day_step',day_step)
	112
	113	!Config Key = iperiod
	114	!Config Desc = periode pour le pas Matsuno
	115	!Config Def = 5
	116	!Config Help = periode pour le pas Matsuno (en pas de temps)
	117	iperiod = 5
	118	CALL getin('iperiod',iperiod)
	119
	120
	121	!Config Key = iconser
	122	!Config Desc = periode de sortie des variables de controle
	123	!Config Def = 240
	124	!Config Help = periode de sortie des variables de controle
	125	!Config (En pas de temps)
	126	iconser = 240
	127	CALL getin('iconser', iconser)
	128
	129	!Config Key = iecri
	130	!Config Desc = periode d'ecriture du fichier histoire
	131	!Config Def = 1
	132	!Config Help = periode d'ecriture du fichier histoire (en jour)
	133	iecri = 1
	134	CALL getin('iecri',iecri)
	135
	136
	137	!Config Key = periodav
	138	!Config Desc = periode de stockage fichier histmoy
	139	!Config Def = 1
	140	!Config Help = periode de stockage fichier histmoy (en jour)
	141	periodav = 1.
	142	CALL getin('periodav',periodav)
	143
	144	!Config Key = idissip
	145	!Config Desc = periode de la dissipation
	146	!Config Def = 10
	147	!Config Help = periode de la dissipation
	148	!Config (en pas) ... a completer !
	149	idissip = 10
	150	CALL getin('idissip',idissip)
	151
	152	ccc .... P. Le Van , modif le 29/04/97 .pour la dissipation ...
	153	ccc
	154
	155	!Config Key = lstardis
	156	!Config Desc = choix de l'operateur de dissipation
	157	!Config Def = y
	158	!Config Help = choix de l'operateur de dissipation
	159	!Config 'y' si on veut star et 'n' si on veut non-start !
	160	!Config Moi y en a pas comprendre !
	161	lstardis = .TRUE.
	162	CALL getin('lstardis',lstardis)
	163
	164
	165	!Config Key = nitergdiv
	166	!Config Desc = Nombre d'iteration de gradiv
	167	!Config Def = 1
	168	!Config Help = nombre d'iterations de l'operateur de dissipation
	169	!Config gradiv
	170	nitergdiv = 1
	171	CALL getin('nitergdiv',nitergdiv)
	172
	173	!Config Key = nitergrot
	174	!Config Desc = nombre d'iterations de nxgradrot
	175	!Config Def = 2
	176	!Config Help = nombre d'iterations de l'operateur de dissipation
	177	!Config nxgradrot
	178	nitergrot = 2
	179	CALL getin('nitergrot',nitergrot)
	180
	181
	182	!Config Key = niterh
	183	!Config Desc = nombre d'iterations de divgrad
	184	!Config Def = 2
	185	!Config Help = nombre d'iterations de l'operateur de dissipation
	186	!Config divgrad
	187	niterh = 2
	188	CALL getin('niterh',niterh)
	189
	190
	191	!Config Key = tetagdiv
	192	!Config Desc = temps de dissipation pour div
	193	!Config Def = 7200
	194	!Config Help = temps de dissipation des plus petites longeur
	195	!Config d'ondes pour u,v (gradiv)
	196	tetagdiv = 7200.
	197	CALL getin('tetagdiv',tetagdiv)
	198
	199	!Config Key = tetagrot
	200	!Config Desc = temps de dissipation pour grad
	201	!Config Def = 7200
	202	!Config Help = temps de dissipation des plus petites longeur
	203	!Config d'ondes pour u,v (nxgradrot)
	204	tetagrot = 7200.
	205	CALL getin('tetagrot',tetagrot)
	206
	207	!Config Key = tetatemp
	208	!Config Desc = temps de dissipation pour h
	209	!Config Def = 7200
	210	!Config Help = temps de dissipation des plus petites longeur
	211	!Config d'ondes pour h (divgrad)
	212	tetatemp = 7200.
	213	CALL getin('tetatemp',tetatemp )
	214
	215	!Config Key = coefdis
	216	!Config Desc = coefficient pour gamdissip
	217	!Config Def = 0
	218	!Config Help = coefficient pour gamdissip
	219	coefdis = 0.
	220	CALL getin('coefdis',coefdis)
	221
	222	!Config Key = purmats
	223	!Config Desc = Schema d'integration
	224	!Config Def = n
	225	!Config Help = Choix du schema d'integration temporel.
	226	!Config y = pure Matsuno sinon c'est du Matsuno-leapfrog
	227	purmats = .FALSE.
	228	CALL getin('purmats',purmats)
	229
	230	c ...............................................................
	231
	232	!Config Key = physic
	233	!Config Desc = Avec ls physique
	234	!Config Def = y
	235	!Config Help = Permet de faire tourner le modele sans
	236	!Config physique.
	237	physic = .TRUE.
	238	CALL getin('physic',physic)
	239
	240
	241	!Config Key = iphysiq
	242	!Config Desc = Periode de la physique
	243	!Config Def = 5
	244	!Config Help = Periode de la physique en pas de temps de la dynamique.
	245	iphysiq = 5
	246	CALL getin('iphysiq', iphysiq)
	247
	248	ccc .... P.Le Van, ajout le 03/01/96 pour l'ecriture phys ...
	249	c
	250
	251	!Config Key = ecritphy
	252	!Config Desc = Frequence d'ecriture de la physique
	253	!Config Def = 1
	254	!Config Help = frequence de l'ecriture du fichier histphy
	255	!Config en jours.
	256	ecritphy = 1
	257	CALL getin('ecritphy',ecritphy)
	258
	259	!Config Key = cycle_diurne
	260	!Config Desc = Cycle ddiurne
	261	!Config Def = y
	262	!Config Help = Cette option permet d'eteidre le cycle diurne.
	263	!Config Peut etre util pour accelerer le code !
	264	cycle_diurne = .TRUE.
	265	CALL getin('cycle_diurne',cycle_diurne)
	266
	267	!Config Key = soil_model
	268	!Config Desc = Modele de sol
	269	!Config Def = y
	270	!Config Help = Choix du modele de sol (Thermique ?)
	271	!Config Option qui pourait un string afin de pouvoir
	272	!Config plus de choix ! Ou meme une liste d'options !
	273	soil_model = .TRUE.
	274	CALL getin('soil_model',soil_model)
	275
	276	!Config Key = new_oliq
	277	!Config Desc = Nouvelle eau liquide
	278	!Config Def = y
	279	!Config Help = Permet de mettre en route la
	280	!Config nouvelle parametrisation de l'eau liquide !
	281	new_oliq = .TRUE.
	282	CALL getin('new_oliq',new_oliq)
	283
	284	!Config Key = ok_orodr
	285	!Config Desc = Orodr ???
	286	!Config Def = y
	287	!Config Help = Y en a pas comprendre !
	288	!Config
	289	ok_orodr = .TRUE.
	290	CALL getin('ok_orodr',ok_orodr)
	291
	292	!Config Key = ok_orolf
	293	!Config Desc = Orolf ??
	294	!Config Def = y
	295	!Config Help = Connais pas !
	296	ok_orolf = .TRUE.
	297	CALL getin('ok_orolf', ok_orolf)
	298
	299	!Config Key = ok_limitvrai
	300	!Config Desc = Force la lecture de la bonne annee
	301	!Config Def = n
	302	!Config Help = On peut forcer le modele a lire le
	303	!Config fichier SST de la bonne annee. C'est une tres bonne
	304	!Config idee, pourquoi ne pas mettre toujours a y ???
	305	ok_limitvrai = .FALSE.
	306	CALL getin('ok_limitvrai',ok_limitvrai)
	307
	308	!Config Key = nbapp_rad
	309	!Config Desc = Frequence d'appel au rayonnement
	310	!Config Def = 12
	311	!Config Help = Nombre d'appels des routines de rayonnements
	312	!Config par jour.
	313	nbapp_rad = 12
	314	CALL getin('nbapp_rad',nbapp_rad)
	315
	316	!Config Key = iflag_con
	317	!Config Desc = Flag de convection
	318	!Config Def = 2
	319	!Config Help = Flag pour la convection les options suivantes existent :
	320	!Config 1 pour LMD,
	321	!Config 2 pour Tiedtke,
	322	!Config 3 pour CCM(NCAR)
	323	iflag_con = 2
	324	CALL getin('iflag_con',iflag_con)
	325
	326	DO i = 1, longcles
	327	clesphy0(i) = 0.
	328	ENDDO
	329	clesphy0(1) = FLOAT( iflag_con )
	330	clesphy0(2) = FLOAT( nbapp_rad )
	331
	332	IF( cycle_diurne ) clesphy0(3) = 1.
	333	IF( soil_model ) clesphy0(4) = 1.
	334	IF( new_oliq ) clesphy0(5) = 1.
	335	IF( ok_orodr ) clesphy0(6) = 1.
	336	IF( ok_orolf ) clesphy0(7) = 1.
	337	IF( ok_limitvrai ) clesphy0(8) = 1.
	338
	339
	340	ccc .... P. Le Van , ajout le 7/03/95 .pour le zoom ...
	341	c ......... ( modif le 17/04/96 ) .........
	342	c
	343	IF( etatinit ) GO TO 100
	344
	345	!Config Key = clon
	346	!Config Desc = centre du zoom, longitude
	347	!Config Def = 0
	348	!Config Help = longitude en degres du centre
	349	!Config du zoom
	350	clonn = 0.
	351	CALL getin('clon',clonn)
	352
	353	!Config Key = clat
	354	!Config Desc = centre du zoom, latitude
	355	!Config Def = 0
	356	!Config Help = latitude en degres du centre du zoom
	357	!Config
	358	clatt = 0.
	359	CALL getin('clat',clatt)
	360
	361	c
	362	c
	363	IF( ABS(clat - clatt).GE. 0.001 ) THEN
	364	PRINT *,' La valeur de clat passee par run.def est differente de
	365	* celle lue sur le fichier start '
	366	STOP
	367	ENDIF
	368
	369	!Config Key = grossismx
	370	!Config Desc = zoom en longitude
	371	!Config Def = 1.0
	372	!Config Help = facteur de grossissement du zoom,
	373	!Config selon la longitude
	374	grossismxx = 1.0
	375	CALL getin('grossismx',grossismxx)
	376
	377
	378	IF( ABS(grossismx - grossismxx).GE. 0.001 ) THEN
	379	PRINT *,' La valeur de grossismx passee par run.def est differente
	380	* de celle lue sur le fichier start '
	381	STOP
	382	ENDIF
	383
	384	!Config Key = grossismy
	385	!Config Desc = zoom en latitude
	386	!Config Def = 1.0
	387	!Config Help = facteur de grossissement du zoom,
	388	!Config selon la latitude
	389	grossismyy = 1.0
	390	CALL getin('grossismy',grossismyy)
	391
	392	IF( ABS(grossismy - grossismyy).GE. 0.001 ) THEN
	393	PRINT *,' La valeur de grossismy passee par run.def est differen
	394	* te de celle lue sur le fichier start '
	395	STOP
	396	ENDIF
	397
	398	IF( grossismx.LT.1. ) THEN
	399	PRINT ,' ATTENTION !! grossismx < 1 . * '
	400	STOP
	401	ELSE
	402	alphax = 1. - 1./ grossismx
	403	ENDIF
	404
	405
	406	IF( grossismy.LT.1. ) THEN
	407	PRINT ,' ATTENTION !! grossismy < 1 . * '
	408	STOP
	409	ELSE
	410	alphay = 1. - 1./ grossismy
	411	ENDIF
	412
	413	PRINT *,' alphax alphay defrun ',alphax,alphay
	414	c
	415	c alphax et alphay sont les anciennes formulat. des grossissements
	416	c
	417	c
	418
	419	!Config Key = fxyhypb
	420	!Config Desc = Fonction hyperbolique
	421	!Config Def = y
	422	!Config Help = Fonction f(y) hyperbolique si = .true.
	423	!Config sinon sinusoidale
	424	fxyhypbb = .TRUE.
	425	CALL getin('fxyhypb',fxyhypbb)
	426
	427	IF( .NOT.fxyhypb ) THEN
	428	IF( fxyhypbb ) THEN
	429	PRINT ,' ***** PBS DANS DEFRUN ****** '
	430	PRINT ,' ** fxyhypb lu sur le fichier start est F ',
	431	* 'alors qu il est T sur run.def ***'
	432	STOP
	433	ENDIF
	434	ELSE
	435	IF( .NOT.fxyhypbb ) THEN
	436	PRINT ,' ***** PBS DANS DEFRUN ****** '
	437	PRINT ,' ** fxyhypb lu sur le fichier start est T ',
	438	* 'alors qu il est F sur run.def **** '
	439	STOP
	440	ENDIF
	441	ENDIF
	442	c
	443	!Config Key = dzoomx
	444	!Config Desc = extension en longitude
	445	!Config Def = 0
	446	!Config Help = extension en longitude de la zone du zoom
	447	!Config ( fraction de la zone totale)
	448	dzoomxx = 0.0
	449	CALL getin('dzoomx',dzoomxx)
	450
	451	IF( fxyhypb ) THEN
	452	IF( ABS(dzoomx - dzoomxx).GE. 0.001 ) THEN
	453	PRINT *,' La valeur de dzoomx passee par run.def est differente
	454	* de celle lue sur le fichier start '
	455	STOP
	456	ENDIF
	457	ENDIF
	458
	459	!Config Key = dzoomy
	460	!Config Desc = extension en latitude
	461	!Config Def = 0
	462	!Config Help = extension en latitude de la zone du zoom
	463	!Config ( fraction de la zone totale)
	464	dzoomyy = 0.0
	465	CALL getin('dzoomy',dzoomyy)
	466
	467	IF( fxyhypb ) THEN
	468	IF( ABS(dzoomy - dzoomyy).GE. 0.001 ) THEN
	469	PRINT *,' La valeur de dzoomy passee par run.def est differente
	470	* de celle lue sur le fichier start '
	471	STOP
	472	ENDIF
	473	ENDIF
	474
[240]	475	!Config Key = taux
	476	!Config Desc = raideur du zoom en X
	477	!Config Def = 3
	478	!Config Help = raideur du zoom en X
	479	tauxx = 3.0
	480	CALL getin('taux',tauxx)
	481
	482	IF( fxyhypb ) THEN
	483	IF( ABS(taux - tauxx).GE. 0.001 ) THEN
	484	PRINT *,' La valeur de taux passee par run.def est differente
	485	* de celle lue sur le fichier start '
	486	STOP
	487	ENDIF
	488	ENDIF
	489
	490	!Config Key = tauyy
	491	!Config Desc = raideur du zoom en Y
	492	!Config Def = 3
	493	!Config Help = raideur du zoom en Y
	494	tauyy = 3.0
	495	CALL getin('tauy',tauyy)
	496
	497	IF( fxyhypb ) THEN
	498	IF( ABS(tauy - tauyy).GE. 0.001 ) THEN
	499	PRINT *,' La valeur de tauy passee par run.def est differente
	500	* de celle lue sur le fichier start '
	501	STOP
	502	ENDIF
	503	ENDIF
	504
[237]	505	cc
	506	IF( .NOT.fxyhypb ) THEN
	507
	508	!Config Key = ysinus
	509	!Config IF = !fxyhypb
	510	!Config Desc = Fonction en Sinus
	511	!Config Def = y
	512	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
	513	!Config sinon y = latit.
	514	ysinuss = .TRUE.
	515	CALL getin('ysinus',ysinuss)
	516
	517	IF( .NOT.ysinus ) THEN
	518	IF( ysinuss ) THEN
	519	PRINT ,' ***** PBS DANS DEFRUN ****** '
	520	PRINT ,' ** ysinus lu sur le fichier start est F ',
	521	* 'alors qu il est T sur run.def ***'
	522	STOP
	523	ENDIF
	524	ELSE
	525	IF( .NOT.ysinuss ) THEN
	526	PRINT ,' ***** PBS DANS DEFRUN ****** '
	527	PRINT ,' ** ysinus lu sur le fichier start est T ',
	528	* 'alors qu il est F sur run.def **** '
	529	STOP
	530	ENDIF
	531	ENDIF
	532	ENDIF
	533	c
	534
	535	write(tapeout,*)' #########################################'
	536	write(tapeout,*)' Configuration des parametres du gcm: '
	537	write(tapeout,*)' dayref = ', dayref
	538	write(tapeout,*)' anneeref = ', anneeref
	539	write(tapeout,*)' nday = ', nday
	540	write(tapeout,*)' day_step = ', day_step
	541	write(tapeout,*)' iperiod = ', iperiod
	542	write(tapeout,*)' iconser = ', iconser
	543	write(tapeout,*)' iecri = ', iecri
	544	write(tapeout,*)' periodav = ', periodav
	545	write(tapeout,*)' idissip = ', idissip
	546	write(tapeout,*)' lstardis = ', lstardis
	547	write(tapeout,*)' nitergdiv = ', nitergdiv
	548	write(tapeout,*)' nitergrot = ', nitergrot
	549	write(tapeout,*)' niterh = ', niterh
	550	write(tapeout,*)' tetagdiv = ', tetagdiv
	551	write(tapeout,*)' tetagrot = ', tetagrot
	552	write(tapeout,*)' tetatemp = ', tetatemp
	553	write(tapeout,*)' coefdis = ', coefdis
	554	write(tapeout,*)' purmats = ', purmats
	555	write(tapeout,*)' physic = ', physic
	556	write(tapeout,*)' iphysiq = ', iphysiq
	557	write(tapeout,*)' ecritphy = ', ecritphy
	558	write(tapeout,*)' cycle_diurne = ', cycle_diurne
	559	write(tapeout,*)' soil_model = ', soil_model
	560	write(tapeout,*)' new_oliq = ', new_oliq
	561	write(tapeout,*)' ok_orodr = ', ok_orodr
	562	write(tapeout,*)' ok_orolf = ', ok_orolf
	563	write(tapeout,*)' ok_limitvrai = ', ok_limitvrai
	564	write(tapeout,*)' nbapp_rad = ', nbapp_rad
	565	write(tapeout,*)' iflag_con = ', iflag_con
	566	write(tapeout,*)' clonn = ', clonn
	567	write(tapeout,*)' clatt = ', clatt
	568	write(tapeout,*)' grossismx = ', grossismx
	569	write(tapeout,*)' grossismy = ', grossismy
	570	write(tapeout,*)' fxyhypbb = ', fxyhypbb
	571	write(tapeout,*)' dzoomxx = ', dzoomxx
[240]	572	write(tapeout,*)' dzoomy = ', dzoomyy
	573	write(tapeout,*)' tauxx = ', tauxx
	574	write(tapeout,*)' tauyy = ', tauyy
[237]	575
	576	RETURN
	577	c ...............................................
	578	c
	579	100 CONTINUE
	580	!Config Key = clon
	581	!Config Desc = centre du zoom, longitude
	582	!Config Def = 0
	583	!Config Help = longitude en degres du centre
	584	!Config du zoom
	585	clon = 0.
	586	CALL getin('clon',clon)
	587
	588	!Config Key = clat
	589	!Config Desc = centre du zoom, latitude
	590	!Config Def = 0
	591	!Config Help = latitude en degres du centre du zoom
	592	!Config
	593	clat = 0.
	594	CALL getin('clat',clat)
	595
	596	!Config Key = grossismx
	597	!Config Desc = zoom en longitude
	598	!Config Def = 1.0
	599	!Config Help = facteur de grossissement du zoom,
	600	!Config selon la longitude
	601	grossismx = 1.0
	602	CALL getin('grossismx',grossismx)
	603
	604	!Config Key = grossismy
	605	!Config Desc = zoom en latitude
	606	!Config Def = 1.0
	607	!Config Help = facteur de grossissement du zoom,
	608	!Config selon la latitude
	609	grossismy = 1.0
	610	CALL getin('grossismy',grossismy)
	611
	612	IF( grossismx.LT.1. ) THEN
	613	PRINT ,' ATTENTION !! grossismx < 1 . * '
	614	STOP
	615	ELSE
	616	alphax = 1. - 1./ grossismx
	617	ENDIF
	618
	619
	620	IF( grossismy.LT.1. ) THEN
	621	PRINT ,' ATTENTION !! grossismy < 1 . * '
	622	STOP
	623	ELSE
	624	alphay = 1. - 1./ grossismy
	625	ENDIF
	626
	627	PRINT *,' alphax alphay defrun ',alphax,alphay
	628	c
	629	c alphax et alphay sont les anciennes formulat. des grossissements
	630	c
	631	c
	632
	633	!Config Key = fxyhypb
	634	!Config Desc = Fonction hyperbolique
	635	!Config Def = y
	636	!Config Help = Fonction f(y) hyperbolique si = .true.
	637	!Config sinon sinusoidale
	638	fxyhypb = .TRUE.
	639	CALL getin('fxyhypb',fxyhypb)
	640
	641	!Config Key = dzoomx
	642	!Config Desc = extension en longitude
	643	!Config Def = 0
	644	!Config Help = extension en longitude de la zone du zoom
	645	!Config ( fraction de la zone totale)
	646	dzoomx = 0.0
	647	CALL getin('dzoomx',dzoomx)
	648
	649	!Config Key = dzoomy
	650	!Config Desc = extension en latitude
	651	!Config Def = 0
	652	!Config Help = extension en latitude de la zone du zoom
	653	!Config ( fraction de la zone totale)
	654	dzoomy = 0.0
	655	CALL getin('dzoomy',dzoomy)
	656
[240]	657	!Config Key = taux
	658	!Config Desc = raideur du zoom en X
	659	!Config Def = 3
	660	!Config Help = raideur du zoom en X
	661	taux = 3.0
	662	CALL getin('taux',taux)
[237]	663
[240]	664	!Config Key = tauy
	665	!Config Desc = raideur du zoom en Y
	666	!Config Def = 3
	667	!Config Help = raideur du zoom en Y
	668	tauy = 3.0
	669	CALL getin('tauy',tauy)
	670
[237]	671	!Config Key = ysinus
	672	!Config IF = !fxyhypb
	673	!Config Desc = Fonction en Sinus
	674	!Config Def = y
	675	!Config Help = Fonction f(y) avec y = Sin(latit.) si = .true.
	676	!Config sinon y = latit.
	677	ysinus = .TRUE.
	678	CALL getin('ysinus',ysinus)
	679	c
	680
	681	write(tapeout,*)' #########################################'
	682	write(tapeout,*)' Configuration des parametres du gcm: '
	683	write(tapeout,*)' dayref = ', dayref
	684	write(tapeout,*)' anneeref = ', anneeref
	685	write(tapeout,*)' nday = ', nday
	686	write(tapeout,*)' day_step = ', day_step
	687	write(tapeout,*)' iperiod = ', iperiod
	688	write(tapeout,*)' iconser = ', iconser
	689	write(tapeout,*)' iecri = ', iecri
	690	write(tapeout,*)' periodav = ', periodav
	691	write(tapeout,*)' idissip = ', idissip
	692	write(tapeout,*)' lstardis = ', lstardis
	693	write(tapeout,*)' nitergdiv = ', nitergdiv
	694	write(tapeout,*)' nitergrot = ', nitergrot
	695	write(tapeout,*)' niterh = ', niterh
	696	write(tapeout,*)' tetagdiv = ', tetagdiv
	697	write(tapeout,*)' tetagrot = ', tetagrot
	698	write(tapeout,*)' tetatemp = ', tetatemp
	699	write(tapeout,*)' coefdis = ', coefdis
	700	write(tapeout,*)' purmats = ', purmats
	701	write(tapeout,*)' physic = ', physic
	702	write(tapeout,*)' iphysiq = ', iphysiq
	703	write(tapeout,*)' ecritphy = ', ecritphy
	704	write(tapeout,*)' cycle_diurne = ', cycle_diurne
	705	write(tapeout,*)' soil_model = ', soil_model
	706	write(tapeout,*)' new_oliq = ', new_oliq
	707	write(tapeout,*)' ok_orodr = ', ok_orodr
	708	write(tapeout,*)' ok_orolf = ', ok_orolf
	709	write(tapeout,*)' ok_limitvrai = ', ok_limitvrai
	710	write(tapeout,*)' nbapp_rad = ', nbapp_rad
	711	write(tapeout,*)' iflag_con = ', iflag_con
	712	write(tapeout,*)' clonn = ', clonn
	713	write(tapeout,*)' clatt = ', clatt
	714	write(tapeout,*)' grossismx = ', grossismx
	715	write(tapeout,*)' grossismy = ', grossismy
	716	write(tapeout,*)' fxyhypbb = ', fxyhypbb
[240]	717	write(tapeout,*)' dzoomx = ', dzoomx
[237]	718	write(tapeout,*)' dzoomy = ', dzoomy
[240]	719	write(tapeout,*)' taux = ', taux
	720	write(tapeout,*)' tauy = ', tauy
[237]	721	c
	722	RETURN
	723	END

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

Download in other formats:

Original Format