source: LMDZ4/branches/LMDZ4-dev-20091210/libf/dyn3d/conf_gcm.F @ 5000

Last change on this file since 5000 was 1201, checked in by Laurent Fairhead, 15 years ago

Modifications nécessaires a l'inclusion d'un calendrier réaliste.
La date courante est calculée dans leapfrog.F et exprimée en Jour Julien
(modifié). On en a profité pour faire un peu de ménage dans la gestion des dates
du modèle.
Dans la physique, on utilise les routines de passages entre calendrier Julien et
Gregorien incluses dans IOIPSL pour calculer le nombre de jours écoulés depuis le
1er janvier (pour les conditions aux limites) ou l'equinoxe (pour le calcul de
la longitude solaire). Le calcul de l'orbite reprend celui du gcm planétaire
(codé par FH)
On décide du calendrier à utiliser à l'aide du paramètre calend du run.def. Par
défaut celui-ci est à earth_360d
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 24.9 KB
Line 
1!
2! $Id: conf_gcm.F 1201 2009-07-07 14:01:00Z abarral $
3!
4c
5c
6      SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
7c
8#ifdef CPP_IOIPSL
9      use IOIPSL
10#else
11! if not using IOIPSL, we still need to use (a local version of) getin
12      use ioipsl_getincom
13#endif
14      IMPLICIT NONE
15c-----------------------------------------------------------------------
16c     Auteurs :   L. Fairhead , P. Le Van  .
17c
18c     Arguments :
19c
20c     tapedef   :
21c     etatinit  :     = TRUE   , on ne  compare pas les valeurs des para-
22c     -metres  du zoom  avec  celles lues sur le fichier start .
23c      clesphy0 :  sortie  .
24c
25       LOGICAL etatinit
26       INTEGER tapedef
27
28       INTEGER        longcles
29       PARAMETER(     longcles = 20 )
30       REAL clesphy0( longcles )
31c
32c   Declarations :
33c   --------------
34#include "dimensions.h"
35#include "paramet.h"
36#include "control.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     
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  = iperiod
166!Config  Desc = periode pour le pas Matsuno
167!Config  Def  = 5
168!Config  Help = periode pour le pas Matsuno (en pas de temps)
169       iperiod = 5
170       CALL getin('iperiod',iperiod)
171
172!Config  Key  = iapp_tracvl
173!Config  Desc = frequence du groupement des flux
174!Config  Def  = iperiod
175!Config  Help = frequence du groupement des flux (en pas de temps)
176       iapp_tracvl = iperiod
177       CALL getin('iapp_tracvl',iapp_tracvl)
178
179!Config  Key  = iconser
180!Config  Desc = periode de sortie des variables de controle
181!Config  Def  = 240 
182!Config  Help = periode de sortie des variables de controle
183!Config         (En pas de temps)
184       iconser = 240 
185       CALL getin('iconser', iconser)
186
187!Config  Key  = iecri
188!Config  Desc = periode d'ecriture du fichier histoire
189!Config  Def  = 1
190!Config  Help = periode d'ecriture du fichier histoire (en jour)
191       iecri = 1
192       CALL getin('iecri',iecri)
193
194
195!Config  Key  = periodav
196!Config  Desc = periode de stockage fichier histmoy
197!Config  Def  = 1
198!Config  Help = periode de stockage fichier histmoy (en jour)
199       periodav = 1.
200       CALL getin('periodav',periodav)
201
202!Config  Key  = output_grads_dyn
203!Config  Desc = output dynamics diagnostics in 'dyn.dat' file
204!Config  Def  = n
205!Config  Help = output dynamics diagnostics in Grads-readable 'dyn.dat' file
206       output_grads_dyn=.false.
207       CALL getin('output_grads_dyn',output_grads_dyn)
208
209!Config  Key  = idissip
210!Config  Desc = periode de la dissipation
211!Config  Def  = 10
212!Config  Help = periode de la dissipation
213!Config         (en pas) ... a completer !
214       idissip = 10
215       CALL getin('idissip',idissip)
216
217ccc  ....   P. Le Van , modif le 29/04/97 .pour la dissipation  ...
218ccc
219
220!Config  Key  = lstardis
221!Config  Desc = choix de l'operateur de dissipation
222!Config  Def  = y
223!Config  Help = choix de l'operateur de dissipation
224!Config         'y' si on veut star et 'n' si on veut non-start !
225!Config         Moi y en a pas comprendre !
226       lstardis = .TRUE.
227       CALL getin('lstardis',lstardis)
228
229
230!Config  Key  = nitergdiv
231!Config  Desc = Nombre d'iteration de gradiv
232!Config  Def  = 1
233!Config  Help = nombre d'iterations de l'operateur de dissipation
234!Config         gradiv
235       nitergdiv = 1
236       CALL getin('nitergdiv',nitergdiv)
237
238!Config  Key  = nitergrot
239!Config  Desc = nombre d'iterations de nxgradrot
240!Config  Def  = 2
241!Config  Help = nombre d'iterations de l'operateur de dissipation 
242!Config         nxgradrot
243       nitergrot = 2
244       CALL getin('nitergrot',nitergrot)
245
246
247!Config  Key  = niterh
248!Config  Desc = nombre d'iterations de divgrad
249!Config  Def  = 2
250!Config  Help = nombre d'iterations de l'operateur de dissipation
251!Config         divgrad
252       niterh = 2
253       CALL getin('niterh',niterh)
254
255
256!Config  Key  = tetagdiv
257!Config  Desc = temps de dissipation pour div
258!Config  Def  = 7200
259!Config  Help = temps de dissipation des plus petites longeur
260!Config         d'ondes pour u,v (gradiv)
261       tetagdiv = 7200.
262       CALL getin('tetagdiv',tetagdiv)
263
264!Config  Key  = tetagrot
265!Config  Desc = temps de dissipation pour grad
266!Config  Def  = 7200
267!Config  Help = temps de dissipation des plus petites longeur
268!Config         d'ondes pour u,v (nxgradrot)
269       tetagrot = 7200.
270       CALL getin('tetagrot',tetagrot)
271
272!Config  Key  = tetatemp
273!Config  Desc = temps de dissipation pour h
274!Config  Def  = 7200
275!Config  Help =  temps de dissipation des plus petites longeur
276!Config         d'ondes pour h (divgrad)   
277       tetatemp  = 7200.
278       CALL getin('tetatemp',tetatemp )
279
280! Parametres controlant la variation sur la verticale des constantes de
281! dissipation.
282! Pour le moment actifs uniquement dans la version a 39 niveaux
283! avec ok_strato=y
284
285       dissip_factz=4.
286       dissip_deltaz=10.
287       dissip_zref=30.
288       CALL getin('dissip_factz',dissip_factz )
289       CALL getin('dissip_deltaz',dissip_deltaz )
290       CALL getin('dissip_zref',dissip_zref )
291
292       iflag_top_bound=1
293       tau_top_bound=1.e-5
294       CALL getin('iflag_top_bound',iflag_top_bound)
295       CALL getin('tau_top_bound',tau_top_bound)
296
297!Config  Key  = coefdis
298!Config  Desc = coefficient pour gamdissip
299!Config  Def  = 0
300!Config  Help = coefficient pour gamdissip 
301       coefdis = 0.
302       CALL getin('coefdis',coefdis)
303
304!Config  Key  = purmats
305!Config  Desc = Schema d'integration
306!Config  Def  = n
307!Config  Help = Choix du schema d'integration temporel.
308!Config         y = pure Matsuno sinon c'est du Matsuno-leapfrog
309       purmats = .FALSE.
310       CALL getin('purmats',purmats)
311
312!Config  Key  = ok_guide
313!Config  Desc = Guidage
314!Config  Def  = n
315!Config  Help = Guidage
316       ok_guide = .FALSE.
317       CALL getin('ok_guide',ok_guide)
318
319c    ...............................................................
320
321!Config  Key  =  read_start
322!Config  Desc = Initialize model using a 'start.nc' file
323!Config  Def  = y
324!Config  Help = y: intialize dynamical fields using a 'start.nc' file
325!               n: fields are initialized by 'iniacademic' routine
326       read_start= .true.
327       CALL getin('read_start',read_start)
328
329!Config  Key  = iflag_phys
330!Config  Desc = Avec ls physique
331!Config  Def  = 1
332!Config  Help = Permet de faire tourner le modele sans
333!Config         physique.
334       iflag_phys = 1
335       CALL getin('iflag_phys',iflag_phys)
336
337
338!Config  Key  =  iphysiq
339!Config  Desc = Periode de la physique
340!Config  Def  = 5
341!Config  Help = Periode de la physique en pas de temps de la dynamique.
342       iphysiq = 5
343       CALL getin('iphysiq', iphysiq)
344
345!Config  Key  = ip_ebil_dyn
346!Config  Desc = PRINT level for energy conserv. diag.
347!Config  Def  = 0
348!Config  Help = PRINT level for energy conservation diag. ;
349!               les options suivantes existent :
350!Config         0 pas de print
351!Config         1 pas de print
352!Config         2 print,
353       ip_ebil_dyn = 0
354       CALL getin('ip_ebil_dyn',ip_ebil_dyn)
355!
356
357      DO i = 1, longcles
358       clesphy0(i) = 0.
359      ENDDO
360
361ccc  ....   P. Le Van , ajout  le 7/03/95 .pour le zoom ...
362c     .........   (  modif  le 17/04/96 )   .........
363c
364      IF( etatinit ) GO TO 100
365
366!Config  Key  = clon
367!Config  Desc = centre du zoom, longitude
368!Config  Def  = 0
369!Config  Help = longitude en degres du centre
370!Config         du zoom
371       clonn = 0.
372       CALL getin('clon',clonn)
373
374!Config  Key  = clat
375!Config  Desc = centre du zoom, latitude
376!Config  Def  = 0
377!Config  Help = latitude en degres du centre du zoom
378!Config         
379       clatt = 0.
380       CALL getin('clat',clatt)
381
382c
383c
384      IF( ABS(clat - clatt).GE. 0.001 )  THEN
385        write(lunout,*)'conf_gcm: La valeur de clat passee par run.def',
386     &    ' est differente de celle lue sur le fichier  start '
387        STOP
388      ENDIF
389
390!Config  Key  = grossismx
391!Config  Desc = zoom en longitude
392!Config  Def  = 1.0
393!Config  Help = facteur de grossissement du zoom,
394!Config         selon la longitude
395       grossismxx = 1.0
396       CALL getin('grossismx',grossismxx)
397
398
399      IF( ABS(grossismx - grossismxx).GE. 0.001 )  THEN
400        write(lunout,*)'conf_gcm: La valeur de grossismx passee par ',
401     &  'run.def est differente de celle lue sur le fichier  start '
402        STOP
403      ENDIF
404
405!Config  Key  = grossismy
406!Config  Desc = zoom en latitude
407!Config  Def  = 1.0
408!Config  Help = facteur de grossissement du zoom,
409!Config         selon la latitude
410       grossismyy = 1.0
411       CALL getin('grossismy',grossismyy)
412
413      IF( ABS(grossismy - grossismyy).GE. 0.001 )  THEN
414        write(lunout,*)'conf_gcm: La valeur de grossismy passee par ',
415     & 'run.def est differente de celle lue sur le fichier  start '
416        STOP
417      ENDIF
418     
419      IF( grossismx.LT.1. )  THEN
420        write(lunout,*)
421     &       'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
422         STOP
423      ELSE
424         alphax = 1. - 1./ grossismx
425      ENDIF
426
427
428      IF( grossismy.LT.1. )  THEN
429        write(lunout,*)
430     &       'conf_gcm: ***  ATTENTION !! grossismy < 1 .   *** '
431         STOP
432      ELSE
433         alphay = 1. - 1./ grossismy
434      ENDIF
435
436      write(lunout,*)'conf_gcm: alphax alphay',alphax,alphay
437c
438c    alphax et alphay sont les anciennes formulat. des grossissements
439c
440c
441
442!Config  Key  = fxyhypb
443!Config  Desc = Fonction  hyperbolique
444!Config  Def  = y
445!Config  Help = Fonction  f(y)  hyperbolique  si = .true. 
446!Config         sinon  sinusoidale
447       fxyhypbb = .TRUE.
448       CALL getin('fxyhypb',fxyhypbb)
449
450      IF( .NOT.fxyhypb )  THEN
451         IF( fxyhypbb )     THEN
452            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
453            write(lunout,*)' *** fxyhypb lu sur le fichier start est ',
454     *       'F alors  qu il est  T  sur  run.def  ***'
455              STOP
456         ENDIF
457      ELSE
458         IF( .NOT.fxyhypbb )   THEN
459            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
460            write(lunout,*)' ***  fxyhypb lu sur le fichier start est ',
461     *        'T alors  qu il est  F  sur  run.def  ****  '
462              STOP
463         ENDIF
464      ENDIF
465c
466!Config  Key  = dzoomx
467!Config  Desc = extension en longitude
468!Config  Def  = 0
469!Config  Help = extension en longitude  de la zone du zoom 
470!Config         ( fraction de la zone totale)
471       dzoomxx = 0.0
472       CALL getin('dzoomx',dzoomxx)
473
474      IF( fxyhypb )  THEN
475       IF( ABS(dzoomx - dzoomxx).GE. 0.001 )  THEN
476        write(lunout,*)'conf_gcm: La valeur de dzoomx passee par ',
477     *  'run.def est differente de celle lue sur le fichier  start '
478        STOP
479       ENDIF
480      ENDIF
481
482!Config  Key  = dzoomy
483!Config  Desc = extension en latitude
484!Config  Def  = 0
485!Config  Help = extension en latitude de la zone  du zoom 
486!Config         ( fraction de la zone totale)
487       dzoomyy = 0.0
488       CALL getin('dzoomy',dzoomyy)
489
490      IF( fxyhypb )  THEN
491       IF( ABS(dzoomy - dzoomyy).GE. 0.001 )  THEN
492        write(lunout,*)'conf_gcm: La valeur de dzoomy passee par ',
493     * 'run.def est differente de celle lue sur le fichier  start '
494        STOP
495       ENDIF
496      ENDIF
497     
498!Config  Key  = taux
499!Config  Desc = raideur du zoom en  X
500!Config  Def  = 3
501!Config  Help = raideur du zoom en  X
502       tauxx = 3.0
503       CALL getin('taux',tauxx)
504
505      IF( fxyhypb )  THEN
506       IF( ABS(taux - tauxx).GE. 0.001 )  THEN
507        write(lunout,*)'conf_gcm: La valeur de taux passee par ',
508     * 'run.def est differente de celle lue sur le fichier  start '
509        STOP
510       ENDIF
511      ENDIF
512
513!Config  Key  = tauyy
514!Config  Desc = raideur du zoom en  Y
515!Config  Def  = 3
516!Config  Help = raideur du zoom en  Y
517       tauyy = 3.0
518       CALL getin('tauy',tauyy)
519
520      IF( fxyhypb )  THEN
521       IF( ABS(tauy - tauyy).GE. 0.001 )  THEN
522        write(lunout,*)'conf_gcm: La valeur de tauy passee par ',
523     * 'run.def est differente de celle lue sur le fichier  start '
524        STOP
525       ENDIF
526      ENDIF
527
528cc
529      IF( .NOT.fxyhypb  )  THEN
530
531!Config  Key  = ysinus
532!Config  IF   = !fxyhypb
533!Config  Desc = Fonction en Sinus
534!Config  Def  = y
535!Config  Help = Fonction  f(y) avec y = Sin(latit.) si = .true.
536!Config         sinon y = latit.
537       ysinuss = .TRUE.
538       CALL getin('ysinus',ysinuss)
539
540        IF( .NOT.ysinus )  THEN
541          IF( ysinuss )     THEN
542            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
543            write(lunout,*)' *** ysinus lu sur le fichier start est F',
544     *       ' alors  qu il est  T  sur  run.def  ***'
545            STOP
546          ENDIF
547        ELSE
548          IF( .NOT.ysinuss )   THEN
549            write(lunout,*)' ********  PBS DANS  CONF_GCM  ******** '
550            write(lunout,*)' *** ysinus lu sur le fichier start est T',
551     *        ' alors  qu il est  F  sur  run.def  ****  '
552              STOP
553          ENDIF
554        ENDIF
555      ENDIF ! of IF( .NOT.fxyhypb  )
556c
557!Config  Key  = offline
558!Config  Desc = Nouvelle eau liquide
559!Config  Def  = n
560!Config  Help = Permet de mettre en route la
561!Config         nouvelle parametrisation de l'eau liquide !
562       offline = .FALSE.
563       CALL getin('offline',offline)
564
565!Config  Key  = config_inca
566!Config  Desc = Choix de configuration de INCA
567!Config  Def  = none
568!Config  Help = Choix de configuration de INCA :
569!Config         'none' = sans INCA
570!Config         'chem' = INCA avec calcul de chemie
571!Config         'aero' = INCA avec calcul des aerosols
572      config_inca = 'none'
573      CALL getin('config_inca',config_inca)
574
575
576!Config  Key  = ok_dynzon
577!Config  Desc = calcul et sortie des transports
578!Config  Def  = n
579!Config  Help = Permet de mettre en route le calcul des transports
580!Config         
581      ok_dynzon = .FALSE.
582      CALL getin('ok_dynzon',ok_dynzon)
583
584      write(lunout,*)' #########################################'
585      write(lunout,*)' Configuration des parametres du gcm: '
586      write(lunout,*)' planet_type = ', planet_type
587      write(lunout,*)' calend = ', calend
588      write(lunout,*)' dayref = ', dayref
589      write(lunout,*)' anneeref = ', anneeref
590      write(lunout,*)' nday = ', nday
591      write(lunout,*)' day_step = ', day_step
592      write(lunout,*)' iperiod = ', iperiod
593      write(lunout,*)' iconser = ', iconser
594      write(lunout,*)' iecri = ', iecri
595      write(lunout,*)' periodav = ', periodav
596      write(lunout,*)' output_grads_dyn = ', output_grads_dyn
597      write(lunout,*)' idissip = ', idissip
598      write(lunout,*)' lstardis = ', lstardis
599      write(lunout,*)' nitergdiv = ', nitergdiv
600      write(lunout,*)' nitergrot = ', nitergrot
601      write(lunout,*)' niterh = ', niterh
602      write(lunout,*)' tetagdiv = ', tetagdiv
603      write(lunout,*)' tetagrot = ', tetagrot
604      write(lunout,*)' tetatemp = ', tetatemp
605      write(lunout,*)' coefdis = ', coefdis
606      write(lunout,*)' purmats = ', purmats
607      write(lunout,*)' read_start = ', read_start
608      write(lunout,*)' iflag_phys = ', iflag_phys
609      write(lunout,*)' iphysiq = ', iphysiq
610      write(lunout,*)' clonn = ', clonn
611      write(lunout,*)' clatt = ', clatt
612      write(lunout,*)' grossismx = ', grossismx
613      write(lunout,*)' grossismy = ', grossismy
614      write(lunout,*)' fxyhypbb = ', fxyhypbb
615      write(lunout,*)' dzoomxx = ', dzoomxx
616      write(lunout,*)' dzoomy = ', dzoomyy
617      write(lunout,*)' tauxx = ', tauxx
618      write(lunout,*)' tauyy = ', tauyy
619      write(lunout,*)' offline = ', offline
620      write(lunout,*)' config_inca = ', config_inca
621      write(lunout,*)' ok_dynzon = ', ok_dynzon
622
623      RETURN
624c   ...............................................
625c
626100   CONTINUE
627!Config  Key  = clon
628!Config  Desc = centre du zoom, longitude
629!Config  Def  = 0
630!Config  Help = longitude en degres du centre
631!Config         du zoom
632       clon = 0.
633       CALL getin('clon',clon)
634
635!Config  Key  = clat
636!Config  Desc = centre du zoom, latitude
637!Config  Def  = 0
638!Config  Help = latitude en degres du centre du zoom
639!Config         
640       clat = 0.
641       CALL getin('clat',clat)
642
643!Config  Key  = grossismx
644!Config  Desc = zoom en longitude
645!Config  Def  = 1.0
646!Config  Help = facteur de grossissement du zoom,
647!Config         selon la longitude
648       grossismx = 1.0
649       CALL getin('grossismx',grossismx)
650
651!Config  Key  = grossismy
652!Config  Desc = zoom en latitude
653!Config  Def  = 1.0
654!Config  Help = facteur de grossissement du zoom,
655!Config         selon la latitude
656       grossismy = 1.0
657       CALL getin('grossismy',grossismy)
658
659      IF( grossismx.LT.1. )  THEN
660        write(lunout,*)
661     &   'conf_gcm: ***  ATTENTION !! grossismx < 1 .   *** '
662         STOP
663      ELSE
664         alphax = 1. - 1./ grossismx
665      ENDIF
666
667
668      IF( grossismy.LT.1. )  THEN
669        write(lunout,*)
670     &  'conf_gcm: ***  ATTENTION !! grossismy < 1 .   *** '
671         STOP
672      ELSE
673         alphay = 1. - 1./ grossismy
674      ENDIF
675
676      write(lunout,*)'conf_gcm: alphax alphay ',alphax,alphay
677c
678c    alphax et alphay sont les anciennes formulat. des grossissements
679c
680c
681
682!Config  Key  = fxyhypb
683!Config  Desc = Fonction  hyperbolique
684!Config  Def  = y
685!Config  Help = Fonction  f(y)  hyperbolique  si = .true. 
686!Config         sinon  sinusoidale
687       fxyhypb = .TRUE.
688       CALL getin('fxyhypb',fxyhypb)
689
690!Config  Key  = dzoomx
691!Config  Desc = extension en longitude
692!Config  Def  = 0
693!Config  Help = extension en longitude  de la zone du zoom 
694!Config         ( fraction de la zone totale)
695       dzoomx = 0.0
696       CALL getin('dzoomx',dzoomx)
697
698!Config  Key  = dzoomy
699!Config  Desc = extension en latitude
700!Config  Def  = 0
701!Config  Help = extension en latitude de la zone  du zoom 
702!Config         ( fraction de la zone totale)
703       dzoomy = 0.0
704       CALL getin('dzoomy',dzoomy)
705
706!Config  Key  = taux
707!Config  Desc = raideur du zoom en  X
708!Config  Def  = 3
709!Config  Help = raideur du zoom en  X
710       taux = 3.0
711       CALL getin('taux',taux)
712
713!Config  Key  = tauy
714!Config  Desc = raideur du zoom en  Y
715!Config  Def  = 3
716!Config  Help = raideur du zoom en  Y
717       tauy = 3.0
718       CALL getin('tauy',tauy)
719
720!Config  Key  = ysinus
721!Config  IF   = !fxyhypb
722!Config  Desc = Fonction en Sinus
723!Config  Def  = y
724!Config  Help = Fonction  f(y) avec y = Sin(latit.) si = .true.
725!Config         sinon y = latit.
726       ysinus = .TRUE.
727       CALL getin('ysinus',ysinus)
728c
729!Config  Key  = offline
730!Config  Desc = Nouvelle eau liquide
731!Config  Def  = n
732!Config  Help = Permet de mettre en route la
733!Config         nouvelle parametrisation de l'eau liquide !
734       offline = .FALSE.
735       CALL getin('offline',offline)
736
737!Config  Key  = config_inca
738!Config  Desc = Choix de configuration de INCA
739!Config  Def  = none
740!Config  Help = Choix de configuration de INCA :
741!Config         'none' = sans INCA
742!Config         'chem' = INCA avec calcul de chemie
743!Config         'aero' = INCA avec calcul des aerosols
744      config_inca = 'none'
745      CALL getin('config_inca',config_inca)
746
747!Config  Key  = ok_dynzon
748!Config  Desc = calcul et sortie des transports
749!Config  Def  = n
750!Config  Help = Permet de mettre en route le calcul des transports
751!Config         
752       ok_dynzon = .FALSE.
753       CALL getin('ok_dynzon',ok_dynzon)
754
755!Config key = ok_strato
756!Config  Desc = activation de la version strato
757!Config  Def  = .FALSE.
758!Config  Help = active la version stratosphérique de LMDZ de F. Lott
759
760      ok_strato=.FALSE.
761      CALL getin('ok_strato',ok_strato)
762
763!Config  Key  = ok_gradsfile
764!Config  Desc = activation des sorties grads du guidage
765!Config  Def  = n
766!Config  Help = active les sorties grads du guidage
767
768       ok_gradsfile = .FALSE.
769       CALL getin('ok_gradsfile',ok_gradsfile)
770
771      write(lunout,*)' #########################################'
772      write(lunout,*)' Configuration des parametres du gcm: '
773      write(lunout,*)' planet_type = ', planet_type
774      write(lunout,*)' calend = ', calend
775      write(lunout,*)' dayref = ', dayref
776      write(lunout,*)' anneeref = ', anneeref
777      write(lunout,*)' nday = ', nday
778      write(lunout,*)' day_step = ', day_step
779      write(lunout,*)' iperiod = ', iperiod
780      write(lunout,*)' iconser = ', iconser
781      write(lunout,*)' iecri = ', iecri
782      write(lunout,*)' periodav = ', periodav
783      write(lunout,*)' output_grads_dyn = ', output_grads_dyn
784      write(lunout,*)' idissip = ', idissip
785      write(lunout,*)' lstardis = ', lstardis
786      write(lunout,*)' nitergdiv = ', nitergdiv
787      write(lunout,*)' nitergrot = ', nitergrot
788      write(lunout,*)' niterh = ', niterh
789      write(lunout,*)' tetagdiv = ', tetagdiv
790      write(lunout,*)' tetagrot = ', tetagrot
791      write(lunout,*)' tetatemp = ', tetatemp
792      write(lunout,*)' coefdis = ', coefdis
793      write(lunout,*)' purmats = ', purmats
794      write(lunout,*)' read_start = ', read_start
795      write(lunout,*)' iflag_phys = ', iflag_phys
796      write(lunout,*)' iphysiq = ', iphysiq
797      write(lunout,*)' clon = ', clon
798      write(lunout,*)' clat = ', clat
799      write(lunout,*)' grossismx = ', grossismx
800      write(lunout,*)' grossismy = ', grossismy
801      write(lunout,*)' fxyhypb = ', fxyhypb
802      write(lunout,*)' dzoomx = ', dzoomx
803      write(lunout,*)' dzoomy = ', dzoomy
804      write(lunout,*)' taux = ', taux
805      write(lunout,*)' tauy = ', tauy
806      write(lunout,*)' offline = ', offline
807      write(lunout,*)' config_inca = ', config_inca
808      write(lunout,*)' ok_dynzon = ', ok_dynzon
809      write(lunout,*)' ok_strato = ', ok_strato
810      write(lunout,*)' ok_gradsfile = ', ok_gradsfile
811c
812      RETURN
813      END
Note: See TracBrowser for help on using the repository browser.