source: trunk/MESOSCALE/LMDZ.MARS.new/myGCM/log @ 1097

Last change on this file since 1097 was 1097, checked in by aslmd, 12 years ago

MESOSCALE. cleaning bounding GCM stuff 2/2.
File size: 90.9 KB
Line 
1 iniadvtrac: Reading file traceur.def
2  Choix du shema  Van-Leer  pourle traceur no            1
3  Choix du shema  Van-Leer  pourle traceur no            2
4  Choix du shema  Van-Leer  pourle traceur no            3
5  Choix du shema  Van-Leer  pourle traceur no            4
6  Choix du shema  Van-Leer  pourle traceur no            5
7  Choix du shema  Van-Leer  pourle traceur no            6
8  Choix du shema  Van-Leer  pour la vapeur d'eau
9 dynetat0: rad,omeg,g,cpp,kappa   3397200.      7.0776514E-05   3.720000   
10   744.4990      0.2567930   
11 dynetat0: loading tracers
12 DEFRUN_NEW: reading file run.def
13 
14 Number of days to run:
15 COUNT :           53
16 getin_checkcohe : Found a problem on key iaervar
17 getin_checkcohe : The following values were encoutered :
18                 iaervar ==
19 24                                                                             
20                     
21                 iaervar ==
22 26                                                                             
23                     
24 getin_checkcohe : We will keep only the last value
25 COUNT :           59
26 getin_checkcohe : Found a problem on key callnlte
27 getin_checkcohe : The following values were encoutered :
28                 callnlte ==
29 .true.                                                                         
30                     
31                 callnlte ==
32 .false.                                                                       
33                     
34 getin_checkcohe : We will keep only the last value
35 COUNT :           61
36 getin_checkcohe : Found a problem on key nltemodel
37 getin_checkcohe : The following values were encoutered :
38                 nltemodel ==
39 2                                                                             
40                     
41                 nltemodel ==
42 0                                                                             
43                     
44 getin_checkcohe : We will keep only the last value
45 COUNT :           64
46 getin_checkcohe : Found a problem on key nircorr
47 getin_checkcohe : The following values were encoutered :
48                 nircorr ==
49 1                                                                             
50                     
51                 nircorr ==
52 0                                                                             
53                     
54 getin_checkcohe : We will keep only the last value
55 COUNT :           95
56 getin_checkcohe : Found a problem on key photochem
57 getin_checkcohe : The following values were encoutered :
58                 photochem ==
59 .true.                                                                         
60                     
61                 photochem ==
62 .false.                                                                       
63                     
64 getin_checkcohe : We will keep only the last value
65 COUNT :           97
66 getin_checkcohe : Found a problem on key callthermos
67 getin_checkcohe : The following values were encoutered :
68                 callthermos ==
69 .true.                                                                         
70                     
71                 callthermos ==
72 .false.                                                                       
73                     
74 getin_checkcohe : We will keep only the last value
75 COUNT :          100
76 getin_checkcohe : Found a problem on key callconduct
77 getin_checkcohe : The following values were encoutered :
78                 callconduct ==
79 .true.                                                                         
80                     
81                 callconduct ==
82 .false.                                                                       
83                     
84 getin_checkcohe : We will keep only the last value
85 COUNT :          102
86 getin_checkcohe : Found a problem on key callmolvis
87 getin_checkcohe : The following values were encoutered :
88                 callmolvis ==
89 .true.                                                                         
90                     
91                 callmolvis ==
92 .false.                                                                       
93                     
94 getin_checkcohe : We will keep only the last value
95 COUNT :          104
96 getin_checkcohe : Found a problem on key callmoldiff
97 getin_checkcohe : The following values were encoutered :
98                 callmoldiff ==
99 .true.                                                                         
100                     
101                 callmoldiff ==
102 .false.                                                                       
103                     
104 getin_checkcohe : We will keep only the last value
105 COUNT :          106
106 getin_checkcohe : Found a problem on key thermochem
107 getin_checkcohe : The following values were encoutered :
108                 thermochem ==
109 .true.                                                                         
110                     
111                 thermochem ==
112 .false.                                                                       
113                     
114 getin_checkcohe : We will keep only the last value
115 COUNT :          108
116 getin_checkcohe : Found a problem on key calleuv
117 getin_checkcohe : The following values were encoutered :
118                 calleuv ==
119 .true.                                                                         
120                     
121                 calleuv ==
122 .false.                                                                       
123                     
124 getin_checkcohe : We will keep only the last value
125  nday =    3.000000   
126 USING DEFAULTS : ndynstep =       -9999
127 
128 Number of dynamical steps per day:(should be a multiple of iperiod)
129  day_step =          960
130 
131 periode pour le pas Matsuno (en pas)
132  iperiod =            5
133 
134 periode de sortie des variables de controle (en pas)
135  iconser =          120
136 
137 periode de la dissipation (en pas)
138  idissip =            1
139 
140 choix de l'operateur de dissipation (star ou  non star )
141  lstardis =  T
142 
143 avec ou sans coordonnee hybrides
144  hybrid =  T
145 
146 nombre d'iterations de l'operateur de dissipation   gradiv
147  nitergdiv =            1
148 
149 nombre d'iterations de l'operateur de dissipation  nxgradrot
150  nitergrot =            2
151 
152 nombre d'iterations de l'operateur de dissipation  divgrad
153  niterh =            2
154 
155 temps de dissipation des plus petites long.d ondes pour u,v (gradiv)
156  tetagdiv =    2000.000   
157 
158 temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
159  tetagrot =    6000.000   
160 
161 temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
162  tetatemp =    6000.000   
163 
164 coefficient pour gamdissip
165  coefdis =   0.0000000E+00
166 
167 choix du shema d'integration temporelle
168 (true==Matsuno ou false==Matsuno-leapfrog)
169  purmats =  F
170 
171 avec ou sans physique
172  physic =  T
173 
174 periode de la physique (en pas)
175  iphysiq =           10
176 
177 choix d'une grille reguliere
178  grireg =  T
179 
180 frequence (en pas) de l'ecriture du fichier diagfi.nc
181  ecritphy =    40.00000   
182 
183 longitude en degres du centre du zoom
184  clon =    63.00000   
185 
186 latitude en degres du centre du zoom
187  clat =   0.0000000E+00
188 
189 facteur de grossissement du zoom, selon longitude
190  grossismx =    1.000000   
191 
192 facteur de grossissement du zoom , selon latitude
193  grossismy =    1.000000   
194  Defrun  alphax alphay    0.0000000E+00  0.0000000E+00
195 
196 Fonction  f(y)  hyperbolique  si = .true., sinon  sinusoidale
197  fxyhypb =  F
198 
199 extension en longitude de la zone du zoom (fraction de la zone totale)
200  dzoomx =   0.0000000E+00
201 
202 extension en latitude de la zone du zoom (fraction de la zone totale)
203  dzoomy =   0.0000000E+00
204 
205 raideur du zoom en  X
206  taux =    2.000000   
207 
208 raideur du zoom en  Y
209  tauy =    2.000000   
210 
211 Avec sponge layer
212  callsponge =  T
213 
214 Sponge: number of layers over which sponge extends
215 USING DEFAULTS : nsponge =           3
216  nsponge =            3
217 
218 Sponge mode: (forcing is towards ...
219   over upper nsponge layers)
220   0: (h=hmean,u=v=0)
221   1: (h=hmean,u=umean,v=0)
222   2: (h=hmean,u=umean,v=vmean)
223  mode_sponge =            2
224 
225 Sponge: characteristice time scale tetasponge
226 (seconds) at topmost layer (time scale then
227  doubles with decreasing layer index).
228  tetasponge =    50000.00   
229 -----------------------------------------------
230 
231 
232 -----------------------------------------------
233 
234 
235 No zoom ? -> fxyhypb set to False           (It runs better that way)
236  R  CPP  Kappa    191.1821       744.4990      0.2567930   
237 ****************************
238 Lecture de z2sig.def
239 ****************************
240 *******************************
241 Systeme en coordonnees hybrides
242 
243  BP
244   1.000000      0.9990893      0.9961299      0.9887227      0.9736167   
245  0.9470957      0.9053387      0.8451603      0.7657363      0.6698347   
246  0.5632617      0.4541319      0.3506572      0.2588113      0.1820027   
247  0.1212465      7.5687021E-02  4.3286249E-02  2.1493321E-02  7.8194505E-03
248  6.6300534E-04  3.9641478E-11  0.0000000E+00  0.0000000E+00  0.0000000E+00
249  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
250  AP
251  0.0000000E+00  9.1087818E-03  3.8738251E-02  0.1130879      0.2655351   
252  0.5357349      0.9672797       1.600810       2.453982       3.502020   
253   4.676154       5.868705       6.962464       7.865741       8.523151   
254   8.916158       9.054893       8.963833       8.661031       8.111788   
255   6.921022       4.003914       2.106964       1.072680      0.5344583   
256  0.2654041      0.1317958      6.5447852E-02  3.2500453E-02  0.0000000E+00
257  BPs
258  0.9995446      0.9976096      0.9924263      0.9811697      0.9603562   
259  0.9262172      0.8752495      0.8054483      0.7177855      0.6165482   
260  0.5086969      0.4023946      0.3047343      0.2204070      0.1516246   
261  9.8466769E-02  5.9486635E-02  3.2389786E-02  1.4656385E-02  4.2412281E-03
262  3.3150270E-04  1.9820739E-11  0.0000000E+00  0.0000000E+00  0.0000000E+00
263  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
264  APs
265  4.5543909E-03  2.3923516E-02  7.5913072E-02  0.1893115      0.4006350   
266  0.7515073       1.284045       2.027396       2.978001       4.089087   
267   5.272429       6.415585       7.414103       8.194447       8.719654   
268   8.985525       9.009363       8.812431       8.386410       7.516405   
269   5.462468       3.055439       1.589822      0.8035691      0.3999312   
270  0.1985999      9.8621808E-02  4.8974153E-02  2.4319850E-02
271  PRESNIVS
272   609.7268       608.5658       605.4559       598.7028       586.2180   
273   565.7440       535.1863       493.3509       440.8271       380.1835   
274   315.5775       251.8763       193.3020       142.6427       101.2106   
275   69.05025       45.29621       28.57020       17.32681       10.10355   
276   5.664685       3.055439       1.589822      0.8035691      0.3999312   
277  0.1985999      9.8621808E-02  4.8974153E-02  2.4319850E-02
278 Pseudo altitude des Presnivs :
279  4.4802926E-03  2.3539927E-02  7.4772738E-02  0.1869370      0.3976738   
280  0.7531729       1.308441       2.122384       3.248062       4.728050   
281   6.590548       8.845211       11.49205       14.53116       17.96255   
282   21.78624       26.00236       30.61095       35.61204       41.00571   
283   46.79208       52.96536       59.49837       66.32151       73.29922   
284   80.29922       87.29922       94.29922       101.2992   
285
286
287           ....  INIGEOM  date du 01/06/98   .....
288
289        Calcul des elongations cu et cv  comme sommes des 4
290      elong. cuij1, .. 4  , cvij1,.. 4  qui les entourent , aux 
291      memes endroits que les aires aireij1,...j4   .
292
293  gamdi_gd   0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
294           1           2           2
295
296
297
298 *** Inigeom ,  Y = Latitude  , der. sinusoid . ***
299    ***  Coordonnees de la grille  ***
300
301
302    LONGITUDES  aux pts.   V  ( degres ) 
303
304
305  -180.00 -174.38 -168.75 -163.12 -157.50 -151.88 -146.25 -140.62
306  -135.00 -129.38 -123.75 -118.12 -112.50 -106.88 -101.25  -95.63
307   -90.00  -84.37  -78.75  -73.12  -67.50  -61.87  -56.25  -50.63
308   -45.00  -39.38  -33.75  -28.12  -22.50  -16.88  -11.25   -5.63
309     0.00    5.63   11.25   16.88   22.50   28.12   33.75   39.38
310    45.00   50.62   56.25   61.87   67.50   73.12   78.75   84.38
311    90.00   95.63  101.25  106.88  112.50  118.12  123.75  129.38
312   135.00  140.62  146.25  151.88  157.50  163.13  168.75  174.38
313   180.00
314
315
316    LATITUDES   aux pts.   V  ( degres ) 
317
318
319    88.12   84.37   80.62   76.87   73.12   69.37   65.62   61.87
320    58.12   54.37   50.62   46.87   43.12   39.37   35.62   31.87
321    28.12   24.37   20.62   16.87   13.12    9.37    5.62    1.87
322    -1.88   -5.63   -9.38  -13.13  -16.88  -20.63  -24.38  -28.13
323   -31.88  -35.63  -39.38  -43.13  -46.88  -50.62  -54.38  -58.13
324   -61.88  -65.62  -69.38  -73.12  -76.88  -80.62  -84.37  -88.12
325
326
327    LONGITUDES  aux pts.   U  ( degres ) 
328
329
330  -177.19 -171.56 -165.94 -160.31 -154.69 -149.06 -143.44 -137.81
331  -132.19 -126.56 -120.94 -115.31 -109.69 -104.06  -98.44  -92.81
332   -87.19  -81.56  -75.94  -70.31  -64.69  -59.06  -53.44  -47.81
333   -42.19  -36.56  -30.94  -25.31  -19.69  -14.06   -8.44   -2.81
334     2.81    8.44   14.06   19.69   25.31   30.94   36.56   42.19
335    47.81   53.44   59.06   64.69   70.31   75.94   81.56   87.19
336    92.81   98.44  104.06  109.69  115.31  120.94  126.56  132.19
337   137.81  143.44  149.06  154.69  160.31  165.94  171.56  177.19
338   182.81
339
340
341    LATITUDES   aux pts.   U  ( degres ) 
342
343
344    90.00   86.25   82.50   78.75   75.00   71.25   67.50   63.75
345    60.00   56.25   52.50   48.75   45.00   41.25   37.50   33.75
346    30.00   26.25   22.50   18.75   15.00   11.25    7.50    3.75
347     0.00   -3.75   -7.50  -11.25  -15.00  -18.75  -22.50  -26.25
348   -30.00  -33.75  -37.50  -41.25  -45.00  -48.75  -52.50  -56.25
349   -60.00  -63.75  -67.50  -71.25  -75.00  -78.75  -82.50  -86.25
350   -90.00
351
352
353  EIGNVL
354 -0.940309E-05-0.999202E+00-0.999209E+00-0.398716E+01-0.398718E+01
355 -0.893513E+01-0.893514E+01-0.157954E+02-0.157955E+02-0.245020E+02
356 -0.245020E+02-0.349710E+02-0.349710E+02-0.471016E+02-0.471017E+02
357 -0.607770E+02-0.607770E+02-0.758655E+02-0.758655E+02-0.922217E+02
358 -0.922218E+02-0.109688E+03-0.109688E+03-0.128097E+03-0.128097E+03
359 -0.147270E+03-0.147270E+03-0.167023E+03-0.167023E+03-0.187167E+03
360 -0.187167E+03-0.207506E+03-0.207506E+03-0.227845E+03-0.227845E+03
361 -0.247988E+03-0.247988E+03-0.267742E+03-0.267742E+03-0.286915E+03
362 -0.286915E+03-0.305323E+03-0.305323E+03-0.322790E+03-0.322790E+03
363 -0.339146E+03-0.339146E+03-0.354234E+03-0.354235E+03-0.367910E+03
364 -0.367910E+03-0.380040E+03-0.380041E+03-0.390509E+03-0.390510E+03
365 -0.399216E+03-0.399216E+03-0.406076E+03-0.406076E+03-0.411024E+03
366 -0.411024E+03-0.414012E+03-0.414012E+03-0.415012E+03
367
368                Inifilr colat0 alphax    0.5000000E+00   0.0000000E+00
369  TRUNCATION AT           64
370  jfiltnv jfiltsv jfiltnu jfiltsu            8          41           8
371          41
372    Modes premiers  v 
373   4 10 14 20 26 32 40 52 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64
374  64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 52 40 32 26 20 14 10  4
375    Modes premiers  u 
376  64  6 12 18 24 30 36 46 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64
377  64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 46 36 30 24 18 12  6
378 Dissipation :
379 Multiplication de la dissipation en altitude : fac_mid, fac_up =   3.000000   
380   30.00000   
381 Transition mid /up : startalt, delta =   70.00000       20.00000     (km)
382 calcul des valeurs propres
383 Constantes de temps de la diffusion horizontale
384  INIDI tetamin dtvr    66.66684       92.47396               5
385  INIDI tetamin idissip    66.66684               1
386  INIDI idissip dtdiss            1   92.47396   
387 pseudoZ(km)  zvert  dt(tetagdiv) dt(tetagrot) dt(divgrad)
388  4.4802926E-03   1.000001      4.6237011E-02  1.5412334E-02  1.5412334E-02
389  2.3539927E-02   1.000011      4.6237506E-02  1.5412500E-02  1.5412500E-02
390  7.4772738E-02   1.000113      4.6242204E-02  1.5414066E-02  1.5414066E-02
391  0.1869370       1.000712      4.6269897E-02  1.5423298E-02  1.5423298E-02
392  0.3976738       1.003286      4.6388932E-02  1.5462974E-02  1.5462974E-02
393  0.7531729       1.012164      4.6799421E-02  1.5599805E-02  1.5599805E-02
394   1.308441       1.038334      4.8009411E-02  1.6003134E-02  1.6003134E-02
395   2.122384       1.105890      5.1133037E-02  1.7044343E-02  1.7044343E-02
396   3.248062       1.256737      5.8107711E-02  1.9369235E-02  1.9369235E-02
397   4.728050       1.535234      7.0984595E-02  2.3661530E-02  2.3661530E-02
398   6.590548       1.930723      8.9270793E-02  2.9756924E-02  2.9756924E-02
399   8.845211       2.338095      0.1081065      3.6035489E-02  3.6035489E-02
400   11.49205       2.645828      0.1223351      4.0778361E-02  4.0778361E-02
401   14.53116       2.829569      0.1308307      4.3610230E-02  4.3610230E-02
402   17.96255       2.923871      0.1351909      4.5063645E-02  4.5063645E-02
403   21.78624       2.967935      0.1372284      4.5742787E-02  4.5742787E-02
404   26.00236       2.987214      0.1381198      4.6039917E-02  4.6039917E-02
405   30.61095       2.995183      0.1384882      4.6162728E-02  4.6162728E-02
406   35.61204       2.998292      0.1386320      4.6210654E-02  4.6210654E-02
407   41.00571       2.999433      0.1386847      4.6228237E-02  4.6228237E-02
408   46.79208       2.999824      0.1387028      4.6234269E-02  4.6234269E-02
409   52.96536       2.999952      0.1387087      4.6236239E-02  4.6236239E-02
410   59.49837       3.000110      0.1387160      4.6238672E-02  4.6238672E-02
411   66.32151       3.007357      0.1390511      4.6350367E-02  4.6350367E-02
412   73.29922       3.475960      0.1607179      5.3572625E-02  5.3572625E-02
413   80.29922       17.70859      0.8187916      0.2729305      0.2729305   
414   87.29922       29.66585       1.371659      0.4572197      0.4572197   
415   94.29922       29.99493       1.386875      0.4622915      0.4622915   
416   101.2992       29.99992       1.387106      0.4623686      0.4623686   
417 PRESSION SURFACE
418 MIN:    162.2194   
419 MAX:    1119.796   
42055555555555555555555555555555555555555555555555555555555555555555
42166666666566655555566666666677766666666666666666666666666666666666
42266666667777766677777777888777777776666677777777777777777777777766
42377777777777777777888888888888888888888888888888887777777777777777
42477777777777777777778888888888888888888888888888877777777777777777
42577777777776766677777888888888888888888888888888777777777777777777
42677777766666666666777888888888888888888888888887777777777777777777
42777776666666666666667778888888888888888888888777777777777777777777
42877776666666556556666778888887777777778888887777777777777777777777
42977776666655555555666777888877777777777777777777777777777777777777
43077777666555555555566667788877777777777777777777777777777777777777
43177777766555544444555556777777777777777777777777777777777777777777
43277777766554333444444556777777777776767766666777777777787777776777
43377777776554222344444455777777777666666655566667777777787777776777
43477777776544322344444556777777776666665554456666777777777777766777
43577777776544333444444555777777766665555555555556677777777776666777
43677777765554333344445566777777666666555555544455677777777754556777
43777777765554333334445566777777665666555555544445667777777643456677
43877777765244333334445555777776656556555554444445666677777643456667
43966677764024333334454555677766555555555544444446666666777644556666
44066666665233333233454455667666565555555544444446776666666655566666
44166666665443331233454455567666555555554544443346776555666666666666
44266666665443332233454444556766655555544444443445776554566666666666
44366665655432222223443444455666655555544444443344665554456666666666
44466665655432212223333344455655655555444444343344444454445666666666
44555665554332112223324344445565655555544444333334344444444655666665
44655555544331111123345333555665555554443343333334334444434455556655
44755454443321112122235433445555555544443333334344333444333444455555
44844444443321122122233334664555555444433333344444333444334343445554
44943444443332222222333333444555555444443333344444333444333333344444
45044444443332222222333333444455554444443334445554443334433334444444
45143333433333222223333334444544444444344334455565544444433334444444
45243333443333222223332334444444444444344444468986655444443334444444
453443334433333222223233334444444444444334446***98766544443344444444
454444334433333332222233334455444444444334459***99876544444443444444
455444443333333333222333344555444444443334469***99986544444334443444
456334433333333333333333445775444444443334459**999986544444334444333
45733443333334333333333445677654444443334445799999975444444434444333
45844443433334443333444445666544444443333445678888654444444444443344
45944433333333343333444345565444443333333444455565554444444444333334
46043333334333343334444444444443333333333344444445443433344444333334
46133344443333333334444444444443333333333333444444443333344444333333
46233444433333333334444444443333333333333333344444333334444444433333
46333443333333333334344444333333333333333333333333333334444444434333
46433333333333333333333333333333333333333333333333333333444333333333
46533333333333333333333333333333333332333333333333333333333333333333
46633333333333333333333333333333333332223333333333333333333333333333
46733332222223333332222222222222222222222222222222223322222222232223
46822222222222222222222222222222222222222222222222222222222222222222
4691
470               run du pas      0  au pas   2880  c'est a dire du jour  13380   au jour  13383
471
472
473 dynredem0: iim,jjm,llm,idayref          64          48          29       13383
474 dynredem0: rad,omeg,g,cpp,kappa   3397200.      7.0776514E-05   3.720000   
475   744.4990      0.2567930   
476 gcm: callgroupeun set to: F
477  GCM ---- Test_period apres continue   OK ! ----- itau:            0
4780**********    pas      0     jour*****heure  0.0    date   0.00000    **********
479          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
480        0.947077 0.976960E+00     0.951     1.085     0.966    18.390     0.825
481
482 
483 sponge mode           2
484 nsponge tetasponge            3   50000.00   
485 Coeffs for the sponge layer
486 Z (km)     tau      cst
487   87.29922       200000.0      4.6236982E-04
488   94.29921       100000.0      9.2473964E-04
489   105.3311       50000.00      1.8494793E-03
490 
491 in groupeun.F, ngroup=           2
492 Shema  Amont nouveau  appele dans  Vanleer   
493 
494 
495 --------------------------------------------
496  inifis: Parameters for the physics (callphys.def)
497 --------------------------------------------
498 Directory where external input files are:
499  datafile = /planeto/emlmd/run_mcd5.v0_new/datadir
500 Run with or without tracer transport ?
501  tracer =  T
502 Diurnal cycle ?
503 (if diurnal=False, diurnal averaged solar heating)
504  diurnal =  T
505 Seasonal cycle ?
506 (if season=False, Ls stays constant, to value set in 'start'
507  season =  T
508 Write some extra output to the screen ?
509  lwrite =  F
510 Save statistics in file stats.nc ?
511  callstats =  T
512 Save EOF profiles in file 'profiles' for Climate Database?
513  calleofdump =  F
514 Dust scenario: 1=constant dust (read from startfi
515  or set as tauvis); 2=Viking scenario; =3 MGS scenario,
516 =6 cold (low dust) scenario; =7 warm (high dust) scenario
517 =24,25 ... 30 :Mars Year 24, ... or 30 from TES assimilation
518  iaervar =           26
519 Reference (visible) dust opacity at 610 Pa (matters only if iaervar=1)
520  tauvis =   0.1000000   
521 Dust vertical distribution:
522 (=1 top set by topdustref parameter; =2 Viking scenario; =3 MGS scenario)
523  iddist =            3
524 Dust top altitude (km). (Matters only if iddist=1)
525  topdustref =    55.00000   
526 Prescribed surface thermal flux (H/(rho*cp),K m/s)
527 USING DEFAULTS : tke_heat_flux =  0.0000000E+00
528  tke_heat_flux =   0.0000000E+00
529  0 means the usual schemes are computing
530 call radiative transfer ?
531  callrad =  T
532 call slope insolation scheme ?(matters only if callrad=T)
533 USING DEFAULTS : callslope = F
534  callslope =  F
535 call NLTE radiative schemes ?(matters only if callrad=T)
536  callnlte =  F
537 NLTE model?
538 0 -> old model, static O
539 1 -> old model, dynamic O
540 2 -> new model
541 (matters only if callnlte=T)
542  nltemodel =            0
543 call CO2 NIR absorption ?(matters only if callrad=T)
544  callnirco2 =  T
545 New NIR NLTE correction ?0-> old model (no correction)1-> new correction
546 (matters only if callnirco2=T)
547  nircorr =            0
548 call turbulent vertical diffusion ?
549  calldifv =  T
550 call thermals ?
551  calltherm =  T
552 output thermal diagnostics ?
553 USING DEFAULTS : outptherm = F
554  outptherm =  F
555 call convective adjustment ?
556  calladj =  T
557 call Richardson-based surface layer ?
558  callrichsl =  T
559 call CO2 condensation ?
560  callcond =  T
561 call thermal conduction in the soil ?
562  callsoil =  T
563 call Lott's gravity wave/subgrid topography scheme ?
564  calllott =  T
565 rad.transfer is computed every iradia physical timestep
566  iradia =            1
567 Output of the exchange coefficient mattrix ?(for diagnostics only)
568  callg2d =  F
569 Rayleigh scattering : (should be .false. for now)
570  rayleigh =  F
571 Transported dust ? (if >0, use 'dustbin' dust bins)
572  dustbin =            2
573 Radiatively active dust ? (matters if dustbin>0)
574  active =  T
575 use mass and number mixing ratios to predict dust size ?
576  doubleq =  T
577 USING DEFAULTS : submicron = F
578  submicron =  F
579 dust lifted by GCM surface winds ?
580  lifting =  T
581  dust lifted by dust devils ?
582  callddevil =  F
583 Gravitationnal sedimentation ?
584  sedimentation =  T
585 Radiatively active transported atmospheric water ice ?
586  activice =  T
587 Compute water cycle ?
588  water =  T
589 Activate the thermal inertia feedback ?
590 USING DEFAULTS : tifeedback = F
591  tifeedback =  F
592 effective variance for water ice clouds ?
593 water_param nueff Sedimentation:  0.1000000   
594 water param CCN reduc. factor ?  0.0000000E+00
595 USING DEFAULTS : ccn_factor =   4.500000   
596  ccn_factor =    4.500000   
597 Careful: only used when microphys=F, otherwise
598 the contact parameter is used instead;
599 Microphysical scheme for water-ice clouds?
600  microphys =  T
601 water contact parameter ?
602 mteta =   0.9500000   
603 Dust scavenging by H2O/CO2 snowfall ?
604  scavenging =  T
605 Permanent water caps at poles ? .true. is RECOMMENDED
606 (with .true., North cap is a source of water and South pole is a cold trap)
607  caps =  T
608 water ice albedo ?
609  albedo_h2o_ice =   0.3500000   
610 water ice thermal inertia ?
611  inert_h2o_ice =    800.0000   
612 frost thickness threshold for albedo ?
613  frost_albedo_threshold =   4.9999999E-03
614 Titus crocus line ?
615 USING DEFAULTS : tituscap = F
616 tituscap F
617 photochemistry: include chemical species
618  photochem =  F
619 call thermosphere ?
620  callthermos =  F
621  water included without cycle (only if water=.false.)
622  thermoswater =  F
623 call thermal conduction ? (only if callthermos=.true.)
624  callconduct =  F
625 call EUV heating ? (only if callthermos=.true.)
626  calleuv =  F
627 call molecular viscosity ? (only if callthermos=.true.)
628  callmolvis =  F
629 call molecular diffusion ? (only if callthermos=.true.)
630  callmoldiff =  F
631 call thermospheric photochemistry ? (only if callthermos=.true.)
632  thermochem =  F
633 Method to include solar variability
634 0-> old method (using solarcondate); 1-> variability wit E10.7
635  solvarmod =            0
636 date for solar flux calculation: (1985 < date < 2002)
637  (Only used if solvarmod=0)
638 (Solar min=1996.4 ave=1993.4 max=1990.6)
639  solarcondate =    1993.400   
640 Solar variability as observed for MY:
641 Only if solvarmod=1
642  solvaryear =           24
643 UV heating efficiency:measured values between 0.19 and 0.23 (Fox et al. 1996)
644 lower values may be used to compensate low 15 um cooling
645  euveff =   0.2100000   
646 --------------------------------------------
647 
648 
649 
650 inifis: daysec   88775.00   
651 
652 inifis: The radiative transfer is computed:
653            each            1  physical time-step
654         or each    924.7396      seconds
655 
656 
657 TABFI de phyeta0           0           0
658 tabfi: tab_cntrl   3010.000       29.00000       13380.00      0.0000000E+00
659   3397200.      7.0776514E-05   3.720000       43.49000      0.2567930   
660   88775.00       924.7396      0.0000000E+00  0.0000000E+00   669.0000   
661   206.6600       249.2200       485.0000       25.19000      9.9999998E-03
662   30.00000      1.0000000E-06  0.5800000      0.5600000      0.9000000   
663  0.9000000      0.9500000      0.1000000      0.0000000E+00  0.0000000E+00
664  0.0000000E+00  9.9999997E-05  9.9999997E-05  0.5000000      0.5000000   
665   1000000.      0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
666  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
667  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
668  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
669  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
670  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
671  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
672  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
673  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
674  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
675  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
676  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
677  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
678  0.0000000E+00
679 IN tabfi day_ini=       13380
680 *****************************************************
681 Reading tab_cntrl when calling tabfi before changes
682 *****************************************************
683 (1)      = ngridmx?     3010.00     3010.00
684 (2)            lmax       29.00       29.00
685 (3)         day_ini    13380.00    13380.00
686 (5)             rad  3397200.00  3397200.00
687 (10)         daysec    88775.00    88775.00
688 (6)            omeg   0.707765E-04   0.707765E-04
689 (7)               g        3.72        3.72
690 (8)           mugaz       43.49       43.49
691 (9)             rcp        0.26        0.26
692 (11)        dtphys?   0.924740E+03   0.924740E+03
693 (14)       year_day      669.00      669.00
694 (15)       periheli      206.66      206.66
695 (16)        aphelie      249.22      249.22
696 (17)       peri_day      485.00      485.00
697 (18)       obliquit       25.19       25.19
698 (19)     z0_default   0.100000E-01   0.100000E-01
699 (21)      emin_turb   0.100000E-05   0.100000E-05
700 (20)        lmixmin       30.00       30.00
701 (26)        emissiv        0.95        0.95
702 (24)     emisice(1)        0.90        0.90
703 (25)     emisice(2)        0.90        0.90
704 (22)    albedice(1)        0.58        0.58
705 (23)    albedice(2)        0.56        0.56
706 (31)   iceradius(1)   0.100000E-03   0.100000E-03
707 (32)   iceradius(2)   0.100000E-03   0.100000E-03
708 (33)   dtemisice(1)        0.50        0.50
709 (34)   dtemisice(2)        0.50        0.50
710 (27)         tauvis        0.10        0.10
711 (35)        volcapa  1000000.00  1000000.00
712 
713 Lmodif in tabfi!!!!!!!           0
714 Geopotentiel au sol <phisfi>:  -26336.29       52485.20   
715 Albedo du sol nu <albedodat>:  0.1003705      0.3208752   
716 <zmea>:  -7079.710       14109.67   
717 <zstd>:   8.602325       5754.301   
718 <zsig>:  1.4299675E-04  6.9831863E-02
719 <zgam>:  1.0000000E-07   1.000000   
720 <zthe>:  -89.90435       89.88921   
721 CO2 ice cover <co2ice>:  0.0000000E+00   958.8163   
722 phyetat0: Le champ <tsurf> est present
723           J ignore donc les autres temperatures TS**
724 Temperature du sol <tsurf>   143.1560       291.0728   
725 Surface emissivity <emis>:  0.5347768      0.9500000   
726 Surface roughness <z0>:  7.9248937E-05  1.7229829E-02
727 pbl wind variance <q2>:  8.3298417E-08   157.1757   
728 tracer on surface <co2>:  0.0000000E+00  0.0000000E+00
729 tracer on surface <dust_number>: -2.1804593E+16  3.8433894E+14
730 tracer on surface <dust_mass>:  -1790.340       634.2337   
731 tracer on surface <ccn_number>:  1.5558199E+09  2.8858525E+16
732 tracer on surface <ccn_mass>:  4.8576188E-05   910.5092   
733 tracer on surface <h2o_ice>:  -40.61012       13.81997   
734 phyetat0: loading surface tracer h2o_ice instead of h2o_vap
735 tracer on surface <h2o_ice>:  -40.61012       13.81997   
736 
737 Soil volumetric heat capacity:   1000000.   
738 Soil thermal inertia <inertiedat>:   40.57156       1200.000   
739 Soil temperature <tsoil>:   143.1567       289.7982   
740 In physiq day_ini =       13380
741 initracer: found all expected tracers, namely:
742                  1  co2
743                  2  dust_number
744                  3  dust_mass
745                  4  ccn_number
746                  5  ccn_mass
747                  6  h2o_ice
748                  7  h2o_vap
749 initracer: doubleq_param reff_lift:  3.0000001E-06
750 initracer: doubleq_param nueff_lift:  0.5000000   
751 initracer: doubleq_param alpha_lift:  1.0000000E-06
752 
753 ******** initracer : dust transport parameters :
754 alpha_lift =   0.0000000E+00  1.1936620E+07  1.0000000E-06  1.0000000E-30
755  1.0000000E-30  0.0000000E+00  0.0000000E+00
756 alpha_devil =   0.0000000E+00   107429.6      8.9999999E-09  1.0000000E-30
757  1.0000000E-30  0.0000000E+00  0.0000000E+00
758 radius  =   0.0000000E+00  3.0000001E-06  3.0000001E-06  3.0000001E-06
759  3.0000001E-06  3.0000001E-06  0.0000000E+00
760 reff_lift (um) =    3.0000001E-06
761 size distribution variance  =   0.6367614   
762 r3n_q , ref_r0 :   1.5401418E-05   2.755676   
763 Ice dryness ?
764 USING DEFAULTS : icedryness =   1.000000   
765  icedryness =    1.000000   
766 Surfini: predefined ice caps
767 water ice caps distribution for 64x48 resolution
768 latitude | longitude | ig
769 ice water cap   90.00000      0.0000000E+00           1
770 ice water cap   86.24999      -180.0000               2
771 ice water cap   86.24999      -174.3750               3
772 ice water cap   86.24999      -168.7500               4
773 ice water cap   86.24999      -163.1250               5
774 ice water cap   86.24999      -157.5000               6
775 ice water cap   86.24999      -151.8750               7
776 ice water cap   86.24999      -146.2500               8
777 ice water cap   86.24999      -140.6250               9
778 ice water cap   86.24999      -135.0000              10
779 ice water cap   86.24999      -129.3750              11
780 ice water cap   86.24999      -123.7500              12
781 ice water cap   86.24999      -118.1250              13
782 ice water cap   86.24999      -112.5000              14
783 ice water cap   86.24999      -106.8750              15
784 ice water cap   86.24999      -101.2500              16
785 ice water cap   86.24999      -95.62501              17
786 ice water cap   86.24999      -90.00000              18
787 ice water cap   86.24999      -84.37499              19
788 ice water cap   86.24999      -78.75000              20
789 ice water cap   86.24999      -73.12500              21
790 ice water cap   86.24999      -67.50000              22
791 ice water cap   86.24999      -61.87499              23
792 ice water cap   86.24999      -56.24999              24
793 ice water cap   86.24999      -50.62501              25
794 ice water cap   86.24999      -45.00000              26
795 ice water cap   86.24999      -39.37500              27
796 ice water cap   86.24999      -33.75000              28
797 ice water cap   86.24999      -28.12500              29
798 ice water cap   86.24999      -22.49999              30
799 ice water cap   86.24999      -16.87501              31
800 ice water cap   86.24999      -11.25000              32
801 ice water cap   86.24999      -5.625002              33
802 ice water cap   86.24999      0.0000000E+00          34
803 ice water cap   86.24999       5.625002              35
804 ice water cap   86.24999       11.25000              36
805 ice water cap   86.24999       16.87501              37
806 ice water cap   86.24999       22.49999              38
807 ice water cap   86.24999       28.12500              39
808 ice water cap   86.24999       33.75000              40
809 ice water cap   86.24999       39.37500              41
810 ice water cap   86.24999       45.00000              42
811 ice water cap   86.24999       50.62499              43
812 ice water cap   86.24999       56.25001              44
813 ice water cap   86.24999       61.87499              45
814 ice water cap   86.24999       67.50001              46
815 ice water cap   86.24999       73.12500              47
816 ice water cap   86.24999       78.74998              48
817 ice water cap   86.24999       84.37500              49
818 ice water cap   86.24999       89.99999              50
819 ice water cap   86.24999       95.62501              51
820 ice water cap   86.24999       101.2500              52
821 ice water cap   86.24999       106.8750              53
822 ice water cap   86.24999       112.5000              54
823 ice water cap   86.24999       118.1250              55
824 ice water cap   86.24999       123.7500              56
825 ice water cap   86.24999       129.3750              57
826 ice water cap   86.24999       135.0000              58
827 ice water cap   86.24999       140.6250              59
828 ice water cap   86.24999       146.2500              60
829 ice water cap   86.24999       151.8750              61
830 ice water cap   86.24999       157.5000              62
831 ice water cap   86.24999       163.1250              63
832 ice water cap   86.24999       168.7500              64
833 ice water cap   86.24999       174.3750              65
834 ice water cap   82.49999      -129.3750              75
835 ice water cap   82.49999      -123.7500              76
836 ice water cap   82.49999      -118.1250              77
837 ice water cap   82.49999      -112.5000              78
838 ice water cap   82.49999      -106.8750              79
839 ice water cap   82.49999      -101.2500              80
840 ice water cap   82.49999      -95.62501              81
841 ice water cap   82.49999      -90.00000              82
842 ice water cap   82.49999      -84.37499              83
843 ice water cap   82.49999      -78.75000              84
844 ice water cap   82.49999      -73.12500              85
845 ice water cap   82.49999      -28.12500              93
846 ice water cap   82.49999      -22.49999              94
847 ice water cap   82.49999      -16.87501              95
848 ice water cap   82.49999      -11.25000              96
849 ice water cap   82.49999      -5.625002              97
850 ice water cap   82.49999      0.0000000E+00          98
851 ice water cap   82.49999       5.625002              99
852 ice water cap   82.49999       11.25000             100
853 ice water cap   82.49999       16.87501             101
854 ice water cap   82.49999       22.49999             102
855 ice water cap   82.49999       28.12500             103
856 ice water cap   82.49999       33.75000             104
857 ice water cap   82.49999       39.37500             105
858 ice water cap   82.49999       45.00000             106
859 ice water cap   82.49999       50.62499             107
860 ice water cap   82.49999       56.25001             108
861 ice water cap   82.49999       61.87499             109
862 ice water cap   82.49999       67.50001             110
863 ice water cap   82.49999       73.12500             111
864 ice water cap   82.49999       78.74998             112
865 ice water cap   82.49999       84.37500             113
866 ice water cap   82.49999       89.99999             114
867 ice water cap   78.75000      -146.2500             136
868 ice water cap   78.75000      -135.0000             138
869 ice water cap   78.75000      -123.7500             140
870 ice water cap   78.75000      -112.5000             142
871 ice water cap   78.75000      -5.625002             161
872 ice water cap   78.75000       5.625002             163
873 ice water cap   78.75000       16.87501             165
874 ice water cap   78.75000       118.1250             183
875 ice water cap   78.75000       129.3750             185
876 ice water cap   78.75000       140.6250             187
877 ice water cap   78.75000       151.8750             189
878 ice water cap   78.75000       163.1250             191
879 ice water cap   78.75000       174.3750             193
880 ice water cap   75.00000      -180.0000             194
881 ice water cap   75.00000      -129.3750             203
882 ice water cap   75.00000      -106.8750             207
883 ice water cap   75.00000       101.2500             244
884 ice water cap   75.00000       112.5000             246
885 ice water cap   75.00000       135.0000             250
886 ice water cap   75.00000       146.2500             252
887 ice water cap   75.00000       157.5000             254
888 ice water cap   75.00000       168.7500             256
889 surfini: water frost thickness  4.9999999E-03
890 surfini: water ice albedo:  0.3500000   
891 surfini: water ice TI:   800.0000   
892 minimum albedo sans water caps  0.1003705   
893 maximum albedo sans water caps  0.3208752   
894 minimum albedo avec water caps  0.1003705   
895 maximum albedo avec water caps  0.3500000   
896 minimum albedo avec givre et co2  0.1008788   
897 maximum albedo avec givre et co2  0.5800000   
898 Perihelie en Mkm     206.6600   
899 Aphelise  en Mkm     249.2200   
900 obliquite en degres  :   25.19000   
901 e_elips  9.3357898E-02
902 p_elips   1.510404   
903 zanom     1.728111   
904 zx0      1.818617   
905 longitude solaire du perihelie timeperi =    1.908231   
906 physiq: water_param Surface water ice albedo:  0.3500000   
907 nlthermeq: transition to NLTE centred at   0.10Pa
908               half-width (scale heights)   0.50
909          suggested LTE coverage at least   0.04Pa
910 LTE rad. calculations up to layer           29
911 callradite: found all scatterers, namely:
912                  1  dust_doubleq
913                  2  h2o_ice
914 Scatterer: dust_doubleq
915   corresponding files:
916 VIS: optprop_dustvis_TM_n50.dat
917 IR : optprop_dustir_n50.dat
918 Scatterer: h2o_ice
919   corresponding files:
920 VIS: optprop_icevis_n30.dat
921 IR : optprop_iceir_n30.dat
922 Splitting radiative calculations:  ngridmx,ngrid,ndomainsz,ndomain        3010
923        3010         151          20
924 aeroptproperties: interpolation grid
925 Scatterer:            1
926 Domain:            1
927 Min radius (m):   3.5353878E-08
928 Max radius (m):   9.9999997E-06
929 Scatterer:            1
930 Domain:            2
931 Min radius (m):   3.5353878E-08
932 Max radius (m):   9.9999997E-06
933 Scatterer:            2
934 Domain:            1
935 Min radius (m):   1.1303746E-07
936 Max radius (m):   2.9999999E-05
937 Scatterer:            2
938 Domain:            2
939 Min radius (m):   1.1303746E-07
940 Max radius (m):   2.9999999E-05
941 aeropacity: i_ice=           6
942 Typical profiles of Qext(vis)/Qext(IR)
943   and Qext(IR)/Qabs(IR):
944 Aerosol #            1
945 solsir:    1.600077       1.615189       1.632739       1.656156   
946   1.673179       1.688357       1.700979       1.712606       1.732043   
947   1.767485       1.846160       1.956535       2.126256       2.416159   
948   2.924124       3.627217       3.967699       4.619432       4.999025   
949   5.531231       5.867186       6.121830       6.161941       5.507154   
950   3.473106       1.283761      0.4759978      0.3981150      0.3957708   
951 Qext/Qabs(IR):    1.397849       1.393870       1.389535       1.383844   
952   1.379724       1.376071       1.373089       1.370418       1.366129   
953   1.359017       1.344959       1.327367       1.303859       1.271339   
954   1.229796       1.184386       1.146510       1.112741       1.088421   
955   1.071170       1.051863       1.034979       1.018531       1.007443   
956   1.002272       1.000459       1.000086       1.000062       1.000061   
957 Aerosol #            2
958 solsir:    1.576552       1.570121       1.565980       1.568253   
959   1.568665       1.579205       1.586318       1.594997       1.588596   
960   1.593094       1.581340       1.569868       1.526461       1.506043   
961   1.604984       1.702579       1.836655       1.946286       1.991414   
962   2.106776       2.331118       2.288174       2.249288       2.215015   
963   2.107166       2.065983       2.065983       2.065983       2.065983   
964 Qext/Qabs(IR):    1.302334       1.306394       1.309110       1.314693   
965   1.311184       1.305743       1.300053       1.297721       1.303020   
966   1.308674       1.323397       1.357420       1.410284       1.449922   
967   1.444335       1.432657       1.395102       1.358613       1.356317   
968   1.338657       1.233698       1.142038       1.084853       1.030944   
969   1.005017       1.000078       1.000078       1.000078       1.000078   
970Ls =   0.004706 check lat = -0.000007 lon =   0.000000
971 tauref(610. Pa) = 0.593544 tau(610. Pa) = 0.593544
972 In sugwd NKTOPG=           7
973 In vdifc: Tcond(P=1mb)=   136.2700      Lcond=   590000.0   
974           acond,bcond,ccond  3.2403751E-04  7.3383721E-03  3.3921041E-04
975 USING DEFAULTS : dtke_thermals = F
976 USING DEFAULTS : nsplit_thermals =          10
977 USING DEFAULTS : r_aspect_thermals =   1.000000   
978 Newcondens: improved_ztcond= T
979 Newcondens: bound_qco2= F
980 In newcondens:Tcond(P=1mb)=   136.2700      Lcond=   590000.0   
981 acond,bcond,ccond  3.2403751E-04  7.3383721E-03  3.3921041E-04
982  albedocaps: Use TES Cap albedoes ?
983  albedocaps: TESicealbedo =  T
984  albedocaps: Coefficient for Northern Cap ?
985  albedocaps: TESice_Ncoef =    1.600000   
986  albedocaps: Coefficient for Southern Cap ?
987  albedocaps: TESice_Scoef =    1.600000   
988 watercloud: igcm_h2o_vap=           7
989             igcm_h2o_ice=           6
990 time subsampling for microphysic ?
991 USING DEFAULTS : imicro =          15
992 imicro =           15
993 Physical timestep is   924.7396   
994 Microphysics timestep is   61.64931   
995 vrat_cld   31.6227757253798     
996 
997 Microphysics: size bin information:
998 i,rb_cld(i), rad_cld(i),dr_cld(i)
999 -----------------------------------
1000 1   0.100000E-09    0.100000E-06    0.124591E-06
1001 2   0.124691E-06    0.316228E-06    0.269617E-06
1002 3   0.394309E-06    0.100000E-05    0.852605E-06
1003 4   0.124691E-05    0.316228E-05    0.269617E-05
1004 5   0.394309E-05    0.100000E-04    0.999606E-02
1005 6   0.100000E-01
1006 -----------------------------------
1007 water_param contact parameter:  0.9500000   
1008 Variance of ice & CCN distribs :  0.3087235   
1009 nuice for sedimentation:  0.1000000   
1010 Volume of a water molecule:  1.956521739130435E-005
1011 
1012 dear user, please keep in mind that
1013 contact parameter IS constant
1014 
1015 correction for the shape of the ice particles ?
1016 USING DEFAULTS : ice_shape =  0.7500000   
1017  ice_shape =   0.7500000   
1018 water_param nueff Sedimentation:  0.1000000   
1019 water_param nueff Radiative:  0.1000000   
1020 
1021                         || STATS ||
1022 
1023 daysec   88775.00   
1024 dtphys   924.7396   
1025 nsteppd=          96
1026 istats=           8
1027 Storing           12 times per day
1028 thus every            8 physical timestep
1029 
1030 iniwrite: iim,jjm,llm,idayref          64          48          29       13380
1031 iniwrite: rad,omeg,g,mugaz,rcp   3397200.      7.0776514E-05   3.720000   
1032   43.49000      0.2567930   
1033 iniwrite: daysec,dtphys   88775.00       924.7396   
1034Ls =   0.009924 check lat = -0.000007 lon =   0.000000
1035 tauref(610. Pa) = 0.592217 tau(610. Pa) = 0.592217
1036Ls =   0.015805 check lat = -0.000007 lon =   0.000000
1037 tauref(610. Pa) = 0.590757 tau(610. Pa) = 0.590757
1038Ls =   0.021037 check lat = -0.000007 lon =   0.000000
1039 tauref(610. Pa) = 0.589297 tau(610. Pa) = 0.589297
1040 WRITEDIAGFI: date=   4.1666668E-02
1041 ==========================
1042 DIAGFI: creating variable tsurf
1043 ==========================
1044 DIAGFI: creating variable ps
1045 ==========================
1046 DIAGFI: creating variable co2ice
1047 ==========================
1048 DIAGFI: creating variable temp
1049 ==========================
1050 DIAGFI: creating variable u
1051 ==========================
1052 DIAGFI: creating variable v
1053 ==========================
1054 DIAGFI: creating variable emis
1055 ==========================
1056 DIAGFI: creating variable tsoil
1057 ==========================
1058 DIAGFI: creating variable inertiedat
1059 ==========================
1060 DIAGFI: creating variable co2
1061 ==========================
1062 DIAGFI: creating variable co2col
1063 ==========================
1064 DIAGFI: creating variable q01
1065 ==========================
1066 DIAGFI: creating variable q02
1067 ==========================
1068 DIAGFI: creating variable qsurf02
1069 ==========================
1070 DIAGFI: creating variable mtot
1071 ==========================
1072 DIAGFI: creating variable icetot
1073 ==========================
1074 DIAGFI: creating variable vmr_h2oice
1075 ==========================
1076 DIAGFI: creating variable vmr_h2ovap
1077 ==========================
1078 DIAGFI: creating variable reffice
1079 ==========================
1080 DIAGFI: creating variable Nccntot
1081 ==========================
1082 DIAGFI: creating variable Mccntot
1083 ==========================
1084 DIAGFI: creating variable rice
1085 ==========================
1086 DIAGFI: creating variable h2o_ice_s
1087 ==========================
1088 DIAGFI: creating variable albedo
1089 ==========================
1090 DIAGFI: creating variable tauref
1091 ==========================
1092 DIAGFI: creating variable tau
1093 ==========================
1094 DIAGFI: creating variable dqsdust
1095 ==========================
1096 DIAGFI: creating variable dqndust
1097 ==========================
1098 DIAGFI: creating variable reffdust
1099 ==========================
1100 DIAGFI: creating variable dustq
1101 ==========================
1102 DIAGFI: creating variable dustN
1103 ==========================
1104 DIAGFI: creating variable ccn
1105 ==========================
1106 DIAGFI: creating variable ccnN
1107 ==========================
1108 DIAGFI: creating variable ccnq
1109 ==========================
1110 DIAGFI: creating variable zmax_th
1111 ==========================
1112 DIAGFI: creating variable hfmax_th
1113 ==========================
1114 DIAGFI: creating variable wstar
1115Ls =   0.026262 check lat = -0.000007 lon =   0.000000
1116 tauref(610. Pa) = 0.587969 tau(610. Pa) = 0.587969
1117Ls =   0.031480 check lat = -0.000007 lon =   0.000000
1118 tauref(610. Pa) = 0.586509 tau(610. Pa) = 0.586509
1119Ls =   0.036719 check lat = -0.000007 lon =   0.000000
1120 tauref(610. Pa) = 0.585049 tau(610. Pa) = 0.585049
1121Ls =   0.041944 check lat = -0.000007 lon =   0.000000
1122 tauref(610. Pa) = 0.583722 tau(610. Pa) = 0.583722
1123 =====================
1124 STATS: creation de ps
1125 =====================
1126 STATS: creation de tsurf
1127 =====================
1128 STATS: creation de co2ice
1129 =====================
1130 STATS: creation de tauref
1131 =====================
1132 STATS: creation de fluxsurf_lw
1133 =====================
1134 STATS: creation de fluxsurf_sw
1135 =====================
1136 STATS: creation de fluxtop_lw
1137 =====================
1138 STATS: creation de fluxtop_sw
1139 =====================
1140 STATS: creation de temp
1141 =====================
1142 STATS: creation de u
1143 =====================
1144 STATS: creation de v
1145 =====================
1146 STATS: creation de w
1147 =====================
1148 STATS: creation de rho
1149 =====================
1150 STATS: creation de pressure
1151 =====================
1152 STATS: creation de q2
1153 =====================
1154 STATS: creation de emis
1155 =====================
1156 STATS: creation de zmax_th
1157 =====================
1158 STATS: creation de hfmax_th
1159 =====================
1160 STATS: creation de wstar
1161 =====================
1162 STATS: creation de vmr_h2ovap
1163 =====================
1164 STATS: creation de vmr_h2oice
1165 =====================
1166 STATS: creation de vmr_h2osat
1167 =====================
1168 STATS: creation de h2o_ice_s
1169 =====================
1170 STATS: creation de albedo
1171 =====================
1172 STATS: creation de mtot
1173 =====================
1174 STATS: creation de icetot
1175 =====================
1176 STATS: creation de reffice
1177 =====================
1178 STATS: creation de Nccntot
1179 =====================
1180 STATS: creation de Mccntot
1181 =====================
1182 STATS: creation de rice
1183 =====================
1184 STATS: creation de tauTES
1185 =====================
1186 STATS: creation de tau
1187 =====================
1188 STATS: creation de dqsdust
1189 =====================
1190 STATS: creation de dqndust
1191 =====================
1192 STATS: creation de reffdust
1193 =====================
1194 STATS: creation de dustq
1195 =====================
1196 STATS: creation de dustN
1197 =====================
1198 STATS: creation de ccnq
1199 =====================
1200 STATS: creation de ccnN
1201 WRITEDIAGFI: date=   8.3333336E-02
1202Ls =   0.047169 check lat = -0.000007 lon =   0.000000
1203 tauref(610. Pa) = 0.582262 tau(610. Pa) = 0.582262
1204Ls =   0.053043 check lat = -0.000007 lon =   0.000000
1205 tauref(610. Pa) = 0.580801 tau(610. Pa) = 0.580801
1206Ls =   0.058268 check lat = -0.000007 lon =   0.000000
1207 tauref(610. Pa) = 0.579474 tau(610. Pa) = 0.579474
1208Ls =   0.063500 check lat = -0.000007 lon =   0.000000
1209 tauref(610. Pa) = 0.578014 tau(610. Pa) = 0.578014
1210 WRITEDIAGFI: date=   0.1250000   
12110**********    pas    120     jour*****heure  3.0    date   0.12500    **********
1212          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1213        0.947168 0.845509E+00     0.951     1.076     0.966    18.279     0.828
1214
1215Ls =   0.068732 check lat = -0.000007 lon =   0.000000
1216 tauref(610. Pa) = 0.576554 tau(610. Pa) = 0.576554
1217Ls =   0.073950 check lat = -0.000007 lon =   0.000000
1218 tauref(610. Pa) = 0.575226 tau(610. Pa) = 0.575226
1219Ls =   0.079176 check lat = -0.000007 lon =   0.000000
1220 tauref(610. Pa) = 0.573766 tau(610. Pa) = 0.573766
1221Ls =   0.085056 check lat = -0.000007 lon =   0.000000
1222 tauref(610. Pa) = 0.572306 tau(610. Pa) = 0.572306
1223 WRITEDIAGFI: date=   0.1666667   
1224Ls =   0.090275 check lat = -0.000007 lon =   0.000000
1225 tauref(610. Pa) = 0.570979 tau(610. Pa) = 0.570979
1226Ls =   0.095507 check lat = -0.000007 lon =   0.000000
1227 tauref(610. Pa) = 0.569519 tau(610. Pa) = 0.569519
1228Ls =   0.100732 check lat = -0.000007 lon =   0.000000
1229 tauref(610. Pa) = 0.568058 tau(610. Pa) = 0.568058
1230Ls =   0.105957 check lat = -0.000007 lon =   0.000000
1231 tauref(610. Pa) = 0.566731 tau(610. Pa) = 0.566731
1232 WRITEDIAGFI: date=   0.2083333   
1233Ls =   0.111182 check lat = -0.000007 lon =   0.000000
1234 tauref(610. Pa) = 0.565271 tau(610. Pa) = 0.565271
1235Ls =   0.117056 check lat = -0.000007 lon =   0.000000
1236 tauref(610. Pa) = 0.563811 tau(610. Pa) = 0.563811
1237Ls =   0.122281 check lat = -0.000007 lon =   0.000000
1238 tauref(610. Pa) = 0.562483 tau(610. Pa) = 0.562483
1239Ls =   0.127513 check lat = -0.000007 lon =   0.000000
1240 tauref(610. Pa) = 0.561023 tau(610. Pa) = 0.561023
1241 WRITEDIAGFI: date=   0.2500000   
12420**********    pas    240     jour*****heure  6.0    date   0.25000    **********
1243          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1244        0.947275 0.894936E+00     0.951     1.116     0.966    18.153     0.829
1245
1246Ls =   0.132738 check lat = -0.000007 lon =   0.000000
1247 tauref(610. Pa) = 0.559563 tau(610. Pa) = 0.559563
1248Ls =   0.137956 check lat = -0.000007 lon =   0.000000
1249 tauref(610. Pa) = 0.558236 tau(610. Pa) = 0.558236
1250Ls =   0.143188 check lat = -0.000007 lon =   0.000000
1251 tauref(610. Pa) = 0.556776 tau(610. Pa) = 0.556776
1252Ls =   0.149062 check lat = -0.000007 lon =   0.000000
1253 tauref(610. Pa) = 0.555315 tau(610. Pa) = 0.555315
1254 WRITEDIAGFI: date=   0.2916667   
1255Ls =   0.154287 check lat = -0.000007 lon =   0.000000
1256 tauref(610. Pa) = 0.553988 tau(610. Pa) = 0.553988
1257Ls =   0.159519 check lat = -0.000007 lon =   0.000000
1258 tauref(610. Pa) = 0.552528 tau(610. Pa) = 0.552528
1259Ls =   0.164731 check lat = -0.000007 lon =   0.000000
1260 tauref(610. Pa) = 0.551068 tau(610. Pa) = 0.551068
1261Ls =   0.169956 check lat = -0.000007 lon =   0.000000
1262 tauref(610. Pa) = 0.549740 tau(610. Pa) = 0.549740
1263 WRITEDIAGFI: date=   0.3333333   
1264Ls =   0.175181 check lat = -0.000007 lon =   0.000000
1265 tauref(610. Pa) = 0.548280 tau(610. Pa) = 0.548280
1266Ls =   0.181061 check lat = -0.000007 lon =   0.000000
1267 tauref(610. Pa) = 0.546820 tau(610. Pa) = 0.546820
1268Ls =   0.185638 check lat = -0.000007 lon =   0.000000
1269 tauref(610. Pa) = 0.545493 tau(610. Pa) = 0.545493
1270Ls =   0.191512 check lat = -0.000007 lon =   0.000000
1271 tauref(610. Pa) = 0.544033 tau(610. Pa) = 0.544032
1272 WRITEDIAGFI: date=   0.3750000   
12730**********    pas    360     jour*****heure  9.0    date   0.37500    **********
1274          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1275        0.947374 0.906323E+00     0.952     1.101     0.966    18.186     0.827
1276
1277Ls =   0.196730 check lat = -0.000007 lon =   0.000000
1278 tauref(610. Pa) = 0.542572 tau(610. Pa) = 0.542572
1279Ls =   0.201962 check lat = -0.000007 lon =   0.000000
1280 tauref(610. Pa) = 0.541245 tau(610. Pa) = 0.541245
1281Ls =   0.207173 check lat = -0.000007 lon =   0.000000
1282 tauref(610. Pa) = 0.539785 tau(610. Pa) = 0.539785
1283Ls =   0.213054 check lat = -0.000007 lon =   0.000000
1284 tauref(610. Pa) = 0.538325 tau(610. Pa) = 0.538325
1285 WRITEDIAGFI: date=   0.4166667   
1286Ls =   0.217630 check lat = -0.000007 lon =   0.000000
1287 tauref(610. Pa) = 0.536997 tau(610. Pa) = 0.536997
1288Ls =   0.223504 check lat = -0.000007 lon =   0.000000
1289 tauref(610. Pa) = 0.535537 tau(610. Pa) = 0.535537
1290Ls =   0.228729 check lat = -0.000007 lon =   0.000000
1291 tauref(610. Pa) = 0.534077 tau(610. Pa) = 0.534077
1292Ls =   0.233954 check lat = -0.000007 lon =   0.000000
1293 tauref(610. Pa) = 0.532750 tau(610. Pa) = 0.532750
1294 WRITEDIAGFI: date=   0.4583333   
1295Ls =   0.239166 check lat = -0.000007 lon =   0.000000
1296 tauref(610. Pa) = 0.531290 tau(610. Pa) = 0.531289
1297Ls =   0.245054 check lat = -0.000007 lon =   0.000000
1298 tauref(610. Pa) = 0.529829 tau(610. Pa) = 0.529829
1299Ls =   0.249616 check lat = -0.000007 lon =   0.000000
1300 tauref(610. Pa) = 0.528502 tau(610. Pa) = 0.528502
1301Ls =   0.255497 check lat = -0.000007 lon =   0.000000
1302 tauref(610. Pa) = 0.527042 tau(610. Pa) = 0.527042
1303 WRITEDIAGFI: date=   0.5000000   
13040**********    pas    480     jour*****heure 12.0    date   0.50000    **********
1305          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1306        0.947464 0.101995E+01     0.952     1.069     0.966    18.359     0.826
1307
1308Ls =   0.260722 check lat = -0.000007 lon =   0.000000
1309 tauref(610. Pa) = 0.525582 tau(610. Pa) = 0.525582
1310Ls =   0.265933 check lat = -0.000007 lon =   0.000000
1311 tauref(610. Pa) = 0.524254 tau(610. Pa) = 0.524254
1312Ls =   0.271159 check lat = -0.000007 lon =   0.000000
1313 tauref(610. Pa) = 0.522794 tau(610. Pa) = 0.522794
1314Ls =   0.277032 check lat = -0.000007 lon =   0.000000
1315 tauref(610. Pa) = 0.521334 tau(610. Pa) = 0.521334
1316 WRITEDIAGFI: date=   0.5416667   
1317Ls =   0.281609 check lat = -0.000007 lon =   0.000000
1318 tauref(610. Pa) = 0.520007 tau(610. Pa) = 0.520007
1319Ls =   0.287483 check lat = -0.000007 lon =   0.000000
1320 tauref(610. Pa) = 0.518546 tau(610. Pa) = 0.518546
1321Ls =   0.292701 check lat = -0.000007 lon =   0.000000
1322 tauref(610. Pa) = 0.517086 tau(610. Pa) = 0.517086
1323Ls =   0.297926 check lat = -0.000007 lon =   0.000000
1324 tauref(610. Pa) = 0.515759 tau(610. Pa) = 0.515759
1325 WRITEDIAGFI: date=   0.5833333   
1326Ls =   0.303151 check lat = -0.000007 lon =   0.000000
1327 tauref(610. Pa) = 0.514299 tau(610. Pa) = 0.514299
1328Ls =   0.309018 check lat = -0.000007 lon =   0.000000
1329 tauref(610. Pa) = 0.512839 tau(610. Pa) = 0.512839
1330Ls =   0.313601 check lat = -0.000007 lon =   0.000000
1331 tauref(610. Pa) = 0.511511 tau(610. Pa) = 0.511511
1332Ls =   0.319468 check lat = -0.000007 lon =   0.000000
1333 tauref(610. Pa) = 0.510051 tau(610. Pa) = 0.510051
1334 WRITEDIAGFI: date=   0.6250000   
13350**********    pas    600     jour*****heure 15.0    date   0.62500    **********
1336          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1337        0.947557 0.882835E+00     0.952     1.055     0.967    18.467     0.826
1338
1339Ls =   0.324694 check lat = -0.000007 lon =   0.000000
1340 tauref(610. Pa) = 0.508591 tau(610. Pa) = 0.508591
1341Ls =   0.329905 check lat = -0.000007 lon =   0.000000
1342 tauref(610. Pa) = 0.507264 tau(610. Pa) = 0.507263
1343Ls =   0.335130 check lat = -0.000007 lon =   0.000000
1344 tauref(610. Pa) = 0.505803 tau(610. Pa) = 0.505803
1345Ls =   0.341004 check lat = -0.000007 lon =   0.000000
1346 tauref(610. Pa) = 0.504343 tau(610. Pa) = 0.504343
1347 WRITEDIAGFI: date=   0.6666667   
1348Ls =   0.345573 check lat = -0.000007 lon =   0.000000
1349 tauref(610. Pa) = 0.503016 tau(610. Pa) = 0.503016
1350Ls =   0.351447 check lat = -0.000007 lon =   0.000000
1351 tauref(610. Pa) = 0.501556 tau(610. Pa) = 0.501556
1352Ls =   0.356666 check lat = -0.000007 lon =   0.000000
1353 tauref(610. Pa) = 0.500096 tau(610. Pa) = 0.500095
1354Ls =   0.361884 check lat = -0.000007 lon =   0.000000
1355 tauref(610. Pa) = 0.498768 tau(610. Pa) = 0.498768
1356 WRITEDIAGFI: date=   0.7083333   
1357Ls =   0.367109 check lat = -0.000007 lon =   0.000000
1358 tauref(610. Pa) = 0.497308 tau(610. Pa) = 0.497308
1359Ls =   0.372983 check lat = -0.000007 lon =   0.000000
1360 tauref(610. Pa) = 0.495848 tau(610. Pa) = 0.495848
1361Ls =   0.377552 check lat = -0.000007 lon =   0.000000
1362 tauref(610. Pa) = 0.494521 tau(610. Pa) = 0.494521
1363Ls =   0.383426 check lat = -0.000007 lon =   0.000000
1364 tauref(610. Pa) = 0.493060 tau(610. Pa) = 0.493060
1365 WRITEDIAGFI: date=   0.7500000   
13660**********    pas    720     jour*****heure 18.0    date   0.75000    **********
1367          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1368        0.947652 0.947625E+00     0.952     1.055     0.966    18.448     0.826
1369
1370Ls =   0.388651 check lat = -0.000007 lon =   0.000000
1371 tauref(610. Pa) = 0.491600 tau(610. Pa) = 0.491600
1372Ls =   0.393870 check lat = -0.000007 lon =   0.000000
1373 tauref(610. Pa) = 0.490273 tau(610. Pa) = 0.490273
1374Ls =   0.399088 check lat = -0.000007 lon =   0.000000
1375 tauref(610. Pa) = 0.488813 tau(610. Pa) = 0.488813
1376Ls =   0.404962 check lat = -0.000007 lon =   0.000000
1377 tauref(610. Pa) = 0.487353 tau(610. Pa) = 0.487352
1378 WRITEDIAGFI: date=   0.7916667   
1379Ls =   0.409531 check lat = -0.000007 lon =   0.000000
1380 tauref(610. Pa) = 0.486025 tau(610. Pa) = 0.486025
1381Ls =   0.415405 check lat = -0.000007 lon =   0.000000
1382 tauref(610. Pa) = 0.484565 tau(610. Pa) = 0.484565
1383Ls =   0.420617 check lat = -0.000007 lon =   0.000000
1384 tauref(610. Pa) = 0.483105 tau(610. Pa) = 0.483105
1385Ls =   0.425849 check lat = -0.000007 lon =   0.000000
1386 tauref(610. Pa) = 0.481778 tau(610. Pa) = 0.481777
1387 WRITEDIAGFI: date=   0.8333333   
1388Ls =   0.431060 check lat = -0.000007 lon =   0.000000
1389 tauref(610. Pa) = 0.480317 tau(610. Pa) = 0.480317
1390Ls =   0.436934 check lat = -0.000007 lon =   0.000000
1391 tauref(610. Pa) = 0.478857 tau(610. Pa) = 0.478857
1392Ls =   0.441503 check lat = -0.000007 lon =   0.000000
1393 tauref(610. Pa) = 0.477530 tau(610. Pa) = 0.477530
1394Ls =   0.447377 check lat = -0.000007 lon =   0.000000
1395 tauref(610. Pa) = 0.476070 tau(610. Pa) = 0.476070
1396 WRITEDIAGFI: date=   0.8750000   
13970**********    pas    840     jour*****heure 21.0    date   0.87500    **********
1398          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1399        0.947726 0.970058E+00     0.951     1.098     0.965    18.392     0.825
1400
1401Ls =   0.452596 check lat = -0.000007 lon =   0.000000
1402 tauref(610. Pa) = 0.474610 tau(610. Pa) = 0.474610
1403Ls =   0.457814 check lat = -0.000007 lon =   0.000000
1404 tauref(610. Pa) = 0.473282 tau(610. Pa) = 0.473282
1405Ls =   0.463032 check lat = -0.000007 lon =   0.000000
1406 tauref(610. Pa) = 0.471822 tau(610. Pa) = 0.471822
1407Ls =   0.468906 check lat = -0.000007 lon =   0.000000
1408 tauref(610. Pa) = 0.470362 tau(610. Pa) = 0.470362
1409 WRITEDIAGFI: date=   0.9166667   
1410Ls =   0.473469 check lat = -0.000007 lon =   0.000000
1411 tauref(610. Pa) = 0.469034 tau(610. Pa) = 0.469034
1412Ls =   0.479343 check lat = -0.000007 lon =   0.000000
1413 tauref(610. Pa) = 0.467574 tau(610. Pa) = 0.467574
1414Ls =   0.484561 check lat = -0.000007 lon =   0.000000
1415 tauref(610. Pa) = 0.466114 tau(610. Pa) = 0.466114
1416Ls =   0.489772 check lat = -0.000007 lon =   0.000000
1417 tauref(610. Pa) = 0.464787 tau(610. Pa) = 0.464787
1418 WRITEDIAGFI: date=   0.9583333   
1419Ls =   0.494997 check lat = -0.000007 lon =   0.000000
1420 tauref(610. Pa) = 0.463327 tau(610. Pa) = 0.463327
1421Ls =   0.500871 check lat = -0.000007 lon =   0.000000
1422 tauref(610. Pa) = 0.461867 tau(610. Pa) = 0.461867
1423Ls =   0.505441 check lat = -0.000007 lon =   0.000000
1424 tauref(610. Pa) = 0.460539 tau(610. Pa) = 0.460539
1425Ls =   0.511308 check lat = -0.000007 lon =   0.000000
1426 tauref(610. Pa) = 0.459079 tau(610. Pa) = 0.459079
1427 WRITEDIAGFI: date=    1.000000   
14280**********    pas    960     jour*****heure  0.0    date   1.00000    **********
1429          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1430        0.947798 0.103170E+01     0.951     1.081     0.965    18.338     0.825
1431
1432 LTE rad. calculations up to layer           29
1433Ls =   0.516526 check lat = -0.000007 lon =   0.000000
1434 tauref(610. Pa) = 0.458992 tau(610. Pa) = 0.458992
1435Ls =   0.521744 check lat = -0.000007 lon =   0.000000
1436 tauref(610. Pa) = 0.459039 tau(610. Pa) = 0.459039
1437Ls =   0.526963 check lat = -0.000007 lon =   0.000000
1438 tauref(610. Pa) = 0.459089 tau(610. Pa) = 0.459089
1439Ls =   0.532830 check lat = -0.000007 lon =   0.000000
1440 tauref(610. Pa) = 0.459140 tau(610. Pa) = 0.459140
1441 WRITEDIAGFI: date=    1.041667   
1442Ls =   0.537399 check lat = -0.000007 lon =   0.000000
1443 tauref(610. Pa) = 0.459186 tau(610. Pa) = 0.459186
1444Ls =   0.543266 check lat = -0.000007 lon =   0.000000
1445 tauref(610. Pa) = 0.459237 tau(610. Pa) = 0.459237
1446Ls =   0.548485 check lat = -0.000007 lon =   0.000000
1447 tauref(610. Pa) = 0.459288 tau(610. Pa) = 0.459288
1448Ls =   0.553703 check lat = -0.000007 lon =   0.000000
1449 tauref(610. Pa) = 0.459334 tau(610. Pa) = 0.459334
1450 WRITEDIAGFI: date=    1.083333   
1451Ls =   0.558928 check lat = -0.000007 lon =   0.000000
1452 tauref(610. Pa) = 0.459385 tau(610. Pa) = 0.459385
1453Ls =   0.564795 check lat = -0.000007 lon =   0.000000
1454 tauref(610. Pa) = 0.459436 tau(610. Pa) = 0.459436
1455Ls =   0.569365 check lat = -0.000007 lon =   0.000000
1456 tauref(610. Pa) = 0.459482 tau(610. Pa) = 0.459482
1457Ls =   0.575232 check lat = -0.000007 lon =   0.000000
1458 tauref(610. Pa) = 0.459533 tau(610. Pa) = 0.459533
1459 WRITEDIAGFI: date=    1.125000   
14600**********    pas   1080     jour*****heure  3.0    date   1.12500    **********
1461          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1462        0.947890 0.843030E+00     0.951     1.039     0.965    18.230     0.827
1463
1464Ls =   0.580450 check lat = -0.000007 lon =   0.000000
1465 tauref(610. Pa) = 0.459583 tau(610. Pa) = 0.459583
1466Ls =   0.585661 check lat = -0.000007 lon =   0.000000
1467 tauref(610. Pa) = 0.459629 tau(610. Pa) = 0.459629
1468Ls =   0.590880 check lat = -0.000007 lon =   0.000000
1469 tauref(610. Pa) = 0.459680 tau(610. Pa) = 0.459680
1470Ls =   0.596747 check lat = -0.000007 lon =   0.000000
1471 tauref(610. Pa) = 0.459731 tau(610. Pa) = 0.459731
1472 WRITEDIAGFI: date=    1.166667   
1473Ls =   0.601316 check lat = -0.000007 lon =   0.000000
1474 tauref(610. Pa) = 0.459777 tau(610. Pa) = 0.459777
1475Ls =   0.607190 check lat = -0.000007 lon =   0.000000
1476 tauref(610. Pa) = 0.459828 tau(610. Pa) = 0.459828
1477Ls =   0.612402 check lat = -0.000007 lon =   0.000000
1478 tauref(610. Pa) = 0.459879 tau(610. Pa) = 0.459879
1479Ls =   0.617620 check lat = -0.000007 lon =   0.000000
1480 tauref(610. Pa) = 0.459925 tau(610. Pa) = 0.459925
1481 WRITEDIAGFI: date=    1.208333   
1482Ls =   0.622838 check lat = -0.000007 lon =   0.000000
1483 tauref(610. Pa) = 0.459976 tau(610. Pa) = 0.459976
1484Ls =   0.628698 check lat = -0.000007 lon =   0.000000
1485 tauref(610. Pa) = 0.460026 tau(610. Pa) = 0.460026
1486Ls =   0.633268 check lat = -0.000007 lon =   0.000000
1487 tauref(610. Pa) = 0.460073 tau(610. Pa) = 0.460073
1488Ls =   0.639142 check lat = -0.000007 lon =   0.000000
1489 tauref(610. Pa) = 0.460123 tau(610. Pa) = 0.460123
1490 WRITEDIAGFI: date=    1.250000   
14910**********    pas   1200     jour*****heure  6.0    date   1.25000    **********
1492          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1493        0.947987 0.935672E+00     0.951     1.057     0.965    18.134     0.828
1494
1495Ls =   0.644353 check lat = -0.000007 lon =   0.000000
1496 tauref(610. Pa) = 0.460174 tau(610. Pa) = 0.460174
1497Ls =   0.649571 check lat = -0.000007 lon =   0.000000
1498 tauref(610. Pa) = 0.460220 tau(610. Pa) = 0.460220
1499Ls =   0.654790 check lat = -0.000007 lon =   0.000000
1500 tauref(610. Pa) = 0.460271 tau(610. Pa) = 0.460271
1501Ls =   0.660650 check lat = -0.000007 lon =   0.000000
1502 tauref(610. Pa) = 0.460322 tau(610. Pa) = 0.460322
1503 WRITEDIAGFI: date=    1.291667   
1504Ls =   0.665219 check lat = -0.000007 lon =   0.000000
1505 tauref(610. Pa) = 0.460368 tau(610. Pa) = 0.460368
1506Ls =   0.671087 check lat = -0.000007 lon =   0.000000
1507 tauref(610. Pa) = 0.460419 tau(610. Pa) = 0.460419
1508Ls =   0.676305 check lat = -0.000007 lon =   0.000000
1509 tauref(610. Pa) = 0.460470 tau(610. Pa) = 0.460470
1510Ls =   0.681516 check lat = -0.000007 lon =   0.000000
1511 tauref(610. Pa) = 0.460516 tau(610. Pa) = 0.460516
1512 WRITEDIAGFI: date=    1.333333   
1513Ls =   0.686741 check lat = -0.000007 lon =   0.000000
1514 tauref(610. Pa) = 0.460567 tau(610. Pa) = 0.460567
1515Ls =   0.692602 check lat = -0.000007 lon =   0.000000
1516 tauref(610. Pa) = 0.460617 tau(610. Pa) = 0.460617
1517Ls =   0.697164 check lat = -0.000007 lon =   0.000000
1518 tauref(610. Pa) = 0.460663 tau(610. Pa) = 0.460664
1519Ls =   0.703031 check lat = -0.000007 lon =   0.000000
1520 tauref(610. Pa) = 0.460714 tau(610. Pa) = 0.460714
1521 WRITEDIAGFI: date=    1.375000   
15220**********    pas   1320     jour*****heure  9.0    date   1.37500    **********
1523          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1524        0.948083 0.947274E+00     0.952     1.060     0.966    18.040     0.827
1525
1526Ls =   0.708243 check lat = -0.000007 lon =   0.000000
1527 tauref(610. Pa) = 0.460765 tau(610. Pa) = 0.460765
1528Ls =   0.713461 check lat = -0.000007 lon =   0.000000
1529 tauref(610. Pa) = 0.460811 tau(610. Pa) = 0.460811
1530Ls =   0.718679 check lat = -0.000007 lon =   0.000000
1531 tauref(610. Pa) = 0.460862 tau(610. Pa) = 0.460862
1532Ls =   0.724540 check lat = -0.000007 lon =   0.000000
1533 tauref(610. Pa) = 0.460913 tau(610. Pa) = 0.460913
1534 WRITEDIAGFI: date=    1.416667   
1535Ls =   0.729116 check lat = -0.000007 lon =   0.000000
1536 tauref(610. Pa) = 0.460959 tau(610. Pa) = 0.460959
1537Ls =   0.734976 check lat = -0.000007 lon =   0.000000
1538 tauref(610. Pa) = 0.461010 tau(610. Pa) = 0.461010
1539Ls =   0.740188 check lat = -0.000007 lon =   0.000000
1540 tauref(610. Pa) = 0.461060 tau(610. Pa) = 0.461060
1541Ls =   0.745399 check lat = -0.000007 lon =   0.000000
1542 tauref(610. Pa) = 0.461107 tau(610. Pa) = 0.461107
1543 WRITEDIAGFI: date=    1.458333   
1544Ls =   0.750617 check lat = -0.000007 lon =   0.000000
1545 tauref(610. Pa) = 0.461157 tau(610. Pa) = 0.461157
1546Ls =   0.755829 check lat = -0.000007 lon =   0.000000
1547 tauref(610. Pa) = 0.461208 tau(610. Pa) = 0.461208
1548Ls =   0.761047 check lat = -0.000007 lon =   0.000000
1549 tauref(610. Pa) = 0.461254 tau(610. Pa) = 0.461254
1550Ls =   0.766914 check lat = -0.000007 lon =   0.000000
1551 tauref(610. Pa) = 0.461305 tau(610. Pa) = 0.461305
1552 WRITEDIAGFI: date=    1.500000   
15530**********    pas   1440     jour*****heure 12.0    date   1.50000    **********
1554          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1555        0.948176 0.934652E+00     0.952     1.057     0.966    18.081     0.826
1556
1557Ls =   0.772126 check lat = -0.000007 lon =   0.000000
1558 tauref(610. Pa) = 0.461356 tau(610. Pa) = 0.461356
1559Ls =   0.777337 check lat = -0.000007 lon =   0.000000
1560 tauref(610. Pa) = 0.461402 tau(610. Pa) = 0.461402
1561Ls =   0.782555 check lat = -0.000007 lon =   0.000000
1562 tauref(610. Pa) = 0.461453 tau(610. Pa) = 0.461453
1563Ls =   0.787773 check lat = -0.000007 lon =   0.000000
1564 tauref(610. Pa) = 0.461504 tau(610. Pa) = 0.461504
1565 WRITEDIAGFI: date=    1.541667   
1566Ls =   0.792985 check lat = -0.000007 lon =   0.000000
1567 tauref(610. Pa) = 0.461550 tau(610. Pa) = 0.461550
1568Ls =   0.798852 check lat = -0.000007 lon =   0.000000
1569 tauref(610. Pa) = 0.461601 tau(610. Pa) = 0.461601
1570Ls =   0.804063 check lat = -0.000007 lon =   0.000000
1571 tauref(610. Pa) = 0.461651 tau(610. Pa) = 0.461651
1572Ls =   0.809275 check lat = -0.000007 lon =   0.000000
1573 tauref(610. Pa) = 0.461698 tau(610. Pa) = 0.461698
1574 WRITEDIAGFI: date=    1.583333   
1575Ls =   0.814493 check lat = -0.000007 lon =   0.000000
1576 tauref(610. Pa) = 0.461748 tau(610. Pa) = 0.461748
1577Ls =   0.819705 check lat = -0.000007 lon =   0.000000
1578 tauref(610. Pa) = 0.461799 tau(610. Pa) = 0.461799
1579Ls =   0.824916 check lat = -0.000007 lon =   0.000000
1580 tauref(610. Pa) = 0.461845 tau(610. Pa) = 0.461845
1581Ls =   0.830783 check lat = -0.000007 lon =   0.000000
1582 tauref(610. Pa) = 0.461896 tau(610. Pa) = 0.461896
1583 WRITEDIAGFI: date=    1.625000   
15840**********    pas   1560     jour*****heure 15.0    date   1.62500    **********
1585          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1586        0.948276 0.830403E+00     0.952     1.070     0.965    18.139     0.827
1587
1588Ls =   0.835995 check lat = -0.000007 lon =   0.000000
1589 tauref(610. Pa) = 0.461947 tau(610. Pa) = 0.461947
1590Ls =   0.841206 check lat = -0.000007 lon =   0.000000
1591 tauref(610. Pa) = 0.461993 tau(610. Pa) = 0.461993
1592Ls =   0.846424 check lat = -0.000007 lon =   0.000000
1593 tauref(610. Pa) = 0.462044 tau(610. Pa) = 0.462044
1594Ls =   0.851636 check lat = -0.000007 lon =   0.000000
1595 tauref(610. Pa) = 0.462095 tau(610. Pa) = 0.462095
1596 WRITEDIAGFI: date=    1.666667   
1597Ls =   0.856847 check lat = -0.000007 lon =   0.000000
1598 tauref(610. Pa) = 0.462141 tau(610. Pa) = 0.462141
1599Ls =   0.862721 check lat = -0.000007 lon =   0.000000
1600 tauref(610. Pa) = 0.462191 tau(610. Pa) = 0.462192
1601Ls =   0.867926 check lat = -0.000007 lon =   0.000000
1602 tauref(610. Pa) = 0.462242 tau(610. Pa) = 0.462242
1603Ls =   0.873137 check lat = -0.000007 lon =   0.000000
1604 tauref(610. Pa) = 0.462288 tau(610. Pa) = 0.462288
1605 WRITEDIAGFI: date=    1.708333   
1606Ls =   0.878349 check lat = -0.000007 lon =   0.000000
1607 tauref(610. Pa) = 0.462339 tau(610. Pa) = 0.462339
1608Ls =   0.883560 check lat = -0.000007 lon =   0.000000
1609 tauref(610. Pa) = 0.462390 tau(610. Pa) = 0.462390
1610Ls =   0.888772 check lat = -0.000007 lon =   0.000000
1611 tauref(610. Pa) = 0.462436 tau(610. Pa) = 0.462436
1612Ls =   0.893983 check lat = -0.000007 lon =   0.000000
1613 tauref(610. Pa) = 0.462487 tau(610. Pa) = 0.462487
1614 WRITEDIAGFI: date=    1.750000   
16150**********    pas   1680     jour*****heure 18.0    date   1.75000    **********
1616          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1617        0.948370 0.854859E+00     0.951     1.075     0.965    18.091     0.827
1618
1619Ls =   0.899850 check lat = -0.000007 lon =   0.000000
1620 tauref(610. Pa) = 0.462538 tau(610. Pa) = 0.462538
1621Ls =   0.905061 check lat = -0.000007 lon =   0.000000
1622 tauref(610. Pa) = 0.462584 tau(610. Pa) = 0.462584
1623Ls =   0.910280 check lat = -0.000007 lon =   0.000000
1624 tauref(610. Pa) = 0.462635 tau(610. Pa) = 0.462635
1625Ls =   0.915484 check lat = -0.000007 lon =   0.000000
1626 tauref(610. Pa) = 0.462685 tau(610. Pa) = 0.462685
1627 WRITEDIAGFI: date=    1.791667   
1628Ls =   0.920696 check lat = -0.000007 lon =   0.000000
1629 tauref(610. Pa) = 0.462732 tau(610. Pa) = 0.462731
1630Ls =   0.925914 check lat = -0.000007 lon =   0.000000
1631 tauref(610. Pa) = 0.462782 tau(610. Pa) = 0.462782
1632Ls =   0.931774 check lat = -0.000007 lon =   0.000000
1633 tauref(610. Pa) = 0.462833 tau(610. Pa) = 0.462833
1634Ls =   0.936979 check lat = -0.000007 lon =   0.000000
1635 tauref(610. Pa) = 0.462879 tau(610. Pa) = 0.462879
1636 WRITEDIAGFI: date=    1.833333   
1637Ls =   0.942197 check lat = -0.000007 lon =   0.000000
1638 tauref(610. Pa) = 0.462930 tau(610. Pa) = 0.462930
1639Ls =   0.947416 check lat = -0.000007 lon =   0.000000
1640 tauref(610. Pa) = 0.462981 tau(610. Pa) = 0.462981
1641Ls =   0.952620 check lat = -0.000007 lon =   0.000000
1642 tauref(610. Pa) = 0.463027 tau(610. Pa) = 0.463027
1643Ls =   0.957825 check lat = -0.000007 lon =   0.000000
1644 tauref(610. Pa) = 0.463078 tau(610. Pa) = 0.463078
1645 WRITEDIAGFI: date=    1.875000   
16460**********    pas   1800     jour*****heure 21.0    date   1.87500    **********
1647          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1648        0.948442 0.914820E+00     0.951     1.076     0.965    17.987     0.826
1649
1650Ls =   0.963692 check lat = -0.000007 lon =   0.000000
1651 tauref(610. Pa) = 0.463129 tau(610. Pa) = 0.463128
1652Ls =   0.968910 check lat = -0.000007 lon =   0.000000
1653 tauref(610. Pa) = 0.463175 tau(610. Pa) = 0.463175
1654Ls =   0.974115 check lat = -0.000007 lon =   0.000000
1655 tauref(610. Pa) = 0.463225 tau(610. Pa) = 0.463226
1656Ls =   0.979333 check lat = -0.000007 lon =   0.000000
1657 tauref(610. Pa) = 0.463276 tau(610. Pa) = 0.463276
1658 WRITEDIAGFI: date=    1.916667   
1659Ls =   0.984538 check lat = -0.000007 lon =   0.000000
1660 tauref(610. Pa) = 0.463322 tau(610. Pa) = 0.463322
1661Ls =   0.989742 check lat = -0.000007 lon =   0.000000
1662 tauref(610. Pa) = 0.463373 tau(610. Pa) = 0.463373
1663Ls =   0.995609 check lat = -0.000007 lon =   0.000000
1664 tauref(610. Pa) = 0.463424 tau(610. Pa) = 0.463424
1665Ls =   1.000165 check lat = -0.000007 lon =   0.000000
1666 tauref(610. Pa) = 0.463470 tau(610. Pa) = 0.463470
1667 WRITEDIAGFI: date=    1.958333   
1668Ls =   1.006039 check lat = -0.000007 lon =   0.000000
1669 tauref(610. Pa) = 0.463521 tau(610. Pa) = 0.463521
1670Ls =   1.011237 check lat = -0.000007 lon =   0.000000
1671 tauref(610. Pa) = 0.463572 tau(610. Pa) = 0.463572
1672Ls =   1.016455 check lat = -0.000007 lon =   0.000000
1673 tauref(610. Pa) = 0.463618 tau(610. Pa) = 0.463618
1674Ls =   1.021667 check lat = -0.000007 lon =   0.000000
1675 tauref(610. Pa) = 0.463669 tau(610. Pa) = 0.463668
1676 WRITEDIAGFI: date=    2.000000   
16770**********    pas   1920     jour*****heure  0.0    date   2.00000    **********
1678          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1679        0.948507 0.923695E+00     0.950     1.048     0.964    18.032     0.827
1680
1681 LTE rad. calculations up to layer           29
1682Ls =   1.027520 check lat = -0.000007 lon =   0.000000
1683 tauref(610. Pa) = 0.463678 tau(610. Pa) = 0.463678
1684Ls =   1.032083 check lat = -0.000007 lon =   0.000000
1685 tauref(610. Pa) = 0.463683 tau(610. Pa) = 0.463683
1686Ls =   1.037950 check lat = -0.000007 lon =   0.000000
1687 tauref(610. Pa) = 0.463688 tau(610. Pa) = 0.463688
1688Ls =   1.043154 check lat = -0.000007 lon =   0.000000
1689 tauref(610. Pa) = 0.463694 tau(610. Pa) = 0.463694
1690 WRITEDIAGFI: date=    2.041667   
1691Ls =   1.048359 check lat = -0.000007 lon =   0.000000
1692 tauref(610. Pa) = 0.463699 tau(610. Pa) = 0.463699
1693Ls =   1.053577 check lat = -0.000007 lon =   0.000000
1694 tauref(610. Pa) = 0.463704 tau(610. Pa) = 0.463704
1695Ls =   1.059431 check lat = -0.000007 lon =   0.000000
1696 tauref(610. Pa) = 0.463710 tau(610. Pa) = 0.463710
1697Ls =   1.063993 check lat = -0.000007 lon =   0.000000
1698 tauref(610. Pa) = 0.463714 tau(610. Pa) = 0.463714
1699 WRITEDIAGFI: date=    2.083333   
1700Ls =   1.069853 check lat = -0.000007 lon =   0.000000
1701 tauref(610. Pa) = 0.463720 tau(610. Pa) = 0.463720
1702Ls =   1.075065 check lat = -0.000007 lon =   0.000000
1703 tauref(610. Pa) = 0.463725 tau(610. Pa) = 0.463725
1704Ls =   1.080276 check lat = -0.000007 lon =   0.000000
1705 tauref(610. Pa) = 0.463730 tau(610. Pa) = 0.463730
1706Ls =   1.085481 check lat = -0.000007 lon =   0.000000
1707 tauref(610. Pa) = 0.463736 tau(610. Pa) = 0.463736
1708 WRITEDIAGFI: date=    2.125000   
17090**********    pas   2040     jour*****heure  3.0    date   2.12500    **********
1710          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1711        0.948574 0.860545E+00     0.950     1.043     0.964    17.996     0.829
1712
1713Ls =   1.091334 check lat = -0.000007 lon =   0.000000
1714 tauref(610. Pa) = 0.463741 tau(610. Pa) = 0.463741
1715Ls =   1.095904 check lat = -0.000007 lon =   0.000000
1716 tauref(610. Pa) = 0.463746 tau(610. Pa) = 0.463746
1717Ls =   1.101764 check lat = -0.000007 lon =   0.000000
1718 tauref(610. Pa) = 0.463751 tau(610. Pa) = 0.463751
1719Ls =   1.106962 check lat = -0.000007 lon =   0.000000
1720 tauref(610. Pa) = 0.463757 tau(610. Pa) = 0.463757
1721 WRITEDIAGFI: date=    2.166667   
1722Ls =   1.112180 check lat = -0.000007 lon =   0.000000
1723 tauref(610. Pa) = 0.463762 tau(610. Pa) = 0.463762
1724Ls =   1.117392 check lat = -0.000007 lon =   0.000000
1725 tauref(610. Pa) = 0.463767 tau(610. Pa) = 0.463767
1726Ls =   1.123245 check lat = -0.000007 lon =   0.000000
1727 tauref(610. Pa) = 0.463772 tau(610. Pa) = 0.463772
1728Ls =   1.127808 check lat = -0.000007 lon =   0.000000
1729 tauref(610. Pa) = 0.463777 tau(610. Pa) = 0.463777
1730 WRITEDIAGFI: date=    2.208333   
1731Ls =   1.133661 check lat = -0.000007 lon =   0.000000
1732 tauref(610. Pa) = 0.463783 tau(610. Pa) = 0.463783
1733Ls =   1.138873 check lat = -0.000007 lon =   0.000000
1734 tauref(610. Pa) = 0.463788 tau(610. Pa) = 0.463788
1735Ls =   1.144077 check lat = -0.000007 lon =   0.000000
1736 tauref(610. Pa) = 0.463793 tau(610. Pa) = 0.463793
1737Ls =   1.149289 check lat = -0.000007 lon =   0.000000
1738 tauref(610. Pa) = 0.463798 tau(610. Pa) = 0.463799
1739 WRITEDIAGFI: date=    2.250000   
17400**********    pas   2160     jour*****heure  6.0    date   2.25000    **********
1741          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1742        0.948658 0.957657E+00     0.950     1.091     0.964    17.950     0.831
1743
1744Ls =   1.155142 check lat = -0.000007 lon =   0.000000
1745 tauref(610. Pa) = 0.463804 tau(610. Pa) = 0.463804
1746Ls =   1.159705 check lat = -0.000007 lon =   0.000000
1747 tauref(610. Pa) = 0.463809 tau(610. Pa) = 0.463809
1748Ls =   1.165558 check lat = -0.000007 lon =   0.000000
1749 tauref(610. Pa) = 0.463814 tau(610. Pa) = 0.463814
1750Ls =   1.170770 check lat = -0.000007 lon =   0.000000
1751 tauref(610. Pa) = 0.463820 tau(610. Pa) = 0.463820
1752 WRITEDIAGFI: date=    2.291667   
1753Ls =   1.175974 check lat = -0.000007 lon =   0.000000
1754 tauref(610. Pa) = 0.463824 tau(610. Pa) = 0.463824
1755Ls =   1.181192 check lat = -0.000007 lon =   0.000000
1756 tauref(610. Pa) = 0.463830 tau(610. Pa) = 0.463830
1757Ls =   1.187039 check lat = -0.000007 lon =   0.000000
1758 tauref(610. Pa) = 0.463835 tau(610. Pa) = 0.463835
1759Ls =   1.191602 check lat = -0.000007 lon =   0.000000
1760 tauref(610. Pa) = 0.463840 tau(610. Pa) = 0.463840
1761 WRITEDIAGFI: date=    2.333333   
1762Ls =   1.197462 check lat = -0.000007 lon =   0.000000
1763 tauref(610. Pa) = 0.463846 tau(610. Pa) = 0.463846
1764Ls =   1.202667 check lat = -0.000007 lon =   0.000000
1765 tauref(610. Pa) = 0.463851 tau(610. Pa) = 0.463851
1766Ls =   1.207871 check lat = -0.000007 lon =   0.000000
1767 tauref(610. Pa) = 0.463856 tau(610. Pa) = 0.463856
1768Ls =   1.213083 check lat = -0.000007 lon =   0.000000
1769 tauref(610. Pa) = 0.463861 tau(610. Pa) = 0.463861
1770 WRITEDIAGFI: date=    2.375000   
17710**********    pas   2280     jour*****heure  9.0    date   2.37500    **********
1772          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1773        0.948753 0.977855E+00     0.951     1.100     0.964    18.044     0.831
1774
1775Ls =   1.218936 check lat = -0.000007 lon =   0.000000
1776 tauref(610. Pa) = 0.463867 tau(610. Pa) = 0.463867
1777Ls =   1.223492 check lat = -0.000007 lon =   0.000000
1778 tauref(610. Pa) = 0.463872 tau(610. Pa) = 0.463871
1779Ls =   1.229352 check lat = -0.000007 lon =   0.000000
1780 tauref(610. Pa) = 0.463877 tau(610. Pa) = 0.463877
1781Ls =   1.234563 check lat = -0.000007 lon =   0.000000
1782 tauref(610. Pa) = 0.463882 tau(610. Pa) = 0.463882
1783 WRITEDIAGFI: date=    2.416667   
1784Ls =   1.239761 check lat = -0.000007 lon =   0.000000
1785 tauref(610. Pa) = 0.463887 tau(610. Pa) = 0.463887
1786Ls =   1.244973 check lat = -0.000007 lon =   0.000000
1787 tauref(610. Pa) = 0.463893 tau(610. Pa) = 0.463893
1788Ls =   1.250819 check lat = -0.000007 lon =   0.000000
1789 tauref(610. Pa) = 0.463898 tau(610. Pa) = 0.463898
1790Ls =   1.255382 check lat = -0.000007 lon =   0.000000
1791 tauref(610. Pa) = 0.463903 tau(610. Pa) = 0.463903
1792 WRITEDIAGFI: date=    2.458333   
1793Ls =   1.261235 check lat = -0.000007 lon =   0.000000
1794 tauref(610. Pa) = 0.463908 tau(610. Pa) = 0.463908
1795Ls =   1.266447 check lat = -0.000007 lon =   0.000000
1796 tauref(610. Pa) = 0.463914 tau(610. Pa) = 0.463914
1797Ls =   1.271651 check lat = -0.000007 lon =   0.000000
1798 tauref(610. Pa) = 0.463919 tau(610. Pa) = 0.463919
1799Ls =   1.276863 check lat = -0.000007 lon =   0.000000
1800 tauref(610. Pa) = 0.463924 tau(610. Pa) = 0.463924
1801 WRITEDIAGFI: date=    2.500000   
18020**********    pas   2400     jour*****heure 12.0    date   2.50000    **********
1803          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1804        0.948844 0.863164E+00     0.951     1.073     0.964    18.211     0.830
1805
1806Ls =   1.282716 check lat = -0.000007 lon =   0.000000
1807 tauref(610. Pa) = 0.463929 tau(610. Pa) = 0.463930
1808Ls =   1.287265 check lat = -0.000007 lon =   0.000000
1809 tauref(610. Pa) = 0.463934 tau(610. Pa) = 0.463934
1810Ls =   1.293119 check lat = -0.000007 lon =   0.000000
1811 tauref(610. Pa) = 0.463940 tau(610. Pa) = 0.463940
1812Ls =   1.298337 check lat = -0.000007 lon =   0.000000
1813 tauref(610. Pa) = 0.463945 tau(610. Pa) = 0.463945
1814 WRITEDIAGFI: date=    2.541667   
1815Ls =   1.303535 check lat = -0.000007 lon =   0.000000
1816 tauref(610. Pa) = 0.463950 tau(610. Pa) = 0.463950
1817Ls =   1.308746 check lat = -0.000007 lon =   0.000000
1818 tauref(610. Pa) = 0.463955 tau(610. Pa) = 0.463955
1819Ls =   1.314593 check lat = -0.000007 lon =   0.000000
1820 tauref(610. Pa) = 0.463961 tau(610. Pa) = 0.463961
1821Ls =   1.319149 check lat = -0.000007 lon =   0.000000
1822 tauref(610. Pa) = 0.463966 tau(610. Pa) = 0.463966
1823 WRITEDIAGFI: date=    2.583333   
1824Ls =   1.325009 check lat = -0.000007 lon =   0.000000
1825 tauref(610. Pa) = 0.463971 tau(610. Pa) = 0.463971
1826Ls =   1.330213 check lat = -0.000007 lon =   0.000000
1827 tauref(610. Pa) = 0.463977 tau(610. Pa) = 0.463977
1828Ls =   1.335418 check lat = -0.000007 lon =   0.000000
1829 tauref(610. Pa) = 0.463982 tau(610. Pa) = 0.463981
1830Ls =   1.340629 check lat = -0.000007 lon =   0.000000
1831 tauref(610. Pa) = 0.463987 tau(610. Pa) = 0.463987
1832 WRITEDIAGFI: date=    2.625000   
18330**********    pas   2520     jour*****heure 15.0    date   2.62500    **********
1834          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1835        0.948931 0.874282E+00     0.951     1.079     0.964    18.370     0.830
1836
1837Ls =   1.346483 check lat = -0.000007 lon =   0.000000
1838 tauref(610. Pa) = 0.463992 tau(610. Pa) = 0.463992
1839Ls =   1.351032 check lat = -0.000007 lon =   0.000000
1840 tauref(610. Pa) = 0.463997 tau(610. Pa) = 0.463997
1841Ls =   1.356899 check lat = -0.000007 lon =   0.000000
1842 tauref(610. Pa) = 0.464003 tau(610. Pa) = 0.464003
1843Ls =   1.362097 check lat = -0.000007 lon =   0.000000
1844 tauref(610. Pa) = 0.464008 tau(610. Pa) = 0.464008
1845 WRITEDIAGFI: date=    2.666667   
1846Ls =   1.367295 check lat = -0.000007 lon =   0.000000
1847 tauref(610. Pa) = 0.464013 tau(610. Pa) = 0.464013
1848Ls =   1.372513 check lat = -0.000007 lon =   0.000000
1849 tauref(610. Pa) = 0.464018 tau(610. Pa) = 0.464018
1850Ls =   1.378353 check lat = -0.000007 lon =   0.000000
1851 tauref(610. Pa) = 0.464024 tau(610. Pa) = 0.464024
1852Ls =   1.382908 check lat = -0.000007 lon =   0.000000
1853 tauref(610. Pa) = 0.464029 tau(610. Pa) = 0.464028
1854 WRITEDIAGFI: date=    2.708333   
1855Ls =   1.388769 check lat = -0.000007 lon =   0.000000
1856 tauref(610. Pa) = 0.464034 tau(610. Pa) = 0.464034
1857Ls =   1.393966 check lat = -0.000007 lon =   0.000000
1858 tauref(610. Pa) = 0.464039 tau(610. Pa) = 0.464039
1859Ls =   1.399171 check lat = -0.000007 lon =   0.000000
1860 tauref(610. Pa) = 0.464044 tau(610. Pa) = 0.464044
1861Ls =   1.404382 check lat = -0.000007 lon =   0.000000
1862 tauref(610. Pa) = 0.464050 tau(610. Pa) = 0.464050
1863 WRITEDIAGFI: date=    2.750000   
18640**********    pas   2640     jour*****heure 18.0    date   2.75000    **********
1865          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1866        0.949009 0.982761E+00     0.951     1.068     0.964    18.391     0.831
1867
1868Ls =   1.410229 check lat = -0.000007 lon =   0.000000
1869 tauref(610. Pa) = 0.464055 tau(610. Pa) = 0.464055
1870Ls =   1.414785 check lat = -0.000007 lon =   0.000000
1871 tauref(610. Pa) = 0.464060 tau(610. Pa) = 0.464060
1872Ls =   1.420645 check lat = -0.000007 lon =   0.000000
1873 tauref(610. Pa) = 0.464065 tau(610. Pa) = 0.464065
1874Ls =   1.425843 check lat = -0.000007 lon =   0.000000
1875 tauref(610. Pa) = 0.464071 tau(610. Pa) = 0.464071
1876 WRITEDIAGFI: date=    2.791667   
1877Ls =   1.431048 check lat = -0.000007 lon =   0.000000
1878 tauref(610. Pa) = 0.464076 tau(610. Pa) = 0.464076
1879Ls =   1.436252 check lat = -0.000007 lon =   0.000000
1880 tauref(610. Pa) = 0.464081 tau(610. Pa) = 0.464081
1881Ls =   1.442106 check lat = -0.000007 lon =   0.000000
1882 tauref(610. Pa) = 0.464087 tau(610. Pa) = 0.464087
1883Ls =   1.446655 check lat = -0.000007 lon =   0.000000
1884 tauref(610. Pa) = 0.464091 tau(610. Pa) = 0.464091
1885 WRITEDIAGFI: date=    2.833333   
1886Ls =   1.452515 check lat = -0.000007 lon =   0.000000
1887 tauref(610. Pa) = 0.464097 tau(610. Pa) = 0.464097
1888Ls =   1.457713 check lat = -0.000007 lon =   0.000000
1889 tauref(610. Pa) = 0.464102 tau(610. Pa) = 0.464102
1890Ls =   1.462917 check lat = -0.000007 lon =   0.000000
1891 tauref(610. Pa) = 0.464107 tau(610. Pa) = 0.464107
1892Ls =   1.468122 check lat = -0.000007 lon =   0.000000
1893 tauref(610. Pa) = 0.464113 tau(610. Pa) = 0.464112
1894 WRITEDIAGFI: date=    2.875000   
18950**********    pas   2760     jour*****heure 21.0    date   2.87500    **********
1896          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1897        0.949064 0.871051E+00     0.950     1.078     0.964    18.339     0.829
1898
1899Ls =   1.473975 check lat = -0.000007 lon =   0.000000
1900 tauref(610. Pa) = 0.464118 tau(610. Pa) = 0.464118
1901Ls =   1.478531 check lat = -0.000007 lon =   0.000000
1902 tauref(610. Pa) = 0.464123 tau(610. Pa) = 0.464123
1903Ls =   1.484378 check lat = -0.000007 lon =   0.000000
1904 tauref(610. Pa) = 0.464128 tau(610. Pa) = 0.464128
1905Ls =   1.489575 check lat = -0.000007 lon =   0.000000
1906 tauref(610. Pa) = 0.464134 tau(610. Pa) = 0.464134
1907 WRITEDIAGFI: date=    2.916667   
1908Ls =   1.494794 check lat = -0.000007 lon =   0.000000
1909 tauref(610. Pa) = 0.464139 tau(610. Pa) = 0.464139
1910Ls =   1.499991 check lat = -0.000007 lon =   0.000000
1911 tauref(610. Pa) = 0.464144 tau(610. Pa) = 0.464144
1912Ls =   1.505838 check lat = -0.000007 lon =   0.000000
1913 tauref(610. Pa) = 0.464149 tau(610. Pa) = 0.464149
1914Ls =   1.510387 check lat = -0.000007 lon =   0.000000
1915 tauref(610. Pa) = 0.464154 tau(610. Pa) = 0.464154
1916 WRITEDIAGFI: date=    2.958333   
1917Ls =   1.516247 check lat = -0.000007 lon =   0.000000
1918 tauref(610. Pa) = 0.464160 tau(610. Pa) = 0.464160
1919Ls =   1.521445 check lat = -0.000007 lon =   0.000000
1920 tauref(610. Pa) = 0.464165 tau(610. Pa) = 0.464165
1921Ls =   1.526650 check lat = -0.000007 lon =   0.000000
1922 tauref(610. Pa) = 0.464170 tau(610. Pa) = 0.464170
1923Ls =   1.531854 check lat = -0.000007 lon =   0.000000
1924 tauref(610. Pa) = 0.464175 tau(610. Pa) = 0.464175
1925 PHYSIQ: for physdem ztime_fin =   1.000000   
1926 physdem1: ngridmx:         3010
1927 physdem1: skipping water vapour tracer
1928 Writing stats...
1929 OK variable aps                                               
1930 OK variable ps                                               
1931 OK variable tsurf                                             
1932 OK variable co2ice                                           
1933 OK variable tauref                                           
1934 OK variable fluxsurf_lw                                       
1935 OK variable fluxsurf_sw                                       
1936 OK variable fluxtop_lw                                       
1937 OK variable fluxtop_sw                                       
1938 OK variable temp                                             
1939 OK variable u                                                 
1940 OK variable v                                                 
1941 OK variable w                                                 
1942 OK variable rho                                               
1943 OK variable pressure                                         
1944 OK variable q2                                               
1945 OK variable emis                                             
1946 OK variable zmax_th                                           
1947 OK variable hfmax_th                                         
1948 OK variable wstar                                             
1949 OK variable vmr_h2ovap                                       
1950 OK variable vmr_h2oice                                       
1951 OK variable vmr_h2osat                                       
1952 OK variable h2o_ice_s                                         
1953 OK variable albedo                                           
1954 OK variable mtot                                             
1955 OK variable icetot                                           
1956 OK variable reffice                                           
1957 OK variable Nccntot                                           
1958 OK variable Mccntot                                           
1959 OK variable rice                                             
1960 OK variable tauTES                                           
1961 OK variable tau                                               
1962 OK variable dqsdust                                           
1963 OK variable dqndust                                           
1964 OK variable reffdust                                         
1965 OK variable dustq                                             
1966 OK variable dustN                                             
1967 OK variable ccnq                                             
1968 OK variable ccnN                                             
1969 WRITEDIAGFI: date=    3.000000   
1970  GCM: Appel test_period avant redem ; itau=        2880
1971 Enregistrement pour            1  0.0000000E+00
19720**********    pas   2880     jour*****heure  0.0    date   3.00000    **********
1973          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
1974        0.949109 0.103297E+01     0.950     1.096     0.964    18.331     0.829
1975
1976 in abort_gcm
1977 Stopping in gcm                 
1978 Reason =
1979 Simulation finished                                                           
1980 
1981 Everything is cool
Note: See TracBrowser for help on using the repository browser.