source: LMDZ5/trunk/libf/dyn3dmem/conf_gcm.F @ 2119

Last change on this file since 2119 was 2111, checked in by lguez, 10 years ago

Not possible to adjust MPI bands with nudging.

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