source: trunk/libf/dyn3d/conf_gcm.F @ 6

Last change on this file since 6 was 6, checked in by slebonnois, 14 years ago

cf commit_v6.log :

  • manipulation traceurs
  • homogeneisation .def
  • bilan_dyn
  • etats initiaux start.nc
  • appels specifiques pour physique
File size: 27.3 KB
Line 
1!
2! $Id: conf_gcm.F 1418 2010-07-19 15:11:24Z jghattas $
3!
4c
5c
6      SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
7c
8      USE control_mod
9#ifdef CPP_IOIPSL
10      use IOIPSL
11#else
12! if not using IOIPSL, we still need to use (a local version of) getin
13      use ioipsl_getincom
14#endif
15      IMPLICIT NONE
16c-----------------------------------------------------------------------
17c     Auteurs :   L. Fairhead , P. Le Van  .
18c
19c     Arguments :
20c
21c     tapedef   :
22c     etatinit  :     = TRUE   , on ne  compare pas les valeurs des para-
23c     -metres  du zoom  avec  celles lues sur le fichier start .
24c      clesphy0 :  sortie  .
25c
26       LOGICAL etatinit
27       INTEGER tapedef
28
29       INTEGER        longcles
30       PARAMETER(     longcles = 20 )
31       REAL clesphy0( longcles )
32c
33c   Declarations :
34c   --------------
35#include "dimensions.h"
36#include "paramet.h"
37#include "logic.h"
38#include "serre.h"
39#include "comdissnew.h"
40#include "temps.h"
41#include "comconst.h"
42
43! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
44! #include "clesphys.h"
45#include "iniprint.h"
46c
47c
48c   local:
49c   ------
50
51      CHARACTER ch1*72,ch2*72,ch3*72,ch4*12
52      REAL clonn,clatt,grossismxx,grossismyy
53      REAL dzoomxx,dzoomyy, tauxx,tauyy
54      LOGICAL  fxyhypbb, ysinuss
55      INTEGER i
56      LOGICAL use_filtre_fft
57c
58c  -------------------------------------------------------------------
59c
60c       .........     Version  du 29/04/97       ..........
61c
62c   Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot,
63c      tetatemp   ajoutes  pour la dissipation   .
64c
65c   Autre parametre ajoute en fin de liste de tapedef : ** fxyhypb **
66c
67c  Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb.
68c    Sinon , choix de fxynew  , a derivee sinusoidale  ..
69c
70c   ......  etatinit = . TRUE. si defrun  est appele dans ETAT0_LMD  ou
71c         LIMIT_LMD  pour l'initialisation de start.dat (dic) et
72c                de limit.dat ( dic)                        ...........
73c           Sinon  etatinit = . FALSE .
74c
75c   Donc etatinit = .F.  si on veut comparer les valeurs de  grossismx ,
76c    grossismy,clon,clat, fxyhypb  lues sur  le fichier  start  avec
77c   celles passees  par run.def ,  au debut du gcm, apres l'appel a
78c    lectba . 
79c   Ces parmetres definissant entre autres la grille et doivent etre
80c   pareils et coherents , sinon il y aura  divergence du gcm .
81c
82c-----------------------------------------------------------------------
83c   initialisations:
84c   ----------------
85
86!Config  Key  = lunout
87!Config  Desc = unite de fichier pour les impressions
88!Config  Def  = 6
89!Config  Help = unite de fichier pour les impressions
90!Config         (defaut sortie standard = 6)
91      lunout=6
92      CALL getin('lunout', lunout)
93      IF (lunout /= 5 .and. lunout /= 6) THEN
94        OPEN(lunout,FILE='lmdz.out')
95      ENDIF
96
97!Config  Key  = prt_level
98!Config  Desc = niveau d'impressions de débogage
99!Config  Def  = 0
100!Config  Help = Niveau d'impression pour le débogage
101!Config         (0 = minimum d'impression)
102      prt_level = 0
103      CALL getin('prt_level',prt_level)
104
105c-----------------------------------------------------------------------
106c  Parametres de controle du run:
107c-----------------------------------------------------------------------
108!Config  Key  = planet_type
109!Config  Desc = planet type ("earth", "mars", "venus", ...)
110!Config  Def  = earth
111!Config  Help = this flag sets the type of atymosphere that is considered
112      planet_type="earth"
113      CALL getin('planet_type',planet_type)
114
115!Config  Key  = calend
116!Config  Desc = type de calendrier utilise
117!Config  Def  = earth_360d
118!Config  Help = valeur possible: earth_360d, earth_365d, earth_366d
119!Config         
120      calend = 'earth_360d'
121      CALL getin('calend', calend)
122
123!Config  Key  = dayref
124!Config  Desc = Jour de l'etat initial
125!Config  Def  = 1
126!Config  Help = Jour de l'etat initial ( = 350  si 20 Decembre ,
127!Config         par expl. ,comme ici ) ... A completer
128      dayref=1
129      CALL getin('dayref', dayref)
130
131!Config  Key  = anneeref
132!Config  Desc = Annee de l'etat initial
133!Config  Def  = 1998
134!Config  Help = Annee de l'etat  initial
135!Config         (   avec  4  chiffres   ) ... A completer
136      anneeref = 1998
137      CALL getin('anneeref',anneeref)
138
139!Config  Key  = raz_date
140!Config  Desc = Remise a zero de la date initiale
141!Config  Def  = 0 (pas de remise a zero)
142!Config  Help = Remise a zero de la date initiale
143!Config         0 pas de remise a zero, on garde la date du fichier restart
144!Config         1 prise en compte de la date de gcm.def avec remise a zero
145!Config         des compteurs de pas de temps
146      raz_date = 0
147      CALL getin('raz_date', raz_date)
148
149!Config  Key  = nday
150!Config  Desc = Nombre de jours d'integration
151!Config  Def  = 10
152!Config  Help = Nombre de jours d'integration
153!Config         ... On pourait aussi permettre des mois ou des annees !
154      nday = 10
155      CALL getin('nday',nday)
156
157!Config  Key  = day_step
158!Config  Desc = nombre de pas par jour
159!Config  Def  = 240
160!Config  Help = nombre de pas par jour (multiple de iperiod) (
161!Config          ici pour  dt = 1 min )
162       day_step = 240
163       CALL getin('day_step',day_step)
164
165!Config  Key  = nsplit_phys
166!Config  Desc = nombre de subdivisions par pas physique
167!Config  Def  = 1
168!Config  Help = nombre de subdivisions par pas physique
169       nsplit_phys = 1
170       CALL getin('nsplit_phys',nsplit_phys)
171
172!Config  Key  = iperiod
173!Config  Desc = periode pour le pas Matsuno
174!Config  Def  = 5
175!Config  Help = periode pour le pas Matsuno (en pas de temps)
176       iperiod = 5
177       CALL getin('iperiod',iperiod)
178
179!Config  Key  = iapp_tracvl
180!Config  Desc = frequence du groupement des flux
181!Config  Def  = iperiod
182!Config  Help = frequence du groupement des flux (en pas de temps)
183       iapp_tracvl = iperiod
184       CALL getin('iapp_tracvl',iapp_tracvl)
185
186!Config  Key  = iconser
187!Config  Desc = periode de sortie des variables de controle
188!Config  Def  = 240 
189!Config  Help = periode de sortie des variables de controle
190!Config         (En pas de temps)
191       iconser = 240 
192       CALL getin('iconser', iconser)
193
194!Config  Key  = iecri
195!Config  Desc = periode d'ecriture du fichier histoire
196!Config  Def  = 1
197!Config  Help = periode d'ecriture du fichier histoire (en jour)
198       iecri = 1
199       CALL getin('iecri',iecri)
200
201
202!Config  Key  = periodav
203!Config  Desc = periode de stockage fichier histmoy
204!Config  Def  = 1
205!Config  Help = periode de stockage fichier histmoy (en jour)
206       periodav = 1.
207       CALL getin('periodav',periodav)
208
209!Config  Key  = output_grads_dyn
210!Config  Desc = output dynamics diagnostics in 'dyn.dat' file
211!Config  Def  = n
212!Config  Help = output dynamics diagnostics in Grads-readable 'dyn.dat' file
213       output_grads_dyn=.false.
214       CALL getin('output_grads_dyn',output_grads_dyn)
215
216!Config  Key  = idissip
217!Config  Desc = periode de la dissipation
218!Config  Def  = 10
219!Config  Help = periode de la dissipation
220!Config         (en pas) ... a completer !
221       idissip = 10
222       CALL getin('idissip',idissip)
223
224ccc  ....   P. Le Van , modif le 29/04/97 .pour la dissipation  ...
225ccc
226
227!Config  Key  = lstardis
228!Config  Desc = choix de l'operateur de dissipation
229!Config  Def  = y
230!Config  Help = choix de l'operateur de dissipation
231!Config         'y' si on veut star et 'n' si on veut non-start !
232!Config         Moi y en a pas comprendre !
233       lstardis = .TRUE.
234       CALL getin('lstardis',lstardis)
235
236
237!Config  Key  = nitergdiv
238!Config  Desc = Nombre d'iteration de gradiv
239!Config  Def  = 1
240!Config  Help = nombre d'iterations de l'operateur de dissipation
241!Config         gradiv
242       nitergdiv = 1
243       CALL getin('nitergdiv',nitergdiv)
244
245!Config  Key  = nitergrot
246!Config  Desc = nombre d'iterations de nxgradrot
247!Config  Def  = 2
248!Config  Help = nombre d'iterations de l'operateur de dissipation 
249!Config         nxgradrot
250       nitergrot = 2
251       CALL getin('nitergrot',nitergrot)
252
253
254!Config  Key  = niterh
255!Config  Desc = nombre d'iterations de divgrad
256!Config  Def  = 2
257!Config  Help = nombre d'iterations de l'operateur de dissipation
258!Config         divgrad
259       niterh = 2
260       CALL getin('niterh',niterh)
261
262
263!Config  Key  = tetagdiv
264!Config  Desc = temps de dissipation pour div
265!Config  Def  = 7200
266!Config  Help = temps de dissipation des plus petites longeur
267!Config         d'ondes pour u,v (gradiv)
268       tetagdiv = 7200.
269       CALL getin('tetagdiv',tetagdiv)
270
271!Config  Key  = tetagrot
272!Config  Desc = temps de dissipation pour grad
273!Config  Def  = 7200
274!Config  Help = temps de dissipation des plus petites longeur
275!Config         d'ondes pour u,v (nxgradrot)
276       tetagrot = 7200.
277       CALL getin('tetagrot',tetagrot)
278
279!Config  Key  = tetatemp
280!Config  Desc = temps de dissipation pour h
281!Config  Def  = 7200
282!Config  Help =  temps de dissipation des plus petites longeur
283!Config         d'ondes pour h (divgrad)   
284       tetatemp  = 7200.
285       CALL getin('tetatemp',tetatemp )
286
287! Parametres controlant la variation sur la verticale des constantes de
288! dissipation.
289! Pour le moment actifs uniquement dans la version a 39 niveaux
290! avec ok_strato=y
291
292       dissip_factz=4.
293       dissip_deltaz=10.
294       dissip_zref=30.
295       CALL getin('dissip_factz',dissip_factz )
296       CALL getin('dissip_deltaz',dissip_deltaz )
297       CALL getin('dissip_zref',dissip_zref )
298
299       iflag_top_bound=1
300       tau_top_bound=1.e-5
301       CALL getin('iflag_top_bound',iflag_top_bound)
302       CALL getin('tau_top_bound',tau_top_bound)
303
304!Config  Key  = coefdis
305!Config  Desc = coefficient pour gamdissip
306!Config  Def  = 0
307!Config  Help = coefficient pour gamdissip 
308       coefdis = 0.
309       CALL getin('coefdis',coefdis)
310
311!Config  Key  = purmats
312!Config  Desc = Schema d'integration
313!Config  Def  = n
314!Config  Help = Choix du schema d'integration temporel.
315!Config         y = pure Matsuno sinon c'est du Matsuno-leapfrog
316       purmats = .FALSE.
317       CALL getin('purmats',purmats)
318
319!Config  Key  = ok_guide
320!Config  Desc = Guidage
321!Config  Def  = n
322!Config  Help = Guidage
323       ok_guide = .FALSE.
324       CALL getin('ok_guide',ok_guide)
325
326c    ...............................................................
327
328!Config  Key  =  read_start
329!Config  Desc = Initialize model using a 'start.nc' file
330!Config  Def  = y
331!Config  Help = y: intialize dynamical fields using a 'start.nc' file
332!               n: fields are initialized by 'iniacademic' routine
333       read_start= .true.
334       CALL getin('read_start',read_start)
335
336!Config  Key  = iflag_phys
337!Config  Desc = Avec ls physique
338!Config  Def  = 1
339!Config  Help = Permet de faire tourner le modele sans
340!Config         physique.
341       iflag_phys = 1
342       CALL getin('iflag_phys',iflag_phys)
343
344
345!Config  Key  =  iphysiq
346!Config  Desc = Periode de la physique
347!Config  Def  = 5
348!Config  Help = Periode de la physique en pas de temps de la dynamique.
349       iphysiq = 5
350       CALL getin('iphysiq', iphysiq)
351
352!Config  Key  = ip_ebil_dyn
353!Config  Desc = PRINT level for energy conserv. diag.
354!Config  Def  = 0
355!Config  Help = PRINT level for energy conservation diag. ;
356!               les options suivantes existent :
357!Config         0 pas de print
358!Config         1 pas de print
359!Config         2 print,
360       ip_ebil_dyn = 0
361       CALL getin('ip_ebil_dyn',ip_ebil_dyn)
362
363      DO i = 1, longcles
364       clesphy0(i) = 0.
365      ENDDO
366
367ccc  ....   P. Le Van , ajout  le 7/03/95 .pour le zoom ...
368c     .........   (  modif  le 17/04/96 )   .........
369c
370      IF( etatinit ) GO TO 100
371
372!Config  Key  = clon
373!Config  Desc = centre du zoom, longitude
374!Config  Def  = 0
375!Config  Help = longitude en degres du centre
376!Config         du zoom
377       clonn = 0.
378       CALL getin('clon',clonn)
379
380!Config  Key  = clat
381!Config  Desc = centre du zoom, latitude
382!Config  Def  = 0
383!Config  Help = latitude en degres du centre du zoom
384!Config         
385       clatt = 0.
386       CALL getin('clat',clatt)
387
388c
389c
390      IF( ABS(clat - clatt).GE. 0.001 )  THEN
391        write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
392     &    ' est differente de celle lue sur le fichier  start '
393        STOP
394      ENDIF
395
396!Config  Key  = grossismx
397!Config  Desc = zoom en longitude
398!Config  Def  = 1.0
399!Config  Help = facteur de grossissement du zoom,
400!Config         selon la longitude
401       grossismxx = 1.0
402       CALL getin('grossismx',grossismxx)
403
404
405      IF( ABS(grossismx - grossismxx).GE. 0.001 )  THEN
406        write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
407     &  'run.def est differente de celle lue sur le fichier  start '
408        STOP
409      ENDIF
410
411!Config  Key  = grossismy
412!Config  Desc = zoom en latitude
413!Config  Def  = 1.0
414!Config  Help = facteur de grossissement du zoom,
415!Config         selon la latitude
416       grossismyy = 1.0
417       CALL getin('grossismy',grossismyy)
418
419      IF( ABS(grossismy - grossismyy).GE. 0.001 )  THEN
420        write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
421     & 'run.def est differente de celle lue sur le fichier  start '
422        STOP
423      ENDIF
424     
425      IF( grossismx.LT.1. )  THEN
426        write(lunout,*)
427     &       'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
428         STOP
429      ELSE
430         alphax = 1. - 1./ grossismx
431      ENDIF
432
433
434      IF( grossismy.LT.1. )  THEN
435        write(lunout,*)
436     &       'conf_gcm: ***  ATTENTION !! grossismy < 1 .   *** '
437         STOP
438      ELSE
439         alphay = 1. - 1./ grossismy
440      ENDIF
441
442      write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
443c
444c    alphax et alphay sont les anciennes formulat. des grossissements
445c
446c
447
448!Config  Key  = fxyhypb
449!Config  Desc = Fonction  hyperbolique
450!Config  Def  = y
451!Config  Help = Fonction  f(y)  hyperbolique  si = .true. 
452!Config         sinon  sinusoidale
453       fxyhypbb = .TRUE.
454       CALL getin('fxyhypb',fxyhypbb)
455
456      IF( .NOT.fxyhypb )  THEN
457         IF( fxyhypbb )     THEN
458            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
459            write(lunout,*)' *** fxyhypb lu sur le fichier start est ',
460     *       'F alors  qu il est  T  sur  run.def  ***'
461              STOP
462         ENDIF
463      ELSE
464         IF( .NOT.fxyhypbb )   THEN
465            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
466            write(lunout,*)' ***  fxyhypb lu sur le fichier start est ',
467     *        'T alors  qu il est  F  sur  run.def  ****  '
468              STOP
469         ENDIF
470      ENDIF
471c
472!Config  Key  = dzoomx
473!Config  Desc = extension en longitude
474!Config  Def  = 0
475!Config  Help = extension en longitude  de la zone du zoom 
476!Config         ( fraction de la zone totale)
477       dzoomxx = 0.0
478       CALL getin('dzoomx',dzoomxx)
479
480      IF( fxyhypb )  THEN
481       IF( ABS(dzoomx - dzoomxx).GE. 0.001 )  THEN
482        write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
483     *  'run.def est differente de celle lue sur le fichier  start '
484        STOP
485       ENDIF
486      ENDIF
487
488!Config  Key  = dzoomy
489!Config  Desc = extension en latitude
490!Config  Def  = 0
491!Config  Help = extension en latitude de la zone  du zoom 
492!Config         ( fraction de la zone totale)
493       dzoomyy = 0.0
494       CALL getin('dzoomy',dzoomyy)
495
496      IF( fxyhypb )  THEN
497       IF( ABS(dzoomy - dzoomyy).GE. 0.001 )  THEN
498        write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
499     * 'run.def est differente de celle lue sur le fichier  start '
500        STOP
501       ENDIF
502      ENDIF
503     
504!Config  Key  = taux
505!Config  Desc = raideur du zoom en  X
506!Config  Def  = 3
507!Config  Help = raideur du zoom en  X
508       tauxx = 3.0
509       CALL getin('taux',tauxx)
510
511      IF( fxyhypb )  THEN
512       IF( ABS(taux - tauxx).GE. 0.001 )  THEN
513        write(lunout,*)'conf_gcm: La valeur de taux passee par ',
514     * 'run.def est differente de celle lue sur le fichier  start '
515        STOP
516       ENDIF
517      ENDIF
518
519!Config  Key  = tauyy
520!Config  Desc = raideur du zoom en  Y
521!Config  Def  = 3
522!Config  Help = raideur du zoom en  Y
523       tauyy = 3.0
524       CALL getin('tauy',tauyy)
525
526      IF( fxyhypb )  THEN
527       IF( ABS(tauy - tauyy).GE. 0.001 )  THEN
528        write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
529     * 'run.def est differente de celle lue sur le fichier  start '
530        STOP
531       ENDIF
532      ENDIF
533
534cc
535      IF( .NOT.fxyhypb  )  THEN
536
537!Config  Key  = ysinus
538!Config  IF   = !fxyhypb
539!Config  Desc = Fonction en Sinus
540!Config  Def  = y
541!Config  Help = Fonction  f(y) avec y = Sin(latit.) si = .true.
542!Config         sinon y = latit.
543       ysinuss = .TRUE.
544       CALL getin('ysinus',ysinuss)
545
546        IF( .NOT.ysinus )  THEN
547          IF( ysinuss )     THEN
548            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
549            write(lunout,*)' *** ysinus lu sur le fichier start est F',
550     *       ' alors  qu il est  T  sur  run.def  ***'
551            STOP
552          ENDIF
553        ELSE
554          IF( .NOT.ysinuss )   THEN
555            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
556            write(lunout,*)' *** ysinus lu sur le fichier start est T',
557     *        ' alors  qu il est  F  sur  run.def  ****  '
558              STOP
559          ENDIF
560        ENDIF
561      ENDIF ! of IF( .NOT.fxyhypb  )
562c
563!Config  Key  = offline
564!Config  Desc = Nouvelle eau liquide
565!Config  Def  = n
566!Config  Help = Permet de mettre en route la
567!Config         nouvelle parametrisation de l'eau liquide !
568       offline = .FALSE.
569       CALL getin('offline',offline)
570
571!Config  Key  = config_inca
572!Config  Desc = Choix de configuration de INCA
573!Config  Def  = none
574!Config  Help = Choix de configuration de INCA :
575!Config         'none' = sans INCA
576!Config         'chem' = INCA avec calcul de chemie
577!Config         'aero' = INCA avec calcul des aerosols
578      config_inca = 'none'
579      CALL getin('config_inca',config_inca)
580
581!Config  Key  = ok_dynzon
582!Config  Desc = calcul et sortie des transports
583!Config  Def  = n
584!Config  Help = Permet de mettre en route le calcul des transports
585!Config         
586      ok_dynzon = .FALSE.
587      CALL getin('ok_dynzon',ok_dynzon)
588
589!Config  Key  = ok_dyn_ins
590!Config  Desc = sorties instantanees dans la dynamique
591!Config  Def  = n
592!Config  Help =
593!Config         
594      ok_dyn_ins = .FALSE.
595      CALL getin('ok_dyn_ins',ok_dyn_ins)
596
597!Config  Key  = ok_dyn_ave
598!Config  Desc = sorties moyennes dans la dynamique
599!Config  Def  = n
600!Config  Help =
601!Config         
602      ok_dyn_ave = .FALSE.
603      CALL getin('ok_dyn_ave',ok_dyn_ave)
604
605
606      write(lunout,*)' #########################################'
607      write(lunout,*)' Configuration des parametres du gcm: '
608      write(lunout,*)' planet_type = ', planet_type
609      write(lunout,*)' calend = ', calend
610      write(lunout,*)' dayref = ', dayref
611      write(lunout,*)' anneeref = ', anneeref
612      write(lunout,*)' nday = ', nday
613      write(lunout,*)' day_step = ', day_step
614      write(lunout,*)' iperiod = ', iperiod
615      write(lunout,*)' iconser = ', iconser
616      write(lunout,*)' iecri = ', iecri
617      write(lunout,*)' periodav = ', periodav
618      write(lunout,*)' output_grads_dyn = ', output_grads_dyn
619      write(lunout,*)' idissip = ', idissip
620      write(lunout,*)' lstardis = ', lstardis
621      write(lunout,*)' nitergdiv = ', nitergdiv
622      write(lunout,*)' nitergrot = ', nitergrot
623      write(lunout,*)' niterh = ', niterh
624      write(lunout,*)' tetagdiv = ', tetagdiv
625      write(lunout,*)' tetagrot = ', tetagrot
626      write(lunout,*)' tetatemp = ', tetatemp
627      write(lunout,*)' coefdis = ', coefdis
628      write(lunout,*)' purmats = ', purmats
629      write(lunout,*)' read_start = ', read_start
630      write(lunout,*)' iflag_phys = ', iflag_phys
631      write(lunout,*)' iphysiq = ', iphysiq
632      write(lunout,*)' clonn = ', clonn
633      write(lunout,*)' clatt = ', clatt
634      write(lunout,*)' grossismx = ', grossismx
635      write(lunout,*)' grossismy = ', grossismy
636      write(lunout,*)' fxyhypbb = ', fxyhypbb
637      write(lunout,*)' dzoomxx = ', dzoomxx
638      write(lunout,*)' dzoomy = ', dzoomyy
639      write(lunout,*)' tauxx = ', tauxx
640      write(lunout,*)' tauyy = ', tauyy
641      write(lunout,*)' offline = ', offline
642      write(lunout,*)' config_inca = ', config_inca
643      write(lunout,*)' ok_dynzon = ', ok_dynzon
644      write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
645      write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
646
647      RETURN
648c   ...............................................
649c
650100   CONTINUE
651!Config  Key  = clon
652!Config  Desc = centre du zoom, longitude
653!Config  Def  = 0
654!Config  Help = longitude en degres du centre
655!Config         du zoom
656       clon = 0.
657       CALL getin('clon',clon)
658
659!Config  Key  = clat
660!Config  Desc = centre du zoom, latitude
661!Config  Def  = 0
662!Config  Help = latitude en degres du centre du zoom
663!Config         
664       clat = 0.
665       CALL getin('clat',clat)
666
667!Config  Key  = grossismx
668!Config  Desc = zoom en longitude
669!Config  Def  = 1.0
670!Config  Help = facteur de grossissement du zoom,
671!Config         selon la longitude
672       grossismx = 1.0
673       CALL getin('grossismx',grossismx)
674
675!Config  Key  = grossismy
676!Config  Desc = zoom en latitude
677!Config  Def  = 1.0
678!Config  Help = facteur de grossissement du zoom,
679!Config         selon la latitude
680       grossismy = 1.0
681       CALL getin('grossismy',grossismy)
682
683      IF( grossismx.LT.1. )  THEN
684        write(lunout,*)
685     &   'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
686         STOP
687      ELSE
688         alphax = 1. - 1./ grossismx
689      ENDIF
690
691
692      IF( grossismy.LT.1. )  THEN
693        write(lunout,*)
694     &  'conf_gcm: ***  ATTENTION !! grossismy < 1 .   *** '
695         STOP
696      ELSE
697         alphay = 1. - 1./ grossismy
698      ENDIF
699
700      write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
701c
702c    alphax et alphay sont les anciennes formulat. des grossissements
703c
704c
705
706!Config  Key  = fxyhypb
707!Config  Desc = Fonction  hyperbolique
708!Config  Def  = y
709!Config  Help = Fonction  f(y)  hyperbolique  si = .true. 
710!Config         sinon  sinusoidale
711       fxyhypb = .TRUE.
712       CALL getin('fxyhypb',fxyhypb)
713
714!Config  Key  = dzoomx
715!Config  Desc = extension en longitude
716!Config  Def  = 0
717!Config  Help = extension en longitude  de la zone du zoom 
718!Config         ( fraction de la zone totale)
719       dzoomx = 0.0
720       CALL getin('dzoomx',dzoomx)
721
722!Config  Key  = dzoomy
723!Config  Desc = extension en latitude
724!Config  Def  = 0
725!Config  Help = extension en latitude de la zone  du zoom 
726!Config         ( fraction de la zone totale)
727       dzoomy = 0.0
728       CALL getin('dzoomy',dzoomy)
729
730!Config  Key  = taux
731!Config  Desc = raideur du zoom en  X
732!Config  Def  = 3
733!Config  Help = raideur du zoom en  X
734       taux = 3.0
735       CALL getin('taux',taux)
736
737!Config  Key  = tauy
738!Config  Desc = raideur du zoom en  Y
739!Config  Def  = 3
740!Config  Help = raideur du zoom en  Y
741       tauy = 3.0
742       CALL getin('tauy',tauy)
743
744!Config  Key  = ysinus
745!Config  IF   = !fxyhypb
746!Config  Desc = Fonction en Sinus
747!Config  Def  = y
748!Config  Help = Fonction  f(y) avec y = Sin(latit.) si = .true.
749!Config         sinon y = latit.
750       ysinus = .TRUE.
751       CALL getin('ysinus',ysinus)
752c
753!Config  Key  = offline
754!Config  Desc = Nouvelle eau liquide
755!Config  Def  = n
756!Config  Help = Permet de mettre en route la
757!Config         nouvelle parametrisation de l'eau liquide !
758       offline = .FALSE.
759       CALL getin('offline',offline)
760
761!Config  Key  = config_inca
762!Config  Desc = Choix de configuration de INCA
763!Config  Def  = none
764!Config  Help = Choix de configuration de INCA :
765!Config         'none' = sans INCA
766!Config         'chem' = INCA avec calcul de chemie
767!Config         'aero' = INCA avec calcul des aerosols
768      config_inca = 'none'
769      CALL getin('config_inca',config_inca)
770
771!Config  Key  = ok_dynzon
772!Config  Desc = sortie des transports zonaux dans la dynamique
773!Config  Def  = n
774!Config  Help = 
775!Config         
776       ok_dynzon = .FALSE.
777       CALL getin('ok_dynzon',ok_dynzon)
778
779!Config  Key  = ok_dyn_ins
780!Config  Desc = sorties instantanees dans la dynamique
781!Config  Def  = n
782!Config  Help =
783!Config         
784      ok_dyn_ins = .FALSE.
785      CALL getin('ok_dyn_ins',ok_dyn_ins)
786
787!Config  Key  = ok_dyn_ave
788!Config  Desc = sorties moyennes dans la dynamique
789!Config  Def  = n
790!Config  Help =
791!Config         
792      ok_dyn_ave = .FALSE.
793      CALL getin('ok_dyn_ave',ok_dyn_ave)
794
795!Config  Key  = use_filtre_fft
796!Config  Desc = flag d'activation des FFT pour le filtre
797!Config  Def  = false
798!Config  Help = permet d'activer l'utilisation des FFT pour effectuer
799!Config         le filtrage aux poles.
800! Le filtre fft n'est pas implemente dans dyn3d
801      use_filtre_fft=.FALSE.
802      CALL getin('use_filtre_fft',use_filtre_fft)
803
804      IF (use_filtre_fft) THEN
805        write(lunout,*)'STOP !!!'
806        write(lunout,*)'use_filtre_fft n est pas implemente dans dyn3d'
807        STOP
808      ENDIF
809     
810!Config key = ok_strato
811!Config  Desc = activation de la version strato
812!Config  Def  = .FALSE.
813!Config  Help = active la version stratosphérique de LMDZ de F. Lott
814
815      ok_strato=.TRUE.
816      CALL getin('ok_strato',ok_strato)
817
818!Config  Key  = ok_gradsfile
819!Config  Desc = activation des sorties grads du guidage
820!Config  Def  = n
821!Config  Help = active les sorties grads du guidage
822
823       ok_gradsfile = .FALSE.
824       CALL getin('ok_gradsfile',ok_gradsfile)
825
826!Config  Key  = ok_limit
827!Config  Desc = creation des fichiers limit dans create_etat0_limit
828!Config  Def  = y
829!Config  Help = production du fichier limit.nc requise
830
831       ok_limit = .TRUE.
832       CALL getin('ok_limit',ok_limit)
833
834!Config  Key  = ok_etat0
835!Config  Desc = creation des fichiers etat0 dans create_etat0_limit
836!Config  Def  = y
837!Config  Help = production des fichiers start.nc, startphy.nc requise
838
839      ok_etat0 = .TRUE.
840      CALL getin('ok_etat0',ok_etat0)
841
842      write(lunout,*)' #########################################'
843      write(lunout,*)' Configuration des parametres de cel0'
844     &             //'_limit: '
845      write(lunout,*)' planet_type = ', planet_type
846      write(lunout,*)' calend = ', calend
847      write(lunout,*)' dayref = ', dayref
848      write(lunout,*)' anneeref = ', anneeref
849      write(lunout,*)' nday = ', nday
850      write(lunout,*)' day_step = ', day_step
851      write(lunout,*)' iperiod = ', iperiod
852      write(lunout,*)' iconser = ', iconser
853      write(lunout,*)' iecri = ', iecri
854      write(lunout,*)' periodav = ', periodav
855      write(lunout,*)' output_grads_dyn = ', output_grads_dyn
856      write(lunout,*)' idissip = ', idissip
857      write(lunout,*)' lstardis = ', lstardis
858      write(lunout,*)' nitergdiv = ', nitergdiv
859      write(lunout,*)' nitergrot = ', nitergrot
860      write(lunout,*)' niterh = ', niterh
861      write(lunout,*)' tetagdiv = ', tetagdiv
862      write(lunout,*)' tetagrot = ', tetagrot
863      write(lunout,*)' tetatemp = ', tetatemp
864      write(lunout,*)' coefdis = ', coefdis
865      write(lunout,*)' purmats = ', purmats
866      write(lunout,*)' read_start = ', read_start
867      write(lunout,*)' iflag_phys = ', iflag_phys
868      write(lunout,*)' iphysiq = ', iphysiq
869      write(lunout,*)' clon = ', clon
870      write(lunout,*)' clat = ', clat
871      write(lunout,*)' grossismx = ', grossismx
872      write(lunout,*)' grossismy = ', grossismy
873      write(lunout,*)' fxyhypb = ', fxyhypb
874      write(lunout,*)' dzoomx = ', dzoomx
875      write(lunout,*)' dzoomy = ', dzoomy
876      write(lunout,*)' taux = ', taux
877      write(lunout,*)' tauy = ', tauy
878      write(lunout,*)' offline = ', offline
879      write(lunout,*)' config_inca = ', config_inca
880      write(lunout,*)' ok_dynzon = ', ok_dynzon
881      write(lunout,*)' ok_dyn_ins = ', ok_dyn_ins
882      write(lunout,*)' ok_dyn_ave = ', ok_dyn_ave
883      write(lunout,*)' ok_strato = ', ok_strato
884      write(lunout,*)' ok_gradsfile = ', ok_gradsfile
885      write(lunout,*)' ok_limit = ', ok_limit
886      write(lunout,*)' ok_etat0 = ', ok_etat0
887c
888      RETURN
889      END
Note: See TracBrowser for help on using the repository browser.