Context Navigation

conf_gcm.F @ 663

Last change on this file since 663 was 541, checked in by lmdzadmin, 20 years ago

Convergence avec la version d'Olivia Coindreau incluant:

le offline
les thermiques
mellor & yamada dans la couche limite

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 23.7 KB

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

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

Download in other formats:

Original Format