source: trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F @ 1000

Last change on this file since 1000 was 979, checked in by emillour, 12 years ago

Common dynamics:
Updates in the dynamics (seq and ) to keep up with updates
in LMDZ5 (up to LMDZ5 trunk, rev 1760):

  • General stuff (essentially to keep up with Earth model):
  • Updated makelmdz_fcm and makelmdz (more control on dimension.h, added option -mem, although it is not usefull for now)
  • Updated build_gcm with more control over fcm
  • Updated create_make_gcm (enable looking for code in subdirectories)
  • bibio:
  • updates (just renaming the files actually...) new_unit.F90 => new_unit_m.F90, pchsp_95.F90 => pchsp_95_m.F90 and pchfe_95.F90 => pchfe_95_m.F90
  • filtrez:
  • mod_fft.F90: use more baseline CPP directives for preprocessor compatibility
  • mod_filtre_fft_loc.F90: added this new file
  • filtreg_mod.f90: added calls to init_..._loc
  • filtreg.F: fixed calls to DGEMM into SGEMM (preprocessing does the switch)
  • dyn3d:
  • removed obsolete files: etat0_netcdf.F90 limit_netcdf.F90

pres2lev.F90

  • added new file : pres2lev_mod.F90 (module containing "old" pres2lev)
  • gcm.F: changed args to call to inidissip (added arg "vert_prof_dissip")
  • inidissip.F90: added arg "vert_prof_dissip" and the "earth model" discterizations (flagged with "planet_type=="earth")
  • comdissnew.h: added 'vert_prof_dissip' to the common block
  • guide_mod.F90: added the "use pres2lev_mod"
  • conf_gcm.F: cosmetics, and evaluation of vert_prof_dissip, (and also of dissip_* factors, for Earth model)
  • comconst.h : added dissip_factz,dissip_zref variables (for Earth mode dissip)
  • dyn3dpar:
  • removed obsolete files: etat0_netcdf.F90 limit_netcdf.F90

pres2lev.F90 mod_const_para.F90

  • added new files: pres2lev_mod.F90 (module containing "old" pres2lev)

mod_const_mpi

  • abort_gcm : better control of abort in parallel mode
  • gcm.F: changed args to call to inidissip (added arg "vert_prof_dissip")
  • inidissip.F90: added arg "vert_prof_dissip" and the "earth model" discterizations (flagged with "planet_type=="earth")
  • comdissnew.h: added 'vert_prof_dissip' to the common block
  • filtreg_p.F : bug correction (array bounds)
  • guide_p_mod.F90 : added the "use pres2lev_mod"
  • conf_gcm.F : cosmetics (and evaluation of vert_prof_dissip , and

also of dissip_* factors, for Earth model)
plus check if "adjust" is indeed not used in OpenMP

  • comconst.h : add dissip_factz,dissip_zref variables (for Earth mode dissip)

EM

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