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 : -26336.29 52485.20 Albedo du sol nu : 0.1003705 0.3208752 : -7079.710 14109.67 : 8.602325 5754.301 : 1.4299675E-04 6.9831863E-02 : 1.0000000E-07 1.000000 : -89.90435 89.88921 CO2 ice cover : 0.0000000E+00 958.8163 phyetat0: Le champ est present J ignore donc les autres temperatures TS** Temperature du sol 143.1560 291.0728 Surface emissivity : 0.5347768 0.9500000 Surface roughness : 7.9248937E-05 1.7229829E-02 pbl wind variance : 8.3298417E-08 157.1757 tracer on surface : 0.0000000E+00 0.0000000E+00 tracer on surface : -2.1804593E+16 3.8433894E+14 tracer on surface : -1790.340 634.2337 tracer on surface : 1.5558199E+09 2.8858525E+16 tracer on surface : 4.8576188E-05 910.5092 tracer on surface : -40.61012 13.81997 phyetat0: loading surface tracer h2o_ice instead of h2o_vap tracer on surface : -40.61012 13.81997 Soil volumetric heat capacity: 1000000. Soil thermal inertia : 40.57156 1200.000 Soil temperature : 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