source: LMDZ6/branches/DYNAMICO-conv/libf/dyn3d/conf_gcm.F90 @ 4976

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