source: LMDZ.3.3/trunk/libf/dyn3d/conf_gcm.F @ 2793

Last change on this file since 2793 was 237, checked in by lmdzadmin, 23 years ago

Initial release
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 20.6 KB
Line 
1c
2c $Header$
3c
4      SUBROUTINE conf_gcm( tapedef, etatinit, clesphy0 )
5c
6      use IOIPSL
7      IMPLICIT NONE
8c-----------------------------------------------------------------------
9c     Auteurs :   L. Fairhead , P. Le Van  .
10c
11c     Arguments :
12c
13c     tapedef   :
14c     etatinit  :     = TRUE   , on ne  compare pas les valeurs des para-
15c     -metres  du zoom  avec  celles lues sur le fichier start .
16c      clesphy0 :  sortie  .
17c
18       LOGICAL etatinit
19       INTEGER tapedef
20
21       INTEGER        longcles
22       PARAMETER(     longcles = 20 )
23       REAL clesphy0( longcles )
24c
25c   Declarations :
26c   --------------
27#include "dimensions.h"
28#include "paramet.h"
29#include "control.h"
30#include "logic.h"
31#include "serre.h"
32#include "comdissnew.h"
33#include "clesph0.h"
34c
35c
36c   local:
37c   ------
38
39      CHARACTER ch1*72,ch2*72,ch3*72,ch4*12
40      INTEGER   tapeout
41      REAL clonn,clatt,grossismxx,grossismyy
42      REAL dzoomxx,dzoomyy
43      LOGICAL  fxyhypbb, ysinuss
44      INTEGER i
45     
46c
47c  -------------------------------------------------------------------
48c
49c       .........     Version  du 29/04/97       ..........
50c
51c   Nouveaux parametres nitergdiv,nitergrot,niterh,tetagdiv,tetagrot,
52c      tetatemp   ajoutes  pour la dissipation   .
53c
54c   Autre parametre ajoute en fin de liste de tapedef : ** fxyhypb **
55c
56c  Si fxyhypb = .TRUE. , choix de la fonction a derivee tangente hyperb.
57c    Sinon , choix de fxynew  , a derivee sinusoidale  ..
58c
59c   ......  etatinit = . TRUE. si defrun  est appele dans ETAT0_LMD  ou
60c         LIMIT_LMD  pour l'initialisation de start.dat (dic) et
61c                de limit.dat ( dic)                        ...........
62c           Sinon  etatinit = . FALSE .
63c
64c   Donc etatinit = .F.  si on veut comparer les valeurs de  grossismx ,
65c    grossismy,clon,clat, fxyhypb  lues sur  le fichier  start  avec
66c   celles passees  par run.def ,  au debut du gcm, apres l'appel a
67c    lectba . 
68c   Ces parmetres definissant entre autres la grille et doivent etre
69c   pareils et coherents , sinon il y aura  divergence du gcm .
70c
71c-----------------------------------------------------------------------
72c   initialisations:
73c   ----------------
74
75      tapeout = 6
76
77c-----------------------------------------------------------------------
78c  Parametres de controle du run:
79c-----------------------------------------------------------------------
80
81!Config  Key  = dayref
82!Config  Desc = Jour de l'etat initial
83!Config  Def  = 181
84!Config  Help = Jour de l'etat initial ( = 350  si 20 Decembre ,
85!Config         par expl. ,comme ici ) ... A completer
86      dayref=181
87      CALL getin('dayref', dayref)
88
89!Config  Key  = anneeref
90!Config  Desc = Annee de l'etat initial
91!Config  Def  = 1998
92!Config  Help = Annee de l'etat  initial
93!Config         (   avec  4  chiffres   ) ... A completer
94      anneeref = 1998
95      CALL getin('anneeref',anneeref)
96
97!Config  Key  = nday
98!Config  Desc = Nombre de jours d'integration
99!Config  Def  = 10
100!Config  Help = Nombre de jours d'integration
101!Config         ... On pourait aussi permettre des mois ou des annees !
102      nday = 10
103      CALL getin('nday',nday)
104
105!Config  Key  = day_step
106!Config  Desc = nombre de pas par jour
107!Config  Def  = 240
108!Config  Help = nombre de pas par jour (multiple de iperiod) (
109!Config          ici pour  dt = 1 min )
110       day_step = 240
111       CALL getin('day_step',day_step)
112
113!Config  Key  = iperiod
114!Config  Desc = periode pour le pas Matsuno
115!Config  Def  = 5
116!Config  Help = periode pour le pas Matsuno (en pas de temps)
117       iperiod = 5
118       CALL getin('iperiod',iperiod)
119
120
121!Config  Key  = iconser
122!Config  Desc = periode de sortie des variables de controle
123!Config  Def  = 240 
124!Config  Help = periode de sortie des variables de controle
125!Config         (En pas de temps)
126       iconser = 240 
127       CALL getin('iconser', iconser)
128
129!Config  Key  = iecri
130!Config  Desc = periode d'ecriture du fichier histoire
131!Config  Def  = 1
132!Config  Help = periode d'ecriture du fichier histoire (en jour)
133       iecri = 1
134       CALL getin('iecri',iecri)
135
136
137!Config  Key  = periodav
138!Config  Desc = periode de stockage fichier histmoy
139!Config  Def  = 1
140!Config  Help = periode de stockage fichier histmoy (en jour)
141       periodav = 1.
142       CALL getin('periodav',periodav)
143
144!Config  Key  = idissip
145!Config  Desc = periode de la dissipation
146!Config  Def  = 10
147!Config  Help = periode de la dissipation
148!Config         (en pas) ... a completer !
149       idissip = 10
150       CALL getin('idissip',idissip)
151
152ccc  ....   P. Le Van , modif le 29/04/97 .pour la dissipation  ...
153ccc
154
155!Config  Key  = lstardis
156!Config  Desc = choix de l'operateur de dissipation
157!Config  Def  = y
158!Config  Help = choix de l'operateur de dissipation
159!Config         'y' si on veut star et 'n' si on veut non-start !
160!Config         Moi y en a pas comprendre !
161       lstardis = .TRUE.
162       CALL getin('lstardis',lstardis)
163
164
165!Config  Key  = nitergdiv
166!Config  Desc = Nombre d'iteration de gradiv
167!Config  Def  = 1
168!Config  Help = nombre d'iterations de l'operateur de dissipation
169!Config         gradiv
170       nitergdiv = 1
171       CALL getin('nitergdiv',nitergdiv)
172
173!Config  Key  = nitergrot
174!Config  Desc = nombre d'iterations de nxgradrot
175!Config  Def  = 2
176!Config  Help = nombre d'iterations de l'operateur de dissipation 
177!Config         nxgradrot
178       nitergrot = 2
179       CALL getin('nitergrot',nitergrot)
180
181
182!Config  Key  = niterh
183!Config  Desc = nombre d'iterations de divgrad
184!Config  Def  = 2
185!Config  Help = nombre d'iterations de l'operateur de dissipation
186!Config         divgrad
187       niterh = 2
188       CALL getin('niterh',niterh)
189
190
191!Config  Key  = tetagdiv
192!Config  Desc = temps de dissipation pour div
193!Config  Def  = 7200
194!Config  Help = temps de dissipation des plus petites longeur
195!Config         d'ondes pour u,v (gradiv)
196       tetagdiv = 7200.
197       CALL getin('tetagdiv',tetagdiv)
198
199!Config  Key  = tetagrot
200!Config  Desc = temps de dissipation pour grad
201!Config  Def  = 7200
202!Config  Help = temps de dissipation des plus petites longeur
203!Config         d'ondes pour u,v (nxgradrot)
204       tetagrot = 7200.
205       CALL getin('tetagrot',tetagrot)
206
207!Config  Key  = tetatemp
208!Config  Desc = temps de dissipation pour h
209!Config  Def  = 7200
210!Config  Help =  temps de dissipation des plus petites longeur
211!Config         d'ondes pour h (divgrad)   
212       tetatemp  = 7200.
213       CALL getin('tetatemp',tetatemp )
214
215!Config  Key  = coefdis
216!Config  Desc = coefficient pour gamdissip
217!Config  Def  = 0
218!Config  Help = coefficient pour gamdissip 
219       coefdis = 0.
220       CALL getin('coefdis',coefdis)
221
222!Config  Key  = purmats
223!Config  Desc = Schema d'integration
224!Config  Def  = n
225!Config  Help = Choix du schema d'integration temporel.
226!Config         y = pure Matsuno sinon c'est du Matsuno-leapfrog
227       purmats = .FALSE.
228       CALL getin('purmats',purmats)
229
230c    ...............................................................
231
232!Config  Key  = physic
233!Config  Desc = Avec ls physique
234!Config  Def  = y
235!Config  Help = Permet de faire tourner le modele sans
236!Config         physique.
237       physic = .TRUE.
238       CALL getin('physic',physic)
239
240
241!Config  Key  =  iphysiq
242!Config  Desc = Periode de la physique
243!Config  Def  = 5
244!Config  Help = Periode de la physique en pas de temps de la dynamique.
245       iphysiq = 5
246       CALL getin('iphysiq', iphysiq)
247
248ccc   .... P.Le Van, ajout le 03/01/96 pour l'ecriture phys ...
249c
250
251!Config  Key  = ecritphy
252!Config  Desc = Frequence d'ecriture de la physique
253!Config  Def  = 1
254!Config  Help = frequence  de l'ecriture du fichier histphy
255!Config         en jours.
256       ecritphy = 1
257       CALL getin('ecritphy',ecritphy)
258
259!Config  Key  = cycle_diurne
260!Config  Desc = Cycle ddiurne
261!Config  Def  = y
262!Config  Help = Cette option permet d'eteidre le cycle diurne.
263!Config         Peut etre util pour accelerer le code !
264       cycle_diurne = .TRUE.
265       CALL getin('cycle_diurne',cycle_diurne)
266
267!Config  Key  = soil_model
268!Config  Desc = Modele de sol
269!Config  Def  = y
270!Config  Help = Choix du modele de sol (Thermique ?)
271!Config         Option qui pourait un string afin de pouvoir
272!Config         plus de choix ! Ou meme une liste d'options !
273       soil_model = .TRUE.
274       CALL getin('soil_model',soil_model)
275
276!Config  Key  = new_oliq
277!Config  Desc = Nouvelle eau liquide
278!Config  Def  = y
279!Config  Help = Permet de mettre en route la
280!Config         nouvelle parametrisation de l'eau liquide !
281       new_oliq = .TRUE.
282       CALL getin('new_oliq',new_oliq)
283
284!Config  Key  = ok_orodr
285!Config  Desc = Orodr ???
286!Config  Def  = y
287!Config  Help = Y en a pas comprendre !
288!Config         
289       ok_orodr = .TRUE.
290       CALL getin('ok_orodr',ok_orodr)
291
292!Config  Key  =  ok_orolf
293!Config  Desc = Orolf ??
294!Config  Def  = y
295!Config  Help = Connais pas !
296       ok_orolf = .TRUE.
297       CALL getin('ok_orolf', ok_orolf)
298
299!Config  Key  = ok_limitvrai
300!Config  Desc = Force la lecture de la bonne annee
301!Config  Def  = n
302!Config  Help = On peut forcer le modele a lire le
303!Config         fichier SST de la bonne annee. C'est une tres bonne
304!Config         idee, pourquoi ne pas mettre toujours a y ???
305       ok_limitvrai = .FALSE.
306       CALL getin('ok_limitvrai',ok_limitvrai)
307
308!Config  Key  = nbapp_rad
309!Config  Desc = Frequence d'appel au rayonnement
310!Config  Def  = 12
311!Config  Help = Nombre  d'appels des routines de rayonnements
312!Config         par jour.
313       nbapp_rad = 12
314       CALL getin('nbapp_rad',nbapp_rad)
315
316!Config  Key  = iflag_con
317!Config  Desc = Flag de convection
318!Config  Def  = 2
319!Config  Help = Flag  pour la convection les options suivantes existent :
320!Config         1 pour LMD,
321!Config         2 pour Tiedtke,
322!Config         3 pour CCM(NCAR) 
323       iflag_con = 2
324       CALL getin('iflag_con',iflag_con)
325
326      DO i = 1, longcles
327       clesphy0(i) = 0.
328      ENDDO
329                          clesphy0(1) = FLOAT( iflag_con )
330                          clesphy0(2) = FLOAT( nbapp_rad )
331
332       IF( cycle_diurne  ) clesphy0(3) =  1.
333       IF(   soil_model  ) clesphy0(4) =  1.
334       IF(     new_oliq  ) clesphy0(5) =  1.
335       IF(     ok_orodr  ) clesphy0(6) =  1.
336       IF(     ok_orolf  ) clesphy0(7) =  1.
337       IF(  ok_limitvrai ) clesphy0(8) =  1.
338
339
340ccc  ....   P. Le Van , ajout  le 7/03/95 .pour le zoom ...
341c     .........   (  modif  le 17/04/96 )   .........
342c
343      IF( etatinit ) GO TO 100
344
345!Config  Key  = clon
346!Config  Desc = centre du zoom, longitude
347!Config  Def  = 0
348!Config  Help = longitude en degres du centre
349!Config         du zoom
350       clonn = 0.
351       CALL getin('clon',clonn)
352
353!Config  Key  = clat
354!Config  Desc = centre du zoom, latitude
355!Config  Def  = 0
356!Config  Help = latitude en degres du centre du zoom
357!Config         
358       clatt = 0.
359       CALL getin('clat',clatt)
360
361c
362c
363      IF( ABS(clat - clatt).GE. 0.001 )  THEN
364        PRINT *,' La valeur de clat passee par run.def est differente de
365     *  celle lue sur le fichier  start '
366        STOP
367      ENDIF
368
369!Config  Key  = grossismx
370!Config  Desc = zoom en longitude
371!Config  Def  = 1.0
372!Config  Help = facteur de grossissement du zoom,
373!Config         selon la longitude
374       grossismxx = 1.0
375       CALL getin('grossismx',grossismxx)
376
377
378      IF( ABS(grossismx - grossismxx).GE. 0.001 )  THEN
379        PRINT *,' La valeur de grossismx passee par run.def est differente 
380     *  de celle lue sur le fichier  start '
381        STOP
382      ENDIF
383
384!Config  Key  = grossismy
385!Config  Desc = zoom en latitude
386!Config  Def  = 1.0
387!Config  Help = facteur de grossissement du zoom,
388!Config         selon la latitude
389       grossismyy = 1.0
390       CALL getin('grossismy',grossismyy)
391
392      IF( ABS(grossismy - grossismyy).GE. 0.001 )  THEN
393        PRINT *,' La valeur de grossismy passee par run.def est differen
394     * te de celle lue sur le fichier  start '
395        STOP
396      ENDIF
397     
398      IF( grossismx.LT.1. )  THEN
399        PRINT *,' ***  ATTENTION !! grossismx < 1 .   *** '
400         STOP
401      ELSE
402         alphax = 1. - 1./ grossismx
403      ENDIF
404
405
406      IF( grossismy.LT.1. )  THEN
407        PRINT *,' ***  ATTENTION !! grossismy < 1 .   *** '
408         STOP
409      ELSE
410         alphay = 1. - 1./ grossismy
411      ENDIF
412
413      PRINT *,' alphax alphay defrun ',alphax,alphay
414c
415c    alphax et alphay sont les anciennes formulat. des grossissements
416c
417c
418
419!Config  Key  = fxyhypb
420!Config  Desc = Fonction  hyperbolique
421!Config  Def  = y
422!Config  Help = Fonction  f(y)  hyperbolique  si = .true. 
423!Config         sinon  sinusoidale
424       fxyhypbb = .TRUE.
425       CALL getin('fxyhypb',fxyhypbb)
426
427      IF( .NOT.fxyhypb )  THEN
428           IF( fxyhypbb )     THEN
429              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
430              PRINT *,' *** fxyhypb lu sur le fichier start est F ',
431     *       'alors  qu il est  T  sur  run.def  ***'
432              STOP
433           ENDIF
434      ELSE
435           IF( .NOT.fxyhypbb )   THEN
436              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
437              PRINT *,' ***  fxyhypb lu sur le fichier start est T ',
438     *        'alors  qu il est  F  sur  run.def  ****  '
439              STOP
440           ENDIF
441      ENDIF
442c
443!Config  Key  = dzoomx
444!Config  Desc = extension en longitude
445!Config  Def  = 0
446!Config  Help = extension en longitude  de la zone du zoom 
447!Config         ( fraction de la zone totale)
448       dzoomxx = 0.0
449       CALL getin('dzoomx',dzoomxx)
450
451      IF( fxyhypb )  THEN
452       IF( ABS(dzoomx - dzoomxx).GE. 0.001 )  THEN
453        PRINT *,' La valeur de dzoomx passee par run.def est differente
454     *  de celle lue sur le fichier  start '
455        STOP
456       ENDIF
457      ENDIF
458
459!Config  Key  = dzoomy
460!Config  Desc = extension en latitude
461!Config  Def  = 0
462!Config  Help = extension en latitude de la zone  du zoom 
463!Config         ( fraction de la zone totale)
464       dzoomyy = 0.0
465       CALL getin('dzoomy',dzoomyy)
466
467      IF( fxyhypb )  THEN
468       IF( ABS(dzoomy - dzoomyy).GE. 0.001 )  THEN
469        PRINT *,' La valeur de dzoomy passee par run.def est differente
470     * de celle lue sur le fichier  start '
471        STOP
472       ENDIF
473      ENDIF
474     
475cc
476      IF( .NOT.fxyhypb  )  THEN
477
478!Config  Key  = ysinus
479!Config  IF   = !fxyhypb
480!Config  Desc = Fonction en Sinus
481!Config  Def  = y
482!Config  Help = Fonction  f(y) avec y = Sin(latit.) si = .true.
483!Config         sinon y = latit.
484       ysinuss = .TRUE.
485       CALL getin('ysinus',ysinuss)
486
487        IF( .NOT.ysinus )  THEN
488           IF( ysinuss )     THEN
489              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
490              PRINT *,' *** ysinus lu sur le fichier start est F ',
491     *       'alors  qu il est  T  sur  run.def  ***'
492              STOP
493           ENDIF
494        ELSE
495           IF( .NOT.ysinuss )   THEN
496              PRINT *,' ********  PBS DANS  DEFRUN  ******** '
497              PRINT *,' ***  ysinus lu sur le fichier start est T ',
498     *        'alors  qu il est  F  sur  run.def  ****  '
499              STOP
500           ENDIF
501        ENDIF
502      ENDIF
503c
504
505      write(tapeout,*)' #########################################'
506      write(tapeout,*)' Configuration des parametres du gcm: '
507      write(tapeout,*)' dayref = ', dayref
508      write(tapeout,*)' anneeref = ', anneeref
509      write(tapeout,*)' nday = ', nday
510      write(tapeout,*)' day_step = ', day_step
511      write(tapeout,*)' iperiod = ', iperiod
512      write(tapeout,*)' iconser = ', iconser
513      write(tapeout,*)' iecri = ', iecri
514      write(tapeout,*)' periodav = ', periodav
515      write(tapeout,*)' idissip = ', idissip
516      write(tapeout,*)' lstardis = ', lstardis
517      write(tapeout,*)' nitergdiv = ', nitergdiv
518      write(tapeout,*)' nitergrot = ', nitergrot
519      write(tapeout,*)' niterh = ', niterh
520      write(tapeout,*)' tetagdiv = ', tetagdiv
521      write(tapeout,*)' tetagrot = ', tetagrot
522      write(tapeout,*)' tetatemp = ', tetatemp
523      write(tapeout,*)' coefdis = ', coefdis
524      write(tapeout,*)' purmats = ', purmats
525      write(tapeout,*)' physic = ', physic
526      write(tapeout,*)' iphysiq = ', iphysiq
527      write(tapeout,*)' ecritphy = ', ecritphy
528      write(tapeout,*)' cycle_diurne = ', cycle_diurne
529      write(tapeout,*)' soil_model = ', soil_model
530      write(tapeout,*)' new_oliq = ', new_oliq
531      write(tapeout,*)' ok_orodr = ', ok_orodr
532      write(tapeout,*)' ok_orolf = ', ok_orolf
533      write(tapeout,*)' ok_limitvrai = ', ok_limitvrai
534      write(tapeout,*)' nbapp_rad = ', nbapp_rad
535      write(tapeout,*)' iflag_con = ', iflag_con
536      write(tapeout,*)' clonn = ', clonn
537      write(tapeout,*)' clatt = ', clatt
538      write(tapeout,*)' grossismx = ', grossismx
539      write(tapeout,*)' grossismy = ', grossismy
540      write(tapeout,*)' fxyhypbb = ', fxyhypbb
541      write(tapeout,*)' dzoomxx = ', dzoomxx
542      write(tapeout,*)' dzoomy = ', dzoomy
543
544      RETURN
545c   ...............................................
546c
547100   CONTINUE
548!Config  Key  = clon
549!Config  Desc = centre du zoom, longitude
550!Config  Def  = 0
551!Config  Help = longitude en degres du centre
552!Config         du zoom
553       clon = 0.
554       CALL getin('clon',clon)
555
556!Config  Key  = clat
557!Config  Desc = centre du zoom, latitude
558!Config  Def  = 0
559!Config  Help = latitude en degres du centre du zoom
560!Config         
561       clat = 0.
562       CALL getin('clat',clat)
563
564!Config  Key  = grossismx
565!Config  Desc = zoom en longitude
566!Config  Def  = 1.0
567!Config  Help = facteur de grossissement du zoom,
568!Config         selon la longitude
569       grossismx = 1.0
570       CALL getin('grossismx',grossismx)
571
572!Config  Key  = grossismy
573!Config  Desc = zoom en latitude
574!Config  Def  = 1.0
575!Config  Help = facteur de grossissement du zoom,
576!Config         selon la latitude
577       grossismy = 1.0
578       CALL getin('grossismy',grossismy)
579
580      IF( grossismx.LT.1. )  THEN
581        PRINT *,' ***  ATTENTION !! grossismx < 1 .   *** '
582         STOP
583      ELSE
584         alphax = 1. - 1./ grossismx
585      ENDIF
586
587
588      IF( grossismy.LT.1. )  THEN
589        PRINT *,' ***  ATTENTION !! grossismy < 1 .   *** '
590         STOP
591      ELSE
592         alphay = 1. - 1./ grossismy
593      ENDIF
594
595      PRINT *,' alphax alphay defrun ',alphax,alphay
596c
597c    alphax et alphay sont les anciennes formulat. des grossissements
598c
599c
600
601!Config  Key  = fxyhypb
602!Config  Desc = Fonction  hyperbolique
603!Config  Def  = y
604!Config  Help = Fonction  f(y)  hyperbolique  si = .true. 
605!Config         sinon  sinusoidale
606       fxyhypb = .TRUE.
607       CALL getin('fxyhypb',fxyhypb)
608
609!Config  Key  = dzoomx
610!Config  Desc = extension en longitude
611!Config  Def  = 0
612!Config  Help = extension en longitude  de la zone du zoom 
613!Config         ( fraction de la zone totale)
614       dzoomx = 0.0
615       CALL getin('dzoomx',dzoomx)
616
617!Config  Key  = dzoomy
618!Config  Desc = extension en latitude
619!Config  Def  = 0
620!Config  Help = extension en latitude de la zone  du zoom 
621!Config         ( fraction de la zone totale)
622       dzoomy = 0.0
623       CALL getin('dzoomy',dzoomy)
624
625cc
626
627!Config  Key  = ysinus
628!Config  IF   = !fxyhypb
629!Config  Desc = Fonction en Sinus
630!Config  Def  = y
631!Config  Help = Fonction  f(y) avec y = Sin(latit.) si = .true.
632!Config         sinon y = latit.
633       ysinus = .TRUE.
634       CALL getin('ysinus',ysinus)
635c
636
637      write(tapeout,*)' #########################################'
638      write(tapeout,*)' Configuration des parametres du gcm: '
639      write(tapeout,*)' dayref = ', dayref
640      write(tapeout,*)' anneeref = ', anneeref
641      write(tapeout,*)' nday = ', nday
642      write(tapeout,*)' day_step = ', day_step
643      write(tapeout,*)' iperiod = ', iperiod
644      write(tapeout,*)' iconser = ', iconser
645      write(tapeout,*)' iecri = ', iecri
646      write(tapeout,*)' periodav = ', periodav
647      write(tapeout,*)' idissip = ', idissip
648      write(tapeout,*)' lstardis = ', lstardis
649      write(tapeout,*)' nitergdiv = ', nitergdiv
650      write(tapeout,*)' nitergrot = ', nitergrot
651      write(tapeout,*)' niterh = ', niterh
652      write(tapeout,*)' tetagdiv = ', tetagdiv
653      write(tapeout,*)' tetagrot = ', tetagrot
654      write(tapeout,*)' tetatemp = ', tetatemp
655      write(tapeout,*)' coefdis = ', coefdis
656      write(tapeout,*)' purmats = ', purmats
657      write(tapeout,*)' physic = ', physic
658      write(tapeout,*)' iphysiq = ', iphysiq
659      write(tapeout,*)' ecritphy = ', ecritphy
660      write(tapeout,*)' cycle_diurne = ', cycle_diurne
661      write(tapeout,*)' soil_model = ', soil_model
662      write(tapeout,*)' new_oliq = ', new_oliq
663      write(tapeout,*)' ok_orodr = ', ok_orodr
664      write(tapeout,*)' ok_orolf = ', ok_orolf
665      write(tapeout,*)' ok_limitvrai = ', ok_limitvrai
666      write(tapeout,*)' nbapp_rad = ', nbapp_rad
667      write(tapeout,*)' iflag_con = ', iflag_con
668      write(tapeout,*)' clonn = ', clonn
669      write(tapeout,*)' clatt = ', clatt
670      write(tapeout,*)' grossismx = ', grossismx
671      write(tapeout,*)' grossismy = ', grossismy
672      write(tapeout,*)' fxyhypbb = ', fxyhypbb
673      write(tapeout,*)' dzoomxx = ', dzoomxx
674      write(tapeout,*)' dzoomy = ', dzoomy
675c
676      RETURN
677      END
Note: See TracBrowser for help on using the repository browser.