source: LMDZ.3.3/branches/rel-LF/libf/dyn3d/conf_gcm.F @ 406

Last change on this file since 406 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
Line 
1c
2c $Header$
3c
4      SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
5c
6      use IOIPSL
7      IMPLICIT NONE
8c-----------------------------------------------------------------------
9c     Auteurs :   L. Fairhead , P. Le Van  .
10c
11c     Arguments :
12c
13c     tapedef   :
14c     etatinit  :     = TRUE   , on ne  compare pas les valeurs des para-
15c     -metres  du zoom  avec  celles lues sur le fichier start .
16c      clesphy0 :  sortie  .
17c
18       LOGICAL etatinit
19       INTEGER tapedef
20
21       INTEGER        longcles
22       PARAMETER(     longcles = 20 )
23       REAL clesphy0( longcles )
24c
25c   Declarations :
26c   --------------
27#include "dimensions.h"
28#include "paramet.h"
29#include "control.h"
30#include "logic.h"
31#include "serre.h"
32#include "comdissnew.h"
33#include "clesphys.h"
34c
35c
36c   local:
37c   ------
38
39      CHARACTER ch1*72,ch2*72,ch3*72,ch4*12
40      INTEGER   tapeout
41      REAL clonn,clatt,grossismxx,grossismyy
42      REAL dzoomxx,dzoomyy, tauxx,tauyy
43      LOGICAL  fxyhypbb, ysinuss
44      INTEGER i
45     
46c
47c  -------------------------------------------------------------------
48c
49c       .........     Version  du 29/04/97       ..........
50c
51c   Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot,
52c      tetatemp   ajoutes  pour la dissipation   .
53c
54c   Autre parametre ajoute en fin de liste de tapedef : ** fxyhypb **
55c
56c  Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb.
57c    Sinon , choix de fxynew  , a derivee sinusoidale  ..
58c
59c   ......  etatinit = . TRUE. si defrun  est appele dans ETAT0_LMD  ou
60c         LIMIT_LMD  pour l'initialisation de start.dat (dic) et
61c                de limit.dat ( dic)                        ...........
62c           Sinon  etatinit = . FALSE .
63c
64c   Donc etatinit = .F.  si on veut comparer les valeurs de  grossismx ,
65c    grossismy,clon,clat, fxyhypb  lues sur  le fichier  start  avec
66c   celles passees  par run.def ,  au debut du gcm, apres l'appel a
67c    lectba . 
68c   Ces parmetres definissant entre autres la grille et doivent etre
69c   pareils et coherents , sinon il y aura  divergence du gcm .
70c
71c-----------------------------------------------------------------------
72c   initialisations:
73c   ----------------
74
75      tapeout = 6
76
77c-----------------------------------------------------------------------
78c  Parametres de controle du run:
79c-----------------------------------------------------------------------
80
81!Config  Key  = dayref
82!Config  Desc = Jour de l'etat initial
83!Config  Def  = 1
84!Config  Help = Jour de l'etat initial ( = 350  si 20 Decembre ,
85!Config         par expl. ,comme ici ) ... A completer
86      dayref=1
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
152ccc  ....   P. Le Van , modif le 29/04/97 .pour la dissipation  ...
153ccc
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
230c    ...............................................................
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
248ccc   .... P.Le Van, ajout le 03/01/96 pour l'ecriture phys ...
249c
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
340ccc  ....   P. Le Van , ajout  le 7/03/95 .pour le zoom ...
341c     .........   (  modif  le 17/04/96 )   .........
342c
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
361c
362c
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
414c
415c    alphax et alphay sont les anciennes formulat. des grossissements
416c
417c
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
442c
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     
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
505cc
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
533c
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
572      write(tapeout,*)' dzoomy = ', dzoomyy
573      write(tapeout,*)' tauxx = ', tauxx
574      write(tapeout,*)' tauyy = ', tauyy
575
576      RETURN
577c   ...............................................
578c
579100   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
628c
629c    alphax et alphay sont les anciennes formulat. des grossissements
630c
631c
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
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)
663
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
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)
679c
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
717      write(tapeout,*)' dzoomx = ', dzoomx
718      write(tapeout,*)' dzoomy = ', dzoomy
719      write(tapeout,*)' taux = ', taux
720      write(tapeout,*)' tauy = ', tauy
721c
722      RETURN
723      END
Note: See TracBrowser for help on using the repository browser.