Context Navigation

conf_gcm.F @ 558

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

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

Download in other formats:

Original Format