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

Last change on this file since 545 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
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  = 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
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
162ccc  ....   P. Le Van , modif le 29/04/97 .pour la dissipation  ...
163ccc
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
240c    ...............................................................
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
258ccc   .... P.Le Van, ajout le 03/01/96 pour l'ecriture phys ...
259c
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
350ccc  ....   P. Le Van , ajout  le 7/03/95 .pour le zoom ...
351c     .........   (  modif  le 17/04/96 )   .........
352c
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
371c
372c
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
424c
425c    alphax et alphay sont les anciennes formulat. des grossissements
426c
427c
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
452c
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     
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
515cc
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
543c
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
582      write(tapeout,*)' dzoomy = ', dzoomyy
583      write(tapeout,*)' tauxx = ', tauxx
584      write(tapeout,*)' tauyy = ', tauyy
585
586      RETURN
587c   ...............................................
588c
589100   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
638c
639c    alphax et alphay sont les anciennes formulat. des grossissements
640c
641c
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
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)
673
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
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)
689c
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
722      write(tapeout,*)' clonn = ', clon
723      write(tapeout,*)' clatt = ', clat
724      write(tapeout,*)' grossismx = ', grossismx
725      write(tapeout,*)' grossismy = ', grossismy
726      write(tapeout,*)' fxyhypbb = ', fxyhypbb
727      write(tapeout,*)' dzoomx = ', dzoomx
728      write(tapeout,*)' dzoomy = ', dzoomy
729      write(tapeout,*)' taux = ', taux
730      write(tapeout,*)' tauy = ', tauy
731c
732      RETURN
733      END
Note: See TracBrowser for help on using the repository browser.