Index: unk/MESOSCALE/LMDZ.MARS.new/myGCM/log
===================================================================
--- /trunk/MESOSCALE/LMDZ.MARS.new/myGCM/log	(revision 1098)
+++ 	(revision )
@@ -1,1981 +1,0 @@
- iniadvtrac: Reading file traceur.def
-  Choix du shema  Van-Leer  pourle traceur no            1
-  Choix du shema  Van-Leer  pourle traceur no            2
-  Choix du shema  Van-Leer  pourle traceur no            3
-  Choix du shema  Van-Leer  pourle traceur no            4
-  Choix du shema  Van-Leer  pourle traceur no            5
-  Choix du shema  Van-Leer  pourle traceur no            6
-  Choix du shema  Van-Leer  pour la vapeur d'eau
- dynetat0: rad,omeg,g,cpp,kappa   3397200.      7.0776514E-05   3.720000    
-   744.4990      0.2567930    
- dynetat0: loading tracers
- DEFRUN_NEW: reading file run.def
- 
- Number of days to run:
- COUNT :           53
- getin_checkcohe : Found a problem on key iaervar
- getin_checkcohe : The following values were encoutered :
-                 iaervar == 
- 24                                                                             
-                      
-                 iaervar == 
- 26                                                                             
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :           59
- getin_checkcohe : Found a problem on key callnlte
- getin_checkcohe : The following values were encoutered :
-                 callnlte == 
- .true.                                                                         
-                      
-                 callnlte == 
- .false.                                                                        
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :           61
- getin_checkcohe : Found a problem on key nltemodel
- getin_checkcohe : The following values were encoutered :
-                 nltemodel == 
- 2                                                                              
-                      
-                 nltemodel == 
- 0                                                                              
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :           64
- getin_checkcohe : Found a problem on key nircorr
- getin_checkcohe : The following values were encoutered :
-                 nircorr == 
- 1                                                                              
-                      
-                 nircorr == 
- 0                                                                              
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :           95
- getin_checkcohe : Found a problem on key photochem
- getin_checkcohe : The following values were encoutered :
-                 photochem == 
- .true.                                                                         
-                      
-                 photochem == 
- .false.                                                                        
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :           97
- getin_checkcohe : Found a problem on key callthermos
- getin_checkcohe : The following values were encoutered :
-                 callthermos == 
- .true.                                                                         
-                      
-                 callthermos == 
- .false.                                                                        
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :          100
- getin_checkcohe : Found a problem on key callconduct
- getin_checkcohe : The following values were encoutered :
-                 callconduct == 
- .true.                                                                         
-                      
-                 callconduct == 
- .false.                                                                        
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :          102
- getin_checkcohe : Found a problem on key callmolvis
- getin_checkcohe : The following values were encoutered :
-                 callmolvis == 
- .true.                                                                         
-                      
-                 callmolvis == 
- .false.                                                                        
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :          104
- getin_checkcohe : Found a problem on key callmoldiff
- getin_checkcohe : The following values were encoutered :
-                 callmoldiff == 
- .true.                                                                         
-                      
-                 callmoldiff == 
- .false.                                                                        
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :          106
- getin_checkcohe : Found a problem on key thermochem
- getin_checkcohe : The following values were encoutered :
-                 thermochem == 
- .true.                                                                         
-                      
-                 thermochem == 
- .false.                                                                        
-                      
- getin_checkcohe : We will keep only the last value
- COUNT :          108
- getin_checkcohe : Found a problem on key calleuv
- getin_checkcohe : The following values were encoutered :
-                 calleuv == 
- .true.                                                                         
-                      
-                 calleuv == 
- .false.                                                                        
-                      
- getin_checkcohe : We will keep only the last value
-  nday =    3.000000    
- USING DEFAULTS : ndynstep =       -9999
- 
- Number of dynamical steps per day:(should be a multiple of iperiod)
-  day_step =          960
- 
- periode pour le pas Matsuno (en pas)
-  iperiod =            5
- 
- periode de sortie des variables de controle (en pas)
-  iconser =          120
- 
- periode de la dissipation (en pas)
-  idissip =            1
- 
- choix de l'operateur de dissipation (star ou  non star )
-  lstardis =  T
- 
- avec ou sans coordonnee hybrides
-  hybrid =  T
- 
- nombre d'iterations de l'operateur de dissipation   gradiv 
-  nitergdiv =            1
- 
- nombre d'iterations de l'operateur de dissipation  nxgradrot
-  nitergrot =            2
- 
- nombre d'iterations de l'operateur de dissipation  divgrad
-  niterh =            2
- 
- temps de dissipation des plus petites long.d ondes pour u,v (gradiv)
-  tetagdiv =    2000.000    
- 
- temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
-  tetagrot =    6000.000    
- 
- temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
-  tetatemp =    6000.000    
- 
- coefficient pour gamdissip
-  coefdis =   0.0000000E+00
- 
- choix du shema d'integration temporelle 
- (true==Matsuno ou false==Matsuno-leapfrog)
-  purmats =  F
- 
- avec ou sans physique
-  physic =  T
- 
- periode de la physique (en pas)
-  iphysiq =           10
- 
- choix d'une grille reguliere
-  grireg =  T
- 
- frequence (en pas) de l'ecriture du fichier diagfi.nc
-  ecritphy =    40.00000    
- 
- longitude en degres du centre du zoom
-  clon =    63.00000    
- 
- latitude en degres du centre du zoom 
-  clat =   0.0000000E+00
- 
- facteur de grossissement du zoom, selon longitude
-  grossismx =    1.000000    
- 
- facteur de grossissement du zoom , selon latitude
-  grossismy =    1.000000    
-  Defrun  alphax alphay    0.0000000E+00  0.0000000E+00
- 
- Fonction  f(y)  hyperbolique  si = .true., sinon  sinusoidale
-  fxyhypb =  F
- 
- extension en longitude de la zone du zoom (fraction de la zone totale)
-  dzoomx =   0.0000000E+00
- 
- extension en latitude de la zone du zoom (fraction de la zone totale)
-  dzoomy =   0.0000000E+00
- 
- raideur du zoom en  X
-  taux =    2.000000    
- 
- raideur du zoom en  Y
-  tauy =    2.000000    
- 
- Avec sponge layer
-  callsponge =  T
- 
- Sponge: number of layers over which sponge extends
- USING DEFAULTS : nsponge =           3
-  nsponge =            3
- 
- Sponge mode: (forcing is towards ...
-   over upper nsponge layers)
-   0: (h=hmean,u=v=0)
-   1: (h=hmean,u=umean,v=0)
-   2: (h=hmean,u=umean,v=vmean)
-  mode_sponge =            2
- 
- Sponge: characteristice time scale tetasponge
- (seconds) at topmost layer (time scale then 
-  doubles with decreasing layer index).
-  tetasponge =    50000.00    
- -----------------------------------------------
-  
-  
- -----------------------------------------------
-  
-  
- No zoom ? -> fxyhypb set to False           (It runs better that way)
-  R  CPP  Kappa    191.1821       744.4990      0.2567930    
- ****************************
- Lecture de z2sig.def
- ****************************
- *******************************
- Systeme en coordonnees hybrides
- 
-  BP 
-   1.000000      0.9990893      0.9961299      0.9887227      0.9736167    
-  0.9470957      0.9053387      0.8451603      0.7657363      0.6698347    
-  0.5632617      0.4541319      0.3506572      0.2588113      0.1820027    
-  0.1212465      7.5687021E-02  4.3286249E-02  2.1493321E-02  7.8194505E-03
-  6.6300534E-04  3.9641478E-11  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  AP 
-  0.0000000E+00  9.1087818E-03  3.8738251E-02  0.1130879      0.2655351    
-  0.5357349      0.9672797       1.600810       2.453982       3.502020    
-   4.676154       5.868705       6.962464       7.865741       8.523151    
-   8.916158       9.054893       8.963833       8.661031       8.111788    
-   6.921022       4.003914       2.106964       1.072680      0.5344583    
-  0.2654041      0.1317958      6.5447852E-02  3.2500453E-02  0.0000000E+00
-  BPs 
-  0.9995446      0.9976096      0.9924263      0.9811697      0.9603562    
-  0.9262172      0.8752495      0.8054483      0.7177855      0.6165482    
-  0.5086969      0.4023946      0.3047343      0.2204070      0.1516246    
-  9.8466769E-02  5.9486635E-02  3.2389786E-02  1.4656385E-02  4.2412281E-03
-  3.3150270E-04  1.9820739E-11  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  APs
-  4.5543909E-03  2.3923516E-02  7.5913072E-02  0.1893115      0.4006350    
-  0.7515073       1.284045       2.027396       2.978001       4.089087    
-   5.272429       6.415585       7.414103       8.194447       8.719654    
-   8.985525       9.009363       8.812431       8.386410       7.516405    
-   5.462468       3.055439       1.589822      0.8035691      0.3999312    
-  0.1985999      9.8621808E-02  4.8974153E-02  2.4319850E-02
-  PRESNIVS
-   609.7268       608.5658       605.4559       598.7028       586.2180    
-   565.7440       535.1863       493.3509       440.8271       380.1835    
-   315.5775       251.8763       193.3020       142.6427       101.2106    
-   69.05025       45.29621       28.57020       17.32681       10.10355    
-   5.664685       3.055439       1.589822      0.8035691      0.3999312    
-  0.1985999      9.8621808E-02  4.8974153E-02  2.4319850E-02
- Pseudo altitude des Presnivs : 
-  4.4802926E-03  2.3539927E-02  7.4772738E-02  0.1869370      0.3976738    
-  0.7531729       1.308441       2.122384       3.248062       4.728050    
-   6.590548       8.845211       11.49205       14.53116       17.96255    
-   21.78624       26.00236       30.61095       35.61204       41.00571    
-   46.79208       52.96536       59.49837       66.32151       73.29922    
-   80.29922       87.29922       94.29922       101.2992    
-
-
-           ....  INIGEOM  date du 01/06/98   ..... 
-
-        Calcul des elongations cu et cv  comme sommes des 4 
-      elong. cuij1, .. 4  , cvij1,.. 4  qui les entourent , aux  
-      memes endroits que les aires aireij1,...j4   . 
-
-  gamdi_gd   0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-           1           2           2
-
-
-
- *** Inigeom ,  Y = Latitude  , der. sinusoid . ***
-    ***  Coordonnees de la grille  *** 
-
-
-    LONGITUDES  aux pts.   V  ( degres )  
-
-
-  -180.00 -174.38 -168.75 -163.12 -157.50 -151.88 -146.25 -140.62
-  -135.00 -129.38 -123.75 -118.12 -112.50 -106.88 -101.25  -95.63
-   -90.00  -84.37  -78.75  -73.12  -67.50  -61.87  -56.25  -50.63
-   -45.00  -39.38  -33.75  -28.12  -22.50  -16.88  -11.25   -5.63
-     0.00    5.63   11.25   16.88   22.50   28.12   33.75   39.38
-    45.00   50.62   56.25   61.87   67.50   73.12   78.75   84.38
-    90.00   95.63  101.25  106.88  112.50  118.12  123.75  129.38
-   135.00  140.62  146.25  151.88  157.50  163.13  168.75  174.38
-   180.00
-
-
-    LATITUDES   aux pts.   V  ( degres )  
-
-
-    88.12   84.37   80.62   76.87   73.12   69.37   65.62   61.87
-    58.12   54.37   50.62   46.87   43.12   39.37   35.62   31.87
-    28.12   24.37   20.62   16.87   13.12    9.37    5.62    1.87
-    -1.88   -5.63   -9.38  -13.13  -16.88  -20.63  -24.38  -28.13
-   -31.88  -35.63  -39.38  -43.13  -46.88  -50.62  -54.38  -58.13
-   -61.88  -65.62  -69.38  -73.12  -76.88  -80.62  -84.37  -88.12
-
-
-    LONGITUDES  aux pts.   U  ( degres )  
-
-
-  -177.19 -171.56 -165.94 -160.31 -154.69 -149.06 -143.44 -137.81
-  -132.19 -126.56 -120.94 -115.31 -109.69 -104.06  -98.44  -92.81
-   -87.19  -81.56  -75.94  -70.31  -64.69  -59.06  -53.44  -47.81
-   -42.19  -36.56  -30.94  -25.31  -19.69  -14.06   -8.44   -2.81
-     2.81    8.44   14.06   19.69   25.31   30.94   36.56   42.19
-    47.81   53.44   59.06   64.69   70.31   75.94   81.56   87.19
-    92.81   98.44  104.06  109.69  115.31  120.94  126.56  132.19
-   137.81  143.44  149.06  154.69  160.31  165.94  171.56  177.19
-   182.81
-
-
-    LATITUDES   aux pts.   U  ( degres )  
-
-
-    90.00   86.25   82.50   78.75   75.00   71.25   67.50   63.75
-    60.00   56.25   52.50   48.75   45.00   41.25   37.50   33.75
-    30.00   26.25   22.50   18.75   15.00   11.25    7.50    3.75
-     0.00   -3.75   -7.50  -11.25  -15.00  -18.75  -22.50  -26.25
-   -30.00  -33.75  -37.50  -41.25  -45.00  -48.75  -52.50  -56.25
-   -60.00  -63.75  -67.50  -71.25  -75.00  -78.75  -82.50  -86.25
-   -90.00
-
-
-  EIGNVL 
- -0.940309E-05-0.999202E+00-0.999209E+00-0.398716E+01-0.398718E+01
- -0.893513E+01-0.893514E+01-0.157954E+02-0.157955E+02-0.245020E+02
- -0.245020E+02-0.349710E+02-0.349710E+02-0.471016E+02-0.471017E+02
- -0.607770E+02-0.607770E+02-0.758655E+02-0.758655E+02-0.922217E+02
- -0.922218E+02-0.109688E+03-0.109688E+03-0.128097E+03-0.128097E+03
- -0.147270E+03-0.147270E+03-0.167023E+03-0.167023E+03-0.187167E+03
- -0.187167E+03-0.207506E+03-0.207506E+03-0.227845E+03-0.227845E+03
- -0.247988E+03-0.247988E+03-0.267742E+03-0.267742E+03-0.286915E+03
- -0.286915E+03-0.305323E+03-0.305323E+03-0.322790E+03-0.322790E+03
- -0.339146E+03-0.339146E+03-0.354234E+03-0.354235E+03-0.367910E+03
- -0.367910E+03-0.380040E+03-0.380041E+03-0.390509E+03-0.390510E+03
- -0.399216E+03-0.399216E+03-0.406076E+03-0.406076E+03-0.411024E+03
- -0.411024E+03-0.414012E+03-0.414012E+03-0.415012E+03
-
-                Inifilr colat0 alphax    0.5000000E+00   0.0000000E+00
-  TRUNCATION AT           64
-  jfiltnv jfiltsv jfiltnu jfiltsu            8          41           8
-          41
-    Modes premiers  v  
-   4 10 14 20 26 32 40 52 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64
-  64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 52 40 32 26 20 14 10  4
-    Modes premiers  u  
-  64  6 12 18 24 30 36 46 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64
-  64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 64 46 36 30 24 18 12  6
- Dissipation : 
- Multiplication de la dissipation en altitude : fac_mid, fac_up =   3.000000    
-   30.00000    
- Transition mid /up : startalt, delta =   70.00000       20.00000     (km)
- calcul des valeurs propres
- Constantes de temps de la diffusion horizontale
-  INIDI tetamin dtvr    66.66684       92.47396               5
-  INIDI tetamin idissip    66.66684               1
-  INIDI idissip dtdiss            1   92.47396    
- pseudoZ(km)  zvert  dt(tetagdiv) dt(tetagrot) dt(divgrad)
-  4.4802926E-03   1.000001      4.6237011E-02  1.5412334E-02  1.5412334E-02
-  2.3539927E-02   1.000011      4.6237506E-02  1.5412500E-02  1.5412500E-02
-  7.4772738E-02   1.000113      4.6242204E-02  1.5414066E-02  1.5414066E-02
-  0.1869370       1.000712      4.6269897E-02  1.5423298E-02  1.5423298E-02
-  0.3976738       1.003286      4.6388932E-02  1.5462974E-02  1.5462974E-02
-  0.7531729       1.012164      4.6799421E-02  1.5599805E-02  1.5599805E-02
-   1.308441       1.038334      4.8009411E-02  1.6003134E-02  1.6003134E-02
-   2.122384       1.105890      5.1133037E-02  1.7044343E-02  1.7044343E-02
-   3.248062       1.256737      5.8107711E-02  1.9369235E-02  1.9369235E-02
-   4.728050       1.535234      7.0984595E-02  2.3661530E-02  2.3661530E-02
-   6.590548       1.930723      8.9270793E-02  2.9756924E-02  2.9756924E-02
-   8.845211       2.338095      0.1081065      3.6035489E-02  3.6035489E-02
-   11.49205       2.645828      0.1223351      4.0778361E-02  4.0778361E-02
-   14.53116       2.829569      0.1308307      4.3610230E-02  4.3610230E-02
-   17.96255       2.923871      0.1351909      4.5063645E-02  4.5063645E-02
-   21.78624       2.967935      0.1372284      4.5742787E-02  4.5742787E-02
-   26.00236       2.987214      0.1381198      4.6039917E-02  4.6039917E-02
-   30.61095       2.995183      0.1384882      4.6162728E-02  4.6162728E-02
-   35.61204       2.998292      0.1386320      4.6210654E-02  4.6210654E-02
-   41.00571       2.999433      0.1386847      4.6228237E-02  4.6228237E-02
-   46.79208       2.999824      0.1387028      4.6234269E-02  4.6234269E-02
-   52.96536       2.999952      0.1387087      4.6236239E-02  4.6236239E-02
-   59.49837       3.000110      0.1387160      4.6238672E-02  4.6238672E-02
-   66.32151       3.007357      0.1390511      4.6350367E-02  4.6350367E-02
-   73.29922       3.475960      0.1607179      5.3572625E-02  5.3572625E-02
-   80.29922       17.70859      0.8187916      0.2729305      0.2729305    
-   87.29922       29.66585       1.371659      0.4572197      0.4572197    
-   94.29922       29.99493       1.386875      0.4622915      0.4622915    
-   101.2992       29.99992       1.387106      0.4623686      0.4623686    
- PRESSION SURFACE
- MIN:    162.2194    
- MAX:    1119.796    
-55555555555555555555555555555555555555555555555555555555555555555
-66666666566655555566666666677766666666666666666666666666666666666
-66666667777766677777777888777777776666677777777777777777777777766
-77777777777777777888888888888888888888888888888887777777777777777
-77777777777777777778888888888888888888888888888877777777777777777
-77777777776766677777888888888888888888888888888777777777777777777
-77777766666666666777888888888888888888888888887777777777777777777
-77776666666666666667778888888888888888888888777777777777777777777
-77776666666556556666778888887777777778888887777777777777777777777
-77776666655555555666777888877777777777777777777777777777777777777
-77777666555555555566667788877777777777777777777777777777777777777
-77777766555544444555556777777777777777777777777777777777777777777
-77777766554333444444556777777777776767766666777777777787777776777
-77777776554222344444455777777777666666655566667777777787777776777
-77777776544322344444556777777776666665554456666777777777777766777
-77777776544333444444555777777766665555555555556677777777776666777
-77777765554333344445566777777666666555555544455677777777754556777
-77777765554333334445566777777665666555555544445667777777643456677
-77777765244333334445555777776656556555554444445666677777643456667
-66677764024333334454555677766555555555544444446666666777644556666
-66666665233333233454455667666565555555544444446776666666655566666
-66666665443331233454455567666555555554544443346776555666666666666
-66666665443332233454444556766655555544444443445776554566666666666
-66665655432222223443444455666655555544444443344665554456666666666
-66665655432212223333344455655655555444444343344444454445666666666
-55665554332112223324344445565655555544444333334344444444655666665
-55555544331111123345333555665555554443343333334334444434455556655
-55454443321112122235433445555555544443333334344333444333444455555
-44444443321122122233334664555555444433333344444333444334343445554
-43444443332222222333333444555555444443333344444333444333333344444
-44444443332222222333333444455554444443334445554443334433334444444
-43333433333222223333334444544444444344334455565544444433334444444
-43333443333222223332334444444444444344444468986655444443334444444
-443334433333222223233334444444444444334446***98766544443344444444
-444334433333332222233334455444444444334459***99876544444443444444
-444443333333333222333344555444444443334469***99986544444334443444
-334433333333333333333445775444444443334459**999986544444334444333
-33443333334333333333445677654444443334445799999975444444434444333
-44443433334443333444445666544444443333445678888654444444444443344
-44433333333343333444345565444443333333444455565554444444444333334
-43333334333343334444444444443333333333344444445443433344444333334
-33344443333333334444444444443333333333333444444443333344444333333
-33444433333333334444444443333333333333333344444333334444444433333
-33443333333333334344444333333333333333333333333333334444444434333
-33333333333333333333333333333333333333333333333333333444333333333
-33333333333333333333333333333333332333333333333333333333333333333
-33333333333333333333333333333333332223333333333333333333333333333
-33332222223333332222222222222222222222222222222223322222222232223
-22222222222222222222222222222222222222222222222222222222222222222
-1
-               run du pas      0  au pas   2880  c'est a dire du jour  13380   au jour  13383
-
-
- dynredem0: iim,jjm,llm,idayref          64          48          29       13383
- dynredem0: rad,omeg,g,cpp,kappa   3397200.      7.0776514E-05   3.720000    
-   744.4990      0.2567930    
- gcm: callgroupeun set to: F
-  GCM ---- Test_period apres continue   OK ! ----- itau:            0
-0**********    pas      0     jour*****heure  0.0    date   0.00000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947077 0.976960E+00     0.951     1.085     0.966    18.390     0.825
-
- 
- sponge mode           2
- nsponge tetasponge            3   50000.00    
- Coeffs for the sponge layer
- Z (km)     tau      cst
-   87.29922       200000.0      4.6236982E-04
-   94.29921       100000.0      9.2473964E-04
-   105.3311       50000.00      1.8494793E-03
- 
- in groupeun.F, ngroup=           2
- Shema  Amont nouveau  appele dans  Vanleer   
- 
- 
- --------------------------------------------
-  inifis: Parameters for the physics (callphys.def)
- --------------------------------------------
- Directory where external input files are:
-  datafile = /planeto/emlmd/run_mcd5.v0_new/datadir
- Run with or without tracer transport ?
-  tracer =  T
- Diurnal cycle ?
- (if diurnal=False, diurnal averaged solar heating)
-  diurnal =  T
- Seasonal cycle ?
- (if season=False, Ls stays constant, to value set in 'start'
-  season =  T
- Write some extra output to the screen ?
-  lwrite =  F
- Save statistics in file stats.nc ?
-  callstats =  T
- Save EOF profiles in file 'profiles' for Climate Database?
-  calleofdump =  F
- Dust scenario: 1=constant dust (read from startfi
-  or set as tauvis); 2=Viking scenario; =3 MGS scenario,
- =6 cold (low dust) scenario; =7 warm (high dust) scenario 
- =24,25 ... 30 :Mars Year 24, ... or 30 from TES assimilation
-  iaervar =           26
- Reference (visible) dust opacity at 610 Pa (matters only if iaervar=1)
-  tauvis =   0.1000000    
- Dust vertical distribution:
- (=1 top set by topdustref parameter; =2 Viking scenario; =3 MGS scenario)
-  iddist =            3
- Dust top altitude (km). (Matters only if iddist=1)
-  topdustref =    55.00000    
- Prescribed surface thermal flux (H/(rho*cp),K m/s)
- USING DEFAULTS : tke_heat_flux =  0.0000000E+00
-  tke_heat_flux =   0.0000000E+00
-  0 means the usual schemes are computing
- call radiative transfer ?
-  callrad =  T
- call slope insolation scheme ?(matters only if callrad=T)
- USING DEFAULTS : callslope = F
-  callslope =  F
- call NLTE radiative schemes ?(matters only if callrad=T)
-  callnlte =  F
- NLTE model?
- 0 -> old model, static O
- 1 -> old model, dynamic O
- 2 -> new model
- (matters only if callnlte=T)
-  nltemodel =            0
- call CO2 NIR absorption ?(matters only if callrad=T)
-  callnirco2 =  T
- New NIR NLTE correction ?0-> old model (no correction)1-> new correction
- (matters only if callnirco2=T)
-  nircorr =            0
- call turbulent vertical diffusion ?
-  calldifv =  T
- call thermals ?
-  calltherm =  T
- output thermal diagnostics ?
- USING DEFAULTS : outptherm = F
-  outptherm =  F
- call convective adjustment ?
-  calladj =  T
- call Richardson-based surface layer ?
-  callrichsl =  T
- call CO2 condensation ?
-  callcond =  T
- call thermal conduction in the soil ?
-  callsoil =  T
- call Lott's gravity wave/subgrid topography scheme ?
-  calllott =  T
- rad.transfer is computed every iradia physical timestep
-  iradia =            1
- Output of the exchange coefficient mattrix ?(for diagnostics only)
-  callg2d =  F
- Rayleigh scattering : (should be .false. for now)
-  rayleigh =  F
- Transported dust ? (if >0, use 'dustbin' dust bins)
-  dustbin =            2
- Radiatively active dust ? (matters if dustbin>0)
-  active =  T
- use mass and number mixing ratios to predict dust size ?
-  doubleq =  T
- USING DEFAULTS : submicron = F
-  submicron =  F
- dust lifted by GCM surface winds ?
-  lifting =  T
-  dust lifted by dust devils ?
-  callddevil =  F
- Gravitationnal sedimentation ?
-  sedimentation =  T
- Radiatively active transported atmospheric water ice ?
-  activice =  T
- Compute water cycle ?
-  water =  T
- Activate the thermal inertia feedback ?
- USING DEFAULTS : tifeedback = F
-  tifeedback =  F
- effective variance for water ice clouds ?
- water_param nueff Sedimentation:  0.1000000    
- water param CCN reduc. factor ?  0.0000000E+00
- USING DEFAULTS : ccn_factor =   4.500000    
-  ccn_factor =    4.500000    
- Careful: only used when microphys=F, otherwise
- the contact parameter is used instead;
- Microphysical scheme for water-ice clouds?
-  microphys =  T
- water contact parameter ?
- mteta =   0.9500000    
- Dust scavenging by H2O/CO2 snowfall ?
-  scavenging =  T
- Permanent water caps at poles ? .true. is RECOMMENDED
- (with .true., North cap is a source of water and South pole is a cold trap)
-  caps =  T
- water ice albedo ?
-  albedo_h2o_ice =   0.3500000    
- water ice thermal inertia ?
-  inert_h2o_ice =    800.0000    
- frost thickness threshold for albedo ?
-  frost_albedo_threshold =   4.9999999E-03
- Titus crocus line ?
- USING DEFAULTS : tituscap = F
- tituscap F
- photochemistry: include chemical species
-  photochem =  F
- call thermosphere ?
-  callthermos =  F
-  water included without cycle (only if water=.false.)
-  thermoswater =  F
- call thermal conduction ? (only if callthermos=.true.)
-  callconduct =  F
- call EUV heating ? (only if callthermos=.true.)
-  calleuv =  F
- call molecular viscosity ? (only if callthermos=.true.)
-  callmolvis =  F
- call molecular diffusion ? (only if callthermos=.true.)
-  callmoldiff =  F
- call thermospheric photochemistry ? (only if callthermos=.true.)
-  thermochem =  F
- Method to include solar variability
- 0-> old method (using solarcondate); 1-> variability wit E10.7
-  solvarmod =            0
- date for solar flux calculation: (1985 < date < 2002)
-  (Only used if solvarmod=0)
- (Solar min=1996.4 ave=1993.4 max=1990.6)
-  solarcondate =    1993.400    
- Solar variability as observed for MY: 
- Only if solvarmod=1
-  solvaryear =           24
- UV heating efficiency:measured values between 0.19 and 0.23 (Fox et al. 1996)
- lower values may be used to compensate low 15 um cooling
-  euveff =   0.2100000    
- --------------------------------------------
- 
- 
- 
- inifis: daysec   88775.00    
- 
- inifis: The radiative transfer is computed:
-            each            1  physical time-step
-         or each    924.7396      seconds
- 
- 
- TABFI de phyeta0           0           0
- tabfi: tab_cntrl   3010.000       29.00000       13380.00      0.0000000E+00
-   3397200.      7.0776514E-05   3.720000       43.49000      0.2567930    
-   88775.00       924.7396      0.0000000E+00  0.0000000E+00   669.0000    
-   206.6600       249.2200       485.0000       25.19000      9.9999998E-03
-   30.00000      1.0000000E-06  0.5800000      0.5600000      0.9000000    
-  0.9000000      0.9500000      0.1000000      0.0000000E+00  0.0000000E+00
-  0.0000000E+00  9.9999997E-05  9.9999997E-05  0.5000000      0.5000000    
-   1000000.      0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
-  0.0000000E+00
- IN tabfi day_ini=       13380
- *****************************************************
- Reading tab_cntrl when calling tabfi before changes
- *****************************************************
- (1)      = ngridmx?     3010.00     3010.00
- (2)            lmax       29.00       29.00
- (3)         day_ini    13380.00    13380.00
- (5)             rad  3397200.00  3397200.00
- (10)         daysec    88775.00    88775.00
- (6)            omeg   0.707765E-04   0.707765E-04
- (7)               g        3.72        3.72
- (8)           mugaz       43.49       43.49
- (9)             rcp        0.26        0.26
- (11)        dtphys?   0.924740E+03   0.924740E+03
- (14)       year_day      669.00      669.00
- (15)       periheli      206.66      206.66
- (16)        aphelie      249.22      249.22
- (17)       peri_day      485.00      485.00
- (18)       obliquit       25.19       25.19
- (19)     z0_default   0.100000E-01   0.100000E-01
- (21)      emin_turb   0.100000E-05   0.100000E-05
- (20)        lmixmin       30.00       30.00
- (26)        emissiv        0.95        0.95
- (24)     emisice(1)        0.90        0.90
- (25)     emisice(2)        0.90        0.90
- (22)    albedice(1)        0.58        0.58
- (23)    albedice(2)        0.56        0.56
- (31)   iceradius(1)   0.100000E-03   0.100000E-03
- (32)   iceradius(2)   0.100000E-03   0.100000E-03
- (33)   dtemisice(1)        0.50        0.50
- (34)   dtemisice(2)        0.50        0.50
- (27)         tauvis        0.10        0.10
- (35)        volcapa  1000000.00  1000000.00
- 
- Lmodif in tabfi!!!!!!!           0
- Geopotentiel au sol <phisfi>:  -26336.29       52485.20    
- Albedo du sol nu <albedodat>:  0.1003705      0.3208752    
- <zmea>:  -7079.710       14109.67    
- <zstd>:   8.602325       5754.301    
- <zsig>:  1.4299675E-04  6.9831863E-02
- <zgam>:  1.0000000E-07   1.000000    
- <zthe>:  -89.90435       89.88921    
- CO2 ice cover <co2ice>:  0.0000000E+00   958.8163    
- phyetat0: Le champ <tsurf> est present
-           J ignore donc les autres temperatures TS**
- Temperature du sol <tsurf>   143.1560       291.0728    
- Surface emissivity <emis>:  0.5347768      0.9500000    
- Surface roughness <z0>:  7.9248937E-05  1.7229829E-02
- pbl wind variance <q2>:  8.3298417E-08   157.1757    
- tracer on surface <co2>:  0.0000000E+00  0.0000000E+00
- tracer on surface <dust_number>: -2.1804593E+16  3.8433894E+14
- tracer on surface <dust_mass>:  -1790.340       634.2337    
- tracer on surface <ccn_number>:  1.5558199E+09  2.8858525E+16
- tracer on surface <ccn_mass>:  4.8576188E-05   910.5092    
- tracer on surface <h2o_ice>:  -40.61012       13.81997    
- phyetat0: loading surface tracer h2o_ice instead of h2o_vap
- tracer on surface <h2o_ice>:  -40.61012       13.81997    
- 
- Soil volumetric heat capacity:   1000000.    
- Soil thermal inertia <inertiedat>:   40.57156       1200.000    
- Soil temperature <tsoil>:   143.1567       289.7982    
- In physiq day_ini =       13380
- initracer: found all expected tracers, namely:
-                  1  co2
-                  2  dust_number
-                  3  dust_mass
-                  4  ccn_number
-                  5  ccn_mass
-                  6  h2o_ice
-                  7  h2o_vap
- initracer: doubleq_param reff_lift:  3.0000001E-06
- initracer: doubleq_param nueff_lift:  0.5000000    
- initracer: doubleq_param alpha_lift:  1.0000000E-06
- 
- ******** initracer : dust transport parameters :
- alpha_lift =   0.0000000E+00  1.1936620E+07  1.0000000E-06  1.0000000E-30
-  1.0000000E-30  0.0000000E+00  0.0000000E+00
- alpha_devil =   0.0000000E+00   107429.6      8.9999999E-09  1.0000000E-30
-  1.0000000E-30  0.0000000E+00  0.0000000E+00
- radius  =   0.0000000E+00  3.0000001E-06  3.0000001E-06  3.0000001E-06
-  3.0000001E-06  3.0000001E-06  0.0000000E+00
- reff_lift (um) =    3.0000001E-06
- size distribution variance  =   0.6367614    
- r3n_q , ref_r0 :   1.5401418E-05   2.755676    
- Ice dryness ?
- USING DEFAULTS : icedryness =   1.000000    
-  icedryness =    1.000000    
- Surfini: predefined ice caps
- water ice caps distribution for 64x48 resolution
- latitude | longitude | ig
- ice water cap   90.00000      0.0000000E+00           1
- ice water cap   86.24999      -180.0000               2
- ice water cap   86.24999      -174.3750               3
- ice water cap   86.24999      -168.7500               4
- ice water cap   86.24999      -163.1250               5
- ice water cap   86.24999      -157.5000               6
- ice water cap   86.24999      -151.8750               7
- ice water cap   86.24999      -146.2500               8
- ice water cap   86.24999      -140.6250               9
- ice water cap   86.24999      -135.0000              10
- ice water cap   86.24999      -129.3750              11
- ice water cap   86.24999      -123.7500              12
- ice water cap   86.24999      -118.1250              13
- ice water cap   86.24999      -112.5000              14
- ice water cap   86.24999      -106.8750              15
- ice water cap   86.24999      -101.2500              16
- ice water cap   86.24999      -95.62501              17
- ice water cap   86.24999      -90.00000              18
- ice water cap   86.24999      -84.37499              19
- ice water cap   86.24999      -78.75000              20
- ice water cap   86.24999      -73.12500              21
- ice water cap   86.24999      -67.50000              22
- ice water cap   86.24999      -61.87499              23
- ice water cap   86.24999      -56.24999              24
- ice water cap   86.24999      -50.62501              25
- ice water cap   86.24999      -45.00000              26
- ice water cap   86.24999      -39.37500              27
- ice water cap   86.24999      -33.75000              28
- ice water cap   86.24999      -28.12500              29
- ice water cap   86.24999      -22.49999              30
- ice water cap   86.24999      -16.87501              31
- ice water cap   86.24999      -11.25000              32
- ice water cap   86.24999      -5.625002              33
- ice water cap   86.24999      0.0000000E+00          34
- ice water cap   86.24999       5.625002              35
- ice water cap   86.24999       11.25000              36
- ice water cap   86.24999       16.87501              37
- ice water cap   86.24999       22.49999              38
- ice water cap   86.24999       28.12500              39
- ice water cap   86.24999       33.75000              40
- ice water cap   86.24999       39.37500              41
- ice water cap   86.24999       45.00000              42
- ice water cap   86.24999       50.62499              43
- ice water cap   86.24999       56.25001              44
- ice water cap   86.24999       61.87499              45
- ice water cap   86.24999       67.50001              46
- ice water cap   86.24999       73.12500              47
- ice water cap   86.24999       78.74998              48
- ice water cap   86.24999       84.37500              49
- ice water cap   86.24999       89.99999              50
- ice water cap   86.24999       95.62501              51
- ice water cap   86.24999       101.2500              52
- ice water cap   86.24999       106.8750              53
- ice water cap   86.24999       112.5000              54
- ice water cap   86.24999       118.1250              55
- ice water cap   86.24999       123.7500              56
- ice water cap   86.24999       129.3750              57
- ice water cap   86.24999       135.0000              58
- ice water cap   86.24999       140.6250              59
- ice water cap   86.24999       146.2500              60
- ice water cap   86.24999       151.8750              61
- ice water cap   86.24999       157.5000              62
- ice water cap   86.24999       163.1250              63
- ice water cap   86.24999       168.7500              64
- ice water cap   86.24999       174.3750              65
- ice water cap   82.49999      -129.3750              75
- ice water cap   82.49999      -123.7500              76
- ice water cap   82.49999      -118.1250              77
- ice water cap   82.49999      -112.5000              78
- ice water cap   82.49999      -106.8750              79
- ice water cap   82.49999      -101.2500              80
- ice water cap   82.49999      -95.62501              81
- ice water cap   82.49999      -90.00000              82
- ice water cap   82.49999      -84.37499              83
- ice water cap   82.49999      -78.75000              84
- ice water cap   82.49999      -73.12500              85
- ice water cap   82.49999      -28.12500              93
- ice water cap   82.49999      -22.49999              94
- ice water cap   82.49999      -16.87501              95
- ice water cap   82.49999      -11.25000              96
- ice water cap   82.49999      -5.625002              97
- ice water cap   82.49999      0.0000000E+00          98
- ice water cap   82.49999       5.625002              99
- ice water cap   82.49999       11.25000             100
- ice water cap   82.49999       16.87501             101
- ice water cap   82.49999       22.49999             102
- ice water cap   82.49999       28.12500             103
- ice water cap   82.49999       33.75000             104
- ice water cap   82.49999       39.37500             105
- ice water cap   82.49999       45.00000             106
- ice water cap   82.49999       50.62499             107
- ice water cap   82.49999       56.25001             108
- ice water cap   82.49999       61.87499             109
- ice water cap   82.49999       67.50001             110
- ice water cap   82.49999       73.12500             111
- ice water cap   82.49999       78.74998             112
- ice water cap   82.49999       84.37500             113
- ice water cap   82.49999       89.99999             114
- ice water cap   78.75000      -146.2500             136
- ice water cap   78.75000      -135.0000             138
- ice water cap   78.75000      -123.7500             140
- ice water cap   78.75000      -112.5000             142
- ice water cap   78.75000      -5.625002             161
- ice water cap   78.75000       5.625002             163
- ice water cap   78.75000       16.87501             165
- ice water cap   78.75000       118.1250             183
- ice water cap   78.75000       129.3750             185
- ice water cap   78.75000       140.6250             187
- ice water cap   78.75000       151.8750             189
- ice water cap   78.75000       163.1250             191
- ice water cap   78.75000       174.3750             193
- ice water cap   75.00000      -180.0000             194
- ice water cap   75.00000      -129.3750             203
- ice water cap   75.00000      -106.8750             207
- ice water cap   75.00000       101.2500             244
- ice water cap   75.00000       112.5000             246
- ice water cap   75.00000       135.0000             250
- ice water cap   75.00000       146.2500             252
- ice water cap   75.00000       157.5000             254
- ice water cap   75.00000       168.7500             256
- surfini: water frost thickness  4.9999999E-03
- surfini: water ice albedo:  0.3500000    
- surfini: water ice TI:   800.0000    
- minimum albedo sans water caps  0.1003705    
- maximum albedo sans water caps  0.3208752    
- minimum albedo avec water caps  0.1003705    
- maximum albedo avec water caps  0.3500000    
- minimum albedo avec givre et co2  0.1008788    
- maximum albedo avec givre et co2  0.5800000    
- Perihelie en Mkm     206.6600    
- Aphelise  en Mkm     249.2200    
- obliquite en degres  :   25.19000    
- e_elips  9.3357898E-02
- p_elips   1.510404    
- zanom     1.728111    
- zx0      1.818617    
- longitude solaire du perihelie timeperi =    1.908231    
- physiq: water_param Surface water ice albedo:  0.3500000    
- nlthermeq: transition to NLTE centred at   0.10Pa
-               half-width (scale heights)   0.50
-          suggested LTE coverage at least   0.04Pa
- LTE rad. calculations up to layer           29
- callradite: found all scatterers, namely:
-                  1  dust_doubleq
-                  2  h2o_ice
- Scatterer: dust_doubleq
-   corresponding files: 
- VIS: optprop_dustvis_TM_n50.dat
- IR : optprop_dustir_n50.dat
- Scatterer: h2o_ice
-   corresponding files: 
- VIS: optprop_icevis_n30.dat
- IR : optprop_iceir_n30.dat
- Splitting radiative calculations:  ngridmx,ngrid,ndomainsz,ndomain        3010
-        3010         151          20
- aeroptproperties: interpolation grid
- Scatterer:            1
- Domain:            1
- Min radius (m):   3.5353878E-08
- Max radius (m):   9.9999997E-06
- Scatterer:            1
- Domain:            2
- Min radius (m):   3.5353878E-08
- Max radius (m):   9.9999997E-06
- Scatterer:            2
- Domain:            1
- Min radius (m):   1.1303746E-07
- Max radius (m):   2.9999999E-05
- Scatterer:            2
- Domain:            2
- Min radius (m):   1.1303746E-07
- Max radius (m):   2.9999999E-05
- aeropacity: i_ice=           6
- Typical profiles of Qext(vis)/Qext(IR)
-   and Qext(IR)/Qabs(IR):
- Aerosol #            1
- solsir:    1.600077       1.615189       1.632739       1.656156    
-   1.673179       1.688357       1.700979       1.712606       1.732043    
-   1.767485       1.846160       1.956535       2.126256       2.416159    
-   2.924124       3.627217       3.967699       4.619432       4.999025    
-   5.531231       5.867186       6.121830       6.161941       5.507154    
-   3.473106       1.283761      0.4759978      0.3981150      0.3957708    
- Qext/Qabs(IR):    1.397849       1.393870       1.389535       1.383844    
-   1.379724       1.376071       1.373089       1.370418       1.366129    
-   1.359017       1.344959       1.327367       1.303859       1.271339    
-   1.229796       1.184386       1.146510       1.112741       1.088421    
-   1.071170       1.051863       1.034979       1.018531       1.007443    
-   1.002272       1.000459       1.000086       1.000062       1.000061    
- Aerosol #            2
- solsir:    1.576552       1.570121       1.565980       1.568253    
-   1.568665       1.579205       1.586318       1.594997       1.588596    
-   1.593094       1.581340       1.569868       1.526461       1.506043    
-   1.604984       1.702579       1.836655       1.946286       1.991414    
-   2.106776       2.331118       2.288174       2.249288       2.215015    
-   2.107166       2.065983       2.065983       2.065983       2.065983    
- Qext/Qabs(IR):    1.302334       1.306394       1.309110       1.314693    
-   1.311184       1.305743       1.300053       1.297721       1.303020    
-   1.308674       1.323397       1.357420       1.410284       1.449922    
-   1.444335       1.432657       1.395102       1.358613       1.356317    
-   1.338657       1.233698       1.142038       1.084853       1.030944    
-   1.005017       1.000078       1.000078       1.000078       1.000078    
-Ls =   0.004706 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.593544 tau(610. Pa) = 0.593544
- In sugwd NKTOPG=           7
- In vdifc: Tcond(P=1mb)=   136.2700      Lcond=   590000.0    
-           acond,bcond,ccond  3.2403751E-04  7.3383721E-03  3.3921041E-04
- USING DEFAULTS : dtke_thermals = F
- USING DEFAULTS : nsplit_thermals =          10
- USING DEFAULTS : r_aspect_thermals =   1.000000    
- Newcondens: improved_ztcond= T
- Newcondens: bound_qco2= F
- In newcondens:Tcond(P=1mb)=   136.2700      Lcond=   590000.0    
- acond,bcond,ccond  3.2403751E-04  7.3383721E-03  3.3921041E-04
-  albedocaps: Use TES Cap albedoes ?
-  albedocaps: TESicealbedo =  T
-  albedocaps: Coefficient for Northern Cap ?
-  albedocaps: TESice_Ncoef =    1.600000    
-  albedocaps: Coefficient for Southern Cap ?
-  albedocaps: TESice_Scoef =    1.600000    
- watercloud: igcm_h2o_vap=           7
-             igcm_h2o_ice=           6
- time subsampling for microphysic ?
- USING DEFAULTS : imicro =          15
- imicro =           15
- Physical timestep is   924.7396    
- Microphysics timestep is   61.64931    
- vrat_cld   31.6227757253798     
-  
- Microphysics: size bin information:
- i,rb_cld(i), rad_cld(i),dr_cld(i)
- -----------------------------------
- 1   0.100000E-09    0.100000E-06    0.124591E-06
- 2   0.124691E-06    0.316228E-06    0.269617E-06
- 3   0.394309E-06    0.100000E-05    0.852605E-06
- 4   0.124691E-05    0.316228E-05    0.269617E-05
- 5   0.394309E-05    0.100000E-04    0.999606E-02
- 6   0.100000E-01
- -----------------------------------
- water_param contact parameter:  0.9500000    
- Variance of ice & CCN distribs :  0.3087235    
- nuice for sedimentation:  0.1000000    
- Volume of a water molecule:  1.956521739130435E-005
-  
- dear user, please keep in mind that
- contact parameter IS constant
-  
- correction for the shape of the ice particles ?
- USING DEFAULTS : ice_shape =  0.7500000    
-  ice_shape =   0.7500000    
- water_param nueff Sedimentation:  0.1000000    
- water_param nueff Radiative:  0.1000000    
- 
-                         || STATS ||
- 
- daysec   88775.00    
- dtphys   924.7396    
- nsteppd=          96
- istats=           8
- Storing           12 times per day
- thus every            8 physical timestep 
- 
- iniwrite: iim,jjm,llm,idayref          64          48          29       13380
- iniwrite: rad,omeg,g,mugaz,rcp   3397200.      7.0776514E-05   3.720000    
-   43.49000      0.2567930    
- iniwrite: daysec,dtphys   88775.00       924.7396    
-Ls =   0.009924 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.592217 tau(610. Pa) = 0.592217
-Ls =   0.015805 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.590757 tau(610. Pa) = 0.590757
-Ls =   0.021037 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.589297 tau(610. Pa) = 0.589297
- WRITEDIAGFI: date=   4.1666668E-02
- ==========================
- DIAGFI: creating variable tsurf
- ==========================
- DIAGFI: creating variable ps
- ==========================
- DIAGFI: creating variable co2ice
- ==========================
- DIAGFI: creating variable temp
- ==========================
- DIAGFI: creating variable u
- ==========================
- DIAGFI: creating variable v
- ==========================
- DIAGFI: creating variable emis
- ==========================
- DIAGFI: creating variable tsoil
- ==========================
- DIAGFI: creating variable inertiedat
- ==========================
- DIAGFI: creating variable co2
- ==========================
- DIAGFI: creating variable co2col
- ==========================
- DIAGFI: creating variable q01
- ==========================
- DIAGFI: creating variable q02
- ==========================
- DIAGFI: creating variable qsurf02
- ==========================
- DIAGFI: creating variable mtot
- ==========================
- DIAGFI: creating variable icetot
- ==========================
- DIAGFI: creating variable vmr_h2oice
- ==========================
- DIAGFI: creating variable vmr_h2ovap
- ==========================
- DIAGFI: creating variable reffice
- ==========================
- DIAGFI: creating variable Nccntot
- ==========================
- DIAGFI: creating variable Mccntot
- ==========================
- DIAGFI: creating variable rice
- ==========================
- DIAGFI: creating variable h2o_ice_s
- ==========================
- DIAGFI: creating variable albedo
- ==========================
- DIAGFI: creating variable tauref
- ==========================
- DIAGFI: creating variable tau
- ==========================
- DIAGFI: creating variable dqsdust
- ==========================
- DIAGFI: creating variable dqndust
- ==========================
- DIAGFI: creating variable reffdust
- ==========================
- DIAGFI: creating variable dustq
- ==========================
- DIAGFI: creating variable dustN
- ==========================
- DIAGFI: creating variable ccn
- ==========================
- DIAGFI: creating variable ccnN
- ==========================
- DIAGFI: creating variable ccnq
- ==========================
- DIAGFI: creating variable zmax_th
- ==========================
- DIAGFI: creating variable hfmax_th
- ==========================
- DIAGFI: creating variable wstar
-Ls =   0.026262 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.587969 tau(610. Pa) = 0.587969
-Ls =   0.031480 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.586509 tau(610. Pa) = 0.586509
-Ls =   0.036719 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.585049 tau(610. Pa) = 0.585049
-Ls =   0.041944 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.583722 tau(610. Pa) = 0.583722
- =====================
- STATS: creation de ps
- =====================
- STATS: creation de tsurf
- =====================
- STATS: creation de co2ice
- =====================
- STATS: creation de tauref
- =====================
- STATS: creation de fluxsurf_lw
- =====================
- STATS: creation de fluxsurf_sw
- =====================
- STATS: creation de fluxtop_lw
- =====================
- STATS: creation de fluxtop_sw
- =====================
- STATS: creation de temp
- =====================
- STATS: creation de u
- =====================
- STATS: creation de v
- =====================
- STATS: creation de w
- =====================
- STATS: creation de rho
- =====================
- STATS: creation de pressure
- =====================
- STATS: creation de q2
- =====================
- STATS: creation de emis
- =====================
- STATS: creation de zmax_th
- =====================
- STATS: creation de hfmax_th
- =====================
- STATS: creation de wstar
- =====================
- STATS: creation de vmr_h2ovap
- =====================
- STATS: creation de vmr_h2oice
- =====================
- STATS: creation de vmr_h2osat
- =====================
- STATS: creation de h2o_ice_s
- =====================
- STATS: creation de albedo
- =====================
- STATS: creation de mtot
- =====================
- STATS: creation de icetot
- =====================
- STATS: creation de reffice
- =====================
- STATS: creation de Nccntot
- =====================
- STATS: creation de Mccntot
- =====================
- STATS: creation de rice
- =====================
- STATS: creation de tauTES
- =====================
- STATS: creation de tau
- =====================
- STATS: creation de dqsdust
- =====================
- STATS: creation de dqndust
- =====================
- STATS: creation de reffdust
- =====================
- STATS: creation de dustq
- =====================
- STATS: creation de dustN
- =====================
- STATS: creation de ccnq
- =====================
- STATS: creation de ccnN
- WRITEDIAGFI: date=   8.3333336E-02
-Ls =   0.047169 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.582262 tau(610. Pa) = 0.582262
-Ls =   0.053043 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.580801 tau(610. Pa) = 0.580801
-Ls =   0.058268 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.579474 tau(610. Pa) = 0.579474
-Ls =   0.063500 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.578014 tau(610. Pa) = 0.578014
- WRITEDIAGFI: date=   0.1250000    
-0**********    pas    120     jour*****heure  3.0    date   0.12500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947168 0.845509E+00     0.951     1.076     0.966    18.279     0.828
-
-Ls =   0.068732 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.576554 tau(610. Pa) = 0.576554
-Ls =   0.073950 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.575226 tau(610. Pa) = 0.575226
-Ls =   0.079176 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.573766 tau(610. Pa) = 0.573766
-Ls =   0.085056 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.572306 tau(610. Pa) = 0.572306
- WRITEDIAGFI: date=   0.1666667    
-Ls =   0.090275 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.570979 tau(610. Pa) = 0.570979
-Ls =   0.095507 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.569519 tau(610. Pa) = 0.569519
-Ls =   0.100732 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.568058 tau(610. Pa) = 0.568058
-Ls =   0.105957 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.566731 tau(610. Pa) = 0.566731
- WRITEDIAGFI: date=   0.2083333    
-Ls =   0.111182 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.565271 tau(610. Pa) = 0.565271
-Ls =   0.117056 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.563811 tau(610. Pa) = 0.563811
-Ls =   0.122281 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.562483 tau(610. Pa) = 0.562483
-Ls =   0.127513 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.561023 tau(610. Pa) = 0.561023
- WRITEDIAGFI: date=   0.2500000    
-0**********    pas    240     jour*****heure  6.0    date   0.25000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947275 0.894936E+00     0.951     1.116     0.966    18.153     0.829
-
-Ls =   0.132738 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.559563 tau(610. Pa) = 0.559563
-Ls =   0.137956 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.558236 tau(610. Pa) = 0.558236
-Ls =   0.143188 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.556776 tau(610. Pa) = 0.556776
-Ls =   0.149062 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.555315 tau(610. Pa) = 0.555315
- WRITEDIAGFI: date=   0.2916667    
-Ls =   0.154287 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.553988 tau(610. Pa) = 0.553988
-Ls =   0.159519 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.552528 tau(610. Pa) = 0.552528
-Ls =   0.164731 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.551068 tau(610. Pa) = 0.551068
-Ls =   0.169956 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.549740 tau(610. Pa) = 0.549740
- WRITEDIAGFI: date=   0.3333333    
-Ls =   0.175181 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.548280 tau(610. Pa) = 0.548280
-Ls =   0.181061 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.546820 tau(610. Pa) = 0.546820
-Ls =   0.185638 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.545493 tau(610. Pa) = 0.545493
-Ls =   0.191512 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.544033 tau(610. Pa) = 0.544032
- WRITEDIAGFI: date=   0.3750000    
-0**********    pas    360     jour*****heure  9.0    date   0.37500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947374 0.906323E+00     0.952     1.101     0.966    18.186     0.827
-
-Ls =   0.196730 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.542572 tau(610. Pa) = 0.542572
-Ls =   0.201962 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.541245 tau(610. Pa) = 0.541245
-Ls =   0.207173 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.539785 tau(610. Pa) = 0.539785
-Ls =   0.213054 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.538325 tau(610. Pa) = 0.538325
- WRITEDIAGFI: date=   0.4166667    
-Ls =   0.217630 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.536997 tau(610. Pa) = 0.536997
-Ls =   0.223504 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.535537 tau(610. Pa) = 0.535537
-Ls =   0.228729 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.534077 tau(610. Pa) = 0.534077
-Ls =   0.233954 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.532750 tau(610. Pa) = 0.532750
- WRITEDIAGFI: date=   0.4583333    
-Ls =   0.239166 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.531290 tau(610. Pa) = 0.531289
-Ls =   0.245054 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.529829 tau(610. Pa) = 0.529829
-Ls =   0.249616 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.528502 tau(610. Pa) = 0.528502
-Ls =   0.255497 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.527042 tau(610. Pa) = 0.527042
- WRITEDIAGFI: date=   0.5000000    
-0**********    pas    480     jour*****heure 12.0    date   0.50000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947464 0.101995E+01     0.952     1.069     0.966    18.359     0.826
-
-Ls =   0.260722 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.525582 tau(610. Pa) = 0.525582
-Ls =   0.265933 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.524254 tau(610. Pa) = 0.524254
-Ls =   0.271159 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.522794 tau(610. Pa) = 0.522794
-Ls =   0.277032 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.521334 tau(610. Pa) = 0.521334
- WRITEDIAGFI: date=   0.5416667    
-Ls =   0.281609 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.520007 tau(610. Pa) = 0.520007
-Ls =   0.287483 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.518546 tau(610. Pa) = 0.518546
-Ls =   0.292701 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.517086 tau(610. Pa) = 0.517086
-Ls =   0.297926 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.515759 tau(610. Pa) = 0.515759
- WRITEDIAGFI: date=   0.5833333    
-Ls =   0.303151 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.514299 tau(610. Pa) = 0.514299
-Ls =   0.309018 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.512839 tau(610. Pa) = 0.512839
-Ls =   0.313601 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.511511 tau(610. Pa) = 0.511511
-Ls =   0.319468 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.510051 tau(610. Pa) = 0.510051
- WRITEDIAGFI: date=   0.6250000    
-0**********    pas    600     jour*****heure 15.0    date   0.62500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947557 0.882835E+00     0.952     1.055     0.967    18.467     0.826
-
-Ls =   0.324694 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.508591 tau(610. Pa) = 0.508591
-Ls =   0.329905 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.507264 tau(610. Pa) = 0.507263
-Ls =   0.335130 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.505803 tau(610. Pa) = 0.505803
-Ls =   0.341004 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.504343 tau(610. Pa) = 0.504343
- WRITEDIAGFI: date=   0.6666667    
-Ls =   0.345573 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.503016 tau(610. Pa) = 0.503016
-Ls =   0.351447 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.501556 tau(610. Pa) = 0.501556
-Ls =   0.356666 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.500096 tau(610. Pa) = 0.500095
-Ls =   0.361884 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.498768 tau(610. Pa) = 0.498768
- WRITEDIAGFI: date=   0.7083333    
-Ls =   0.367109 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.497308 tau(610. Pa) = 0.497308
-Ls =   0.372983 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.495848 tau(610. Pa) = 0.495848
-Ls =   0.377552 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.494521 tau(610. Pa) = 0.494521
-Ls =   0.383426 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.493060 tau(610. Pa) = 0.493060
- WRITEDIAGFI: date=   0.7500000    
-0**********    pas    720     jour*****heure 18.0    date   0.75000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947652 0.947625E+00     0.952     1.055     0.966    18.448     0.826
-
-Ls =   0.388651 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.491600 tau(610. Pa) = 0.491600
-Ls =   0.393870 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.490273 tau(610. Pa) = 0.490273
-Ls =   0.399088 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.488813 tau(610. Pa) = 0.488813
-Ls =   0.404962 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.487353 tau(610. Pa) = 0.487352
- WRITEDIAGFI: date=   0.7916667    
-Ls =   0.409531 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.486025 tau(610. Pa) = 0.486025
-Ls =   0.415405 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.484565 tau(610. Pa) = 0.484565
-Ls =   0.420617 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.483105 tau(610. Pa) = 0.483105
-Ls =   0.425849 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.481778 tau(610. Pa) = 0.481777
- WRITEDIAGFI: date=   0.8333333    
-Ls =   0.431060 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.480317 tau(610. Pa) = 0.480317
-Ls =   0.436934 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.478857 tau(610. Pa) = 0.478857
-Ls =   0.441503 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.477530 tau(610. Pa) = 0.477530
-Ls =   0.447377 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.476070 tau(610. Pa) = 0.476070
- WRITEDIAGFI: date=   0.8750000    
-0**********    pas    840     jour*****heure 21.0    date   0.87500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947726 0.970058E+00     0.951     1.098     0.965    18.392     0.825
-
-Ls =   0.452596 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.474610 tau(610. Pa) = 0.474610
-Ls =   0.457814 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.473282 tau(610. Pa) = 0.473282
-Ls =   0.463032 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.471822 tau(610. Pa) = 0.471822
-Ls =   0.468906 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.470362 tau(610. Pa) = 0.470362
- WRITEDIAGFI: date=   0.9166667    
-Ls =   0.473469 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.469034 tau(610. Pa) = 0.469034
-Ls =   0.479343 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.467574 tau(610. Pa) = 0.467574
-Ls =   0.484561 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.466114 tau(610. Pa) = 0.466114
-Ls =   0.489772 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464787 tau(610. Pa) = 0.464787
- WRITEDIAGFI: date=   0.9583333    
-Ls =   0.494997 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463327 tau(610. Pa) = 0.463327
-Ls =   0.500871 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461867 tau(610. Pa) = 0.461867
-Ls =   0.505441 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460539 tau(610. Pa) = 0.460539
-Ls =   0.511308 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459079 tau(610. Pa) = 0.459079
- WRITEDIAGFI: date=    1.000000    
-0**********    pas    960     jour*****heure  0.0    date   1.00000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947798 0.103170E+01     0.951     1.081     0.965    18.338     0.825
-
- LTE rad. calculations up to layer           29
-Ls =   0.516526 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.458992 tau(610. Pa) = 0.458992
-Ls =   0.521744 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459039 tau(610. Pa) = 0.459039
-Ls =   0.526963 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459089 tau(610. Pa) = 0.459089
-Ls =   0.532830 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459140 tau(610. Pa) = 0.459140
- WRITEDIAGFI: date=    1.041667    
-Ls =   0.537399 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459186 tau(610. Pa) = 0.459186
-Ls =   0.543266 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459237 tau(610. Pa) = 0.459237
-Ls =   0.548485 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459288 tau(610. Pa) = 0.459288
-Ls =   0.553703 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459334 tau(610. Pa) = 0.459334
- WRITEDIAGFI: date=    1.083333    
-Ls =   0.558928 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459385 tau(610. Pa) = 0.459385
-Ls =   0.564795 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459436 tau(610. Pa) = 0.459436
-Ls =   0.569365 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459482 tau(610. Pa) = 0.459482
-Ls =   0.575232 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459533 tau(610. Pa) = 0.459533
- WRITEDIAGFI: date=    1.125000    
-0**********    pas   1080     jour*****heure  3.0    date   1.12500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947890 0.843030E+00     0.951     1.039     0.965    18.230     0.827
-
-Ls =   0.580450 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459583 tau(610. Pa) = 0.459583
-Ls =   0.585661 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459629 tau(610. Pa) = 0.459629
-Ls =   0.590880 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459680 tau(610. Pa) = 0.459680
-Ls =   0.596747 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459731 tau(610. Pa) = 0.459731
- WRITEDIAGFI: date=    1.166667    
-Ls =   0.601316 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459777 tau(610. Pa) = 0.459777
-Ls =   0.607190 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459828 tau(610. Pa) = 0.459828
-Ls =   0.612402 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459879 tau(610. Pa) = 0.459879
-Ls =   0.617620 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459925 tau(610. Pa) = 0.459925
- WRITEDIAGFI: date=    1.208333    
-Ls =   0.622838 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.459976 tau(610. Pa) = 0.459976
-Ls =   0.628698 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460026 tau(610. Pa) = 0.460026
-Ls =   0.633268 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460073 tau(610. Pa) = 0.460073
-Ls =   0.639142 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460123 tau(610. Pa) = 0.460123
- WRITEDIAGFI: date=    1.250000    
-0**********    pas   1200     jour*****heure  6.0    date   1.25000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.947987 0.935672E+00     0.951     1.057     0.965    18.134     0.828
-
-Ls =   0.644353 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460174 tau(610. Pa) = 0.460174
-Ls =   0.649571 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460220 tau(610. Pa) = 0.460220
-Ls =   0.654790 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460271 tau(610. Pa) = 0.460271
-Ls =   0.660650 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460322 tau(610. Pa) = 0.460322
- WRITEDIAGFI: date=    1.291667    
-Ls =   0.665219 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460368 tau(610. Pa) = 0.460368
-Ls =   0.671087 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460419 tau(610. Pa) = 0.460419
-Ls =   0.676305 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460470 tau(610. Pa) = 0.460470
-Ls =   0.681516 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460516 tau(610. Pa) = 0.460516
- WRITEDIAGFI: date=    1.333333    
-Ls =   0.686741 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460567 tau(610. Pa) = 0.460567
-Ls =   0.692602 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460617 tau(610. Pa) = 0.460617
-Ls =   0.697164 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460663 tau(610. Pa) = 0.460664
-Ls =   0.703031 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460714 tau(610. Pa) = 0.460714
- WRITEDIAGFI: date=    1.375000    
-0**********    pas   1320     jour*****heure  9.0    date   1.37500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948083 0.947274E+00     0.952     1.060     0.966    18.040     0.827
-
-Ls =   0.708243 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460765 tau(610. Pa) = 0.460765
-Ls =   0.713461 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460811 tau(610. Pa) = 0.460811
-Ls =   0.718679 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460862 tau(610. Pa) = 0.460862
-Ls =   0.724540 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460913 tau(610. Pa) = 0.460913
- WRITEDIAGFI: date=    1.416667    
-Ls =   0.729116 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.460959 tau(610. Pa) = 0.460959
-Ls =   0.734976 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461010 tau(610. Pa) = 0.461010
-Ls =   0.740188 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461060 tau(610. Pa) = 0.461060
-Ls =   0.745399 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461107 tau(610. Pa) = 0.461107
- WRITEDIAGFI: date=    1.458333    
-Ls =   0.750617 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461157 tau(610. Pa) = 0.461157
-Ls =   0.755829 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461208 tau(610. Pa) = 0.461208
-Ls =   0.761047 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461254 tau(610. Pa) = 0.461254
-Ls =   0.766914 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461305 tau(610. Pa) = 0.461305
- WRITEDIAGFI: date=    1.500000    
-0**********    pas   1440     jour*****heure 12.0    date   1.50000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948176 0.934652E+00     0.952     1.057     0.966    18.081     0.826
-
-Ls =   0.772126 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461356 tau(610. Pa) = 0.461356
-Ls =   0.777337 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461402 tau(610. Pa) = 0.461402
-Ls =   0.782555 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461453 tau(610. Pa) = 0.461453
-Ls =   0.787773 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461504 tau(610. Pa) = 0.461504
- WRITEDIAGFI: date=    1.541667    
-Ls =   0.792985 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461550 tau(610. Pa) = 0.461550
-Ls =   0.798852 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461601 tau(610. Pa) = 0.461601
-Ls =   0.804063 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461651 tau(610. Pa) = 0.461651
-Ls =   0.809275 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461698 tau(610. Pa) = 0.461698
- WRITEDIAGFI: date=    1.583333    
-Ls =   0.814493 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461748 tau(610. Pa) = 0.461748
-Ls =   0.819705 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461799 tau(610. Pa) = 0.461799
-Ls =   0.824916 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461845 tau(610. Pa) = 0.461845
-Ls =   0.830783 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461896 tau(610. Pa) = 0.461896
- WRITEDIAGFI: date=    1.625000    
-0**********    pas   1560     jour*****heure 15.0    date   1.62500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948276 0.830403E+00     0.952     1.070     0.965    18.139     0.827
-
-Ls =   0.835995 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461947 tau(610. Pa) = 0.461947
-Ls =   0.841206 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.461993 tau(610. Pa) = 0.461993
-Ls =   0.846424 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462044 tau(610. Pa) = 0.462044
-Ls =   0.851636 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462095 tau(610. Pa) = 0.462095
- WRITEDIAGFI: date=    1.666667    
-Ls =   0.856847 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462141 tau(610. Pa) = 0.462141
-Ls =   0.862721 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462191 tau(610. Pa) = 0.462192
-Ls =   0.867926 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462242 tau(610. Pa) = 0.462242
-Ls =   0.873137 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462288 tau(610. Pa) = 0.462288
- WRITEDIAGFI: date=    1.708333    
-Ls =   0.878349 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462339 tau(610. Pa) = 0.462339
-Ls =   0.883560 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462390 tau(610. Pa) = 0.462390
-Ls =   0.888772 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462436 tau(610. Pa) = 0.462436
-Ls =   0.893983 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462487 tau(610. Pa) = 0.462487
- WRITEDIAGFI: date=    1.750000    
-0**********    pas   1680     jour*****heure 18.0    date   1.75000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948370 0.854859E+00     0.951     1.075     0.965    18.091     0.827
-
-Ls =   0.899850 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462538 tau(610. Pa) = 0.462538
-Ls =   0.905061 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462584 tau(610. Pa) = 0.462584
-Ls =   0.910280 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462635 tau(610. Pa) = 0.462635
-Ls =   0.915484 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462685 tau(610. Pa) = 0.462685
- WRITEDIAGFI: date=    1.791667    
-Ls =   0.920696 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462732 tau(610. Pa) = 0.462731
-Ls =   0.925914 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462782 tau(610. Pa) = 0.462782
-Ls =   0.931774 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462833 tau(610. Pa) = 0.462833
-Ls =   0.936979 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462879 tau(610. Pa) = 0.462879
- WRITEDIAGFI: date=    1.833333    
-Ls =   0.942197 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462930 tau(610. Pa) = 0.462930
-Ls =   0.947416 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.462981 tau(610. Pa) = 0.462981
-Ls =   0.952620 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463027 tau(610. Pa) = 0.463027
-Ls =   0.957825 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463078 tau(610. Pa) = 0.463078
- WRITEDIAGFI: date=    1.875000    
-0**********    pas   1800     jour*****heure 21.0    date   1.87500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948442 0.914820E+00     0.951     1.076     0.965    17.987     0.826
-
-Ls =   0.963692 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463129 tau(610. Pa) = 0.463128
-Ls =   0.968910 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463175 tau(610. Pa) = 0.463175
-Ls =   0.974115 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463225 tau(610. Pa) = 0.463226
-Ls =   0.979333 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463276 tau(610. Pa) = 0.463276
- WRITEDIAGFI: date=    1.916667    
-Ls =   0.984538 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463322 tau(610. Pa) = 0.463322
-Ls =   0.989742 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463373 tau(610. Pa) = 0.463373
-Ls =   0.995609 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463424 tau(610. Pa) = 0.463424
-Ls =   1.000165 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463470 tau(610. Pa) = 0.463470
- WRITEDIAGFI: date=    1.958333    
-Ls =   1.006039 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463521 tau(610. Pa) = 0.463521
-Ls =   1.011237 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463572 tau(610. Pa) = 0.463572
-Ls =   1.016455 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463618 tau(610. Pa) = 0.463618
-Ls =   1.021667 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463669 tau(610. Pa) = 0.463668
- WRITEDIAGFI: date=    2.000000    
-0**********    pas   1920     jour*****heure  0.0    date   2.00000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948507 0.923695E+00     0.950     1.048     0.964    18.032     0.827
-
- LTE rad. calculations up to layer           29
-Ls =   1.027520 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463678 tau(610. Pa) = 0.463678
-Ls =   1.032083 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463683 tau(610. Pa) = 0.463683
-Ls =   1.037950 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463688 tau(610. Pa) = 0.463688
-Ls =   1.043154 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463694 tau(610. Pa) = 0.463694
- WRITEDIAGFI: date=    2.041667    
-Ls =   1.048359 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463699 tau(610. Pa) = 0.463699
-Ls =   1.053577 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463704 tau(610. Pa) = 0.463704
-Ls =   1.059431 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463710 tau(610. Pa) = 0.463710
-Ls =   1.063993 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463714 tau(610. Pa) = 0.463714
- WRITEDIAGFI: date=    2.083333    
-Ls =   1.069853 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463720 tau(610. Pa) = 0.463720
-Ls =   1.075065 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463725 tau(610. Pa) = 0.463725
-Ls =   1.080276 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463730 tau(610. Pa) = 0.463730
-Ls =   1.085481 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463736 tau(610. Pa) = 0.463736
- WRITEDIAGFI: date=    2.125000    
-0**********    pas   2040     jour*****heure  3.0    date   2.12500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948574 0.860545E+00     0.950     1.043     0.964    17.996     0.829
-
-Ls =   1.091334 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463741 tau(610. Pa) = 0.463741
-Ls =   1.095904 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463746 tau(610. Pa) = 0.463746
-Ls =   1.101764 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463751 tau(610. Pa) = 0.463751
-Ls =   1.106962 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463757 tau(610. Pa) = 0.463757
- WRITEDIAGFI: date=    2.166667    
-Ls =   1.112180 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463762 tau(610. Pa) = 0.463762
-Ls =   1.117392 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463767 tau(610. Pa) = 0.463767
-Ls =   1.123245 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463772 tau(610. Pa) = 0.463772
-Ls =   1.127808 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463777 tau(610. Pa) = 0.463777
- WRITEDIAGFI: date=    2.208333    
-Ls =   1.133661 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463783 tau(610. Pa) = 0.463783
-Ls =   1.138873 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463788 tau(610. Pa) = 0.463788
-Ls =   1.144077 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463793 tau(610. Pa) = 0.463793
-Ls =   1.149289 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463798 tau(610. Pa) = 0.463799
- WRITEDIAGFI: date=    2.250000    
-0**********    pas   2160     jour*****heure  6.0    date   2.25000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948658 0.957657E+00     0.950     1.091     0.964    17.950     0.831
-
-Ls =   1.155142 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463804 tau(610. Pa) = 0.463804
-Ls =   1.159705 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463809 tau(610. Pa) = 0.463809
-Ls =   1.165558 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463814 tau(610. Pa) = 0.463814
-Ls =   1.170770 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463820 tau(610. Pa) = 0.463820
- WRITEDIAGFI: date=    2.291667    
-Ls =   1.175974 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463824 tau(610. Pa) = 0.463824
-Ls =   1.181192 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463830 tau(610. Pa) = 0.463830
-Ls =   1.187039 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463835 tau(610. Pa) = 0.463835
-Ls =   1.191602 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463840 tau(610. Pa) = 0.463840
- WRITEDIAGFI: date=    2.333333    
-Ls =   1.197462 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463846 tau(610. Pa) = 0.463846
-Ls =   1.202667 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463851 tau(610. Pa) = 0.463851
-Ls =   1.207871 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463856 tau(610. Pa) = 0.463856
-Ls =   1.213083 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463861 tau(610. Pa) = 0.463861
- WRITEDIAGFI: date=    2.375000    
-0**********    pas   2280     jour*****heure  9.0    date   2.37500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948753 0.977855E+00     0.951     1.100     0.964    18.044     0.831
-
-Ls =   1.218936 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463867 tau(610. Pa) = 0.463867
-Ls =   1.223492 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463872 tau(610. Pa) = 0.463871
-Ls =   1.229352 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463877 tau(610. Pa) = 0.463877
-Ls =   1.234563 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463882 tau(610. Pa) = 0.463882
- WRITEDIAGFI: date=    2.416667    
-Ls =   1.239761 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463887 tau(610. Pa) = 0.463887
-Ls =   1.244973 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463893 tau(610. Pa) = 0.463893
-Ls =   1.250819 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463898 tau(610. Pa) = 0.463898
-Ls =   1.255382 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463903 tau(610. Pa) = 0.463903
- WRITEDIAGFI: date=    2.458333    
-Ls =   1.261235 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463908 tau(610. Pa) = 0.463908
-Ls =   1.266447 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463914 tau(610. Pa) = 0.463914
-Ls =   1.271651 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463919 tau(610. Pa) = 0.463919
-Ls =   1.276863 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463924 tau(610. Pa) = 0.463924
- WRITEDIAGFI: date=    2.500000    
-0**********    pas   2400     jour*****heure 12.0    date   2.50000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948844 0.863164E+00     0.951     1.073     0.964    18.211     0.830
-
-Ls =   1.282716 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463929 tau(610. Pa) = 0.463930
-Ls =   1.287265 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463934 tau(610. Pa) = 0.463934
-Ls =   1.293119 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463940 tau(610. Pa) = 0.463940
-Ls =   1.298337 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463945 tau(610. Pa) = 0.463945
- WRITEDIAGFI: date=    2.541667    
-Ls =   1.303535 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463950 tau(610. Pa) = 0.463950
-Ls =   1.308746 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463955 tau(610. Pa) = 0.463955
-Ls =   1.314593 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463961 tau(610. Pa) = 0.463961
-Ls =   1.319149 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463966 tau(610. Pa) = 0.463966
- WRITEDIAGFI: date=    2.583333    
-Ls =   1.325009 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463971 tau(610. Pa) = 0.463971
-Ls =   1.330213 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463977 tau(610. Pa) = 0.463977
-Ls =   1.335418 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463982 tau(610. Pa) = 0.463981
-Ls =   1.340629 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463987 tau(610. Pa) = 0.463987
- WRITEDIAGFI: date=    2.625000    
-0**********    pas   2520     jour*****heure 15.0    date   2.62500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.948931 0.874282E+00     0.951     1.079     0.964    18.370     0.830
-
-Ls =   1.346483 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463992 tau(610. Pa) = 0.463992
-Ls =   1.351032 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.463997 tau(610. Pa) = 0.463997
-Ls =   1.356899 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464003 tau(610. Pa) = 0.464003
-Ls =   1.362097 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464008 tau(610. Pa) = 0.464008
- WRITEDIAGFI: date=    2.666667    
-Ls =   1.367295 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464013 tau(610. Pa) = 0.464013
-Ls =   1.372513 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464018 tau(610. Pa) = 0.464018
-Ls =   1.378353 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464024 tau(610. Pa) = 0.464024
-Ls =   1.382908 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464029 tau(610. Pa) = 0.464028
- WRITEDIAGFI: date=    2.708333    
-Ls =   1.388769 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464034 tau(610. Pa) = 0.464034
-Ls =   1.393966 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464039 tau(610. Pa) = 0.464039
-Ls =   1.399171 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464044 tau(610. Pa) = 0.464044
-Ls =   1.404382 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464050 tau(610. Pa) = 0.464050
- WRITEDIAGFI: date=    2.750000    
-0**********    pas   2640     jour*****heure 18.0    date   2.75000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.949009 0.982761E+00     0.951     1.068     0.964    18.391     0.831
-
-Ls =   1.410229 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464055 tau(610. Pa) = 0.464055
-Ls =   1.414785 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464060 tau(610. Pa) = 0.464060
-Ls =   1.420645 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464065 tau(610. Pa) = 0.464065
-Ls =   1.425843 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464071 tau(610. Pa) = 0.464071
- WRITEDIAGFI: date=    2.791667    
-Ls =   1.431048 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464076 tau(610. Pa) = 0.464076
-Ls =   1.436252 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464081 tau(610. Pa) = 0.464081
-Ls =   1.442106 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464087 tau(610. Pa) = 0.464087
-Ls =   1.446655 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464091 tau(610. Pa) = 0.464091
- WRITEDIAGFI: date=    2.833333    
-Ls =   1.452515 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464097 tau(610. Pa) = 0.464097
-Ls =   1.457713 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464102 tau(610. Pa) = 0.464102
-Ls =   1.462917 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464107 tau(610. Pa) = 0.464107
-Ls =   1.468122 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464113 tau(610. Pa) = 0.464112
- WRITEDIAGFI: date=    2.875000    
-0**********    pas   2760     jour*****heure 21.0    date   2.87500    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.949064 0.871051E+00     0.950     1.078     0.964    18.339     0.829
-
-Ls =   1.473975 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464118 tau(610. Pa) = 0.464118
-Ls =   1.478531 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464123 tau(610. Pa) = 0.464123
-Ls =   1.484378 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464128 tau(610. Pa) = 0.464128
-Ls =   1.489575 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464134 tau(610. Pa) = 0.464134
- WRITEDIAGFI: date=    2.916667    
-Ls =   1.494794 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464139 tau(610. Pa) = 0.464139
-Ls =   1.499991 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464144 tau(610. Pa) = 0.464144
-Ls =   1.505838 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464149 tau(610. Pa) = 0.464149
-Ls =   1.510387 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464154 tau(610. Pa) = 0.464154
- WRITEDIAGFI: date=    2.958333    
-Ls =   1.516247 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464160 tau(610. Pa) = 0.464160
-Ls =   1.521445 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464165 tau(610. Pa) = 0.464165
-Ls =   1.526650 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464170 tau(610. Pa) = 0.464170
-Ls =   1.531854 check lat = -0.000007 lon =   0.000000
- tauref(610. Pa) = 0.464175 tau(610. Pa) = 0.464175
- PHYSIQ: for physdem ztime_fin =   1.000000    
- physdem1: ngridmx:         3010
- physdem1: skipping water vapour tracer
- Writing stats...
- OK variable aps                                               
- OK variable ps                                                
- OK variable tsurf                                             
- OK variable co2ice                                            
- OK variable tauref                                            
- OK variable fluxsurf_lw                                       
- OK variable fluxsurf_sw                                       
- OK variable fluxtop_lw                                        
- OK variable fluxtop_sw                                        
- OK variable temp                                              
- OK variable u                                                 
- OK variable v                                                 
- OK variable w                                                 
- OK variable rho                                               
- OK variable pressure                                          
- OK variable q2                                                
- OK variable emis                                              
- OK variable zmax_th                                           
- OK variable hfmax_th                                          
- OK variable wstar                                             
- OK variable vmr_h2ovap                                        
- OK variable vmr_h2oice                                        
- OK variable vmr_h2osat                                        
- OK variable h2o_ice_s                                         
- OK variable albedo                                            
- OK variable mtot                                              
- OK variable icetot                                            
- OK variable reffice                                           
- OK variable Nccntot                                           
- OK variable Mccntot                                           
- OK variable rice                                              
- OK variable tauTES                                            
- OK variable tau                                               
- OK variable dqsdust                                           
- OK variable dqndust                                           
- OK variable reffdust                                          
- OK variable dustq                                             
- OK variable dustN                                             
- OK variable ccnq                                              
- OK variable ccnN                                              
- WRITEDIAGFI: date=    3.000000    
-  GCM: Appel test_period avant redem ; itau=        2880
- Enregistrement pour            1  0.0000000E+00
-0**********    pas   2880     jour*****heure  0.0    date   3.00000    **********
-          masse    rmsdpdt       energie  enstrophie  entropie   rmsv    mt.ang
-        0.949109 0.103297E+01     0.950     1.096     0.964    18.331     0.829
-
- in abort_gcm
- Stopping in gcm                 
- Reason = 
- Simulation finished                                                            
-  
- Everything is cool
