source: trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F90 @ 3553

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

SL, Venus: new keys for flexibility cp0/cp(T) and Held-Suarez type physics

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