| Version 58 (modified by , 3 years ago) (diff) |
|---|
Archive of POULDP reports
9th September 2021
Duration:1h
Participants: Jiandong, Christophe, Lucas, Francois, Romain, Antony, Ehouarn, Antoine, Siteng, Sandrine, Margaux,
Joseph, Aymeric
News:
====
- Jiandong Liu arrived at LMD for his 1 year long postdoc.
- Science fair on 8-10 October with Mars stuff. Volunteers needed!
Recent commits: https://trac.lmd.jussieu.fr/Planeto/log
==============
- r2558: MCS observables simulator, now computes column-integrated dod
- r2559: Cleanup of the wstats routines to generate the GCM's "stats.nc" file
(a first step towards fixing ticket #81)
Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1
==============
- None
Hot topics:
==========
* Improving the dust cycle (Antoine)
- cdod ratio between MCS column and GCM column normalized
looking at median valies (means can be quite large)
much more dust in GCM than MCS in the clear season (in daytime)
more comparable ratios during nighttime
* Dynamico (Romain)
- Nothing new (Occigen off-line this week)
Romain:
- Planetary Evolution Model: stop substep criterion on surface ice change
seems to work fine
- tests using lower resolution GCM to run faster (and to compare to full GCM-only run)
Christophe:
- 8 year simulation seems converged. To be further analyzed.
Sandrine:
- ACS Tirvim observations: thermal tides, looking at zonal
temperatures => diurnal & semi-diurnal modes,
compares well to GCM, increasing amplitude with altitude
phase difference between GCM and Tirvim (by 2-3 hours)
also compares during Global storm => larger amplitudes (up to 40k)
Joseph:
- will commit soon his last changes
- started writing his paper
2nd September 2021
Duration: 1h Participants: Lucas, Ehouarn, Margaux, Antoine, Romain, Thomas, Francois, Siteng, Joseph, Aymeric News: ==== - NASA Ames will organize a training session on its legacy GCM Recent commits: https://trac.lmd.jussieu.fr/Planeto/log ============== - None Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1 ============== - #81: some "call writediagfi" (vmr_xxx, num_xxx) end up skipped because in an IF(callstats) section in physiq. Hot topics: =========== * Improving the dust cycle (Antoine) - removed the "brighter dust" modification - test simulations with reduced mountain fraction for topflows - test simulation with lower maximum for dustradadjust (set to 3) In both cases gets warmer temperatures than before and the water cycle is better (due to changes in water albedo parameters) - thinking about issues due to the extrapolation of MCS profiles when building column values in the dusty season. One could think of adding a correction. - working on the MCs-normalization post-processing in the simulator (to match integrated dust between MCS and GCM outputs) * Dynamico (Romain) - Tracer initialization in the start files now using the start-archive files Still the model crashes after 0.5° of Ls Romain: - PEM: tests having substeps (of 1 year) in the extrapolation. with a stop criterion (surface ice amount change) when change is too large Joseph: - preparing EPSC. Water cycle tuning with frost albedo. Aymeric: - work on improvement of the online MCD interface following the brainstorm session last Monday. Lucas: - writing his internship report Margaux: - EPSC presentation about HDO modelling with comparisons with ACS. GCM too wet overall. Too dry when supersaturation is removed. Siteng - Thermal tides EPSC presentation, comparisons between TIRVIM and GCM Mode decomposition of the modes Thomas - bug in "seedout" in the MCD, at least with python interface Ehouarn - writing up a report for TAS-I on the LMD work done for ExoMars 2022 EDL
26th August 2021
Duration: 1h
Participants: Romain, Ehouarn, Joseph, Antoine, Siteng, Thomas, Francois, Lucas
News:
====
- Jiandong Liu will start his year long postdoc at LMD on Wednesday September 8th.
- We will have to start organizing MAMO (and EMM meeting in the following week)
Recent commits: https://trac.lmd.jussieu.fr/Planeto/log
==============
- r2554: Make NO and O2 Nightglow emissions computations more flexible
(indexes no longer hard coded)
- r2555: added arch files and IOIPSL install script for Irene-Rome
Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1
==============
- None.
But from Gitlab Issues #23 by Gabriella:
some "call writediagfi" (vmr_xxx, num_xxx) end up skipped because
in an IF(callstats) section in physiq.
Hot topics:
===========
* Dynamico (Romain)
- Test with low resolution work well. Now testing higher resolution (nbp=40) but crashes after a couple of days
Most likely linked to the initialization of tracers.
* Improving the dust cycle (Antoine)
- comparisons with ACS dust (and ice) profiles at low latitudes.
Looks at the impact of the injected dust radius.
ACS (quite noisy looking profiles...) sees high radius near the surface ~ 2um ,
and detached layers at ~20 km
Injection in the GCM of 3um dust works less well than using 2.5um
Will test this configuration for other Mars years / scenarios.
Siteng:
- Thermal tides study with TIRVIM. Focused on MY35 Ls~100, LT 9am/9pm
Looks at the dayly anomaly in data and GCM+averaging kernels
Decomposition in various wave modes. reasonable agreement between data and GCM
Impact of the sampling on the decomposition can be significant
Thomas:
- looking with Antoine at the MCD web interface: let's meet to discuss about some improvement
next Monday 4pm or Wednesday 4pm? (if Aymeric is also available)
- finished the MCD6.1 beta ready to be delivered to ESA.
Joseph:
- just back from vacation. Has lots of commits to make in the near future.
Lucas:
- working on EPSC presentation and has a first draft of a paper on CO2 caps evolution
Romain:
- Planetary Evolution Model: in-between GCM simulations interpolation with a
time step of 100 years. Finds local accumulation (with strongly unequilibrated
initial conditions). Tests using 3 GCM simulated years rather than 2.
Also needs to investigate how to constrain physically the extrapolation step.
Francois:
- Will be looking for 2 postdocs for "Mars through Time" in the coming months.
15th July 2021
Duration: 1h Participants: Antony, Siteng, Ehouarn, Francois, Margaux, Romain, Lucas, Thomas News: ==== - Summer break: no more POULDP Mars until 26th of August! - GENCI A11 call (deadline sept 3rd), we will request access to TGCC for the extra computation time, and test GCMs run fine on Irene-Rome in the Summer. Recent commits: https://trac.lmd.jussieu.fr/Planeto/log ============== - r2544: fix calendar when outputing multiple time steps in restart - r2545: output tab_cntrl with XIOS - r2546: adapt utilities concat/localtime/zrecast to XIOS files - r2547: update deftank reference xml files - r2551: undo unwanted changes in physiq from 2545 Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1 ============== - #79: check if run.def is there and stop if not - #80: zrecast should spot if "pressure" (or other field) is requested twice >> Add "area" to zrecast outputs (Francois) Hot topics: =========== * Time keeping in the start files Should be solved (waiting for Roland's feedback, just in case) * Dynamico (Romain) - the water cycle is now reasonable. - comparisons with MCS temperatures are promising (better in polar regions than LMDZ). - winds are quite different between LMDZ and dynamico. * Improving the Dust Cycle (Antoine) - TBC Thomas: - slope winds parametrization for MCD: comparisons with HR mesoscale simulations around Korolev: fails to well reproduce slope winds during the day. Siteng: - correcting for the seasonal change using GCM outputs doesn't work well - at 10 Pa the model is overall cooler (~5-10K) than obs. Lucas: - surface pressure comparisons with landers. MSL deficit around Ls 180-270. mesocale prediction better than GCM (local circulation effects?); but differences also at Insight landing site. Margaux: - comparing GCM and ACS data (H2O) from Anna, during the storm, more high altitude H2O in the GCm than in the observations.
8th July 2021
Duration: 1h Participants: Lucas, Ehouarn, Siteng, Christophe, Romain, Thomas, Antoine, Antony Recent commits: https://trac.lmd.jussieu.fr/Planeto/log ============== - None Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1 ============== - None => suggestion (Antoine): improve zrecast when user requests "pressure" which is also computed. => in some cases (Romain) the model crashes/stops but job continues until killed => when run.def is missing, the model should stop with a clean error message. Hot topics: =========== * Time keeping in the start files Romain made the modifications & test and will commit. * Dynamico Romain: Too dry dynamico => the problem was in the initial conditions in startfi.nc Fixes the problem for the beginning of the year. But way too much water is second half of year. Checked that water is globally conserved until month #9, where watercaptag decreases in the south * Improving the Dust Cycle Antoine: looking at ACS data around global dust storm of MY34 (paper by M. Luginin and Anna Fedorava ) and will compare mass loading and effective radius with those computed by the GCM. Plans to zrecast (high vertical resolution) or maybe use "profiles" utility Siteng: Sees large seasonal changes in the tides Sharp change at the rim of Hellas in the waves from model. Looks at substracting the seasonal effect from the tides. Christophe: Negative values of h2o_ice at he surface: re-ran 5 years from scratch, no negative values But water cycle clearly not converged yet.
1st July 2021
Duration: 43 min
Participants: Ehouarn, Siteng, Romain, Thomas, Christophe, Antoine, Margaux, Antony, Lucas, Sandrine, François
News:
=====
Roland Young is looking for a postdoc https://rolandyoung.bitbucket.io/601471.pdf
Recent commits: https://trac.lmd.jussieu.fr/Planeto/log
==============
- r2550: Make simu_MCS.F90 compatible with DYNAMICO lon-lat output files, with a check on the latitude array order when doing the interpolation
- r2551: Add possibility of additional tests (NaNs?, but also of unrealistic values) of fields at the begining of physics (i.e. coming from the dynamics) and at the
end of the physics integration
Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1
==============
- None
Wiki updates: https://trac.lmd.jussieu.fr/Planeto/wiki
============
- None
Ongoing discussions:
===================
- Romain : did some comparison between DYNAMICO-LMDZ-MCS : tsurf seems similar between DYNAMICO-LMDZ
some difference in wind between DYNAMICO and LMDZ, could it be the sponge layer?
big difference in water cycle: dryer in DYNAMICO than LMDZ, why?
- Siteng: continue to work on diurnal cycle in the observations, difference at 15 h local time between GCM and obs
- Antoine : work on dust cycle to match with MCS temperature, added the Joseph water, improvements are in the good way
- Thomas : new version of MCD 6.1.alpha: slop winds in high resolution mode.
Good agreements between mesoscale and MCD on winds
Delivery soon
- Christophe : new work on co2clouds => debugging step
17th June 2021
Duration:1h30mn
Participants: Antony, Ehouarn, Siteng, Christophe, Francois, Joseph, Adrien, Lucas, Margaux,
Romain, Thomas, Loic
News:
====
- TGO SWT#23 next week (Wednesday-Thursday), ask Francois if you want to attend.
Recent commits: https://trac.lmd.jussieu.fr/Planeto/log
==============
- None
Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1
==============
- #78: add some checking the validity of input fields coming from the dynamics
at the beginning of the physics. => Ehouarn
Hot topics:
===========
* Time keeping in the start files
- Roland reports some problems when storing multiple time steps in restart/restartfi:
I ran the GCM from sol 422 h 02 to sol 422 h 04, saving the start dump every hour.
startfi.nc has day_ini = 422, hour_ini = 0, and Time = 0.083333 (=2h)
start.nc is the same. OK I think.
After running for two hours,
restartfi.nc has day_ini = 422, hour_ini = 0.0833333, and Time = [0.04166667 0.08333333 0.125 ]
restart.nc has day_ini = 422, hour_ini = 0, and Time = [0.04166667 0.08333333 0.125 ]
- we should probably just get rid of hour_ini altogether? => Romain
* Dynamico
- Crazy behaviour if optim -O3 or -O2 ... (-O1 seems to work)
- strange unphysical patches of surface temperature in test run (seems to depend on physics call frequency)
=> check again in debug?
Christophe:
- "solved" his problem (crash at month#9) re-running the code (updated version)
=> move on to iphysiq=5 (instead of iphysiq=10)
- convergence of results over 8 years of simulations with CO2 clouds, looking at the last 3 years
Margaux:
- Simulations with the "old" (32 level + iphysiq=10) model with HDO, move on to "new" (54 levels+iphysiq=5)
or not?
Joseph:
- Water cycle tuning: influence of the contact parameter => temperature-dependent parameter
Anni recommends the tanh() formulation rather than linear dependence
But for now results are not better...
- Wonders about the possibility of having long term evolution of the aldebo of the long-term water cap
(when frost becomes perennial cap)
Antoine:
- IR-vis opacty conversion coefficient computed online => ranges from 0.1-8 ! (instead of 2.6)
Better (colder) temperatures (compared to MCS) but still not cold enough...
- looks at the brightness of dust in the visible (test with higher +1% omegaVIS wrt reference)
=> perhaps to test with constant dust ?
10th June 2021
Duration: 1h25mn
Participants: Francois, Christophe, Siteng, Thomas, Ehouarn, Margaux, Lucas, Romain, Antoine, Antony
Joseph
News:
====
- EnVision!
- MAMO dates: June 14-17 2022 (confirmed); ice-breaker on Monday 13 evening in Zamansky tower.
- Jiandong Liu will come to LMD on september 1st 2021 (1 year postdoc)
Recent commits: https://trac.lmd.jussieu.fr/Planeto/log
==============
- 2530: re-commit of r2526 for testphys1d
- 2531: fix bug in vdifc that was introduced in r2515
- 2532-33: newstart2icosa update
Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1
==============
- None
=> make a ticket about checking the validity of input fields coming from the dynamics
at the beginning of the physics.
Hot topics:
===========
* Crash with unexpected unrealistic HDO ratio: (Jean-Yves, Margaux, Ehouarn)
HDO_ice ratio of 4000 from the dynamics
Problem identified (with isotopes) to comes from the dynamics.
To be continued...
* 1=1 : (seq == parallel)
- Now OK, as of rev 2531
Christophe:
- some crashes with negative geopotential when entering physics (after 6 years of simulation).
- nighttime CO2 clouds occur ~lon=-45 lat=0, but not when removing topography and diurnal cycles
Romain:
- comparing Dynamico & LMDZ with water/dust/co2 cycles on. => strangely the initial
conditions mismatch
- dissipation fac_up =30 needed for dynamico whilst =10 works for LMDZ
Margaux:
- small study on the MY34 dust storm simulations, compared to MCS.
vertical dust distribution => does not go high enough compared
to observations (worsen by the microphysics), so tried injecting
more dust in altitude to try to match MCS, but doesn't go the right
way in 3D.
- option to try: use prescribed (Conrath) profile with microphysics
Antoine:
- conversion coeff from IRAbs to VisibleExtinction is usually taken to 2.6 inthe GCM.
investigates other values (3.4 ==> 1 um dust, or 1.8 ==> 2 um dust)
Significant effect on the temperatures ~10K; depending on the season
Is running computations where the coeff is computed online based on
GCM column opacities vs dust scenarios columns
- will do a sensitivity study on the brightness of the dust (omega)
Joseph:
- impact of temperature-dependent contact parameter: not much impact on the saturation
=> better water cycle, closer to observations, thins the polar hood.
- tested linear temperature dependence, but will check with hyperbolic formulation.
Thomas:
- compared Insight opacity and dust scenarios (updated to lastest versions);
the decrease after the stormsis steeper in the scenarios than observed.
And overall observed opacities are much higher than in the scenario most of the time.
3rd June 2021
Duration:1h20mn
Participants: Ehouarn, Joseph, Christophe, Siteng, Antoine, Romain, Lucas, Thomas, Francois, Margaux
Loic, Sandrine
News:
====
- selection of VERTIAS and Da-VINCI+ ; Venus rules!
Recent commits: https://trac.lmd.jussieu.fr/Planeto/log
==============
- None
Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1
==============
- None
Hot topics:
===========
* Crash with unexpected unrealistic HDO ratio: (Jean-Yves, Margaux, Ehouarn)
HDO_ice ratio of 4000 from the dynamics
Problem identified (with isotopes) to comes from the dynamics.
To be continued...
* 1=1 : (seq == parallel)
- doesn't work with r2519 , but was OK with r2511 (Ehouarn)
Romain:
- Dynamico: added tab-control in the XIOS outputs + adapted concatnc , zrecast & localtime
=> to be committed
Siteng:
- TIRVIM retrievals good between 5Pa and 200Pa
Antoine:
- registered to IDRIS MPI & OpenMP trainings (september - october)
- dust scaling story: "dustradadjust" & tau_GCM reaches huge values locally
dustradadjustis computed at 14LT and evolves linearly
and at some locations tau_GCM changes by orders of magnitude over a day in the
polar region where extreme values are obtained
Also find the (lesser) problem in the tropics on the volcanoes
=> importance of topflows on volcanoes
Tests with an upper limit of 5 for dustradadjust
Margaux:
- will have an intern work on HDO observable simulator of ACS or plain GCM/observations
comparison (as follow up to Loic's current work), or some "tagging" of water
Joseph:
- writing paper on present day water cycle. Temperature dependent contact parameter
improves things; will look at how this changes saturation profiles.
Thomas:
- slope winds in the MCD: comparisons with mesoscale work on the full wind
or on the anomaly?
Christophe:
- works on improvedco2cloud routine to separate meteoritic dust from the rest
Ehouarn:
- 1=1 ?
- News from Chao Wang : has implemented the dust lifting (Mulholland) scheme
27th May 2021
Duration: 1h Participants: Joseph, Francois, Siteng, Thomas, Romain, Margaux, Antoine, Lucas, Aymeric, Christophe, Antony News: ==== MAMO workshop now tentatively planned on June 14-17. Francois is waiting for confirmation from Sorbone University for the Auditorium. Recent commits: https://trac.lmd.jussieu.fr/Planeto/log ============== @2525 abierjon Mars GCM: Make a cleaner initialization of some variables in some logial variables in MCS simulator @2526 mvals Mars GCM: testphys1d.F: Fix for the 1D, follow up from previous modifications of phys_state_var_init subroutine. @2527 jnaar Bug correction in the 1d GCM, not compiling because of the of the "phys_state_var_init" routine @2528 cmathe MARS: update zrecast utility Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1 ============== - None Ongoing discussions: =================== - Antoine: Fix erroneous dust opacities appearing in the southern polar night when using "dustradadjust". These opacities follow the baroclinic waves structure. Suggestion : work at the fortran code level from a appropriate start files, output with write(*,*) when and where there is the bug, end with "STOP" and try again and again... - Thomas : Slope winds. Struggle with the coriolis forces. Need to remove the sensitivity of friction to stability (Richardson). The fit is not yet perfect but not bad. Francois suggests to use anomaly (diurnal mean ave - instantaneous) even for the computed slope winds for comparison - Siteng : analysis of the role or a priori in the TIRVIM inversion and comparison with GCM - Antony and Romain : fixed the problems with tracers in Dynamico - Mars. Soon will do realistic Mars simulations ! Work on XIOS output - Antony : test on the effect of the dissipation on warming early Mars. - Margaux Vals : new results on D/H cycle and the model. 2 papers in preparation... - Christophe : Find a sensitivity of mesospheric CO2 ice clouds to dust scenario (MY28 vs clim scenario). See CO2 clouds in the afternoon at last!
20th May 2021
Duration: 1h Participants: Francois, Siteng, Francisco, Christophe, Lucas, Thomas, Romain, Antoine, Ehouarn, Margaux Antony News: ==== - Gabriella will join IAA, 1st June (escape of water on Mars). - EPSC deadline is next week! Recent commits: https://trac.lmd.jussieu.fr/Planeto/log ============== - r2516: set reasonable minimal threshold value for water ice/vap in watercloud - r2517: update of Ciclad arch files - r2519: Follow-up of r2500, longitude boundaries were offset wrt cell center - r2522: Adding the possibility of a temperature-dependent water contact parameter - r2524: CO2 clouds; minor fix on diagnostics. Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1 ============== - None Ongoing discussions: =================== * Crash with unexpected unrealistic HDO ratio: (Jean-Yves, Margaux, Ehouarn) HDO_ice ratio of 4000 from the dynamics Problem identified (with isotopes) to comes from the dynamics. To be continued... * Finalizing work on the dust cycle parametrizations (Antoine) - tuning of topflows "19 mountains" (using MY29), still no detached layers during the clear season. - Temperatures are better (wrt GCM 5.3) during clear season but too warm during the dusty season (compared to MCS) in the lower tropo - finds weird high "dustadjust" values in the polar night - tests simulations assuming injection relative to 1um opacities => test MY26 case vs TES dust & temperature * Dynamico Mars (Romain, Antony) - run dynamico + co2 traceur; NaNs appear in prod mode (maybe from the dynamics?) for Romain failure to output tracers from within the dynamics - Antony will commit the latest startarchive2icosa tool Francois: - works on estimation of the 32um CO2 ice obs based on the radiative transfer in the GCM for Christophe Siteng: - comparisons between GCM outputs and TIRVIM, including effect of averaging kernels the smoothing due to kernel will affect the possibly to check amplitude of the tides. Thomas: - Insight solar panel data update and comparison to the 1D model of dust deposition. not so good agreement over the last 100 sols... - Draft of a paper by Ralf Lorenz on this topic is circulating - MCD slope winds: still some tuning of the parameters to do. Lucas: - calibration of Insight pressure sensor. Francisco: - runs simulations to evaluate evolution of water, will submit an abstract to EPSC
6th May 2021
Duration:1h30mn
Participants: Ehouarn, Lucas, Margaux, Christophe, Siteng, Romain, Joseph, Antoine, Thomas, Antony,
Francois
News:
====
- Luca downloaded the latest MCS products and will make an updated MY35 dust scenario.
- Ciclad has been updated with very recent versions of gfortran and ifort.
Recent commits: https://trac.lmd.jussieu.fr/Planeto/log
==============
- r2509: scripts to sync git&svn repositories
- r2510: fix example run_month1 script for Ciclad
- r2511: follow-up of r2507 which broke physics/dynamics separation
- r2512: adding "cap_albedo" flag (default .false.)
- r2513: streamfunction utility, look for cu,cv in start.nc => NB by Francois: it could be
coded/computed based on lon/lat coordinate values. Ticket?
- r2514: update of r2507 to write the correct ztime in restartfi.nc files (dynamico case)
avoid writing diagfi with dynamico + XIOS calendar
- r2515: vdifc update. (latent heat of sublimation + implicit scheme +
adaptative sub-timestep + handling hdo as well)
Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1
==============
- #76: make a utility to compute 32-um surface emissivity to compare to MCS data
Ongoing discussions:
===================
* Brutal drying up in the thermosphere issue: (Francisco, Margaux, Francois)
From Francois's checks, lowering the threshold in watercloud_mod to 1.e-16 is OK.
It should not be a hard coded number...e.g. real,parameter :: qpere_min = 1.e-16 which
is in traceur_mod.F
Margaux ran some simulations and confirms 1.e-16 works well. => Margaux will commit the update
* Crash with unexpected unrealistic HDO ratio: (Jean-Yves, Margaux, Ehouarn)
HDO_ice ratio of 4000 from the dynamics
Problem identified (with isotopes) to comes from the dynamics.
To be continued...
* Dynamico-Mars:
- XIOS output calendar: To keep it "clean" we tell it a timestep is 1s and that the
number of physics time steps is the length of the day. And impose an output time axis
labelled in days and fraction thereof.
- Loading tracers from start file now OK; requires reading a (new generation) tracer.def file
Lucas:
- MSL/REM pressure data & Insight surface pressure investigations
will compare with GCM runs provided by Francois
Francois:
- compiled outputs from various simulations for Lucas.
- updated zrecast sent to Ehouarn & Christophe
Siteng:
- ran GCM for MY35, comparing with TIRVIM temperatures, larger diurnal variations in the GCM
compared to TIRVIM. runs with iphysiq=10; will test with iphysiq=5
Antoine:
- topflow runs. Ran a second year; seems quite converged wrt first year.
- MCS-GCM ratio when comparing DSO: does the extrapolation/column recomputation of
GCM outputs to test versus MCS product to find if their is a constant ratio,
but there doesn't seem to be...
Romain:
- is putting tab_control in XIOS outputs
- will now test Dynamico with tracers (i.e CO2, water, dust cycles included)
- PEM (Planetary Evolution Model) reads/writes start files.
Joseph:
- present day water cycle tuning: work in progress.
- testing taking a water contact parameter temperature dependence into account
29th April 2021
Duration: 1h30mn Participants: Lucas, Joseph, Siteng, Ehouarn, Romain, Thomas, Francisco, Antoine, Francois, Antony Margaux, Aymeric, Christophe Recent commits: https://trac.lmd.jussieu.fr/Planeto/log ============== - r2502: set icelocationmode=4 as default mode (set watercaptags as before) in surfini - r2505: code clean-up (remove obsolete files in LMDZ.MARS/libf/dynphy_lonlat) - r2507: modifications in the calendar values stored in restart files for 1+1=2 - r2508: separating water ice deposit albedo cases between "cap" and "frost" Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1 ============== - None Wiki updates: https://trac.lmd.jussieu.fr/Planeto/wiki ============ - None Ongoing discussions: =================== * Brutal drying up in the thermosphere issue: (Francisco, Margaux, Francois) Francisco found the problem comes from a drastic 1.e-8 threshold correction in waterclouds Strong improvement in the upper atmosphere when setting this threshold to 1.e-12 Francois did some tests with 1.e-16 => does influence the polar hood, but the later is quite variable from a week to the next It also makes sens for HDO computations to lower this threshold. => All agree we should lower the threshold to 1.e-16 => Francois will check saturation at 90K to target the value * Crash with unexpected unrealistic HDO ratio: (Jean-Yves, Margaux, Ehouarn) HDO_ice ratio of 4000 from the dynamics Problem identified (with isotopes) to comes from the dynamics * Comparing GCM runs with TIRVIM (& MCS): (Antoine, Sandrine, Siteng) * Dynamico-Mars: - what time to set when calling the physics from dynamico? - what calendar setup to use in XIOS? Still needs to be checked due to restarts, etc. Christophe: - someone should make a utility to compute surface emissivity at 32um => open a ticket Joseph: - work on the present day water cycle: strong impact on the physics time step (as found by Margaux) which makes the water cycle drier. Some retuning needed frost albedo, contact parameter, variance of the distribution (won't play with that one) - new vdifc for Margaux to test. Siteng: - started with the GCM on Ciclad. Antoine: - mountain slope flows parametrization => only the biggest mountains are now used (6) also did a test run with 19 biggest mountains (converges to similar results) still not very good in terms of detached layers at the beginning of the year... - NB: for computations of DSO there is some latitude on choice of underlying factors for comparisons with MCS - run more years to check convergence. Lucas: - Insight surface pressure measurement drift. Add a correction for sensor temperature and then 1st and 2nd years match - some divergence with the model around Ls=300-360 to further investigate - also looks at MSL pressure data over 4-5 Mars years - checks the VL data and its interpolated derivation that is used to rescale the MCD Thomas: - MCD slope winds parametrization: putting it all together for all times of day and location. Still some problems in the night.
22th April 2021
Duration: 1h35mn
Participants: Francois, Lucas, Ehouarn, Antony, Antoine, Romain, Christophe, Margaux, Thomas, jean-Yves,
Joseph
Guest star: Sandrine
Recent commits: https://trac.lmd.jussieu.fr/Planeto/log
==============
- r2496-99 : Mesoscale work flow
- r2500: Mesh boundaries definitions fix
- r2501: change definition/value of ztime_fin written to restartfi.nc (when using dynamico)
Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1
==============
- None
Wiki updates: https://trac.lmd.jussieu.fr/Planeto/wiki
============
- None
Ongoing discussions:
===================
* Comparing GCM runs with TIRVIM (& MCS): (Antoine, Sandrine)
Sandrine: obs TIRVIM vs "old" GCM (32 layers) and "new" GCM (54 layers) at 30 Pa, worse agreement in the latter.
NB comparison is on observations grid (up to 1Pa); evolution of the diurnal cycle worse with "new" GCM
=> will use "old" (MCDv5.3 like) simulations for now. => rerun and output temp, wind, water, dust 24 times/day
Antoine: Comparions with MCS; "new" run seems slightly warmer (looking at 3am+3pm diurnal ave). Overall a difference
in the phases seems slightly better in the "new" run.
* Brutal drying up in the thermosphere issue: (Francisco, Margaux)
Francisco is focusing on the simulation where it happens
Margaux is checking in various past revisions if the problem is still there. So far the problem always arises (even in 4 year old models). Has a test with 32 layers showing that the problem seems correlated to using the photochemistry...
* Crash with unexpected unrealistic HDO ratio: (Jean-Yves, Margaux)
HDO_ice ratio of 4000 from the dynamics
Problem identified (with isotopes) to comes from the dynamics
* 1+1=2:
Correction to the time axis / housekeeping needed (Romain)
* Dynamico-Mars:
- what time to set when calling the physics from dynamico?
* Calendar:
- ptime sent to the physics in LMDZ == end of physics time step minus1 dynamical step...
e.g. first time the physics is called ptime=dt_phys-dt_dyn (andnot 0 or dt_phys)
why not sent the beginning or the end date of the physical time step date?
Joseph:
- committed at noon today (r2502)the icelocation mode=4 case.
- the new vdifc doesn't reproduce old results exactly. Check if this is linked to the sub-time stepping
- investing the water cycle in current Mars changing mtot and albedo, always too dry during the second half of the year
Lucas:
- pb de calibration of pressure measurements on Insight; also checking 4 MSL data over 4 Mars years.
Thomas:
slope winds in the MCD. had some problems (projection issues) with the Coriolis parameter. Fixed now.
Now uses the PBL heigh (with a min threshold of 300m) to compute reference parameter zfree0 for slope winds to
account for diurnal evolution, but results not that good. Perhaps add a multiplicative factor to the buoyancy?
15th April 2021
Duration: 1h25mn
Participants: Franck L., Christophe, Francisco, Francois, Romain, Siteng, Loic, Thomas, Ehouarn, Lucas,
Antony, Antoine, Joseph
News:
Recent commits: https://trac.lmd.jussieu.fr/Planeto/log
==============
-None
Recent tickets: https://trac.lmd.jussieu.fr/Planeto/report/1
==============
- None
Wiki updates: https://trac.lmd.jussieu.fr/Planeto/wiki
============
- Debug page initiated by Christophe: https://trac.lmd.jussieu.fr/Planeto/wiki/DebugMode
Francisco:
- investigates crash about dried-up atmosphere: Appearance of a cloud at 120km that spreads.
no huge saturation there, nor weird ccns...
=> Check if seting a minimum CCN number helps. Maybe check alsoif problem is reproduced in 1D
Antoine:
- toplows: setting up presence of "big mountains"
- looking at simulations (run_Insight 32 layers, old version of the model) and (run_TGO, 73layers, more recent)
which give different temperatures. Differences in setup (apart altitudes): non_lte param is different,
dissipation is different.
run_TGO calls physics more often, sponge layer is also different...
Largest difference between the two is the water ice column (more in run_Insight)
Tides are a bit better in the run_TGO, compared to MCS.
=> Check the vertical distribution of clouds. Look more closelyat 30 Pa, 10 Pa.
Romain:
- Dynamico runs : modifications for Dynamico (calendar) to implement and commit.
Dynamico (pseudo-random) crashes; happens for Ehouarn too (when asking for averaged outputs).
Antony makes some new initial conditions files with 32 layers + albedo, etc.
=> to be used for future tests
Antony:
- XIOS: tests outputing every 2hours, but has repetitive issues for the first step of the day.
Joseph:
- making a reference water cycle simulation (+latent heat of water) for the present day case.
Checking if the water albedo frost value can be improved.
- implementation: vdifc , icelocationmode=4 => to be committed!
8th April 2021
Durée: 30mn + ... En visio: Christophe, Francois, Lucas, Thomas, Romain, Ehouarn, Antoine, Aymeric, Joseph News: - Debut de stage de Lucas Lange avec Francois (cycle du CO2, Insight) Commissions récentes: https://trac.lmd.jussieu.fr/Planeto/log ==================== - r2494-95: improved CO2clouds: (effet radiatif des nuages de CO2) Tickets récents: https://trac.lmd.jussieu.fr/Planeto/report/1 =============== - RAS Christophe: - besoin de travailler sur co2snow. Thomas: - vents de pente pour la MCD pb: doit changer u* pour le jour ou la nuit pour mieux coller au mésoéchelle. Faire varier le zfree en fonction de l'heure (et la latitude)? - a ajouté le scénario MY35 dans le modèle 1D pour comparer le dépot de poussière sur Insight Antoine: - discute avec Luca Montabone a propos des facteurs de conversion abs IR / ext. visible - entrainement sur les flancs de montagnes, en cours. Joseph: - cycle de l'eau: réglages du icelocationmode=4 . Va committer. - effet de la résolution temporelle de la dynamique sur le cycle de l'eau augmenter day_step (à iphysiq fixé) assèche encore le cycle de l'eau. ...
1st April 2021
Durée: 1h30mn
En visio: Antoine, Christophe, Antony, Thomas, Ehouarn, Romain, Margaux, Francois, Joseph
Franck L.
News:
====
-
Commissions récentes: https://trac.lmd.jussieu.fr/Planeto/log
====================
- r2488: interface dynamico pour les coeff multiplicatif de la dissipation suivant la
verticale à la LMDZ (flag: vert_prof_dissip_type=plugin et co.)
- r2489-90: photochimie: cosmétique + coefficients plus à jour pour
ho2 +ho2 -> h2o2 + o2
- r2491: ajout d'initialisation pour le mésoéchelle
- r2492-93: ajouts workflow mésoéchelle
Tickets récents: https://trac.lmd.jussieu.fr/Planeto/report/1
===============
- RAS
=> fermer le #75?
Antoine:
- ajout d'une page Wiki pour raconter comment faire du profiling.
https://trac.lmd.jussieu.fr/Planeto/wiki/ProfilingMethod
- NB: wiki git vs wiki trac? Faire des refs de l'un vers l'autre.
Chantier sur l'assèchement du GCM:
- Margaux/Francisco regardent; Margaux remonte les versions svn pour tenter d'isoler le basculement
La transition arrive à différents moments de la simu.
Il y a clairement un gros effet de seuil (lorsque la satu > 1000?)
Faut-il ajouter un minimum de CCN (météoritiques) pour éviter d'en avoir trop peu? A tester.
- Franck regarde aussi en remontant les versions.
Chantier HDO:
- Jean-Yves a un pb avec HDO, dans advtrac de la dynamique?
Joseph:
- surfini: détermination des zones watercaptag => option icelocationmode=4 pour
indique les zones (bornes lon et lat, basés sur le 64x48 de référence)
Pour le diagnostique : donner la surface cumulée en "distance au pole (latitude)"
Va dans le bon sens pour des simus en 32x32
=> mettre icelocationmode dans callphys.def?
- haute obliquité (en 32x24), le système se stabilise en ~30 ans, avec chaleur latente
tourne aussi en 64x48 => va comparer avec TES sur des simus actuelles
- vdifc... à commettre (dernier tests hdo pour validation)
Romain:
- pb sur Occigen avec Dynamico simulations pas stable (alors que OK sur 10 sols pour Ehouarn)
- calendrier: pour 1+1=2 OK avec quelques modifs. A commettre.
25th March 2021
Durée: 1h35mn
En visio: Franck L., Francois, Francisco, Joseph, Ehouarn, Franck M., Romain, Christophe,
Jean-Yves, Margaux, Antoine, Aymeric
News:
====
- Penser à faire le ménage sur Occigen
- Trouver une nouvelle date pour le prochain MAMO, mai-juin 2022 ?
Commissions récentes: https://trac.lmd.jussieu.fr/Planeto/log
====================
- RAS
Tickets récents: https://trac.lmd.jussieu.fr/Planeto/report/1
===============
- #75: interaction entre utilitaires localtime et lslin. Si lslin repart d'une sortie localtime,
la variable "time" n'est pas adaptée, a priori. Mais au final, pas de soucis si ca
ne reste vrai qu'à lon=0.
Cycle de l'eau & chimie:
- Pb remarqué par Franciso perte de H2O vap dans la haute atmosphère entre sols 150-300,
lié à l'apparition de glace d'eau (et températures basses); à toutes les latitudes
- Margaux: simulations 32 niveaux à 120km : ne retrouve pas les mêmes structures (le
minimum en V) => lié à la thermosphère?
- Franck L: n'avait pas ce pb en mai 2019 (49 niveaux svn ???)
et déjà le pb en septembre 2019 (svn 2157),
comme pour en 73 niveaux (novembre 2019 svn 2313)
- Francois suggère de regarder la sursaturation et les quelques sols clés ou le pb s'installe.
regarder aussi les ccn
Antoine:
- comparaisons opacités MCS: réunion dédiée cet après-midi à 15h
- mountain flows: repère les "grosses montagnes" pour isoler où on appliquera la paramétrisation
- NB: le surface.nc du datadir continet des valeurs hsmons datant de Chao et pas la
version de Margaux...
Joseph:
- watercaptag: va coder une version ou les domaines des réservoirs sont donné dans le code
pour que ça marche à toutes les résolutions (et avec dynamico)
- vdifc: à venir.
Christophe:
- va commettre des corrections pour les nuages de CO2 radiativement actifs.
updatemicroco2 utilise la température en entrée pour calculer la densité de co2
- reste peut-être des pb avec "co2useh2o"
Margaux:
- possibilité d'injecter la poussière en altitude (e.g. forcée par les obs MCS)? Pour regarder
le hdo en haute altitude.
suggestion de Francois: plutôt imposer le profil de poussière (à la Conrath) pour explorer
la sensibilité.
11th March 2021
Durée: 1h30mn
En visio: Francois, Christophe, Jean-Yves, Margaux, Romain, Joseph, Ehouarn, Thomas, Antony,
Antoine, Aymeric
News:
====
- EPSC 2021 sera virtuel
- Bonne anniversaire Christophe!
Commissions récentes: https://trac.lmd.jussieu.fr/Planeto/log
====================
- RAS
Tickets récents: https://trac.lmd.jussieu.fr/Planeto/report/1
===============
- #73: 1+1=2 => binairement faut pour des raisons de calendrier (Time/controle())
Joseph:
- emballement du cycle de l'eau à basse résolution même en basse résolution
avec Mars actuel. Teste changement du paramètre de contact (m) pour voir
si ça peut se corriger et en réduisant le pas de temps physique.
=> très grosse influence tu pas de temps physiue à basse résolution!
Penser aussi au sous-pas de temps de la microphysique
- vdifc... à commettre (suites aux modifs/validation de Margaux)
Antoine:
- données MCS et comment bien les représenter dans le GCM.
Les fichiers binnés de Lucas donnent bien la même colonne que les scénarios
Arrive à un désaccord entre colonne recalculé dans le GCM via le mmr/reff du GCM.
=> regarde l'évolution des rayons effectifs dans le GCM
- chantier "topdust" qui génère de la poussière en début d'année vers -30° de lat.
tests ponctuels sans dynamique => l'origine de cette bande semble être Hellas.
Antony:
- Dynamico Mars: ça tourne avec qqs corrections pour tourner en prod.
doit faire une commission pour le calendrier
- orographie sous-maille: besoin d'un fichier d'entrée dédié à la "bonne" résolution
(celle de la simulation) qui seront ré-interpolé par XIOS.
Jean-Yves:
- pb à 73 niveaux: plantages à cause du hdo, alors que avec 49 niveaux
(ancienne résolution) pas de soucis... reste à voir d'où vient le pb.
Christophe:
- zrecast : il y a bien extrapolation dans le cas zs... (mais pas dans le cas za)
=> à corriger dans zrecast, par exemple sur la version à commettre de Francois.
- 1+1=2: a des pb sur les coordonnées (lon, lat, aire) mais aussi avec les nuages radiativement
actifs.
Aymeric:
- se fait une base d'états initiaux pour MY34. Attention, si il n'y a pas de traceur.def lors
d'un start2archive on perd les traceurs.... => ticket pour qu'on teste l’existence du
traceur.def lorsque start2archive s'exécute.
4th March 2021
Durée:1h35
En visio: Antony, Christophe, Francisco, Thomas, Margaux, Joseph, Ehouarn, Romain,
Jean-Yves, Francois
News:
====
- Arrivée de Romain (CDD IR sur l'ERC Mars Through Time).
- EGU 2021: le programme est disponible.
Commissions récentes:
====================
- r2467: bug fix sur diagnostiques d'échappements de H, H2 et D
Tickets récents: https://trac.lmd.jussieu.fr/Planeto
===============
- RAS
Christophe:
- bug avec nuages CO2 rad. actifs. => plantages mais pas facile à reproduire car
n'a pas 1+1=2 (même sans les nuages de co2!)
- zrecast: en coordonnées pression est limité par la aps du GCM. En altitude au dessus de
la surface (e.g. 120km) des valeurs sont encore présentes... pb dans zrecast
ou sur le pipeline de calcul?
Romain:
- Obtiens bien 1+1=2 dans le cas test du bench (nuages d'h2o+poussière+traceur co2)
ou presque (? quelques différences très localisées et minimes...)
Francois:
- nouveau zrecast (qui peut travailler sans avoir la température en entrée +
interpolation des espèces en log + ...)
- Rob Lillis trouve des différences en estimant la température dans la thermosphère
pour différentes espèces par rapport à la température GCM...
Francois ne retrouve pas de telles différences + penser au fait que la dynamique agit aussi
Francisco confirme qu'on récupère le bon profil de Température via les espèces
Jean-Yves mentionne que la grille de la diffusion peut jouer un peu.
- rappel "zlay" dans le modèle (physique) devient vite faux à haute altitude
(par contre est mieux estimé dans aeronomars)
Francisco:
- teste l'effect de la chimie ionosphérique sur l’échappement de H et H2.
par moments la différence (par rapport avec simu sans chimie iono) peut
être d'un facteur 4 (parfois bien moins, par ex en fin d'année)
Jean-Yves:
- simulations avec D et hdo. plantages dans la diffusion moléculaire lors de
la reconstruction de l'altitude lorsque hdo_vap devient non-physique en bas
du modèle (pb de conversion hdo_vap et hdo_ice ?) peut-être lié au fait de
tourner avec 73 niveaux?
Thomas:
- toujours sur les vents de pente pour la MCD.
regarde la sensibilité aux paramètres (dz/zlim, pentes, U*) et
compare à des simulations mésoéchelle.
Reste à trouver la bonne échelles des pentes à considérer suivant
l'horizontale (à mettre dans une carte prédéfinie à interroger).
Margaux:
- article hdo est accepté!
- cycle de l'eau: stabilité/conservation sur plusieurs années.
avec 18 années: variabilité de la colonne de glace d'un an à l'autre
surtout entre sols 200 et 300
- conservation de l'eau d'une année à l'autre: légère dérive de 1e-5 relatif par an
tests sans dynamique: meilleur conservation; ce qui nous rappelle que la dynamique
ne conserve pas bien les traceurs... les thermiques sont apparemment aussi source
d'un peu de non-conservation des traceurs...
25th Ferbuary 2021
Durée: 1h05
En visio:Thomas, Antony, Antoine, Margaux, Joseph
News:
====
- Nouveau filet de ping pong!
Commissions récentes:
====================
- v2466 : transmission de watercap au start_archive
Tickets récents:
===============
- RAS
Antony:
- RAS / travaille sur amélioration de la convergence du modèle 1D
Joseph: cycle eau à basse résolution : problème lié à la résolution ?
avec nuages radiativement actifs à LR: zones où les nuages créent un refroidissement localement, ce qu'on ne voit pas en 64x48
Thomas:
- Travail en cours sur les vents de pente: grande annonce la semaine prochaine
Antoine:
- Conservation top dust: bande d'opacité avant LS=120° dans GCM qu'on ne voit pas dans scénario.
- MCS: comparaison dso intégrée aux tau scénarios: sensibilité au rayon effectif ?
- Part se dorer la pilule
Margaux:
- Attends d'avoir 20 ans de simulations
18th February 2021
Durée:1h30 En visio: Antony, Antoine, Joseph, Francisco, Margaux, Ehouarn, Thomas, Francois, André News: ==== - Atterrissage de Perseverance ce soir Commissions récentes: ==================== - r2454-55: pour/dans le mésoéchelle - r2456: modifs pour imposer le type "double" pour les nuages de co2 (?) - r2457-58: scripts mésoéchelle - r2459: suite du r2456 (imposer des variables en double même si GCM compilé en simple precision) - r2460: "fix" mésoéchelle - r2461-63: Ajout de la chimie du deuterium Tickets récents: =============== - RAS Joseph: - vdifc nouveau: se concerte avec Margaux pour validation finale. - cycle de l'eau instable a basse résolution: regarde aussi sans les nuages radiativement actifs, et là les différences sont plus subtiles. Etude de l'emballement à basse résolution, en regardant ce qui se passe au pole (taux de sublimation et transport par l'écoulement moyen) Margaux: - regarder le cycle de l'eau et son évolution sur plusieurs années (32 niveaux ou 54 niveaux). Pas de convergence nette, mais pas de divergence. Variabilité interannuelle? Regarder les ondes baroclines sur les années assez différentes - remarque sur la poussière des simus MY34 (tempête globale); ne s'améliore pas avec les nuages actifs Ehouarn: - regarde l'évolution du cycle du CO2 sur plusieurs années. Sensibilité à la condition initiale sur les premières années. Antoine: - propriétés optiques des poussières (pour bien comparer à MCS) calculées pour des populations (et pas à un seul reff donné), ça change des choses... influence du nu_eff (largeur de la distribution) Thomas: - vents de pente dans la MCD: comparaisons au simus Mésoéchelle. Pb sur le choix de l'échelle de longueur de la pente de référence... Antony: - dynamico Mars: la simu plante (mais pas en débug).
11th February 2021
Durée: 40mn En visio: Ehouarn, Christophe, Thomas, Antony, Antoine, Loic, Margaux News: ==== - missions Hope (EMM) et Tianwen-1 en orbite autour de Mars - demande DARI planéto déposée (6.37Mh! Ca commence à faire un gros bébé...) Commissions récentes: ==================== - RAS Tickets récents: =============== - RAS Antony: - RAS Christophe: - plantage dans une simulation avec H2O comme ccn + nuages radiativement actifs. tester 1+1=2. - réunion nuages CO2 jeudi 4 mars 10h. Thomas: - RAS Antoine: - L'outil aeroptical utilise les mêmes propriétés optiques que le GCM, mais devrait plutôt utiliser les rayons effectifs et distribution obtenues par le GCM. Loic: - A ressoumis le papier HDO. Va reprendre les comparaisons avec ACS. Margaux: - Regarde les modifs sur vdifc de Joseph - HDO: le code de Francisco (chimie HDO) tourne avec ou sans photochimie et ionosphère. - Convergence du cycle de l'eau: à lancé des sim 32 niveaux ou 54 niveaux sur des décennies pour voir si il y a dérive. Ehouarn: - evolution du cycle du CO2 sur plusieurs années. C'était plus probablement un pb de conditions initiales. mais pense regarder plus longtemps que 5 ans. - doit voir avec Francisco et Jean-Yves si les corrections qui ont circulé sont à commettre.
4th February 2021
Durée:1h30 En visio: Francois, Margaux, Thomas, Francisco, Christophe, Antoine, Antony, Aymeric News: ==== - MY36 commence ce dimanche! Commissions récentes: ==================== - r2452: interface dynamico pour lire un z2sig.def (comme LMDZ) pour la discrétisation verticale Tickets récents: =============== - #71 : changer la dénomination "slope winds" en "mountain top dust flows"? Plutôt "mountaintop flow" Tour d'écrans: ============= Francisco: - travail sur la chimie du D. Tout est en place et Ehouarn/Margaux doivent valider et commettre. Encore un changement à faire "à la main" pour tourner avec/sans le HDO. - effet de la chimie sur le haute atmosphère (échapement de H lié à la chimie des ions d'eau): impact pas énorme mais significatif quand même. - compare T,rho avec les occulations IUVS: le GCM semble ~20K trop chaud vers 40-50km. Suggestion de Francois: Est-ce un pb du aux obs (bias chaud à haute pression vu avec SPICAM) Margaux: - Tests de la version de Francisco. Ca tourne. A lancé une simu 1 an. - revue du papier HDO avec Loic; demande de comparer aux obs. - a regardé le vdifc de Joseph; en cours de test. Thomas: - comparaison des "slope winds MCD" avec des simus mésoéchelle. Besoin d'évaluer les pentes à plus grande échelle et donc avec une méthode un peu plus lourde (3 triangles au lieu de 1) le but est de soustraire les vents de pente grande échelle résolus par le GCM - comparaison des accumulations de poussière sur les panneaux Insight et du modèle 1D. Marche bien. Aymeric: - difficultés à raccrocher le méso-échelle avec les versions récentes de la physique. Antony: - Première simu avec Dynamico Mars+ champs SSO avec start_archive. Pb avec le restartfi. Francois: - il faut préparer la BDD pour la mission des EMM. Quels fichiers et champs, à quelle fréquence de sortie, etc. Ehouarn: - regarde l'évolution pluri-annuelle du cycle du CO2 et on dirait que ça n'est pas bien convergé. Problème lié au start utilisé? A la nouvelle résolution verticale qui rend le système plus sensible? Suggestion de Francois: faire en test sans nuages radiativement actif pour voir...
28th January 2021
Durée:1h25
En visio: Margaux, Antony, Christophe, Thomas, Antoine, Joseph, Francois, Ehouarn
News:
====
- N'oubliez pas la demande DARI A10!
Commissions récentes:
====================
- r2449: Ajout de la possibilité d'utiliser le dust scenario 35 (mais attention, un dust_MY35.nc
définitif n'est pas encore disponible)
- r2450: utilitaire lslin: prise en compte des "missing_values" + nettoyage + attention au changement de format d'entrées (lslin.def) si delta Ls automatique
- r2451: util/compile : ajout en commentaire d'un exemple pour compiler sur Occigen
Tickets récents:
===============
- #68: réouverture du ticket (utilisation de lsin avant ou après zrecast) car pb avec les missing_value
dans le cas interpolé (alors que le cas average était OK).=> réparé avec la r2450 mais
pas forcément nécessaire (dixit Francois) car trompeur dans certaines conditions (cycle diurne)
Penser à supprimer la génération du fichier Grads *.ctl (pas trop utile pour la plupart et
Francois a un utilitaire externe qui fait l'équivalent)
Tour d'écrans:
=============
Francois:
- Modification de zrecast pour les cas où la température n'est pas dispo (+ traitement des espèces num_*) à tester puis commettre.
- Pour Rob Lillis: réinterpolation des sorties GCM; différences avec les obs dans les températures thermosphériques à élucider...
- Pour EMM (insertion à 16h42, le 9 février) : on doit faire une base de données GCM en ligne pour l'équipe. A suivre.
- Siteng Fan (postdoc) commence au 1er février, mais à distance.
Antoine:
- Injection de poussières rds: la sédimentation prend le pas sur la vitesse radiative vers ~80km
questions sur la vitesse de sédimentation différente entre masse et nombre...
- simulation "best guess" avec injection 8h-18h. Toujours des DSO plus élevées que MCS.
en générale trop de poussière, qui semble fortement déplacées des zones d'injection
Thomas:
- vents de pente dans le MCD: tests en Real*8 donnent les mêmes problèmes donc va imposer le profil
en log de vent pour calculer K.
- a contacté Ralph Lorenz pour les mesures de déposition de la poussière sur Insight.
- préparation du progress meeting ESA-CNES
Joseph:
- vdifc nettoyé et commenté envoyé à Margaux et Loic.
- simus du cyle de l'eau avec nuages radiativement actifs trop humide même à résolution "normale"
Antony:
- dynamico-mars: choix/utilisation des calendriers à finaliser. Simulations avec l'orographie sous-maille (ré-interpolée de start_archive) à venir. En attendant un retour de David sur l'outil SSO.
Christophe:
- simulations en cours d'analyse. Flux SW en haut de l'atm plus étendu (en latitude) avec les nuages CO2 radiativement actifs (par rapport au cas non radiativement actif). Les rayons des particules de CO2 sont de ~100 microns (sous la limite imposée à 500).
Ehouarn:
- dynamico-mars: mise en place de la lecture du z2sig.def via un système de plugin dans l'interface en cours; fera de même pour la dissip.
- simulations pluri-annelles : est-ce que le cycle de CO2 est stable sur n années de simulation avec MY34 (et clim?). Ca pourra rejoindre le stage de Luca sur le cycle du CO2 Insight.
- Francisco a envoyé une mise à jour de la chimie incluant le HDO. A valider et commettre.
- Contacté par Wuhu Feng (travaille avec John Plane) qui fait la modélisation de la chimie Mg, Fe, Na et obtient de bon résultats comparé à MAVEN.
21th January 2021
Durée: 1h30
En visio: Antony, Christophe, Jean-Yves, Francois, Thomas, Ehouarn, Aymeric, Antoine
News:
====
- demande DARI A10 planéto à faire
Commissions récentes:
====================
- r2441: Ajout dans makelmdz_fcm d'une vérification de l'existence du script "make_scatterers"
lorsqu'on utilise l'option "-s #" (rappel elle est obsolète pour Mars)
- r2442: Modif pour faire tourner "mesoini": le GCM en mode initialisation du mésoéchelle
- r2443: Utilitaire aeroptical: traitement de plus gros fichiers ("64 bit offset mode" NetCDF)
calcul des opacités séparées des poussières (background,storm,..)
calcul des opacités en absorption
- r2444-45: updates "mesoini" + mésoéchelle
- r2447-48: Updates nuages de CO2 (correction du rayon de sédimentation CO2ice) et effet
radiatif des nuages de CO2
Tickets récents:
===============
- RAS
Tour d'écrans:
=============
Christophe:
- Commissions CO2. Ca tourne 1 an avec les améliorations mais sans activer l'effet radiatif des
nuages de CO2.
Jean-Yves:
- comparaisons densité GCM vs NGIMS (2 deep dips). Simu MY32 (avec correction dist Mars-Soleil)
bon accord vers 130km, puis surestimation à plus haute alt. (N2, CO2 et Ar)
sous-estimation légère pour He.
rapport Ar/CO2 et N2/CO2 en bon accord
va poursuivre avec comparaisons avec d'autres deep-dips en MY33 et MY34, et regarder He et O
plus en détail
Antoine:
- comparaisons Nomad-MCS avec Mike: graphes des comparaisons à analyser. Regarder si il y a de
la poussière très haut (~80km). Pour Nomad l'erreur devient de plus en plus grande près de la surface
mais pas trop grande à ~20-30km (dans les couches).
- simulation de référence lancée avec injection 8h-18h pour MY33 en cours.
Antony:
- RAS
Thomas:
- les vents de pente dans la MC: schémas itératif amélioré (résolution simultané de u_l et v_l)
marche avec Coriolis et K fixé => besoin de fixer un profil de K, mais solution très sensible
à K. Mais oscillations entre 2 solutions et besoin d'ajouter un rappel vers l'une.
Peut-être se contenter de K constant?
Francois:
- finalise article Mars primitif avec Martin
- article Insight en finalisation
Aymeric:
- travaille sur la mise à jour du mésoéchelle
- a créé de "tags" sous git pour référencer des rev précises.
Ehouarn:
- RAS
14th January 2021
Durée: 1h10mn
En visio: Margaux, Antony, Thomas, Ehouarn, Francois, Antoine, Joseph, Aymeric, Christophe
News:
====
- arrivée des stagiaires de L3
- Romain Vandemeulebrouck (IR CDD "Mars Through Time") commencera au 1er mars
Commissions récentes:
====================
- r2440: Ajouts sur les utilitaires concatnc et localtime: homogénéisation
(et tests sur le "long_name" de Time pour savoir si sol/Ls)
Tickets récents:
===============
- #68: Ouvert puis fermé, utilisation enchainée des utilitaires zrecast -> lslin avec
gestion des "missing values". C'est en fait OK depuis la r2432
- #69: Pb à la compilation, mais en fait lié à "makelmdz_fcm" lorsqu'on
utilise (à tort) l'option "-s ..." alors qu'elle est obsolète pour Mars
=> corriger makelmdz_fcm
Tour d'écrans:
=============
Francois:
- version nouvelle de zrecast à transmettre (gestion des traceurs num_**) pour commission.
- à propos des Rocket Dust Storm: est ce que dans le mésoéchelle on a de la poussière qui
monte jusqu'à 80km comme dans le GCM?
- a mis a jour le solar occultation database pour TGO + documentation (en ligne sur
www-mars.lmd.jussieu.fr)
Aymeric:
- commence la mise à niveau du mésoéchelle sur les révisions récentes
- mise en place de tags sur git pour identifier les versions "utilisable" mesoscale
Antoine:
- aeroptical: calcul des opacités d'absorptions fait. A commettre.
- occultations solaires sur TGO: A priori pas de pb sur le sza.
- comparaisons Nomad-MCS: à reprendre.
Antony:
- pb avec la création des paramètres d'orographies sous-maille. Vois avec David Cugnet.
- teste en attendant en utilisant l'interpolation des paramètres mis dans le start_archive.nc
Christophe:
- le rayon de sédimentation du CO2: le pb vient avec un flux négatif dans le sédimentation à la
surface. A élucider.
- a lancé un run pour tester l'ajout de l’effet radiatif des nuages de CO2
Joseph:
- regarde l'emballement du cycle de l'eau à basse résolution. La calotte de CO2 se sublime
d'année en année avec l'accumulation de l'eau dans l'atmosphère.
- vdifc nouveau à commettre....
Margaux:
- finalise les livrables pour TAS-I pour l'EDL ExoMars
- travaille sur les retours sur le papier HDO; ajouter des comparaisons aux obs?
Thomas:
- les vents de pente dans la MCD: recherche d'une nouvelle façon d'approximer le vent
pour initialiser le processus itératif.
Ehouarn:
- a fond sur TAS-I EDL ExoMars 2022. Progress meeting Vendredi 15h.
7th January 2021
Durée: 1h35
En visio: Antony, Aymeric, Ehouarn, Antoine, Margaux, Thomas, Christophe, Francois
News:
====
- Relance des "bocal" des jeudis après-midi sur le gather.town de la formation LMDZ:
https://gather.town/app/dounBToM8YBQ0e5I/labLMDZ
- Dépot svn: petit (re-)changement: https n'est plus accessible que par login/mdp
(lmdz-users/lmdz2020 ou vos login/mdp svn)
Mais le dépôt http://svn.lmd.jussieu.fr/Planeto/trunk est à jour, dispo et ouvert à tous
(comme avant; seuls ceux qui ont un compte svn peuvent commettre)
- Les présentations de la formation LMDZ ont été enregistrées et sont en ligne:
https://lmdz.lmd.jussieu.fr/utilisateurs/formation/2020
- CDD IR "Mars Through Time": en attente de la réponse du sélectionné
Commissions récentes:
====================
- r2439: Vérification de la périodicité (longitude ±180) des champs dans le start.nc dynamique
et correction si besoin (c'est ce qui posait problème pour Vénus; pb. dans newstart Vénus
vraisemblablement)
Tickets récents:
===============
- RAS
Tour d'écrans:
=============
Christophe:
- problèmes parfois sur l'ordre des utilitaires à utiliser; ex. localtime => zrecast
=> à préciser (ticket)
- simulation en changeant le rayon de sédimentation CO2cloud (utilisation du "vrai" rayon) et CO2snow
ça tourne, reste à nettoyer et analyser l'impact du changement.
- transfert radiatif dans les nuages de CO2 à venir.
Antoine:
- travaille sur aeroptical pour calculer des opacités d'absorption (plutôt qu'extinction)
- occultations solaires de TGO => le solar zenith angle est il bien calculé? Et y a t'il
un effet en altitude (soleil rasant à la surface, mais plus élevé en altitude)
- regarde l'impact du rayon effectif de la stormdust injectée (e.g. 1 micron au lieu de 3)
Antony:
- newstart dynamico: orographie sous-maille => pb avec le prog de David à élucider
- penser à faire un vrai "newstart" pour dynamico
Margaux:
- préparation pour EDL ExoMars 2022. 3 simus mésoéchelle faites + celle de MY35 à venir
résolution 10km car avec 5-6km ça plante
- simulation LES déjà faite
- reviews pour l'article HDO reçues.
Ehouarn:
- préparation pour TAS-I , scénario "customisé" MY35 => scenario 201 dédié pour la MCD
Thomas:
- sur Vénus pour la VCD => schémas et prototypage en cours
- livrables ESA sur la déposition de poussière. Voudrait comparer à Insight sur 1 an.
- vents de pente dans la MCD. Tenter de rechercher une solution stationnaire via un
pseudo schémas temporel n'est pas idéal; reste à trouver comment améliorer la convergence.
17th December 2020
Durée: 1h43mn En visio: Antoine, Antony, Christophe, Thomas, Margaux, François, Ehouarn, Jospeh News: ==== - "Bocal" cet après-midi sur le gather.town de la formation LMDZ si vous voulez: https://gather.town/app/dounBToM8YBQ0e5I/labLMDZ ? - Formation LMDZ en cours ! - entretiens finis ! (7 candidats) pour l'IR CDD Mars Through Time (recrutement début février) - Stagiaires Emirati : 1 sur simulateur d'observable Hope, 1 sur données TIRVIM avec Sandrine (en présentiel uniquement, sinon ne se fera pas) - Livrables CNES : mise en ligne aujourd'hui - Livrables ESA : à upload sur le serveur CNRS Commissions récentes: ==================== - r2438 : correction d'une typo sur attribut "long name" dans localtime Tickets récents: =============== - RAS Tour de table: ============= Francois: - discussion avec Kevin sur la justesse du solar zenith angle à utiliser pour les solar occultations : le point tangent en haute atmosphère ne correspond pas forcément à un sza de 90° à la surface => vérifier par calcul géométrique que l'erreur est faible par rapport à la maille du GCM. Pb de Mean/True Solar LT ? - Bug Venus avec extension verticale : instabilités énormes aux 2 points de part et d'autre de la jonction ±180° de longitude. Margaux: - simulations mésoéchelle ExoMars : la résolution 5km a planté, nouveaux tests en cours - cycles de l'eau différents avec Christophe : comparaison des versions avec retraçage des modifications - comparaisons GCM-observations TGO de Juan Alde sur le fractionnement (photolyse VS condensation) Thomas : - mise à jour des livrables ESA sur la déposition de poussière à partir des commentaires de Fabrice - Vents de pente dans la MCD : obtient des résultats cohérents avec le schéma implicite (même si un peu lent). Comparaison avec simus mésoéchelles compliquée => voir avec simus mésoéchelles de Margaux? Prendre des pentes plus globales pour la validation avec les simulations (via une baisse de résolution des données MOLA, OU par un calcul de pente à 3 triangles?) Christophe : - A finalisé son dossier JSPS ! - Discussion avec François sur l'effet radiatif des nuages de CO2. Le calcul du rayon de sédimentation par une équation plus fine (comme microphysique de l'eau) qu'une valeur fixée à 5µm fait planter le modèle (Tsurf explose). La prise en compte de meilleurs rayons de sédimentation devrait limiter le grossissement. Antony : - fait tourner le modèle générique 1D - Dynamico Mars à reprendre Ehouarn : - scénarios de poussière pour Thalès - Arrivée du post-doc CSC plutôt en fin d'été finalement Joseph : - nouvelle routine d'interpolation dans newstart, par simple bilinéarité plutôt que par conservation de la masse. Une mauvaise interpolation des propriétés aux pôles peut avoir un gros effet sur le cycle de l'eau. Tester d'atténuer l'effet de l'initialisation en faisant tourner 50 ans? attention à éventuelle rétroaction du cycle de l'eau - nettoie et envoie vdifc chaleur latente à Margaux Antoine : - vérifie le reff de background dust et les DSO comparées à MCS si on change le reffstormdust - convertir DSO d'extinction GCM en absorption via Qext et omega, puis reconvertir en extinction via les hypothèses faites par MCS (reff1.06µm)
10th December 2020
Durée: 1h15mn En visio: Ehouarn, Antony, Thomas, Francois, Joseph, Franck, Antoine, Aymeric News: ==== - on tente un "bocal" cet après-midi sur le gather.town de la formation LMDZ si vous voulez: https://gather.town/app/dounBToM8YBQ0e5I/labLMDZ - entretiens en cours (cette semaine et la prochaine) pour l'IR CDD Mars Through Time Commissions récentes: ==================== - RAS Tickets récents: =============== - RAS Francois: - Stage de 6 mois de Luca (Mars-Septembre) => sur l'analyse de la pression Insight - programme Emiratis: 1 ou 2 stagiaires. Idée de stage: utiliser les images de VMC (et HRSC) de MExpress pour reproduire/illuster les simulations, et donc applicable ensuite aussi pour EMM plus tard. Franck: - enrichissement: impact sur l'ozone. Tests sur la dissipation. Plus d'enrichissement si on dissipe plus fort, mais pas de gros impact sur l'ozone - SWT Nomad: les anglais utilisent beaucoup le modèle Joseph: - RAS Antoine: - rédaction du livrable CNES - tente une simu de référence "best guess" GCMv6 Thomas: - livrable 1D pour l'ESA insolation/dépot de poussière. Fait. - les vents de pente dans la MCD. Quelques corrections pour bien limiter les vents au dessus de zlim. Encore à investiguer comment accélérer/ mieux représenter les choses et sur quelle épaisseur de couche. Antony: - va faire le point sur l'interface Dynamico; reste des trucs à commettre. - bien documenter la création des états initiaux. - encadre un stage L3 mi-janvier Aymeric: - RAS Ehouarn: - RAS - Contacts avec un (futur) post-doc indien.
3rd December 2020
Durée: 1h25mn
En visio: Thomas, Margaux, Christophe, Ehouarn, Loic, Antony, Antoine, Franck, Joseph, Aymeric
News:
====
- gather.town c'est très bien (dit Ehouarn) on pourrait reprendre des "bocal" dans ce cadre
Commissions récentes:
====================
- r2435: mise à jour de l'utilitaire aeroptical (ajout de stormdust)
- r2436: (modèle générique) mise en place de la possibilité d'un nouveau format de traceur.def
=> à propager à l'occasion à la physique martienne?
- r2437: update de la nucléaction dans la micro H2O + writediagmicrophy pour sortir les sous
pas de temps ("diagmicrofi.nc").
Tickets récents:
===============
- RAS
Ehouarn/Francois: Refaire les simus TGO pour sortir aussi la glace d'eau (oubli bête...)
Pour les sorties "writediagfi/stats" des espèces en molec./cm3 le nom est "rho_X"
Ca ne serait pas mieux de mettre "num_X", pour éviter une confusion sur le terme
rho (== densité) ? => plutôt oui dit Franck, voir avec avec Francisco..
Joseph: simulations à basse résol. Le cycle de l'eau (nuages rad. actifs) s’emballe.
peut-être lié à l'interpolation de la haute résol vers basse résol.
voir pour l'inertie thermique en sous sol et quantités en surface.
vdifc : va voir sa version avec Margaux et Loic
Franck: regarde les sorties tauTES pour comparer à TES
Antony: Encadrera un stage de L3 en mode super-TP sur Early Mars
Dynamico: création de l'orographie sous-maille peut se faire via lon-lat réinterpolé
Adapte le programme makesso à Mars.
Antoine: Livrable CNES sur les poussières.
Stormdust à 80km-60km. Toujour une couche de stormdust vers 60km
Question sur l'heure des scénarios; confirmation qu'elle ne compte pas beaucoup sur
la phase du cycle diurne produit.
Comparaison opacité tau_pref_gcm et équivalent reconstruit avec aeroptical. La DSO
modélisation reste très différentes de celle observées par MCS. set-il possible de coller
à la fois aux DSO observées et à la colonne imposée?
Loic: Comparaison HDO avec ACS en cours.
Remarque: peu de nuages d'eau aux poles à la période de l'aphelion belt lorsque sans micro.
Mais pourquoi?
Margaux: Fait les simulations pour TAS-I. LES OK mais Méso plante.
Thomas: Travail sur le 1D avec dépot de poussière. Remis la poussière aérienne à 2 um. Et 6-7um
déposé sur le panneau. A livrer à l'ESA.
Travaux sur l'insolation pour la MCD; quelques corrections supplémentaires. Ca marche
bien. Y compris avec comparaisons au LMD1D.
Vents de pente... toujours des difficultés (convergence ultra lente).
Christophe: Projet Japon en finalisation.
Simulations: la micro du CO2 impacte beaucoup l'émissivité à la surface dans
les régions polaires (le "alpha" dans co2snow est pour un rayon effectif supposé)
Aymeric: RAS
19th November 2020
Durée:1h35mn
En visio: Christophe, Franck, Ehouarn, Margaux, Antony, Thomas, Joseph, Antoine,
Aymeric, Francois
News:
====
- Mise à jour du datadir de référence (ajout des sections efficaces d'absorption
UV de HDO dans datadir/cross_sections) sur LMD Jussieu et Ciclad
Commissions récentes:
====================
- r2431: update du traceur.def.chemistry dans deftank
- r2432: mise à jour utilitaires zrecast + lslin (FF) ; mise à jour README et compile
=> penser à la doc, un jour.
- r2433: ajout de la photodissociation HDO (mis sous commentaire) + modif d'une vitesse
de réaction CO+OH (va changer les résultats). La chimie hétérogène est maintenant
désactivée par défaut (penser à l'activer pour la MCD)
- r2434: mise à jour utilitaires concat + solzenangle + localtime (FF+AB). concat peut
maintenant concaténer des fichiers stats.
Tickets récents:
===============
- #66 : supprimer les "#ifdef NC_DOUBLE " des utilitaires
Franck:
------
- Tout déjà dit (c.f. commentaires sur r2433)
Joseph:
------
- Nouvelle version de vdifc prète. Doit encore voir avec Margaux + Loic si c'est OK pour
le HDO.
- Fait des simulations 32x32 en Mars actuel. Trop d'eau dans le système, probablement lié au
watercaptag basé sur l'inertie thermique
Margaux:
-------
- simulations Mesoéchelle pour TAS-I lancées en 360x360 (~5km de résolution), scénario warm,
Ls=75.7°
- microphysique: a écrit un writediagmicrophy qui gère les sous-pas de temps. A commettre.
travaille toujours sur l'amélioration de la microphysique (avec résolution semi-analytique
pour la condensation) pas très concluant pour l'instant; la version actuelle semble déjà
pas mal. Va commetre l'amélioration sur la nucléation.
- prépare la réunion ACS de la semaine prochaine
Antony:
------
- orographie sous maille avec Dynamico. Le système de pré-traitement terrestre travaille
en lon-lat à résolution fixée. Est-ce qu'il ne suffirait pas simplement de réinterpoler
ce qu'on a dans le startfi? Adapte les scripts de D. cugnet pour surface.nc
- faire un surface.nc à plus haute résolution? Voir avec Aymeric ce qui est fait pour
le mésoéchelle.
Christophe:
----------
- projet pour le Japon en cours de rédaction
- regarde les emissivités qui devienent très basses (avec micoCO2) dans les régions polaires;
emissivités bien trop bases (~0.1). NB: dans stats et diagfi l'unité pour emissivité est
dite "W/m2"; à corriger.
- rayons de glace de CO2 trop grand, même si plafonnés à 500 microns.
Thomas:
------
- les vents de pente dans la MCD: évolution de deltaz entre "zfree" et "zlim".
teste avec des valeurs de "zfree" plus petites. Toujours sur les tests des schémas
implicite/explicite. Ecart entre sol analytique et numériques: en fait les schémas
numériques n'étaient pas convergés (évolution très lente...)
Antoine:
-------
- rocketduststorm: si on n'entraine plus la background dust dans l'updraft on conserve
bien le nombre et la masse des poussières.
- stormdust à 80km : elle est toujours là, même sans entrainement de la background dans
l'updraft
Aymeric:
-------
- RAS. doit mettre le méso à jour.
- le gitlab est a nouveau synchro avec svn.
Francois:
--------
- contrat CNES 2020, c'est parti, il y aura des livrables (calés sur ceux de l'ESA) pour le
15 décembre.
- projet de mise à disposition des données (sorties de modèle) pour EMM
- idem pour TGO.
- André va regarder le givre Omega sur les sorties GCM.
Ehouarn:
-------
- regarde avec Luca pour un scénario dédié pour TAS-I
- simulations TGO MY34 et MY35 relancées avec mise à jour de l'euv
12th November 2020
Durée:1h35mn
En visio: Antoine, Antony, Christophe, Francois, Joseph, Thomas, Francisco, Franck, Ehouarn, Margaux
News:
====
Commissions récentes:
====================
- r2429: Update de jthermcalc_e107.F pour prendre en compte la distance Mars-Soleil
Tickets Récents:
===============
- #66: Vérifier et corriger si besoin option "all" des utilitaires
Francois:
========
- corrections dans les utilitaires (pour la base de donnée TGO)
=> refaire les simulations avec la correction r2429...
- NB: il faudra faire le même type d'exercice pour EMM
=> besoin d'améliorer la doc (README + fichier "compile" + "module load" (action: Ehouarn)
concat,solzenangle, localtime: Antoine
zrecast,lslin: Ehouarn
Change made in the util files :
concat.F90
- rewrite and simplify the handling of time and offset so that any file
can be concatenated, including files from different years or stats file.
- Use modulo to shift starting sols in concat.nc to values between 0 and 669
- add a "adls" option in addition to "sol" or "ls" to add a Ls 1D variable while keeping "sol" as the Time axis
solzenangle.F90
- improve calculation for 1 sol file (stats) to use all local time data
- read the starting sol in the file instead of asking it to the user (except for stats file). Keep the possibility for user to change it.
- ask mean Ls value for stats file instead of sol number
- fix crash when using "all variable" option
zrecast.F90
- Force the vertical interpolation of any variable starting by "rho" to be in log space (as it should for density, molecular concentration in mol.cm-3 etc..)
- set missing value to 1.E20 instead of -9.99E33 for consistancy with other util program.
localtime.F90
- fix crash when using "all variable" option
lslin.F90
- fix various problems in the output.
Francisco:
=========
- la correction de jthermcalc_e107 est importante. Par ex des variation de 40% du CO2 à 150km
ou ~20K à 200km.
- penser à unifier les chauffages UV (Franck+Francisco)
Franck:
======
- implémente la photodissociation de HDO.
- comparaison O3/H2O avec SPICAM
Antoine:
=======
- tests de conservation du nombre stormdust et backgroud dust. Même problème sur les masses.
le problème vient du transport (vertical) dans la sous-colonne qui entraine aussi de la background dust
=> tenter de ne transporter que la stormdust en suposant que la background dust est bien mélangée (et
donc compensation de la subsidance par rapport à l'ascendance)
- Toujours le pb de concervation lors de la sédimentation de stormdust...
Joseph:
======
RAS
Thomas:
======
- livrables sur la déposition de la poussière. Comparaisons sans dust factor corrigé.
- vents de pente: toujours des soucis avec les schémas d'intégration
5th November 2020
Durée:1h30mn
En visio: Joseph, Christophe, Ehouarn, Antoine, Antony, Thomas, Francois
News:
====
- Rappel: La formation aura bien lieu en décembre prochain.
- La formation XIOS a été annulée.
- Deadline bourses JSPS/CNRS attention le deadline est en décembre côté francais
Commissions récentes:
====================
- r2425: Adaptation pour dustinjection/dustscaling_mode. L'axe des temps des scénarios de
poussière est mis à des valeurs entières (car l'heure exacte est définie ailleurs)
- r2426: (dynamique), corrections/typos sans conséquences.
Tickets Récents:
===============
- RAS
Antoine:
-------
- stormdust à 80km: tests de conservation du nombre, sans dynamique.
dans rocketdustorm: trouve des cas où la stormdust augmente... malgré le détrainement...
dans vdifc: RAS
dans sedimentation: pas de sédimentation à la surface de stormdust (dans le code) et
augmentation dans l'atmosphère...
avec sédimentation à la surface de stormdust toujours une augmentation dans
l'atmosphère
Francois:
--------
- prépare les fichiers pour TGO. (runs MY34 et MY35).
Ajout de modif sur zrecast pour réinterpoler (en log) les densitée exprimées en /cm3
Christophe:
----------
- bining des données de M. Giurana => sur Ciclad: /data/cmathe/PFS (températures et Ps, opacité
poussière et de glace)
- sur le dossier JSPS
Thomas:
------
- Dépôt sur les panneaux: est-il possible qu'il y ait des rayons différents à différents endroits?
Antony:
------
- Débute son PhD.
- regarde l'orographie sous-maille.
- Penser à faire un surface.nc à résolution 1/4 de degré.
- va utiliser dynamico-générique pour early Mars
29th October 2020
Durée: 1h30mn En personne: Ehouarn, Francois, Aymeric En visio: Antony, Thomas, Antoine, Margaux News: ==== - Le projet "atmosphères planétaires" a obtenu les 1.4Mh complémentaires demandées sur Occigen (allocation complémentaire A8 pour la période 01/11/2020 => 01/05/2021) - Confinement, le retour... Commissions récentes: ==================== - RAS Tickets Récents: =============== - RAS Antony: ------ - comparaison Dynamico/LMDZ ondes de marées thermiques très similaires - comparaison a même dissipation => retrouve un impact sur le polar warming avec Dynamico => possible de mieux le régler en dissipant moins. Thomas: ------ - vents de pente dans la MCD: cherche à modéliser le terme de flottabilité autrement. Fait diminuer le dz linéairement en fonction de z. Toujours des vents trop intenses et qui convergent pas bien. - corrige les livrables sur la déposition de poussière pour l'ESA Antoine: ------- - les rocket dust storm: montée de la poussière à 80km. les vitesses d’ascension sont rapidement à 10m/s les rayons de stormdust: gros effet de ségrégation par la sédimentation... - A testé une simulation sans dynamique, à analyser. - analyse de la simu dustscaling_mode=2 => températures et marées thermiques du même genre que avec dustscaling_mode=0 Margaux: ------- - L'article HDO est enfin vraiment soumis (avec DOI etc.) - commence les comparaisons avec les obs TGO - continue sur la microphysique (prépare une note technique) avec la tentative d'approche semi-analytique. Christophe: (e-mail) ---------- - j'ai récupéré les données de PFS via Marco Giuaranna, et je suis en train de travailler dessus car elles sont, malheureusement, inutilisables en l'état actuelle. Les variables du fichier .nc sont en 1-Dimensions de longueur N-measurements. Donc, je suis en train de les mettre en 4-Dimensions (Time, altitude, latitude, longitude). Bref, je travaille dessus ^^ - Et je vais vous répondre pour confirmer la réunion concernant la microphysique du co2. Francisco: (e-mail) --------- - teste la version mise à jour de jthermcalc_e107.F avec prise en compte de la distance Terre-Mars. - A des problèmes à l'exec en MPI sur les machines de l'IAA Ehouarn: ------- - axe des temps dans les fichiers scenario de Luca... - run TGO jusqu'à Ls=270° avec update de MY35 de Luca OK. Aymeric: ------- - doit re-tester le mésoéchelle (notamment avec les isotopes) avec les dernières versions
22th October 2020
== 22/10/2020 ==
Durée:1h25mn
En personne: Antoine, Thomas Joseph, Antony
En visio: Ehouarn, Christophe, Francois, Francisco
News:
====
- le postdoc Jiandong Liu compte venir au LMD à partir du 25 Mars 2021 (pour un an).
très intéressé par la haute atmosphère, les obs TGO (HOx).
- Siteng Fan arrivera le 1er février.
- Francisco a aussi une offre de postdoc (début avril) sur l'échappement.
- on lance la recherche d'un(e) IR pour janvier/février sur l'ERC Mars Through Time.
Offre à mettre en ligne ASAP.
- pb avec les diffusion de messages sur la liste planeto-commit => réglé ce matin.
Commissions récentes:
====================
- r2417: Ajout du schéma dustscaling_mode=2 => coeff dust_rad_adjust utilisé
uniquement dans le transfert radiatif
- r2419-20: Mise à jour des arch Ciclad (important pour dynamico)
- r2421: Mise à jour du fichier arch-gfortran.fcm (options -fdefault-double-8
et -finit-real=snan) qui existent depuis gfortran v4.9.3 (au moins)
- r2423: Mise à jour du fichier env (sous ARCH) pour Occigen (important pour dynamico)
- r2424: Mise à jour de l'interface Mars (gestion des isotopes)
Tickets Récents:
===============
- #64: day fraction "fract" uniquement 0 ou 1? Importance pour la chimie?
- #65: Pour l'EUV (impacte la thermosphère), la variation de la distance Mars-Soleil
n'est pas prise en compte...
Francisco:
=========
- teste correction EUV
- travaille sur la chimie du Deutérium.
- Et la CCN3 avec l'ESA?
Christophe:
==========
- continue l'analyse de ses simulations
- possibilité d'ajouter une climatologie des flux météoritiques (cartes? profils?)
en fonction de Ls?
Antoine:
=======
- dust_caling_mode=2: simu d'un an en cours. le dust_rad_adjust se comporte globalement bien
mais valeurs délirantes dans la nuit polaire. Peut-être borner ce coefficient?
regarder les cycles diurnes.
- les rds qui montent à 80km: en cours d'analyse, des gros rayons (~3um) se retrouvent en
altitude ~60km alors qu'il n'y en a pas au début de l’ascension. Bug?
- simulateur d'observables: dso ou mmr pour reconstruire la dso ?
On retrouve des valeurs différentes mais globalement pareil.
Joseph:
======
- runs en 32x32 pour aller plus vite. la chaleur latente marche bien. Regarde pour
un sous pas de temps adaptatif.
Thomas:
======
- accumulation des poussières sur les panneaux solaires. Comparaisons avec
Spirit et Opportunity avec correction des "dust cleaning events"
Pour Phoenix et Pathfinder, pas besoin de cette correction (comme pour Insight)
Pour Phoenix, l’opacité mesurée est parfois assez différente du scénario
- Vents de pente dans la MCD: regarde si on peut chercher à résoudre la
solution stationnaire différemment (procédure itérative), ne converge pas bien.
Antony:
======
- regarde si les utilitaires sont compatible avec les sorties dynamico. Globalement oui,
sauf pour l'aire des mailles... à réfléchir pour le géopotentiel et ps, est-bien conservatif
il faut garder la cohérence avec la topographie.
15th October 2020
Durée:1h30mn En personne: Francois, Thomas, Joseph, Ehouarn, Antoine, Christophe, Aymeric En visio: Loic, Margaux, Antony News: ==== Commissions récentes: ==================== - r2415: Bye bye "tauref" => Séparé en 2 variables plus explictes: "tau_pref_gcm" et "tau_pref_scenario", pour bien séparer la valeur "modèle" et la valeur "cible" - r2416: correction sur l’utilitaire solzenangle (initialisations) Tickets Récents: =============== - #61: Commentaires/documentation de solzenangle fermé. Loic: ==== - valeurs négatives de hdo dans la microphysique suite au bouclage sur les sous-pas de temps => correctif comme pour les ccn - dans watercloud, il y a un seuil à 1.E-8 qui corrige les valeurs; pourquoi pas 1.E-16? - utilise solzenangle pour faire des comparaisons avec ACS. Margaux: ======= - toujours sur la microphysique (optimisation du nombre de pas de temps) Joseph: ====== - chaleur latente: les dernières vérifications (vs. Phoenix: impact de ~4K sur Ts si on inclue les effets de chaleur latente; idem. que ce que trouve Franck Daerden ) - tester des runs avec moins de niveaux verticaux (32x32x.. au lieu du 64x48x54) Antoine: ======= - simulations slope winds +rds, on voit plus de dso dans les simulations que mcs en général, mais le diagnostique est-il bon? Remarque: les dso sont calculés après le transfert radiatif mais avant les autres processus. En utilisant dustq (mmr en fin de pas de temps) pour estimer la dso, on voit des différences (un facteur 2 de moins) - Suite des tests sur la fraction de maille des montagnes au carré => donne des couches détachées plus proches de celles observées par MCS pendant la saison claire Thomas: ====== - les vents de pente dans la MCD. A corrigé des bugs dans son code. Les deux schémas (explicite et implicite) convergent lentement mais vers la même valeur. valeur à convergence par ailleurs trop forte en amplitude (à re-tuner?) Antony: ====== - Pb d'enchainement des calendriers des simulations avec dynamico: lié au stockage de la date finale dans le startfi.nc. A suivre...
8th October 2020
== 08/10/2020 ==
Durée: 1h30mn
En personne: Thomas, Antoine, Antony, Christophe, Francois, Joseph
En visio: Franck
News:
====
Commissions récentes:
====================
- r2414: Corrections sur des commentaires dans le code
Tickets Récents:
===============
- #62: Pb IOIPSL => "Résolu" car il s'agit en fait d'un bug compilo (gfortran 4.8.5)
pas de soucis avec des versions plus récentes (e.g. gfortran 7.2) ou ifort.
- #63: Ne pas oublier le HDO dans le scavenging par le CO2.
Franck:
======
Question de John Wilson sur une intercomparaison des lignes de courant à Ls=90.
Différences entre GCM-LMD et MCD, mais très probablement car sur fenêtres temporelles différentes.
Joseph:
======
Debug ses modifs sur la chaleur latente. Verra ensuite avec Loic et Margaux pour le HDO.
Thomas:
======
- a fait une première version de documents pour l'ESA
- travail sur les vents de pente. Correction sur le schémas explicite mais il faut ~200
itérations. Mise en place d'un schémas semi-implicite; mais converge vers une solution
assez différente du cas explicite.
Antoine:
=======
Test sur 1 an avec les slope winds (et rocket dust storm) mais en changeant l'échelle de
longueur (au carré) => fraction plus petite => effet moindre.
Stormdust à 80km? Run sur 2 sols, regarde quelques profils de rdstorm et background dust.
Le rdstorm monte bien à 80km, avec des gros rayons 1-3 microns... rien ne marche... la
sédimentation donne des flux dans les deux directions...
Antony:
======
Comparaison LMDS - Dynamico en température vont dans le bon sens (par rapport aux diff LMDZ
MCS).
Pb de dates lors d'une reprise de simulations avec Dynamico.
Christophe:
==========
Nuages de CO2 dans un run avec H2O comme noyaux de nucléation
=> beaucoup plus de nuages de CO2 que sans.
Antoine Bierjon
Petite correction par rapport à ce que j'ai dit ce midi : comme Margaux me l'a très justement fait remarquer, la couleur du flux de sédimentation de la poussière stormdust que j'affiche ne donne pas la direction mais dit juste si une couche a gagné ou perdu de la stormdust par sédimentation. Les couches les plus hautes étant bleues (elles perdent de la stormdust par sédimentation) il est normal que celles en-dessous soient rouges (elles reçoivent la stormdust des couches bleues). Tout va bien donc (de ce côté-là en tout cas) !
1st October 2020
Durée: 1h15mn En personne: Ehouarn, Antoine, Antony, Thomas, Joseph, Vincent, Christophe, Aymeric En visio: Francois, Margaux News: ==== - EPSC toujours en cours... - toujours galère la mise en place de visio... - invasion de chinois au LMD en 2021! 3 postdocs (ERC, EMM, CSC)! Commissions récentes: ==================== - r2411: Coefficient de détrainement par défaut des rds mis à 0.05 - r2412: Mise à jour des script d'install de IOIPSL - r2413: Nettoyage aeropacity en préparation du "nouveau" tauscaling Tickets Récents: =============== - RAS Ehouarn: ======= Chantier tauscaling: Revoir si tauref ne devrait pas être renommé (ou scindé suivant les cas) Reste à voir si il suffit de faire un tauscaling_mode=2 pour la nouvelle version, et qu'il n'y a pas de soucis si on appelle plusieurs fois le transfert radiatif. On fait un point mardi matin (Francois, Ehouarn, Margaux, Antoine). Christophe: ========== RAS Francois: ======== RAS Margaux: ======= Le papier HDO a été soumis! Doivent s'occuper de mettre en ligne les données. Travaille sur l'optimisation de la microphysique (H2O), peut-être avec des sous-pas de temps de taille variable? Vincent: ======= Active la condensation du CO2 sur le mésoéchelle, et ça plante... (mais pas grosse surprise) Joseph: ====== Chantier chaleur latente: fait un sous-pas de temps adapté en chaque point de grille. Fait n'importe quoi... mais ça s'arrange... pb des sorties dans ces cas là. Antony: ====== Sur les sorties XIOS: la grille lon-lat de sortie n'est pas tout à fait comme celle des diagfi. A codé une alternative dans le GCM, nécessite un fichier diagfi.nc pour lire coord et aires. Thomas: ====== Fait la doc pour le modèle de déposition de poussière (pour contrat ESA) Les vents de pente: a corrigé les bugs des routines de Léo. Le cas avec friction simple en C_D donne des vents très fort (par rapport aux sorties du mésoéchelle). Avec un modèle Von Karman "K" plus évolué, l'amplitude est moindre, mais besoin d'un sous-pas de temps et convergence assez lente (centaines d'itérations nécessaires). A testé un schémas pseudo-implicite, mais qui diverge... Va tester avec K = cte (là le schémas implicite doit marcher) Antoine: ======= Très pris par l'EPSC. Regarde la couche détachée de stormdust qui apparait à 80km; d'où sort-elle? Fait plusieurs simulations avec différents paramètres d'injection pour comparer à une simulation de référence. Met au propre les comparaisons NOMAD et MCS. Aymeric: ======= RAS
24th September 2020
Durée: 1h25min
En personne: Antoine,Thomas, Antony, Ehouarn, Christophe
En visio:
News:
====
- Les mails automatiques suite à commission ne marchaient plus
(2408: commission Titan et 2409: commission Mars). C'est réparé.
- Formation XIOS à l'IDRIS les 4-6 Novembre:
https://forge.ipsl.jussieu.fr/ioserver/wiki/training_2020
Antony est inscrit.
Commissions récentes:
====================
- 2410: suite de 2409 pour le 1D et newstart
Tickets Récents:
===============
- RAS
Tour de table:
=============
Antony: Comparaisons LMDZ-Dynamico; NB: le dernier niveau est différent (dans
sa définition) sur les deux modèles.
Tests avec LMDZ et sponge "à la dynamico" => impact sur les dernières
couches, comme attendu...
Avec une même dissipation entre LMDZ et Dynamico, les simus se ressemblent plus.
Fait des TPs XIOS
Thomas: Déposition de poussière sur les panneaux, comparaison avec Spirit et Opportunity
l'ordre de grandeur est bon. Moins bon accord avec Pathfinder, ni Phoenix.
Commence à rédiger la doc.
Travaille sur les vents de pente en cours.
Antoine: Tests avec rds sur l'année (plus les slopes winds dans un cas). Ca marche bien
mieux avec une bonne condition initiale. Va mettre un détrainement par défaut à 0.05
dans conf_phys. Regarde pour mieux tuner les slopes winds durant la saison claire.
Fait une doc sur la méthodologie de comparaison Nomad/MCS
Christophe: Tente de faire des simulation, mais Occigen....
Simulations partant des start de Margaux (simu de ref sur le cycle de l'eau)
avec la microphysique du CO2, y compris avec H2O comme CCN.
Le dernier cas (avec H2O comme CCN) est beaucoup plus long (x2 en temps)...
17th September 2020
Durée: 1h15
En personne: Ehouarn, Thomas, Antony, Antoine, Joseph, Francois
En visio: Christophe, Margaux, Loic
News:
====
- l'EPSC, c'est parti, ou presque! (du 21/09 au 9/10)
Commissions récentes:
====================
- 2409: Nettoyage de code. Introduction d'un module "dust_param_mod" pour les flags
et variables propres au cycle des poussières.
Tickets récents:
===============
- RAS
Tour de table:
=============
- Christophe: compare des simulation avec/sans nuages de CO2 mais cycle de H2O.
obtient des rayons de particules trop gros (plafonés) dans les régions polaires
- Margaux: commence à discuter avec Francisco, Jean-Yves, ... pour le D dans la thermosphère
- Loic: Article. Génère des profils pour comparer avec TGO.
- Antony: XIOS: les sorties avec un pas de temps fractionnaire semble marcher.
penser à sortir les valeurs de sol et Ls quoi qu'il arrive.
températures et vents quand même assez différents entre dynamico et LMDZ.
regarder avec une même dissipation (et un même sponge)
- Joseph: La chaleur latente. Ca marche avec des pas implicites. Reste à regaredr la
sensibilité au pas de temps.
Fin du POULDP faute de temps.
10th September 2020
Durée: 1h10mn
En personne: Joseph, Antoine, Thomas, Christophe, Vincent, Francois, Antony, Aymeric
En visio: Ehouarn, Loic, Margaux
Par mail: Mike
News:
====
- Bascule des mails LMD : lmd.jussieu.fr => lmd.ipsl.fr (avec pb connexes vis-à-vis des listes de diffusion...)
Commissions récentes:
====================
- RAS
Tickets récents:
===============
- #62: Peut être des soucis avec les toutes dernières version d'IOIPSL?
Tour de table:
=============
- E-Mike: Dust profiles:
I have almost completed a database of individual MCS profiles (level 2) in order to more completely search for potential TGO-MCS overlaps. I will be adding a search tool using new NAIF capabilities to search for FOV overlaps. However, Armin has indicated that there was 2-month campaign of coordinated MCS-TGO observations, but that they are still working on a tool to best identify the overlaps. So, either way, we will soon be able to better identify the overlaps and I will provide these to Antoine when available. I am also attempting to provide a better error estimate analogous to the information content metrics available with Bayseian techniques (i.e, Rogers).
Frank Daerden(IASB) and Paul Streeter (OU) will be beginning separate analyses of NOMAD opacity profiles in the context of their model results (Frank with GEMS; Paul with his assimilation/reanalysis).
I have continued to add to the NOMAD opacity database as new data become available.
- Loic: comparaison avec les occultations pour HDO en projet. Utilisation de l'outil sza
article HDO en cours de finalisation, soumission imminente.
travaille sur l'utilitaire streamfunction
- Margaux: microphysique avec HDO; c'est fonctionnel.
va rajouter la thermosphère (discussion avec Francisco).
- Joseph: séparation dans vdifc du traitement de l'eau et les autres traceurs
- Antoine: regarde les variations longitudinales de dso GCM vs MCS, avec ou sans les slope winds.
- Vincent: continue la prise en main du modèle mésoéchelle; premières simulations en cours.
- Thomas: reprend les travaux de Léo sur l'insolation sur les pentes dans la MCD; a trouvé quelques bugs à corriger.
va se mettre aux vents de pente dans la MCD.
Modèle d'accumulation de poussière; pas de données des rovers (Spirit, Opportunity) à disposition avant un moment... mais Ralf Lorenz a publié qqchose de probablement récupérable/utilisable.
- Antony: les pas de temps avec dynamico: c'est surtout un soucis XIOS qui demande des pas de temps en secondes entières. Va voir si on peu déclarer les pas de temps en fraction de jours.
regarde comment faire des stats avec XIOS.
- Christophe: continue à dépouiller la simulation avec nuages de CO2.
- Aymeric: RAS (papier Insight)
- Francois: RAS
- Ehouarn: regarde pour mettre un "tauscaling_mode" pour guider la poussière.
Discussion dédiée avec Francois, Antoine, Margaux demain (14h) sur ce sujet.
3rd September 2020
Durée: 1h20
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En personne: Ehouarn, Christophe, Vincent, Antony, Joseph, Francois, Aymeric
En visio:personne!
News:
====
- Nouvelle tentative avec le matériel de visio de l'IPSL... pas d'écran cette fois-ci...mais pas de candidat à la telco non plus...
- Attention le dépot svn du LMD est passé en https; pour mettre en conformité un dépot existant:
svn switch --relocate http://svn.lmd.jussieu.fr/Planeto https://svn.lmd.jussieu.fr/Planeto
Pour un checkout il faut maintenant faire "svn co https://..." (au lieu du "svn co http://...")
- un contact au LMD avec un postdoc potentiel pour travailler sur Mars (modèle + données); à suivre...
Comissions récentes:
===================
- r2406: décommentage des writediagfi dans physiq pour les taux de chauffage SW et LW, vu que les sorties sont pilotées par le diagfi.def
- r2607: implémentation dans la microphysique du HDO
Tickets récents:
===============
- Aucun
Tour de table:
=============
- Joseph: l'utilitaire pour changer l'inertie thermique en 3D finalisé; à tester sur Occigen.
problème de divergence du sous-pas de temps sur la chaeur latente du dépôt de glace
- Antony: Bilan des simus dynamico vs LMDZ (simulations avec poussière fixe sans traceur, puis avec cycle du CO2, puis nuages actifs, etc.)
Besoin d'un pas de temps différent pour dynamico, étude des possibilités en cours
- Vincent: prépare son poster EPSC sur MOLA; comment distinguer entre les différents types de nuages?
s'est lancé sur le modèle méso-échelle.
- Christophe: prépare son poster pour l'EPSC
- Aymeric: RAS (revision du papier Insight en cours)
- Ehouarn: S'est lancé sur le chantier "tauscaling": plein de flags à mettre au propre
- Francois: RAS (travaille sur l'article Insight)
27th August 2020
== 27/08/2020 ==
Durée: 1h30
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En personne: Ehouarn, Joseph, Thomas, Antoine, Jean-Yves, Aymeric, Francois
En visio: Christophe, Antony
News:
====
- c'est la reprise, on essaye d'utiliser la visio en 201, mais c'est très galère...
- il y aura un Progress Meeting MCD ESA/CNES/LMD le mardi 8 septembre après-midi.
Tour de table:
=============
- Christophe: fait des figures pour sa simu d'un an avec nuages de CO2 (pour l'EPSC)
Doit s'occuper de la fusion des routines co2condens.
- Antony: sur dynamico a pu faire tourner 3 simus annuelles avec plus ou moins de paramétrisations enclenchées.
Attention aux pas de temps qui avec Dynamico sont en secondes pour faire des simus comparables.
Passer à 88776s par jour avec dynamico?
- Joseph: Sur la chaleur latente. Crée un utilitaire pour le GCM 3D pour initialiser mieux la température et l'inertie thermique en sous-sol.
- Thomas: cycle du co2 du GCM cherche une meilleur combinaison Inertie/profondeur de la glace.
compare la poussière entre scénario de Luca et mesures Insight.
Différence parfois d'un facteur 2, même en intégrant l'opacité des nuages d'eau. Faire une réunion dédié?
Regarde aussi Pathfinder.
- Jean-Yves: va continuer à comparer avec NGIMS, plus de soucis sur les simus depuis la dernière correction
- Antoine: regarde des comparaisons GCM et MCS de la DSO jour et nuit.
comparaisons MCS-NOMAD: peu de colocalisations (seulement aux hautes latitudes,
ou il y a peu de couches détachées MCS). Va étendre la comparaison à d'autres heures locales même si différentes de MCS.
Va regarder si les couches de dso à 80km dans le GCM ont une raison / un sens physique.
- fin du tour de table faute de temps
9th July 2020
Durée:1h30
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En personne: Ehouarn, Aymeric, Thomas, Léo, Joseph
En visio: Antoine, Margaux, Francois, Antony, Loic,
News:
====
- Bientôt les vacances!
- Margaux a une nouvelle coupe!
- Vacances du POULDP après celui de la semaine prochaine.
Reprise le 27 Aout.
Commissions récentes:
====================
-r2382-93: Nuages de CO2
-r2394-96: Nuages de CO2 + concatnc
-r2397: Correction NIRCO2abs
-r2398-99: Netoyage de code: call_abort_physic() au lieu de "stop"
-r2400: Non-orographic GW: ajout (comme pour venus) d'une limitation sur
lambda_max par la taille de la maille
-r2401: Amélioration zrecast (meilleurs tests des input + possibilité de
transférer les champs 3D (surfaciques)
-r2402: Mêmes améliorations sur le zrecast dans "UTIL"
Tickets récents:
===============
- 57: concatnc. Gérer les stats.nc? => Proposition: stopper proprement si on
tombe sur des stats.nc en expliquant qu'on ne les gère pas.
=> faire un ticket sur la possibilité de stocker plutot que les tendances
les valeurs (zq plutôt pq +dt*pdq)
=> faire une doc "debug" (sur le wiki?)
En cours:
========
- Incorporation de la chaleur latente du givre d'H2O: (Joseph)
problème de propagation de l'incrément dans les sous-pas de temps.
Encore des plantages à résoudre
Peut être un pb avec watersat? (garde-fou sur qsat trop petit)
- Dynamico: (Antony)
Fait une simulation complète poussière+H2O => pb avec surfini ou
on initialise watercaptag en dur => penser à utiliser un getin_p
- Cycle des poussières: (Antoine+Margaux+...)
tauscaling va revenir en version 2.0
- travaux sur la MCD6.0:
- (Thomas) Obtient des beaux résultats qui comparent bien aux obs
avec des poussières de 3um déposées sur les paneaux
- (Léo) Implémentation des vents de pente en cours
--- inchangé ci-dessous---
- Le nom du modèle: LMD-GCM? IPSL-GCM ? LMDZ ?
A revoir? A discuter certainement...
(et en lien avec les autres Vénus, Générique, etc.)
- utilitaire aeroptical (Antoine) comparaisons MCS-GCM
Tourner le bouton du coefficient d'injection permet d'avoir de très bons profils
d'opacité comparé aux obs sans se priver des rocketduststorm
- Cycle des poussières: (Antoine)
Travail sur les start pour démarrer uns simulation avec une distribution
rocket dust storm, mais correctement normalisée par les scénarios.
Lance une simulation avec injection de background dust pour tester.
Mauvaise répartition de la poussière nocturne dans le GCM. Discussion prochaine
sur le retour de l'infâme tauscaling... (semaine prochaine ?)
- Nuages CO2: (Christophe)
La plupart des bugs sont corrigés, les nuages de CO2 se comportent bien :
nuages mésosphériques équatoriaux + nuit polaire. Attention à bien considérer
l'altitude et non la pseudoaltitude (zrecast). Prochaine étape : faire des
jolies figures pour le papier ! Bravo Christophe !
- Cycle du CO2: (Thomas)
La masse de CO2 est bien conservée sur 3 ans, ça colle bien avec les obs Viking1 et 2.
Remonter l'inertie thermique Nord ? A voir
- HDO: Reste à implémenter le hdo dans la microphysique; ça avance. (Margaux/Loic)
Tests sur ccn/tauscaling; faut-il l'appliquer ou utiliser les vrais ccn?
Tombe sur des instabilités...avec une formation énorme de nuages d'eau...
(attention, Loic a des rapports h/d bizarre avec une faible dissipation)
=> penser à mettre à jour les run.def (dissip_fac_*) dans deftank
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
Voir comment on fait pour les cas "dust storm".
Penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée
écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pas le schémas
d'assimilation.
=> faire un ticket.
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder
au Bocal.
=> faire un ticket
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : (Antoine) En stand by .
2nd July 2020
Durée: 1h45
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio:
Franck, Antony, Antoine, Christophe, Margaux, Thomas, Joseph,
Francois, Ehouarn, Léo, Aymeric
News:
====
- Pas de news
Commissions récentes:
====================
-r2380-2381 Pas Mars
Tickets récents:
===============
- Ehouarn : possibilité de run en 2D
En cours:
========
- [inchangé] Le nom du modèle: LMD-GCM? IPSL-GCM ? LMDZ ?
A revoir? A discuter certainement...
(et en lien avec les autres Vénus, Générique, etc.)
- Incorporation de la chaleur latente du givre d'H2O: (Joseph)
Schéma implicite de la chaleur latente est beaucoup trop bourrin !
- plantages 73 niveaux en mode tauscaling: (Franck+Ehouarn)
Le problème semble venir de la dynamique ! Ehouarn fait des tests,
l'oxygène est négatif d'un coup à la latitude 25°, couche 24. La température devient
infinie en sortie de la physique... "Wake up call" sur le fait que d'autres traceurs
peuvent être négatifs à cause des soucis d'arrondi ! Ca demande une inspection
Compliqué de faire un "if" pour reset les rapports de mélange très faibles, il faut
trouver un moyen intelligent de gérer ça -> à suivre
- Dynamico: (Antony)
La conservation du co2 fonctionne bien, ainsi que la répartition de la glace de CO2.
Antony est très content. On repart vers la paramétrisation sous-maille des ondes oro.
Reste à comparer les vents, les nuages... Everything is possible. Stay focused
Sortir les champs zonaux, les ondes et tout ça. Ca va être bien.
- utilitaire aeroptical (Antoine) comparaisons MCS-GCM
Tourner le bouton du coefficient d'injection permet d'avoir de très bons profils
d'opacité comparé aux obs sans se priver des rocketduststorm
- Cycle des poussières: (Antoine)
Travail sur les start pour démarrer uns simulation avec une distribution
rocket dust storm, mais correctement normalisée par les scénarios.
Lance une simulation avec injection de background dust pour tester.
Mauvaise répartition de la poussière nocturne dans le GCM. Discussion prochaine
sur le retour de l'infâme tauscaling... (semaine prochaine ?)
- Nuages CO2: (Christophe)
La plupart des bugs sont corrigés, les nuages de CO2 se comportent bien :
nuages mésosphériques équatoriaux + nuit polaire. Attention à bien considérer
l'altitude et non la pseudoaltitude (zrecast). Prochaine étape : faire des
jolies figures pour le papier ! Bravo Christophe !
- Cycle du CO2: (Thomas)
La masse de CO2 est bien conservée sur 3 ans, ça colle bien avec les obs Viking1 et 2.
Remonter l'inertie thermique Nord ? A voir
- travaux sur la MCD6.0:
- (Thomas)
Surplus d'accumulation de la poussière par rapport aux obs, investigation des
paramètres possibles. Finalement le bouton à tourner semble être Reff, le rayon
effectif de la distribution des particules. Ca aligne les vitesses de sédimentation
mais ça ne corrige pas l'opacité de poussière sur les panneaux. Piste à suivre :
changer le rayon effectif des particules dans le transfert radiatif (agglomération).
- (Léo) Les vents de pente sont implémentés, les tests débutent cet aprem !
--- inchangé ci-dessous ---
- HDO: Reste à implémenter le hdo dans la microphysique; ça avance. (Margaux/Loic)
Tests sur ccn/tauscaling; faut-il l'appliquer ou utiliser les vrais ccn?
Tombe sur des instabilités...avec une formation énorme de nuages d'eau...
(attention, Loic a des rapports h/d bizarre avec une faible dissipation)
=> penser à mettre à jour les run.def (dissip_fac_*) dans deftank
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
Voir comment on fait pour les cas "dust storm".
Penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée
écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pas le schémas
d'assimilation.
=> faire un ticket.
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder
au Bocal.
=> faire un ticket
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : (Antoine) En stand by .
25th June 2020
Durée: 1h30
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio:
Franck, Antony,Antoine, Christophe, Margaux, Thomas, Joseph,
Francois, Ehouarn, Aymeric
News:
====
- EPSC: Les abstracts sont partis
Commissions récentes:
====================
-r2378: newstart (hdo option fix) + new fractionation coefficients
Tickets récents:
===============
- aucun nouveau ticket
En cours:
========
- Le nom du modèle: LMD-GCM? IPSL-GCM ? LMDZ ?
A revoir? A discuter certainement...
(et en lien avec les autres Vénus, Générique, etc.)
- plantages 73 niveaux en mode tauscaling: (Franck+Ehouarn)
Ca varie pas mal en fonction du scenario de poussière & EUV mais
besoin de dissiper plus. Plantages vers Ls=100-110 ... mais aussi
vers Ls=120 (MY34) qui ne passe toujours pas...
bizarrement c'est toujours vers lat~25 que ça part en vrille...
=> penser à regarder aussi si on abaisse l'influence de fac_mid
=> demander à Jean-Yves si la discrétisation de la diffusion peut aider...
- Dynamico: (Antony)
Arrive à tourner un an martien avec co2 avec topo, scenario clim + dust_mass
+ dust_number (nbp=20) => derive du cycle du CO2. Attention à l'"aire"
des mailles dans les sorties XIOS...
Orographie sous-maille: a repris les scripts/routines des terrestres
(de David Cugnet). Que faire pour les hmons, etc... (pour les rds et
slope winds)? Champs à évaluer sur la grille icosa (dans un newstart?)
- utilitaire aeroptical (Antoine) regarde les profils NOMAD de Mike pour comparaison
avec GCM
A des profils concordants avec MCS (année 34, mois 6-9), il y en a 11!
=> Aux hautes latitudes. Va regarder plus en détail l'impact du rayon
effectif des aérosols pour le conversions de longueur d'ondes (avec un
run GCM MY34). Doit regarder les données MCS en altitude (ce qui n'est pas
le cas des données binées de Luca)
- Cycle des poussières: (Antoine)
Travail sur les start pour démarrer uns simulation avec une distribution
rocket dust storm, mais correctement normalisée par les scénarios.
Lance une simulation avec injection de background dust pour tester.
- Nuages CO2: (Christophe)
A ajouter une condition limite dans CO2cloud pour gèrer les cas extrême (peu de
CCN, qui génère des énormes rayons...).
Va voir ce que ça donne sur un run d'un an.
- Cycle du CO2: (Thomas)
A fait une simu test avec les nouvelles valeurs des paramètres sur 3 ans
Doit vérifier si la pression totale est bonne, mais semble OK.
- travaux sur la MCD6.0:
- (Thomas)
les scenarios de Luca sont implémentés dans ESA1D/lmd1d, c'est OK.
Fait des tests d'accumulation , comparé à Insight.
Fit pas super pour l'instant. Trop d'accumulation? Effet radiatif trop fort?
- (Léo) Implémentation des vents de pente en cours....
- Incorporation de la chaleur latente du givre d'H2O: (Joseph)
Mise au propre de vdifc, traitement bien séparé des différents blocs (pour
chaque traceur)
Utilise un sous-pas de temps dans vdifc pour la chaleur latente
--- inchangé ci-dessous ---
- HDO: Reste à implémenter le hdo dans la microphysique; ça avance. (Margaux/Loic)
Tests sur ccn/tauscaling; faut-il l'appliquer ou utiliser les vrais ccn?
Tombe sur des instabilités...avec une formation énorme de nuages d'eau...
(attention, Loic a des rapports h/d bizarre avec une faible dissipation)
=> penser à mettre à jour les run.def (dissip_fac_*) dans deftank
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
Voir comment on fait pour les cas "dust storm".
Penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée
écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pas le schémas
d'assimilation.
=> faire un ticket.
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder
au Bocal.
=> faire un ticket
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : (Antoine) En stand by .
18th June 2020
Durée: 1h
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio:
Lola, Ehouarn, Antoine, Francois, Thomas, Antony, Christophe,
Joseph
News:
====
- Francois présente le projet de simus MY35 et BDD Solar Occultations
à la réunion ACS aujourd'hui => basé sur l'outil d'Antoine
Commissions récentes:
====================
-r2369-70: Légère modifications de makelmdz_fcm (pour Titan)
-r2374: mise sous flag de tests sur la conservation du co2 (adaptés uniquement
lorsque la microphysique du CO2 est activée)
-r2376: mise à jour des xml dans deftank & mise en place d'un
field_def_physics_mars.xml contenant les infos sur les champs disponibles
-r2377: nettoyage des fichiers superflux de l'utilitaire newstart_startarchive2icosa
Tickets récents:
===============
- aucun nouveau ticket
En cours:
========
- Franck: plante toujours à 73 niveaux même iphysiq=5, MY27...
Tests de Ehouarn, avec fac_mid=4, ça passe. A confirmer.
- Cycle du CO2: (Thomas) commits pour fitinertiadepth fait (sur sorties xvik)
Test avec 4 valeurs (run sur 4 ans qui copnvergent bien)
valeurs convergées un peu extrêmes (différentes des "anciennes valeurs")
Fait une simu test avec les nouvelles valeurs des paramètres
- utilitaire aeroptical (Antoine) regarde les profils NOMAD de Mike pour comparaison
avec GCM
Cherche des profils concordants avec MCS (année 34, mois 6-9)
- Cycle des poussières: (Antoine)
Travail sur les start pour démarrer uns simulation avec une distribution
rocket dust storm, mais correctement normalisée par les scénarios.
- Nuages CO2: (Christophe)
Tests en 1D: pb dans improvedCO2cloud qui ne sublime pas lorsqu'il n'y a pas de ccn
=> cas limite à bien comprendre et gérer.
- Dynamico: (Antony) Arrive à tourner un mois martien une fois la bonne dissipation
en place. Attention à la répartition des traceurs dans le start dynamico.
- travaux sur la MCD6.0:
- (Thomas) les scenarios de Luca sont implémentés dans ESA1D/lmd1d,
qui seront à pression de ref 610Pa (vérifier le reste du code)
Va ensuite faire des tests d'accumulation à comparer aux données
- (Léo) Implémentation du flux solaire incident pour des surfaces inclinées
va commencer à regarder des simulations mésoéchelles pour les vents
(simulations de Margaux)
- Incorporation de la chaleur latente du givre d'H2O: (Joseph)
Schéma de Margaux est en explicite, ce quine convient pas trop aux grands pas
de temps. Besoin plutôt de sous-pas de temps (explicite) qui demande de
netoyer/réécrire vdifc et séparer les traceurs plus explicitement (h2o, hdo).
--- inchangé ci-dessous ---
- HDO: Reste à implémenter le hdo dans la microphysique; ça avance. (Margaux/Loic)
Tests sur ccn/tauscaling; faut-il l'appliquer ou utiliser les vrais ccn?
Tombe sur des instabilités...avec une formation énorme de nuages d'eau...
(attention, Loic a des rapports h/d bizarre avec une faible dissipation)
=> penser à mettre à jour les run.def (dissip_fac_*) dans deftank
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
Voir comment on fait pour les cas "dust storm".
Penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée
écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pas le schémas
d'assimilation.
=> faire un ticket.
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder
au Bocal.
=> faire un ticket
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : (Antoine) En stand by .
11th June 2020
Durée: 1h15mn
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio: Ehouarn, Christophe, Léo, Antony, Anni, Lola, Antoine, Francois,
Thomas, Mike, Franck L., Joseph, Margaux, Aymeric
News:
====
- Retour au labo (pour certains)!
- Il faut rédiger un contrat CNES (travaux GCM/MCD - VCD?)
Commissions récentes:
====================
- r2337: Typo dans message d'erreur de dynetat0
- r2338: Utilitaire pour les sorties xvik (tuning du cycle du CO2)
- r2339: Petit nettoyage de makelmdz_fcm
- r2341-53: tentatives de commissions depuis git
- r2355: Update des lectures de profils de traceurs dans le 1D
- r2356-61: Tests commissions depuis git
- r2362: microphysique des nuages de CO2
- r2363: outil pour convertir des start_archive en start dynamico
Tickets récents:
===============
- aucun nouveau ticket
En cours:
========
- Cycle du CO2: (Thomas) commits pour fitinertiadepth fait (sur sorties xvik)
Test avec 4 valeurs (run sur 4 ans qui copnvergent bien)
valeurs convergées un peu extrêmes (différentes des "anciennes valeurs")
Va faire une simu test avec les nouvelles valeurs des paramètres
- utilitaire aeroptical (Antoine) regarde les profils NOMAD de Mike pour comparaison
avec GCM
Cherche des profils concordants avec MCS (année 34, mois 6-9)
- Cycle des poussières: (Antoine) Travaux sur le détrainement.
A fait une simulation avec beta = 1 (impact sur la sédimentation
des poussières) au lien que beta = 0.5 pour voir...
obtient parfois variations importante sur l'altitude des poussières
stormdust, et totale.
Pas terrible côté cycle annuel...
Veux tester rocket dust storm et renormalisation suivant les scénarios
- Nuages CO2: (Christophe) Va committer ses changements (nouveau co2condens_mod4micro), sous
flag "co2cloud" pour l'instant. => va nettoyer co2condens pour faire une version unique
Obtient une formation bizarre de nuages de CO2 à l'équateur. Teste ça en 1D.
- Dynamico: (Antony) Dissipation sous Dynamico, il n'y a que le cas Terrestre...
va coder vert_profile_dissip (comme sous LMDZ.COMMON) dans ICOSAGCM
- Franck: plante toujours à 73 niveaux même avec iphysiq=5, MY27...
- travaux sur la MCD6.0:
- (Thomas) l'outil pour le calcul de la déposition sur les panneaux solaires sera plutôt
inclu dans le modèle ESA1D/lmd1d, avec les scénarios de Luca inclus. Et ça marchera
sur un plan incliné.
- (Léo) Implémentation du flux solaire incident pour des surfaces inclinées
va commencer à regarder des simulations mésoéchelles pour les vents
(simulations de Margaux)
- HDO: Reste à implémenter le hdo dans la microphysique; ça avance. (Margaux/Loic)
Tests sur ccn/tauscaling; faut-il l'appliquer ou utiliser les vrais ccn?
Tombe sur des instabilités...avec une formation énorme de nuages d'eau...
(attention, Loic a des rapports h/d bizarre avec une faible dissipation)
=> penser à mettre à jour les run.def (dissip_fac_*) dans deftank
- Incorporation de la chaleur latente du givre d'H2O: (Joseph)
repartir des travaux de thèse de Margaux. Chantier en cours.
--- inchangé ci-dessous ---
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
Voir comment on fait pour les cas "dust storm".
Penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée
écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pas le schémas
d'assimilation.
=> faire un ticket.
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder
au Bocal.
=> faire un ticket
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : (Antoine) En stand by .
4th June 2020
Durée: 50mn
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio: Antony, Francois, Margaux, Loic, Franck L., Antoine, Aymeric, Thomas
News:
====
-
Commissions récentes:
====================
- r2333-34: ajout de sortie XIOS
- r2335: Mesoscale (follow-up of r2211)
Tickets récents:
===============
- aucun nouveau ticket
En cours:
========
- Dynamico: (Antony) Dissipation sous Dynamico, il n'y a que le cas Terrestre...
paramètre vert_profile_dissip ?
- Franck: toujours des plantages vers Ls=100 avec 73 niveaux. A voir avec Ehouarn.
- travaux sur la MCD6.0: (Thomas)
- l'outil pour le calcul de la déposition sur les panneaux solaires sera plutôt
inclu dans le modèle ESA1D/lmd1d
- (Léo) Implémentation du flux solaire incident pour des surfaces inclinées
va commencer à regarder des simulations mésoéchelles pour les vents
(simulations de Margaux)
- HDO: Reste à implémenter le hdo dans la microphysique; ça avance. (Margaux/Loic)
Question sur ccn/tauscaling; faut-il l'appliquer? Voir l'impact?
Attention, Loic a des rapports h/d bizarre avec une faible dissipation
Fabrication de starts propres. => modifs dans newstart à faire?
=> penser à mettre à jour les run.def (dissip_fac_*) dans deftank
-------------- [en cours mais non discuté] ------------------
- Cycle du CO2: (Thomas) Refais les simus de référence pour xvik.
Penser à faire les commits pour fitinertiadepth
- utilitaire aeroptical (Antoine) regarde les profils NOMAD de Mike
peut-être peut-on trouver des profils concordants avec MCS ?
A fait une simulation avec beta = 1 (impact sur la sédimentation
des poussières) au lien que beta = 0.5 pour voir...
obtient des différences significatives (du même ordre qu'en changeant le détrainement)
- Nuages CO2: (Christophe) Va committer ses changements (nouveau co2condens), sous
flag "co2cloud" pour l'instant.
Obtient une formation bizarre de nuages de CO2 à l'équateur.
- Ticket #51/52 pour lslin (Cristophe) Antoine doit voir pour l'ensemble des utilitaires. Tester avec l'option all.
- simulateur d'observables ACS: (Antoine)
utilitaire solzenangle (pour tout zsa) => A tester.
- simulateur d'observables MCS : (Antoine)
Ajoute la possibilité de sortir d'autres champs du GCM (pas dans MCS)
A améliorer en mettant les autres champs (aerosols, fonctions
poids,...).
- Incorporation de la chaleur latente du givre d'H2O: (Joseph)
repartir des travaux de thèse de Margaux. Chantier à reprendre.
- Détrainement: Antoine regarde les effets de détrainement
(et plages de détrainement 10h-12h ou 10h-14h). Simus les 5è et 6è mois.
Va comparer les profils d'aérosols Nomad avec ceux fourni par Mike Wolff.
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
Voir comment on fait pour les cas "dust storm".
Penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée
écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas
d'assimilation.
=> faire un ticket.
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder
au Bocal.
=> faire un ticket
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : Antoine a fait des simulations annuelles avec
différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 ==
référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer
le système.
En global on a toujours plus de poussière que la cible; taux de
soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la
microphysique
complète.
28th May 2020
Durée: 1h15mn
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio: Léo, Antoine, Ehouarn, Christophe, Thomas, Margaux,
Mike, Antony, Aymeric
News:
====
- EPSC: le call for abstracts est sorti
Commissions récentes:
====================
- r2323-24: isotopes: transport des ratios dans la physique
- r2325: upgrade programme "xvik" => penser à commiter therminertia
dans un sous-répertoire "xvik" de util
- r2327: amélioration outil aeroptical (gestion des "missing values")
- r2329: mésoéchelle (scripts)
- r2332: update/simplification isotopes
Tickets récents:
===============
- #54: Pb opening a dust scenario file with the mesoscale.
En cours:
========
- HDO: Reste à implémenter le hdo dans la microphysique. (Margaux/Loic)
Attention, Loic a des rapports h/d bizarre avec une faible dissipation (?)
=> penser à mettre à jour les run.def (dissip_fac_*) dans deftank
- Nuages CO2: (Christophe) Va committer ses changements (nouveau co2condens), sous
flag "co2cloud" pour l'instant.
Obtient une formation bizarre de nuages de CO2 à l'équateur.
- Cycle du CO2: (Thomas) Refais les simus de référence pour xvik.
Penser à faire les commits pour fitinertiadepth
- travaux sur la MCD6.0: (Thomas)
- travail sur la déposition de poussière: fait un utilitaire pour extraire la
déposition de le poussière sur un panneau solaire et l'atténuation de puissance
- (Léo) Implémentation du flux solaire incident pour des surfaces inclinées
va commencer à regarder des simulations mésoéchelles pour les vents
(simulations de Margaux)
=> prévoir un compte Ciclad pour Léo
- utilitaire aeroptical (Antoine) regarde les profils NOMAD de Mike
peut-être peut-on trouver des profils concordants avec MCS ?
A fait une simulation avec beta = 1 (impact sur la sédimentation
des poussières) au lien que beta = 0.5 pour voir...
obtient des différences significatives (du même ordre qu'en changeant le détrainement)
- Dynamico: (Antony) Fais tourner à partir d'un start issu de starticosa, mais
plate très vite (en haut dumodèle), mais probable problème sur les champs initiaux issu du startarchive2icosa
-------------- [en cours mais non discuté] ------------------
- Ticket #51/52 pour lslin (Cristophe) Antoine doit voir pour l'ensemble des utilitaires. Tester avec l'option all.
- simulateur d'observables ACS: (Antoine)
utilitaire solzenangle (pour tout zsa) => A tester.
- simulateur d'observables MCS : (Antoine)
Ajoute la possibilité de sortir d'autres champs du GCM (pas dans MCS)
A améliorer en mettant les autres champs (aerosols, fonctions
poids,...).
- Incorporation de la chaleur latente du givre d'H2O: (Joseph)
repartir des travaux de thèse de Margaux. Chantier à reprendre.
- Détrainement: Antoine regarde les effets de détrainement
(et plages de détrainement 10h-12h ou 10h-14h). Simus les 5è et 6è mois.
Va comparer les profils d'aérosols Nomad avec ceux fourni par Mike Wolff.
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
Voir comment on fait pour les cas "dust storm".
Penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée
écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas
d'assimilation.
=> faire un ticket.
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder
au Bocal.
=> faire un ticket
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : Antoine a fait des simulations annuelles avec
différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 ==
référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer
le système.
En global on a toujours plus de poussière que la cible; taux de
soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la
microphysique
complète.
14th May 2020
Durée: 1h10
Rappel: Les CR du POULDP sont disponibles sur
http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio: Antoine, Antony, Cristophe, François, Thomas, Margaux, Anni, Lola, Ehouarn, Léo, Frank
News:
====
- EPSC: 100% virtuelle, format à définir. Plus d'infos d'ici la fin du mois.
- Contact avec un étudiant américain pour couplage avec modèle de volcan (Ehouarn)
Commissions récentes:
====================
- r2313: suite concatnc
- r2316: watercap: modifications dans vdifc
- r2317: ajout de l'utilitaire aeroptical (Antoine)
- r2318: suite 2317
- r2321: suite de la chimie ions, presque finalisé
- r2322: transport D/H
Tickets récents:
===============
- #51: lslin: probleme lecture de variables selon leur dimension
- #52: généralisation #51 à tous les utilitaires
- #53: problème compil GCM avec gfortran
En cours:
========
- D/H : (Loic et Margaux)
- commits du transport D/H pour sedimentation et condensation (r2322)
- watercap corrigé
- utilitaire aeroptical (Antoine) calcul des opacités de glace d'eau et de poussière.
Pour comparaisons avec ACS, NOMAD, ...
- Simulateur d'observables (Antoine): ça avance.
- Dynamico: (Antony) Reprise des travaux avec dernière version ICOSAGCM.
La simu basique sans starts tourne bien avec la dernière version (10 jours)
Travail pour newstart_icosa en cours
- Nuages CO2: (Christophe) RAS. Travaille pour sortir des figures les plus pertinentes.
- Ticket #51/52 pour lslin (Cristophe) Antoine doit voir pour l'ensemble des utilitaires. Tester avec l'option all.
- Cycle du CO2: (Thomas) La convergence est bien atteinte avec la microphysique activée
Penser à faire les commits pour xvik, fitinertiadepth
- travaux sur la MCD6.0: (Thomas)
- implémentation de la lecture des scénarios de poussière en journalier
=> interpolation en Ls, 1 valeur par jour. OK.
Implémentation dans la MCD5.3
Prochaine étape: travail sur la déposition de poussière. Biblio.
- (Léo) Implémentation du flux solaire incident pour des surfaces inclinées
-------------- [en cours mais non discuté] ------------------
- simulateur d'observables ACS: (Antoine)
utilitaire solzenangle (pour tout zsa) => A tester.
- simulateur d'observables MCS : (Antoine)
Ajoute la possibilité de sortir d'autres champs du GCM (pas dans MCS)
A améliorer en mettant les autres champs (aerosols, fonctions
poids,...).
- Incorporation de la chaleur latente du givre d'H2O: (Joseph)
repartir des travaux de thèse de Margaux. Chantier à reprendre.
- Détrainement: Antoine regarde les effets de détrainement
(et plages de détrainement 10h-12h ou 10h-14h). Simus les 5è et
6è mois.
Va comparer les profils d'aérosols Nomad avec ceux fourni par Mike
Wolff.
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste
les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
Voir comment on fait pour les cas "dust storm".
Penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est
bien le cas)
- Gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée
écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas
d'assimilation.
=> faire un ticket.
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder
au Bocal.
=> faire un ticket
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : Antoine a fait des simulations annuelles avec
différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 ==
référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer
le système.
En global on a toujours plus de poussière que la cible; taux de
soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la
microphysique
complète.
7th May 2020
Durée:1h20mn
Rappel: Les CR du POULDP sont disponibles sur
http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio: Antony, Franck, Ehouarn, Loic, Christophe,
Léo, Lola, Francois, Francisco, Thomas, Antoine, Aymeric
News:
====
- confinement encore et toujours...
- mise en place d'un Gitlab (IN2P3) autour du modèle générique,
mais pas que... demander à Aymeric pour en être.
Commissions récentes:
====================
- r2306: Les isotopes dans la dynamique (suite)
correction pour lecture traceur.def avec schémas d'advection
- r2307: Les isotopes dans la dynamique; mise au propre des seuils
- r2308: correction dans utilitaire concatnc
+ changement de l'offset de départ
- r2311: nettoyage du code (getin=> getin_p, pas de "stop", etc.)
- r2312: hdo dans la physique
Tickets récents:
===============
- #50: besoin de mettre plus d'info sur la compil/NetCDF quelque part.
En cours:
========
- D/H : (Loic et Margaux)
- en place sans microphysique; une correction sur co2condens encore à faire.
- watercap à corriger (Loic + Francois + Joseph)
sa valeur devrait être toujours négative
=> commettre la correction (Loic)
- Dynamico: (Antony) travail sur outil d'interpolation de start_archive (LMDZ)
vers start dynamico. Attention à la topo sous-maille.
Apprend à communiquer avec XIOS
- Nuages CO2: (Christophe) Quelques figures... simu d'un an
on a bien les nuages de co2 dans les nuits polaires avec
une forte variabilité jour à jour.
Continue de nettoyer les codes.
- Cycle du CO2: (Thomas)
Dérive des simulations sur plusieurs années...
La masse totale de CO2 reste bien constante.
Le système chauffe avec les années.
Problème probablement dû au fait d'avoir mis les
nuages radiativement actif sans microphysique...
Effectivement il y a alors une très forte augmentation
de la quantité de nuages d'une année à la suivante.
- travaux sur la MCD6.0: (Thomas)
- implémentation de la lecture des scénarios de poussière en journalier
=> interpolation en Ls, 1 valeur par jour. OK.
- (Léo) implémentation des vents de pente.
installation/utilisation de la MCD avec python3
-------------- [en cours mais non discuté] ------------------
- simulateur d'observables ACS: (Antoine)
utilitaire solzenangle (pour tout zsa) => A tester.
- simulateur d'observables MCS : (Antoine)
Ajoute la possibilité de sortir d'autres champs du GCM (pas dans MCS)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
- Incorporation de la chaleur latente du givre d'H2O: (Joseph)
repartir des travaux de thèse de Margaux. Chantier à reprendre.
- Détrainement: Antoine regarde les effets de détrainement
(et plages de détrainement 10h-12h ou 10h-14h). Simus les 5è et 6è mois.
Va comparer les profils d'aérosols Nomad avec ceux fourni par Mike Wolff.
Pour les comparaison avec ACS.
Fait un utilitaire pour calculer les opacités
(poussière, glace d'eau, combiné)
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
Voir comment on fait pour les cas "dust storm".
Penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas d'assimilation.
=> faire un ticket.
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder au Bocal.
=> faire un ticket
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
30th April 2020
Durée: 1h30mn
Rappel: Les CR du POULDP sont disponibles sur
http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio: Ehouarn, Antoine, Margaux, Thomas, Aymeric, Christophe,
Antony, Francois, Loic, Francisco, Anni
News:
====
- Fin (au 30/04) de l'allocation A6 sur Occigen, on a dépensé 4.3Mh
(sur 3.8Mh allouées). L'allocation A8 commence le 01/05, on a 3.9Mh.
- EGU (virtuelle) la semaine prochaine
- demande de stage Ingé. ULM (contact F. Forget) 3-4 semaines en Juillet
Commissions récentes:
====================
- r2296: Les isotopes dans la dynamique
- r2298: isotopes dans la dynamique (suite)
- r2302: Ajout de la chimie des ions H2O+
- r2303: Amélioration du concatnc
- r2304: Nettoyage du code
Tickets récents:
===============
- Rien de nouveau...
=> mettre à jour le "compile" des utilitaires
=> Page HowTo / FAQ dédié à NetCDF librairies & compilations ?
A mettre aussi dans le README des utilitaires.
En cours:
========
- Incorporation de la chaleur latente du givre d'H2O: (Joseph) [inchangé]
repartir des travaux de thèse de Margaux. Chantier à reprendre.
- travaux sur la MCD6.0: (Thomas)
- implémentation de la lecture des scénarios de poussière en journalier
(au lieu des moyennes) => interpolation en Ls, 1 valeur par jour.
- Nuages CO2: (Christophe) Ca marche. Regarde pour l'update du code.
(co2condens mis à jour). Regarde les diagnostiques de ses simus.
- Cycle du CO2: (Thomas)
Dérive des simulations sur plusieurs années...
La masse totale de CO2 reste bien constante.
Problème probablement dû au fait d'avoir mis les
nuages radiativement actif sans microphysique...
- D/H : (Loic et Margaux)
- faire la sédimentation conjointe H2O/HDO, c'est mieux.
reste quelquefois des valeurs trop grandes (>100)...
dans le condensation du co2?
- watercap à corriger (Loic + Francois + Joseph)
-------------- [inchangé] ci dessous ------------------
- simulateur d'observables ACS: (Antoine)
utilitaire solzenangle (pour tout zsa) => A tester.
- solliciter Luca pour produire un scénario MY35 (Ls=0-180) (Francois; fait)
- simulateur d'observables MCS : [inchangé] (Antoine)
Ajoute la possibilité de sortir d'autres champs du GCM (pas dans MCS)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
- Dynamico: (Antony) travail sur outil d'interpolation de start_archive (LMDZ)
vers start dynamico. Attention à la topo sous-maille.
- Détrainement: [inchangé] Antoine regarde les effets de détrainement
(et plages de détrainement 10h-12h ou 10h-14h). Simus les 5è et 6è mois.
Va comparer les profils d'aérosols Nomad avec ceux fourni par Mike Wolff.
Pour les comparaison avec ACS.
Fait un utilitaire pour calculer les opacités
(poussière, glace d'eau, combiné)
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
[inchangé] Voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- [inchangé] gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas d'assimilation.
=> faire un ticket.
- pb en local au LMD: [inchangé] on arrive pas à tourner en mode prod, à regarder au Bocal.
=> faire un ticket
- newstart: [inchangé] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
23rd April 2020
Durée:1h15mn
Rappel: Les CR du POULDP sont disponibles sur
http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio: Antony, Francois, Lola, Ehouarn, Antoine, Franck, Thomas,
Christophe, Margaux, Joseph, Aymeric
News:
====
- On a obtenu une rallonge d'espace ~25% sur le scratchdir sur Occigen
(+10 To, +200000 inodes)
Commissions récentes:
====================
- r2289: solzenangle utility (replaces terminator)
- r2290: add possibility to "getin_p" arrays of strings
- r2291-92: Titan/mesoscale
- r2293: Correction pour le 1D, suite à l'ajout de watercap.
- r2294: modif pour cas avec démarrage sans startfi.nc
Tickets récents:
===============
- #45: améliorer zrecast (ajouter les champs surfacique en sortie)
- #46: améliorer concatnc pour concaténer des concat.nc
- #47: modèle générique
- #48: cohérence entres attributs "long_name" et "title" pour les utilitaires
En cours:
========
- travaux sur la MCD6.0: (Thomas)
- implémentation de la lecture des scénarios de poussière
(au lieu des moyennes) soucis d'interpolation à régler.
- simulateur d'observables ACS: (Antoine)
utilitaire solzenangle (pour tout zsa) => A tester.
- solliciter Luca pour produire un scénario MY35 (Ls=0-180) (Francois; fait)
- simulateur d'observables MCS : [inchangé] (Antoine)
Ajoute la possibilité de sortir d'autres champs du GCM (pas dans MCS)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
- D/H : (Loic et Margaux)
- faire une discussion autour de watercap (Loic + Francois + Joseph + Margaux)
- mise en place de l'implémentation dans dyn3dpar "à la Camille" par Margaux
marche bien dans dyn3dpar; à commettre... se pose des questions sur la
sédimentation
- Dynamico: (Antony) travail sur outil d'interpolation de start_archive (LMDZ)
vers start dynamico. Attention à la topo sous-maille.
- Détrainement: [inchangé] Antoine regarde les effets de détrainement
(et plages de détrainement 10h-12h ou 10h-14h). Simus les 5è et 6è mois.
Va comparer les profils d'aérosols Nomad avec ceux fourni par Mike Wolff.
Pour les comparaison avec ACS.
Fait un utilitaire pour calculer les opacités
(poussière, glace d'eau, combiné)
- Gaz non condensables: Franck : Rien de nouveau depuis la semaine dernière.
Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonctiondu % de CO2
C'était mieux en 2010...
NB: Simulation sans nuage radiativement actif en 2010
=> comparer aux simulations de Joseph (Francois+Joseph)
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
[inchangé] Voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Nuages CO2: (Christophe) A simulé un an complet avec bonne conservation de CO2
Penser à inclure tauscaling dans la microphysique
Co2condens mis à jour pour cette version.A voir comment réintégrer dans le trunk.
- Cycle du CO2: (Thomas) [inchangé]
Suite des tests en figeant 2 paramètres à la fois.
Obtiens toujours des dérives assez fortes au cours de années 3-4-5.
Poursuivre les simulations et l'enquête.
- [inchangé] gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas d'assimilation.
=> faire un ticket.
- pb en local au LMD: [inchangé] on arrive pas à tourner en mode prod, à regarder au Bocal.
=> faire un ticket
- newstart: [inchangé] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Joseph : regarde l'incorporation de la chaleur latente du givre d'H2O (travail de thèse
de Margaux). Chantier à reprendre.
16th April 2020
Durée: 1h10
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio: Franck, Ehouarn, Christophe, Antony, Antoine, Margaux, Thomas,
Aymeric, Francois
News:
====
- Le confinement continue...
Commissions récentes:
====================
- r2283: modèle générique
- r2284: ajout de la chimie des ions HCO+
- r2285: initialisation du sous-sol pour le cas sans startfi.nc
- r2286: Vénus
- r2287: Mésoéchelle Vénus
- r2288: Update run.def.MCD5, facteurs de dissipation "mid" et "up"
Tickets récents:
===============
- #44: Penser à nettoyer LMDZ.GENERIC d'anciens scripts et codes (dyn)
- Faire un ticket sur zrecast (ajout des variables de surface) ?
- Faire un ticket: pouvoir concaténer des concat (ajouter controle)
En cours:
========
- simulateur d'observables ACS: (Antoine)
va modifier "terminator" pour s'adapter à d'autres angles solaires
(qui ne dépend pas de l'altitude)
sorties à fournir pour terminator/TGO: T,P,rho,Ps,mmean,espèces (#/cm3 ou kg/m3),...
- solliciter Luca pour produire un scénario MY35 (Francois)
- simulateur d'observables MCS : (Antoine)
Ajoute la possibilité de sortir d'autres champs du GCM (pas dans MCS)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
- Dynamico: (Antony) Premières sorties XIOS d'un run Dynamico. Tout semble correct.
- Détrainement: [inchangé] Antoine regarde les effets de détrainement
(et plages de détrainement 10h-12h ou 10h-14h). Simus les 5è et 6è mois.
Va comparer les profils d'aérosols Nomad avec ceux fourni par Mike Wolff.
Pour les comparaison avec ACS.
- Gaz non condensables: Franck : Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... pas de grosses variation en fonction du % de CO2
C'était mieux en 2010... peut-être sans nuage radiativement actif à l'époque?
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
[inchangé] Voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Nuages CO2: (Christophe) A propos de co2condens: quelques modifs complémentaires à faire.
tests 1D concluant. Test 3D en cours. Peut-être un soucis de conservation
des traceurs dans les thermiques.
A propos des mmr(co2)<1. Pas de soucis en 1D (d'explosion >> 1)
- D/H : (Loic et Margaux) mise en place de l'implémentation dans dyn3dpar "à la Camille"
s'inspire de dyn3dmem (terrestre). Attention au ratio fils/père lorsque père tends vers 0.
Déjà implémenté dans dyn3d.
- Cycle du CO2: (Thomas) [inchangé]
Suite des tests en figeant 2 paramètres à la fois.
Obtiens toujours des dérives assez fortes au cours de années 3-4-5.
Poursuivre les simulations et l'enquête.
- [inchangé] gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + durée écoulée depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas d'assimilation.
- pb en local au LMD: [inchangé] on arrive pas à tourner en mode prod, à regarder au Bocal.
- newstart: [inchangé] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Thomas: travaux sur la MCD:
- utilisation des données Insight comme input (au lieu de VL1)
- implémentation de la lecture des scénarios de poussière (au lieu des moyennes)
9th April 2020
Durée: 1h30
Rappel: Les CR du POULDP sont disponibles sur http://trac.lmd.jussieu.fr/Planeto/wiki/ArchiveReport
En visio: Franck, Christophe, Ehouarn, Lola, Francisco, Anni, Margaux,
Thomas, Antoine, Antony, Francois, Aymeric
News:
====
- changement d’adresse du trac (et d'IP du svn), c'est maintenant sur
http://trac.lmd.jussieu.fr/Planeto
- Lancement de "POULDP Générique" cette semaine
Commissions récentes:
====================
- r2273: ajout (minime) de Francisco pour la chimie ionique
- r2274: nettoyage du code (suppression de miltipl & scatter)
- r2275-76: Mesoéchelle générique - scripts de compilation au propre
- r2277: Mésoéchelle (Venus/générique)
- r2278: nettoyage modèle générique
- r2279: Vénus: mise à jour des .def et xml dans deftank
- r2280: Utilitaire "terminator" pour interpoler diagfi aux terminateurs
et mise à jour dans "localtime" pour lire un stats
- r2281: Ajout de la possibilité de démarrer sans startfi.nc
- r2282: Calcul du dso avec prise en compte de tauscaling (pour cas GCM5, poussière imposée)
Tickets récents:
===============
- #40: fichier used_callphys.def vide en fin de simulation
- #41: sortir dso dans le cas MCDv5.3 (poussière imposée) => r2282
- #42-43: modèle générique
En cours:
========
- modifier "terminator" pour s'adapter à d'autres angles solaires.
- solliciter Luca pour produire un scénario MY35
- sorties à fournir pour terminator/TGO: T,P,rho,Ps,mmean,espèces (#/cm3 ou kg/m3),...
- Dynamico: (Antony) Mise en place d'un démarrage sans startfi.nc.
Prochaine étape, simulation dynamico.
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> Antoine regarde les effets de détrainement (et plages de détrainement
10h-12h ou 10h-14h). Simus les 5è et 6è mois.
Va comparer les profils d'aérosols Nomad avec ceux fourni par Mike Wolff.
Pour les comparaison avec ACS.
=> Gaz non condensables: Franck : Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste les cas
à 94% et à 96% ... en cours
- [inchangé] Reste à voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Nuages CO2: (Christophe) Problèmes de conservation de la masse dans co2cloud
mais mmr>1 parfois après la condensation. A suivre...
- D/H : (Loic et Margaux) Faut-il passer à une implémentation "à la Camille"?
A priori OK pour la dynamique, mais finalement les corrections de Camille
sont dans dyn3dpar/mem (?) terrestre et pas le dyn3dpar planéto... mais peut
-être juste une version intermédiaire de dyn3dpar de Camille à intégrer.
- Cycle du CO2: (Thomas) Suite des tests en figeant 2 paramètres à la fois.
Obtiens toujours des dérives assez fortes au cours de années 3-4-5.
Poursuivre les simulations et l'enquête.
- [inchangé] gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + date éculé depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas d'assimilation.
- pb en local au LMD: [inchangé] on arrive pas à tourner en mode prod, à regarder au Bocal.
- newstart: [inchangé] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- simulateur d'observables MCS : [inchangé] (Antoine)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
- simulateur d'observables ACS: [inchangé] Antoine est dessus pour faire une version dédié "terminateur ACS"
=> faire une sorte de localtime pour interpoler au terminateurs (matin et soir),
température, CO, Ozone, etc.
=> regarde avec Mike pour les opacités (poussière et glace d'eau), réunion prévue semaine pro
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
Antoine Bierjon's mail, 10th April 2020
Bonjour le POULDP !
Je vous fais ce mail pour préciser ce dont on a parlé hier sur l'outil terminator et sa généralisation à plus de solar zenith angle que 90°. Francisco notamment m'avait demandé cela pour étudier les hautes altitudes. Ma question est parle-t-on ici :
- du point de passage dans l'ombre de la planète (que moi j'appelle terminateur, mais peut-être que j'ai tort) : cf PJ 1. Dans ce cas, l'écart de solar zenith angle et de Local Time par rapport à la surface est en effet important quand on monte en altitude, étant à peu près fonction de arccos(Rmars/(Rmars+z)) (soit une déviation de 1.28h (sza=109°) à 200km d'altitude par rapport à la surface)
OU BIEN
- du point tangent à la sphère d'altitude z (confondu à la surface avec le terminateur) : cf PJ 2. Dans ce cas, l'écart de sza et de Local Time est plutôt négligeable (déviation de 3.6x10⁻⁶h à 200km d'altitude par rapport à la surface) et on peut prendre les données interpolées à 90° de sza comme bonne approximation
Merci !
Antoine
2nd April 2020
Durée: 1h30
Rappel: Les CR du POULDP sont disponibles sur http://web.lmd.jussieu.fr/trac-planeto/wiki/ArchiveReport
Visio: Ehouarn, Antoine, Franck, Antony, Christophe, Francois, Margaux, Aymeric, Thomas
News:
====
- MAMO reporté
- Mise en place de POULDP Générique à l'étude
Commissions récentes:
===================
- r2267: Mésoscale
- r2268: Mésoscale
- r2269: Modèle générique: marche à nouveau en OpenMP
- r2270: pour pouvoir compiler avec Dynamico
- r2271: pour Dynamico Mars
- r2272: Nettoyage dans la dynamique commune
=> faire un ticket sur les used_*.def
Rappel: utiliser getin_p au lieu de getin() (ticket #38)
En cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> Antoine regarde les effets de détrainement (et plages de détrainement
10h-12h ou 10h-14h). Simus les 5è et 6è mois.
Va comparer les profils d'aérosols Nomad avec ceux fourni par Mike Wolff.
Pour les comparaison avec ACS.
=> Gas non condensables: Franck : Il y a 10 ans on avait 94% de CO2 et un
meilleur cycle de CO, alors qu'actuellement c'est 95%. Teste une diminution
à 94% ... en cours
- [inchangé] Reste à voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Cycle du CO2: (Thomas) A testé les cas 2D en ne faisant varier que DN,DS ou IN,IS
Va regarder si les prédictions concordent alors.
- Nuages CO2: (Christophe) Regarde quand le mmr(co2)>1, c'est bien après la condensation.
Ce n'est pas qu'une question de température de condensation
Regarde comment limiter mmr(co2) à 1.
- D/H : (Loic et Margaux) a lancé une simu avec fractionnement sur 2 ans.
Résultats comparable à ceux de Franck M, avec ponctuellement de gros rapports
quand il y a peu d'eau. Reste à trouver la meilleure manière de gérer
ces cas.
- [inchangé] gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + date éculé depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas d'assimilation.
- pb en local au LMD: [inchangé] on arrive pas à tourner en mode prod, à regarder au Bocal.
- newstart: [inchangé] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolete car on voudra le faire dans le gcm en interactif.
- simulateur d'observables MCS : [inchangé] (Antoine)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
- simulateur d'observables ACS: [inchangé] Antoine est dessus pour faire une version dédié "terminateur ACS"
=> faire une sorte de localtime pour interpoler au terminateurs (matin et soir),
température, CO, Ozone, etc.
=> regarde avec Mike pour les opacités (poussière et glace d'eau), réunion prévue semaine pro
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico: [inchangé] Ca avance.
26th March 2020
== 26/03/2020 ==
Durée: 1h27
Rappel: Les CR du POULDP sont disponibles sur http://web.lmd.jussieu.fr/trac-planeto/wiki/ArchiveReport
Visio: Ehouarn, Antony, Antoine, Francois, Mike, Margaux, Thomas, Joseph, Christophe, Franck, Loic
News:
====
- Confinement généralisé ! Oskour
- MAMO très certainement annulé... A voir pour annonce prochaine auprès de la communauté.
Commissions récentes:
===================
- r2263: Watercap dans start2archive
- r2264: modification de phyetat0 pour 1+1=2 (qsurf(h20)>=0)
- r2265: dtau poussières pour 1+1=2 GCM 6 --> stockage dans startfi.nc
- r2266: bugfix in the physics for compilation errors
En cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> Antoine regarde les effets de détrainement (et plages de détrainement
10h-12h ou 10h-14h). Comparaison avec les températures MCS ne permet pas de
trouver le bon jeu de paramètres. Antoine va comparer directement les
poussières. François : attention, regarder en moyenne temporelle est peut-être
un peu chaud. Peut-être regarder à différents moments du mois (+ poussière).
- Reste à voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Cycle du CO2: (Thomas) Un peu bloqué. Lié à la convergence des simulations ? Compliqué.
Ehouarn : propose de se ramener à un cas "idéalisé" 2D en ne faisant varier que
les inerties ? Juste pour vérifier que tout marche bien avec l'interpolation et le
reste. Plus méthodique pour attaquer ce problème.
- Nuages CO2: (Christophe) Ca tourne 6 mois, mais le temps de calcul augmente énormément
au 7ème mois (plus de 24h). L'erreur venait d'un random "stop" dans le code,
depuis remplacé par un abort_physics. Pb : VMR(CO2) devient négatif. PB sur le
co2 : conservation de la masse ? Ehouarn a déjà eu ce problème. Schéma de transport
induit des pertes de traceurs. Conservation "à la zob" de la quantité de traceurs
entraine une rétroaction positive du MMR si >1. Lié à la condensation du CO2 ? Cette
équation contient le vmr(co2) qui ne doit pas être >1 --> mettre un test ?
- D/H : (Loic et Margaux) PLusieurs tests en cours:
- certaines différences venaient des optimisations de ifort.
- tests en débug: vdifc semble renvoyer des tendances différentes pour H2O et HDO
même si ces derniers sont identiques en entrée.
Test : cas idéalisé avec H2O et HDO. Margaux & Loic ont remarqué des différences
liées à l'optimisation d'une part, et dans vdifc. HDO est calculé à partir du flux
de H2O.
Rapport D/H initial : vrm(HDO) = vmr(H2O)*2*5*snow ? Conservation de ce rapport ?
Ca a l'air de bien marcher. Sauf quand on rajoute le fractionnement. Transport du
traceur ou du rapport ? En discussion.
- [inchangé] gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + date éculé depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas d'assimilation.
- pb en local au LMD: [inchangé] on arrive pas à tourner en mode prod, à regarder au Bocal.
- newstart: [inchangé] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolete car on voudra le faire dans le gcm en interactif.
- simulateur d'observables MCS : [inchangé] (Antoine)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
- simulateur d'observables ACS: Antoine est dessus pour faire une version dédié "terminateur ACS"
=> faire une sorte de localtime pour interpoler au terminateurs (matin et soir),
température, CO, Ozone, etc.
=> regarde avec Mike pour les opacités (poussière et glace d'eau), réunion prévue semaine pro
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico: Ca avance.
Tour de table :
================
- Antony : Possibilité de démarrer le GCM martien sans start.nc
- Franck : Bug dans newstart problématique pour les gaz non condensables
--> ce n'est pas la cause du problème. De plus, il y a 10 ans, les runs montraient un
meilleur enrichissement. Il semble (Mohammed, Joseph, François) que c'est particulièrement
lié à la quantité de CO2. A creuser et conclure !
- Antoine : Utilitaire d'interpolation de données solar occultation au terminator en cours !
19th March 2020
Durée: 1h20mn
Rappel: Les CR du POULDP sont disponibles sur http://web.lmd.jussieu.fr/trac-planeto/wiki/ArchiveReport
Visio: Aymeric, Ehouarn, Lola, Antoine, Francois, Mike, Margaux, Thomas, Joseph, Christophe, Loic
News:
====
- Confinement généralisé!!
Comissions récentes:
===================
- r2258: Mise à jour de l'option "composition" de newstart
- r2259: Venus
- r2260: watercap: glace d'eau à la surface (remplace le qsurf <0)
- r2261: réintégration des modif de la r2257
- r2262: Vénus
En cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gaz non condensable, Franck regarder ce que ça donne
=> Antoine regarde les effets de détrainement (et plages de détrainement
10h-12h ou 10h-14h).
Compare à présent avec les températures MCS
- Reste à voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Cycle du CO2: (Thomas) Tuning du cycle avec les 4 paramètres de profondeur de glace (h2o)
& inertie thermique Nord/Sud. 4 simulations (on suppose memes paramètres pour
les calottes Nord et Sud) et analyse la 4ème année de simu.
Observe une dérive au bout de 4-5 ans
Bug dans le code des valeurs optimales des paramètre, mais toujours pas de minimum
franc. A suivre...
- Nuages CO2: (Christophe) Ca tourne 6 mois, mais le temps de calcul augmente énormément
au 7ème mois (plus de 24h). SOit il y a eu plantage en cours, soit le temps de
calcul augmente drastiquement... à voir.
- D/H : (Loic et Margaux) PLusieurs tests en cours:
- certaines différences venaient des optimisations de ifort.
- tests en débug: vdifc semble renvoyer des tendances différentes pour H2O et HDO
même si ces derniers sont identiques en entrée.
- 1+1=2 : (Antoine et Joseph) C'est OK en mode GCM5 (poussière+eau).
En mode GCM6 (avec injection, etc.) Antoine confirme qu'on a pas 1+1=2. => Ehouarn.
- [inchangé] gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc
il y a un changement de date car on stoque date de début + date éculé depuis la référence.
A regarder et améliorer, et vérifier que ça ne casse pa le schémas d'assimilation.
- pb en local au LMD: [inchangé] on arrive pas à tourner en mode prod, à regarder au Bocal.
- newstart: [inchangé] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolete car on voudra le faire dans le gcm en interactif.
- simulateur d'observables MCS : [inchangé] (Antoine)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
- simulateur d'observables ACS: Antoine est dessus pour faire une version dédié "terminateur ACS"
=> faire une sorte de localtime pour interpoler au terminateurs (matin et soir),
température, CO, Ozone, etc.
=> regarde avec Mike pour les opacités (pourrière et glac d'eau)
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico: [inchangé] Antony continue à essayer de faire marcher Dynamico en partant du modèle
des planètes géantes... La feuille de route avance bien.
12th March 2020
Rappel: les CR des POULDP sont disponible sur http://lmdz.lmd.jussieu.fr/le-coin-des-developpeurs/cr-de-reunions/pouldp-mars (mais ça va changer, c.f. ci-dessous)
Durée: 1h30mn
Sur place: Francois, Franck L., Ehouarn, Thomas, Antoine, Christophe, Joseph,
Lola, Loic, Margaux
Guets star: Mike W.
News:
====
- ExoMars 2020 est reporté en 2022
- Formation IDRIS MPI la semaine prochaine (Christophe y va).
- Est-ce qu'on mets les CR sur le trac ou le site LMDZ ? => on essaye de les mettre sur Trac les CR
concaténés (sur une même page) (Christophe)
Commissions récentes:
====================
- r2252: bug fix sur les diagnostiques de dso (density-scaled opacity)
- r2253: update pour suivre l'évolution de l'environnement Jean-Zay
- r2254: modèle générique (bug fix sur les aérosols)
- r2255: bug fix dans topmons
- r2256: mise en place environement pour utiliser PGI sur Jean-Zay
- r2257: modifications pour compiler avec PGI; la fonction isnan()
n'est pas standard (utiliser test a/=a plutôt) et derf()
est obsolète, utiliser erf() plutôt.
- Faire un ticket pour mentionner l'utilisation de abort_physic() au lien de stop ou abort
En cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gaz non condensable, faire un coup de
newstart sur le start.nc à Ls=180, mais mettre à jour newstart
avec les valeur de Trainer et al. (prendre les valeurs ACS pour le CO?)
Franck a essayé et trouve un cycle annuel pas terrible...
=> Tests (Ehouarn) de préparation avec des runs:
- problèmes avec le cycle des poussières: derniers tests: ça marche
mieux sans les slope winds... mais il y a le "bug" des rocket dust storm
qui ne détrainent pas par défaut
- Reste à voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5 (=> OK, c'est bien le cas)
- Cycle du CO2: (Thomas) Tuning du cycle avec les 4 paramètres de profondeur de glace (h2o) & inertie thermique Nord/Sud.
4 simulations (on suppose mêmes paramètres pour les calottes Nord et Sud) et analyse la 3ème année de simu.
Quelques différences entre les nouvelles valeurs optimales de profondeur de la glace Sud et les anciennes trouvées par Ehouarn. Peut-être un pb dû au "warm-up" des simulations.
Sinon, le système actuel est déjà pas loin de VL1.
En revanche les valeurs optimales d'inertie thermique correspondent pas trop mal.
- Nuages CO2: (Christophe) Ca marche maintenant sur Occigen! Mais ça plante, dès le 2ème pas de
temps physique dans convadj (pb de zq() de h20_vap)...
- D/H : (Loic et Margaux) Il y avait des pb au niveau de la surface, maintenant corrigés.
Toujours de problèmes de maintient de D/H (même sans fractionnement).
Cherche ou ça dérape dans la physique. et dans watercloud sans sous-pas de temps.
- 1+1=2 : (Antoine et Joseph) C'est OK en mode GCM5 (poussière+eau).
En mode GCM6 (avec injection, etc.) Antoine regarde.
- gestion du temps dans les starts: NB: en regardant 1+1=2, dans start.nc il y a un changement de date
car on stoque date de début + date éculé depuis la référence. A regarder et améliorer, et vérifier que ça ne casse pa le schémas d'assimilation.
- watercap : (Joseph, ticket #20) Ca tourne bien mais les fichiers ne sont pas identiques au bit près...
Toujours d'actualité... suggestion de Francois: un problème lié à l'albédo?
- pb en local au LMD: on arrive pas à tourner en mode prod, à regarder au Bocal.
- newstart: [inchangé] Attention aux VMR proposés en commentaires dans l'option "composition" qui semblent ne plus être d'actualité.
pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolète car on voudra le faire dans le gcm en interactif.
- simulateur d'observables MCS : [inchangé] (Antoine)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
Antoine est en train de résoudre les tickets sur le transfert radiatif des poussières.
- simulateur d'observables ACS: Antoine est dessus pour faire une version dédié "terminateur ACS"
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico: [inchangé] Antony continue à essayer de faire marcher Dynamico en partant du modèle des planètes géantes... La feuille de route avance bien.
5th March 2020
Durée : 1h30
Sur place : Christophe, Joseph, Antony, Antoine, Aymeric, Thomas
Visio : Ehouarn
News:
====
- Conf grand public MAMO le jeudi soir (19h-20h30)
- Annulation du MAMO ? (Coronavirus)
Commissions récentes:
====================
- r2241 - r2245: Modèle Générique/ Titan (jvatant)
- r2246 : resolved ticket #32 --> gros nettoyage des calculs des opacités des différentes poussières
- r2247 : partially resolved ticket #14 --> Amélioration des sorties diagfi pour pour connaitre
les longueurs d'ondes précises des calculs d'opacité (voir attributs de variable).
- r2248 : bugfix slopewinds --> corrections de bugs dans topmons (pb actualisation après rocketduststorms,
les inventaires de poussière n'étaient donc plus consistants).
D'autres bugs sont présents dans topmons en cours de débug.
- r2249 - r2251: Mesoscales
En Cours:
========
- MCD6.0: [inchangé] Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gaz non condensable, faire un coup de
newstart sur le start.nc à Ls=180.
il faut revoir le H2 => à renormaliser... (confirmé par Jean-Yves)
=> Tests (Ehouarn) de préparation avec des runs:
- problèmes avec le cycle des poussières: derniers tests: ça marche
mieux sans les slope winds...
- Reste à voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5
- Cycle du CO2: (Thomas) Tuning du cycle avec les 4 paramètres de profondeur de glace (h2o) & inertie thermique Nord/Sud.
4 simulations (on suppose memes paramètres pour les calottes Nord et Sud) et analyse la 3ème année de simu.
Nettoyage et amélioration de xvik. On le passe ensuite à la moulinette par fit_inertia_depth.
Grosses différences entre les nouvelles valeurs optimales de profondeur de la glace Sud et les anciennes trouvées par Ehouarn.
En revanche les valeurs optimales d'inertie thermique correspondent pas trop mal.
Thomas a lancé des nouvelles simulations avec les valeurs optimales qu'il a trouvé pour comparer avec les obs/anciennes simus.
- Nuages CO2: (Christophe) La microphysique du CO2 tourne en mode débug sur ciclad (40 jours).
Pb de script pour lancer une année entière de simus. Est-il nécessaire de tourner un an en mode débug ?
Ehouarn suggère de passer direct à l'exécutable classique, et d'essayer par acquis de conscience le mode débug sur OCCIGEN.
- D/H : [inchangé] (Loic) Il y a des cas 3D (mais pas en 1D) ou le ratio D/H même décorrélé,
n'est pas préservé, mais dans les cas ou les valeurs de traceur sont trop
faibles. Pour l'instant la correction est: lorsque le flux d'H2O est "trop
petit" on mets le flux d'HDO (et de H2O) à zéro. Et ça marche.
Va poursuivre la validation encourageante.
- watercap : (Joseph, ticket #20) Ca tourne bien mais les fichiers ne sont pas identiques au bit près...
Oscillations de la masse totale d'eau analogues à celles de Christophe pour le CO2.
Pour l'instant la variable tient simplement le compte d'eau sublimée dans les calottes (watercaptag)
- newstart: Attention aux VMR proposés en commentaires dans l'option "composition" qui semblent ne plus être d'actualité.
pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolete car on voudra le faire dans le gcm en interactif.
- simulateur d'observables MCS : [inchangé] (Antoine)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...).
Antoine est en train de résoudre les tickets sur le transfert radiatif des poussières.
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico: (Antony) Antony continue à essayer de faire marcher Dynamico en partant du modèle des planètes géantes...
La feuille de route avance bien.
Tour de table:
=============
- Christophe : le CR est sur le trac planéto !
- Aymeric : Deux règles absolues pour les commits.
1) On ne commit que les changements qui compilent !
2) S'assurer de la rétrocompatibilité !! (flag .false. par défaut)
13th February 2020
Durée: 45 min
Sur place: Christophe, Joseph, Lola, Vincent, Antoine
News:
====
- Meeting TGO la semaine prochaine
- Comité des utilisateurs de l'IDRIS le 5 mars ! Préparer les retours d'ici dans pas longtemps
- RDV conférence Elbereth du 26 au 28 février à l'IAP (Vincent, Antoine, Antony, Thomas, Joseph)
Commissions récentes:
====================
-r 2240 : Initialisation des des variables mu0_int et ztim1 -> résolution des crashs en mode débug 1D
En Cours:
========
- MCD6.0: [inchangé] Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gaz non condensable, faire un coup de
newstart sur le start.nc à Ls=180.
il faut revoir le H2 => à renormaliser... (confirmé par Jean-Yves)
=> Tests (Ehouarn) de préparation avec des runs:
- problèmes avec le cycle des poussières: derniers tests: ça marche
mieux sans les slope winds...
- Reste à voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5
- Cycle du CO2: [inchangé] penser à mettre à jour xvik... (Thomas a fait une simu sur
deux ans pour regarder)
- Nuages CO2: (Christophe)
Test le 3D en debug sur Ciclad, tente de REtourner sur Occigen.
- D/H : [inchangé] (Loic) Il y a des cas 3D (mais pas en 1D) ou le ratio D/H même décorrélé,
n'est pas préservé, mais dans les cas ou les valeurs de traceur sont trop
faibles. Pour l'instant la correction est: lorsque le flux d'H2O est "trop
petit" on mets le flux d'HDO (et de H2O) à zéro. Et ça marche.
Va poursuivre la validation encourageante.
- watercap : [inchangé] En cours de debug... (Joseph, ticket #20) problèmes de reproductibilité
Problèmes en local avec "ulimit -s unlimited", sur Occigen aussi...
- newstart: [inchangé] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolete car on voudra le faire dans le gcm en interactif.
- simulateur d'observables MCS : (Antoine)
A améliorer en mettant les autres champs (aerosols, fonctions poids,...). Antoine est en train de résoudre les tickets sur le transfert radiatif des poussières.
- lessivage par le CO2 : [inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico: [inchangé] (Antony)Feuille de route sur le portage établie (c.f. ticket #26).
Antony fait tourner le cas DYNAMICO-générique en mode "Mars", en cours.
Tour de table:
=============
- Vincent démarre la prise en main du modèle méso-échelle. Félicitations.
- Christophe continue à mettre à jour le wiki -> Compte-rendus du pouldp sur le wiki.
- Joseph : ouverture prochaine d'un ticket pour améliorer la représentation de l'eau liquide
- /!\ A faire pour les japonisants : VACCINS /!\
6th February 2020
Rappel: les CR des POULDP sont disponible sur http://lmdz.lmd.jussieu.fr/le-coin-des-developpeurs/cr-de-reunions/pouldp-mars
Durée: 65mn
Sur place: Antony, Thomas, Christophe, Antoine, Ehouarn, Joseph, Aymeric
News:
====
- Demande DARI pour demain! (11h)
Commissions récentes:
====================
- r2233: Modèle générique, possibilité de compiler/tourner sans IOIPSL
- r2234-37: Mésoéchelle (cf ticket #30: le mésoéchelle ne marche plus
avec callcond=.true.); amélioration de l'install sur Ciclad
- r2238: Simulateur d'observables MCS (petites améliorations) en vue
de l'ajout des aérosols.
NB: une commission de Christophe n'est pas passé... à élucider...
En Cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gaz non condensable, faire un coup de
newstart sur le start.nc à Ls=180.
il faut revoir le H2 => à renormaliser... (confirmé par Jean-Yves)
=> Tests (Ehouarn) de préparation avec des runs:
- problèmes avec le cycle des poussières: derniers tests: ça marche
mieux sans les slope winds...
- Reste à voir comment on fait pour les cas "dust storm".
penser à vérifier que ces simus sont avec iphysiq=5
- Cycle du CO2: penser à mettre à jour xvik... (Thomas a fait une simu sur
deux ans pour regarder)
- Nuages CO2: (Christophe) Simus 3D sur 11 jours.
Problème de conservation de la masse totale de CO2 due à la dynamique.
Simu de 1an en 1D marche.
Test le 3D en debug sur Ciclad
Tente de tourner sur Occigen, ça plante après quelques pas.
- D/H : [inchangé] (Loic) Il y a des cas 3D (mais pas en 1D) ou le ratio D/H même décorrélé,
n'est pas préservé, mais dans les cas ou les valeurs de traceur sont trop
faibles. Pour l'instant la correction est: lorsque le flux d'H2O est "trop
petit" on mets le flux d'HDO (et de H2O) à zéro. Et ça marche.
Va poursuivre la validation encourageante.
- watercap : En cours de debug... (Joseph, ticket #20) problèmes de reproductibilité
Problèmes en local avec "ulimit -s unlimited", sur Occigen aussi...
- newstart: [inchangé] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23),
probablement obsolete car on voudra le faire dans le gcm en interactif.
- simulateur d'observables MCS:
A améliorer en mettant les autres champs (aerosols, fonctions poids,...)
- lessivage par le CO2:[inchangé] Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico: (Antony)Feuille de route sur le portage établie (c.f. ticket #26).
Antony fait tourner le cas DYNAMICO-générique en mode "Mars", en cours.
Tour de table:
=============
- Ehouarn: va lancer les simulations MCDv6, mais problème avec les scenarios...
- Joseph: on pourrait se partager à tour de rôle la sténo du Pouldp?
- Christophe: voudrait organiser une heure ou deux pour mettre à plat le plan du track (Wiki)
peut-être la semaine prochaine?
30th January 2020
== 30/01/2020 ==
Rappel: les CR des POULDP sont disponible sur http://lmdz.lmd.jussieu.fr/le-coin-des-developpeurs/cr-de-reunions/pouldp-mars
Durée: 1h (12h-13h)
En Visio: Personne --> CF News !
Sur place: Jean-Yves, Ehouarn, Vincent, Christophe, Thomas, Antony, Joseph
News:
====
- on laisse tomber la visio?
- Semaine du 18-21 --> Team meeting NOMAD & ACS et François risque d'être absent...
- Soutenance de Jan 31 janvier 2020
Commissions récentes:
====================
- r2224: Bug fix 1D avec les pentes (étaient systématiquement mises à 0)
- r2225: Nettoyage autour des ondes non-orographiques
- r2226: Correction de bugs dans les "rocket dust storms" & "slope winds".
tauscaling utilisé mais pas initialisé, et gradient température
au premier niveau dans topmons non plus...
- r2227: corrections dans l'utilitaire "localtime"
- r2228: ajout pour le 1D (initialisation du jour du péréhélie par la
valeur de Ls de périhélie fournie par l'utilisateur)
- r2229-2232 : modèle générique (modèle des thermiques)
En Cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gaz non condensable, faire un coup de
newstart sur le start.nc à Ls=180. Retuner H2? Jean-Yves a regardé : vmr(h2) trop important. A re-tuner
Jean-Yves doit envoyer le coefficient à corriger à Ehouarn ! Finalement He OK.
Souci de conservation des traceurs dans la dynamique ? A voir...
=> Tests (Ehouarn) de préparation avec des runs:
- sans thermosphère, avec 54 couches et toutes les
nouvelles paramétrisation (lessivage+slope winds+
rocket dust storms+ GW oro et non-oro), OK ça a tourné 2 ans.
- avec en plus la thermosphère et 73 couches, ça tourne avec la r2213
- Reste à voir comment on fait pour les cas "dust storm".
Ca tourne mais le cycle des poussières beaucoup trop lourd... Problèmes liés à la révision r2226 (paramétrisation des slopes winds)
- Super-cycling de la chimie. [INCHANGE] On en a finalement peut-être pas besoin? A priori on part plus sur ça...
- Cycle du CO2 [INCHANGE]: penser à mettre à jour xvik... (Thomas Pierron prend les choses en main)
- Nuages CO2: (Christophe) Simus 3D sur 11 jours.
Problème de conservation de la masse totale de CO2
La physique ne crée une variation que de 1.e-9%, par appel à la physique
le problème vient de la dynamique variation de 1e-3% sur 100 appels
Problème 1D corrigé (fallice), à lancer sur 1 an dans un coin (histoire de vérifier...)
Problème avec ou sans microphysique du CO2...
Va tester sur 1 mois/1 an.
Ca avance, augmentation de moyens (Christophe part sur Occigen). Débug en cours (prometteur).
- D/H : [INCHANGE] (Loic) Il y a des cas 3D (mais pas en 1D) ou le ratio D/H même décorrélé,
n'est pas préservé, mais dans les cas ou les valeurs de traceur sont trop
faibles. Pour l'instant la correction est: lorsque le flux d'H2O est "trop
petit" on mets le flux d'HDO (et de H2O) à zéro. Et ça marche.
Va poursuivre la validation encourageante.
- watercap : Débug en cours (Joseph, ticket #20)
- newstart: [INCHANGE] pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23) (redondant avec le flag r2228 mais à mettre en place
tout de même pour faciliter l'utilisation & les tests)
- simulateur d'observables MCS [INCHANGE]: A tester... (Antoine et Antony) mais avec les simulations...
A améliorer en mettant les autres champs (aerosols, fonctions poids,...)
- lessivage par le CO2 [INCHANGE]: Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico : Feuille de route sur le portage établie (c.f. ticket #26).
Antony se lance...
Sorties XIOS mises en place (1st step !)
Tour de table / Remarques complémentaires:
=========================================
- Ehouarn: test du cycle des poussières --> abandon des slopes winds pour la MCD 6 ?
Création des scénarios de poussière clim/cold/warm pour la MCD6. Réactualisation avec MY32, MY33, MY34 (encore plus clean qu'avant !).
Cold/warm sont créés en prenant les extrema (min/max) d'opacité sur 10 ans et en retirant (respectivement ajoutant) 50% d'opacité.
- Christophe: Ajouter le numéro de révision à l'exec... Le ticket a été fait.
Nettoyage en cours du code de microphysique CO2 (Anni). Il s'est approprié le trac et le wiki ! Merci à toi
- Antony : si les formations FORTRAN intéressent des gens, Patrick part à la retraite ! Il ne faut pas trainer... Niveau avancé Mars/Novembre !
23th January 2020
Rappel: les CR des POULDP sont disponible sur http://lmdz.lmd.jussieu.fr/le-coin-des-developpeurs/cr-de-reunions/pouldp-mars
Durée: 65mn
En Visio:
Sur place: Ehouarn, Vincent, Christophe, Lola, Antoine, Joseph, Aymeric
Guest star: Archibald (stagiaire de 3ème)
News:
====
- Thomas et Antony sont à la formation Fortran à l'IDRIS
- Ne pas oublier la demande DARI A8 planéto
Commissions récentes:
====================
- r2219: Update pour le Mesoéchelle (script de génération des starts du GCM)
- r2220-2221: Update/nettoyage de la paramétrisation des ondes de gravité non-orographiques
- r2222: Modèle générique, update des exemples d'XML pour XIOS
- r2223: Mise en place des sorties XIOS dans la physique Mars
En Cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gaz non condensable, faire un coup de
newstart sur le start.nc à Ls=180. Retuner H2? (Jean-Yves regarde)
=> tuning des ondes orographiques et non-orographiques
(Antoine et Antony). Les tests montrent que les ondes non-orographiques
ont un effet au dessus de 60km tandis que les ondes orographiques ont un
effet plutôt à partir de 80km... On reste avec le tuning de référence
pour l'instant.
=> Tests (Ehouarn) de préparation avec des runs:
- sans thermosphère, avec 54 couches et toutes les
nouvelles paramétrisation (lessivage+slope winds+
rocket dust storms+ GW oro et non-oro), OK ça a tourné 2 ans.
- avec en plus la thermosphère et 73 couches, ça tourne avec la r2213
- Reste à voir comment on fait pour les cas "dust storm".
- Super-cycling de la chimie. On en a finalement peut-être pas besoin?
- Cycle du CO2: penser à mettre à jour xvik...
- Nuages CO2: (Christophe) Simus 3D sur 11 jours.
Problème de conservation de la masse totale de CO2
La physique ne crée une variation que de 1.e-9%, par appel à la physique
le problème vient de la dynamique variation de 1e-3% sur 100 appels
Simu de 1an en 1D, ça plante dans co2snow (fallice négatif!)
Problème avec ou sans microphysique du CO2...
Va tester sur 1 mois/1 an.
- D/H : (Loic) Il y a des cas 3D (mais pas en 1D) ou le ratio D/H même décorrélé,
n'est pas préservé, mais dans les cas ou les valeurs de traceur sont trop
faibles. Pour l'instant la correction est: lorsque le flux d'H2O est "trop
petit" on mets le flux d'HDO (et de H2O) à zéro. Et ça marche.
Va poursuivre la validation encourageante.
- watercap : A mettre en place... (Joseph, ticket #20)
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23)
- simulateur d'observables MCS: A tester... (Antoine et Antony) mais avec les simulations...
A améliorer en mettant les autres champs (aerosols, fonctions poids,...)
- lessivage par le CO2: Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années.
C'est le cas entre année 2 et 3; il faut 1/2 année pour bien lancer le système.
En global on a toujours plus de poussière que la cible; taux de soulèvement à revoir?
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico: Feuille de route sur le portage établie (c.f. ticket #26).
Antony se lance...
Tour de table / Remarques complémentaires:
=========================================
- Ehouarn: quelques problèmes en debug sur un cas test 54 niveau avec toutes les paramétrisations enclenchées...
- Antoine/Joseph: attention: pb avec "slopes" en 1D. A corriger...
- Christophe: Ajouter le numéro de révision à l'exec... Faire un ticket
16th January 2020
Rappel: les CR des POULDP sont disponible sur http://lmdz.lmd.jussieu.fr/le-coin-des-developpeurs/cr-de-reunions/pouldp-mars
Durée: 1h25mn
En Visio:
Sur place: Jean-Yves, Lola, Antony, Christophe, Thomas, Antoine, Ehouarn, Loic
News:
====
- Il faut préparer la demande DARI A8 pour le 7 février. => Ehouarn
- Le projet LMDZ-Transport va peut-être s'arrêter là...
Commissions récentes:
====================
- r2214: Bug fix 1D (il y avait pb si heure de départ différente de 0)
- r2215: correction dans l'utilitaire "localtime" pour gérer les valeurs
indéfinies (e.g. issues de zrecast à certaines pressions ou
altitudes sous la surface)
- r2216: Exemples de scripts de soumission sur Ciclad & Occigen dans
deftank
- r2217: Bug fix sur l'option "composition" (initialisation des valeurs
des non-condensables à celles de MSL)
- r2218: Flag "latentheat" mis à true par défaut, et changé de nom pour
"latentheat_surfwater" pour être plus explicite
En Cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gas non condensable, faire un coup de
newstart sur le start.nc à Ls=180. Retuner H2?
=> tuning des ondes orographiques et non-orographiques
(Antoine et Antony). Les simus tournent mais avec ondes non-orographique
ça crashe sur Ciclad avec la r2216...
Antoine a des simulations avec ondes oro, dont une "qui va bien"...
reste à voir ce qui change avec ondes non-orographiques.
=> Tests (Ehouarn) de préparation avec des runs:
- sans thermosphère, avec 54 couches et toutes les
nouvelles paramétrisation (lessivage+slope winds+
rocket dust storms+ GW oro et non-oro), OK ça a tourné 2 ans.
- avec en plus la thermosphère et 73 couches, ça tourne avec la r2213
- Reste à voir comment on fait pour les cas "dust storm".
- Super-cycling de la chimie. On en a finalement peut-être pas besoin?
- Cycle du CO2: penser à mettre à jour xvik...
- Nuages CO2: (Christophe) finalement il n'y a aps de problème avec les thermiques.
Les tests de conservation sont tous OK.
Simu de 1an en 1D, ça plante dans co2snow (fallice négatif!)
Va faire aussi des tests en 3D sur plusieurs jours.
- D/H : (Loic) Il y a des cas 3D (mais pas en 1D) ou le ratio D/H même décorrélé,
n'est pas préservé, mais dans les cas ou les valeurs de traceur sont trop
faibles. Pour l'instant la correction est: lorsque le flux d'H2O est "trop
petit" on mets le flux d'HDO (et de H2O) à zéro. Et ça marche.
Va poursuivre la validation encourageante.
- watercap : A mettre en place... (Joseph, ticket #20)
- newstart: pouvoir mettre le Ls du périhélie en input (Joseph, ticket #23)
- simulateur d'observables MCS: A tester... (Antoine et Antony) mais avec les simulations...
A améliorer en mettant les autres champs (aerosols, fonctions poids,...)
- lessivage par le CO2: Antoine a fait des simulations annuelles avec différents
paramètres pour voir l'effet du scavenging ratio. Analyse en cours.
Une simulation sur 3 ans tourne (scavenging ratio de 100 == référence) afin de
regarder si on préserve le cycle de poussière sur des années. Résultat mitigé
sur les 2 premières années (différences jusqu'à 50% en relatif, mais plutot dans
les régions polaires.
Sera à valider avec les observations (notamment PFS) et la microphysique
complète.
- Dynamico: Feuille de route sur le portage établie (c.f. ticket #26).
Antony se lance...
Tour de table / Remarques complémentaires:
=========================================
- Jean Yves: Finalement c'est OK pour He, mais peut-être pas pour H2.
=>Ehouarn doit mettre à disposition les stats de la simu de référence 6.0alpha
- Antoine: Pourrait-on stocker le numéro de révision du GCM dans l'executable. et/ou
les sorties?
- Christophe: Les commentaires: est-ce qu'on peut faire le ménage de "vieux trucs"?
Oui! Besoin d'une "doc utilisateur" pour expliquer/désigner des configurations
(en complément de deftank); faire des HowTo? Dans le doc (en LateX) ou pages html?
- Lola: amènera un stagaire de 3ème au POULDP la semaine prochaine
9th January 2020
Rappel: les CR des POULDP sont disponible sur http://lmdz.lmd.jussieu.fr/le-coin-des-developpeurs/cr-de-reunions/pouldp-mars
Durée: 1h10mn
En Visio:
Sur place: Antoine, Joseph, Thomas, Antony, Christophe, Anni, Ehouarn, Jean-Yves
News:
====
- Il faut préparer la demande DARI A8 pour le 7 février.
- Le Kick-Off du nouveau contrat ESA MCD (01/2020->01/2022) a eu lieu le 08/01
Commissions récentes:
====================
- r2203: Vénus
- r2204: Update de l'utilitaire localtime: il y avait un pb si le fichier
d'entrée était sorti d'un zrecast (car mauvais nombre de niveaux
pour aps() et bps())
- r2205: Mésoéchelle...
- r2206: début de mise à jour de la doc par Antony
- r2207: correction sur localtime pour compiler avec ifort
- r2208-2212: Mésoéchelle...
- r2213: Correction de bug dans calchim (problème uniquement avec l'ionosphère)
En Cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gaz non condensable, faire un coup de
newstart sur le start.nc à Ls=180
=> tuning des ondes orographiques et non-orographiques
(Antoine et Antony)
=> Tests (Ehouarn) de préparation avec des runs:
- sans thermosphère, avec 54 couches et toutes les
nouvelles paramétrisation (lessivage+slope winds+
rocket dust storms+ GW oro et non-oro), OK ça a tourné 2 ans.
- avec en plus la thermosphère et 73 couches, ça tourne avec la r2213
- tournera aussi avec chaleur latente de la glace d'eau en surface.
- Super-cycling de la chimie.
- Cycle du CO2: penser à mettre à jour xvik...
- Nuages CO2: co2condens marche bien maintenant, mais pb avec les thermiques
qui ne conservent pas ni le gaz de CO2 ni la glace de CO2.
La hauteur limite des plumes zmax est parfois trop faible
(quelques couches)et le bilan entrainement/détrainement n'est
pas bon. => marche avec une correction: en discute avec Aymeric
=> bug dans thermiques faire un ticket (Christophe)
- D/H : Loic est dessus
- watercap : A mettre en place... (Joseph, ticket #20)
- chaleur latente de la glace d'eau à la surface: Joseph regarde sur 10 ans.
A la dixième année: écart en moyenne de l'ordre de 1% sur les températures.
Donc c'est OK => A enclencher par défaut (dans vdifc_mod.F mettre
latentheat=.true.) => Joseph fait la commission
- newstart: pouvoir mettre le Ls du périhéli en input (Joseph, ticket #23)
- simulateur d'observables MCS: A tester... (Antoine et Antony)
- lessivage par le CO2: Antoine voit avec Anni et Christophe sur le scavenging
ratio. Sera validable avec les observations (e.g. PFS) et la microphysique
complète. => simulations test en cours (avec/sans slope winds). NB les
dernières correction (non_initialisation) changent les choses.
Tour de table / Remarques complémentaires:
=========================================
- Jean-Yves: a fait des comparaisons avec NGIMS (deep dips de MAVEN)
desités de CO2, Ar, N2, He. Marche bien pour N2,
moins bien pour l'Ar
Le He est très trop abondant (facteur 2-3)
dans le GCM dans la haute atmosphère (valeur de ref. à réviser)
- Joseph: a une stagiaire L3 (modèle 1D)
à quand l'outil pour évaluer la dissipation?
- Thomas: utilise la MCD6 et les données Viking 1 et 2, vérifie le cycle du CO2
- Christophe: en 1D la lecture des profils initiaux (fichiers profil_**) n'est
pas prévu pour les CCN de CO2. Attention, il faut que CCN mass et number
soient compatibles car on en déduit le rayon des particules!
- Anni: quel article de référence citer pour le GCM? On attend Francois!
- Ehouarn: La priorité: finaliser un GCM6.0 et une MCD6.0
- Antoine: Il faut des métriques de validation du GCM.
2nd January 2020
Rappel: les CR des POULDP sont disponible sur http://lmdz.lmd.jussieu.fr/le-coin-des-developpeurs/cr-de-reunions/pouldp-mars
Durée: 35 mn
En Visio:
Sur place: Ehouarn, Antony, Thomas
News:
====
- Arrivée de Thomas Pierron (ingé CDD) au LMD aujourd'hui.
Commissions récentes:
====================
- r2202: utilitaire "simu_MCS_temp" (simulateurs d'observables MCS de Antoine)
- r2203: Vénus
En Cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gas non condensable, faire un coup de
newstart sur le start.nc à Ls=180
=> tuning des ondes orographiques et non-orographiques
(Antoine et Antony)
=> Tests (Ehouarn) de préparation avec des runs:
- sans thermosphère, avec 54 couches et toutes les
nouvelles paramétrisation (lessivage+slope winds+
rocket dust storms+ GW oro et non-oro), OK ça a tourné 2 ans.
- avec en plus la thermosphère et 73 couches, ça plante...
- Cycle du CO2: mettre à jour xvik...
- Nuages CO2: co2condens marche bien maintenant, mais pb avec les thermiques
qui ne conservent pas ni le gaz de CO2 ni la glace de CO2
=> bug dans thermiques faire un ticket (Christophe)
- D/H : Loic est dessus
- watercap : A mettre en place... (Joseph, ticket #20)
- chaleur latente de la glace d'eau à la surface: Joseph regarde sur 10 ans
- newstart: pouvoir mettre le LS du périhéli en input (Joseph, ticket #23)
- simulateur d'observables MCS: A tester...
- lessivage par le CO2: Antoine voit avec Anni et Christophe sur le scavenging
ratio. Sera validable avec les observations (e.g. PFS) et la microphysique
complète.
Tour de table / Remarques complémentaires:
=========================================
- Antony: toujours des problèmes avec ses simulations sur Ciclad ou Occigen.
2020
19th December 2019
Durée: 1h05
En Visio:
Sur place: Ehouarn, Christophe, Anni, Antony, Margaux, Antoine, Joseph
News:
====
- Anni, de retour de AGU: Tout le monde attend le MAMO!
- pas de pouldp Mars les deux prochaines semaines, c'est Noël!
Comissions récentes:
===================
- r2186-88 : Vénus
- r2189 : correction de bug avec le scavenging par le CO2
- r2190 : amélioration du script d'auto-install IOIPSL, plus besoin de ksh!
- r2191 : ajout des options C dans l'arch-gfortran
- r2192-94 : Vénus
- r2195 : Vénus, correction du géopot calculé dans la dynamique et passé à
la physique. Pas d'incidence (autre que changements d'arrondis)
pour Mars
- r2196-8 : Vénus
- r2199 : Paramétrisation des slope winds de Margaux
- r2200 : Vénus
- r2201 : Netoyage de la paramétrization des rocket dust storms.
=> transport Van Leer amélioré ou on saute plusieurs mailles
mis de côté; on utilise des sous-pas de temps (plus robuste)
En Cours:
========
- MCD6.0: Configuration GCM6.0: A mettre en place (ticket #24)
=> pour tuner les gas non condensable, faire un coup de
newstart sur le start.nc à Ls=180 et
=> tuning des ondes orographiques et non-orographiques
(Antoine et Antony)
=> Tests (Ehouarn) de préparation avec des runs:
- sans thermosphère, avec 54 couches et toutes les
nouvelles paramétrisation (lessivage+slope winds+
rocket dust storms+ GW oro et non-oro), en cours
- avec en plus la thermosphère et 73 couches, ça plante...
- Nuages CO2: co2condens marche bien maintenant, mais pb avec les thermiques
qui ne conservent pas ni le gaz de CO2 ni la glace de CO2
=> bug dans thermiques
=> faire un ticket (Christophe)
- D/H : Loic est dessus
- watercap : A mettre en place... (Joseph)
- chaleur latente de la glace d'eau à la surface: Joseph regarde sur 10 ans
- newstart: pouvoir mettre le LS du périhéli en input
- simulateur d'observables MCD: Antoine doit commiter sa version (temperature de
jour et nuit), en entrée on peut prendre des diagfi ou stats (journées
perpétuelles).
- lessivage par le CO2: Antoine voit avec Anni et Christophe sur le scavenging
ratio. Sera validable avec les observations (e.g. PFS) et la microphysique
complète.
Tour de table / Remarques complémentaires:
=============
- Christophe: travaille aussi sur le ticket #22 "localtime"
- Anthony: a aussi des problèmes en combinant zrecast et localtime...
Joue avec xarray
12th December 2019
Durée: 1h5 (rappel les CR des POULDP sont en ligne:
http://lmdz.lmd.jussieu.fr/le-coin-des-developpeurs/cr-de-reunions/pouldp-mars )
En visio:
Sur place: Franck L., Antony, Joseph, Ehouarn, Margaux, Aymeric, Antoine
Guest star: Sébastien L.
News:
====
- La formation LMDZ, c'est fait! Avec 7-8 planétos...
Commissions récentes:
====================
- r2185: mise à jour du script d'install IOIPSL sur machine CentOS du LMD
Tour de table:
=============
- Ehouarn: mise en place de scripts "automatiques" d'install des versions
planétos pour la formation LMDZ.
disponibles sur http://www.lmd.jussieu.fr/~lmdz/planets/
Scripts à parfaire (et à mettre également sur svn)
- Aymeric: Deborah travaille sur une version unifiée des ondes non-oro
de Francois Lott.
- Franck: RAS
- Margaux: Slopes winds au propre, va commettre.
- Antony: runs thermospherique TGO sans la iono plante après 5 mois...
joue avec la dissipation pour voir
NB: Aymeric mentionne qu'il y a un article de Jesch et al.
dans MDPI Atmosphere du 15 oct 2019 sur l'aerobraking TGO
- Joseph: - regarde les simulations avec l'effet de la chaleur latente de la
glace d'eau à la surface. Simulations pluriannuelles.
- mise en place de la variable "watercap" en cours
- Veut modifier newstart pour changer le Ls du périhélie (et non pas
le sol du périhélie comme c'est actuellement le cas) => ticket!
- sujet de stage L3: jouer avec le 1D à haute obliquité
- Seb : ravi d'être là (RDEL)
- Antoine: RDEL. Doit toujours voir avec Anni pour les rapports de lessivage
Simulateurs d'observables en stand-by.
Doit voir la semaine prochaine avec Margaux pour la passation.
A faire: => une rubrique "En cours" à la place...
=======
Tickets:
=======
- Fermeture des tickets #2 (réglé avec la r2170) et #19 (redondant avec #12)
5th December 2019
Durée: 50mn
En visio: Ehouarn
Sur place: Christophe, Antoine, Antony, Margaux
News:
====
- Antoine a soutenu son stage.
- Margaux est docteure!
Commissions récentes:
====================
- r2183: mise à jour des scripts d'install IOIPSL sur Ciclad avec NetCDF récent
- r2184: mise en place du lessivage instantané de poussière et H2O par le CO2
=> voir avec Anni pour la valeur du taux d'agrégation.(action Antoine)
Tour de table:
=============
- Antoine: lessivage par le CO2. A tester et retester (par ex. par Christophe),
notamment pour le lessivage de la glace d'eau.
continue sur le simulateur d'observable MCS. Roland et Luca sont
intéressés pour faire un simulateur PFS
- Antony: aérofreinage TGO: va faire des simus thermosphere sans les ions pour
aller plus vite, éventuellement sur Occigen.
va comparer ses résultats sur Ciclad avec ceux obtenus par Ehouarn
sur Occigen.
va mettre à jour la doc du modèle sur le svn.
- Christophe: Va refaire un co2condens à part pour tester dans son cas de figure
pour clarifier les problèmes qu'il rencontre.
- Margaux: cherche à optimiser ses dernières semaines au LMD. Réfléchit à un
article pour les slope winds. Va commiter les slope winds.
- Ehouarn: A relancé les discussions avec Chao Wang qui a des problèmes pour
installer et faire tourner le GCM sur Ciclad.
Prépare la formation LMDZ et la livraison de la MCDv6alpha
A mis un ticket (#22) pour un problème remonté par Katsuyuki
Nogushi (Nara University) avec l'utilitaire localtime: quand le
fichier d'entrées contient des variables indéfinies (-9.99E33),
par exemple suite à un zrecast, alors les résultats interpolés
utilisent cette valeur...
A faire:
=======
Tickets:
=======
- #22 : pb avec localtime si l'entrée contient des valeurs "undefined (aka
-9.99E33) lors des interpolations
28th November 2019
Durée 1h10.
Présents:
========
En visio:
Sur place: Vincent, Antoine, Loic, Ehouarn, Jean-Yves
News:
====
- Francois en route pour l'Antarctique
- Soutenance de thèse de Margaux demain!
Commissions récentes:
===================
- r2175: option "composition" in Newstart
- r2179: option to the sublimation/condensation latent heat
release from the surface water ice
- r2180: script to install IOIPSL on local CentOS machines => scripts
à améliorer en utilisant les arch du GCM =>Ticket #21
- r2181: Enforce using dissip_fac_mid and dissip_fac_up parameters
read from .def files to be used by inidissip.
- r2182: Fixes for the 1D model (CO2 saturation outputs & watercaptag)
Tour de table:
=============
- Vincent: Toujours sur les données, encore du traitement à appliquer
sur l'altitude touchée (surface ou aérienne).
- Antoine: continue l'évaluation des ondes orographiques (en poussière forcée
ou injectée), résultats moins bons avec l'injection par rapport
à MCS.
Chantier lessivage par le CO2/H2O: doit nettoyer le code et
commettre. Prépare sa soutenance (pour lundi prochain)
- Loic: Toujours sur les tests HDO sans fractionnement. Le problème arrive
surtout lorsqu'il y a peu de vapeur (possible pb d'arrondi du calcul
de D sur une valeur très petite de H). => tests de mettre un seuil
pour corriger si valeurs trop petites.
- Jean-Yves: A récupéré les données NGIMS dans la thermo, va faire des
comparaisons avec les simulations sur les espèces lors de deep dips
(~120-130km).
Avec la nouvelle diffusion les rapports entre espèces sont meilleurs
(sauf pour Ls=0-30) et l’échappement de H est plus intense qu'avant.
- Ehouarn: Séjour au Weizmann institute avec Eran Vos (équipe Oded Ahranson)
pour suivre ses travaux. => ajout d'une variable 2D "z_lag" dans le
startfi pour savoir à quelle profondeur est enterrée la glace.
=> modifications (mineurs) dans vdifc_mod, qu'il faut encore valider
et faire converger avec la version de Joseph (variable "watercap").
21th November 2019
Présents: ======== Sur place : Loic, Christophe, Antony, Antoine, François, Joseph, Lola, Margaux, Anni, Aymeric, Frank Visio : Francisco News: ======== - François : Voyage en Antarctique avancé à jeudi 28 novembre (retour début Février ?) - Joseph : admin du site exobiologie.fr - Francisco (rappel) : arrive dans une semaine ! Tour de table: ======== - Joseph : RAS - Lola : RAS - Margaux : préparation de thèse + commit option "latent heat" du givre (.false. par défaut). - Anni : RAS - Loic : HDO -> fin des problèmes dans vdifc. Il reste des problèmes en haute altitude sur le fractionnement (alors qu'il est désactivé). Effet numérique (atm très sèche à cette altitude) ou bien problème de transport ? Plusieurs hypothèses à tester : transport horizontal, vertical ? Commentaire de François : dans le modèle terrestre HDO et H2O ne sont pas considérés comme des traceurs indépendants (transport de D/H au lieu de transporter HDO). Solution possible pour Loic : envisager d'amalgamer HDO et H2O pour s'affranchir de cet effet. - Christophe : perte de masse CO2 toujours problématique. La compréhension générale avance cependant. - Antony : Pb avec co2condens simulation avec thermosphère et chimie (-> aérofreinage TGO) ? Probablement lié à un état initial instable. En fait certainement lié au fait que CICLAD refuse les jobs interactifs trop gros. A voir. - Francisco : fournir un jeu de données centralisé pour les équipes qui pourraient être intéressées (e.g. ACS) ? François : excellente idée, Ehouarn ou Antoine pourraient lancer ce genre de simulation sur OCCIGEN et sortir des profils très faciles à lire (ASCII). Pb : pas de scénario de poussière ni UV pour MY35. Demander à Luca de faire des MaJ de son scénario tous les 3 mois ? - François : utilitaire traceur2column à finaliser. Vérifier qu'il marche bien avec les start.nc. Pas mal de nouvelles enigmes atmosphériques : nouvelles mesures Ar, N2, comportement anormal de l'oxygène. Beaucoup de travail en perspective en Antarctique. Arrivée inopinée de Frank : "Il faut être prudents, comme pour le méthane.". Pour François, il semble compliqué que le comportement de l'oxygène soit lié à la chimie du fait des énormes quantités d'oxygène mises en jeu. Frank répond "Si c'est vrai !". Une autre possibilité : filament d'air polaire avec concentrations différentes ? Bof bof car les autres gaz ne se comportent pas de la même manière. "Il est urgent d'attendre." - Frank : en train d'extraire des profils pour SPICAM. Une première : run de 5 années martiennes consécutives avec scénarios de poussière adaptés. A voir si l'eau se comporte bien (inquiétude de François) et CO. - Aymeric : RAS - Antoine : présentation sur les ondes de gravité orographiques. Paramétrisation du coefficient de proportionnalité "GKDRAG". Exploration avec les valeurs 0.1, 1, 10 (censée être une constante ~1 d'après la littérature - Forget 1999). Utilisation du nouveau cycle des poussières (MCD 6.alpha). Objectif de l'étude : possibilité de rephaser les ondes de marées dans le GCM par rapport à MCS à l'aide des ondes de gravité. Avec poussière libre ça marche pas trop, mais quand on revient sur de la poussière prescrite le GKDRAG=1 est vraiment chouette. Vérifier à d'autres Ls car Ls=0 est un moment compliqué dynamiquement. De manière générale l'étude va dans le bon sens. A faire : ======= - Joseph : test de l'option "latent heat" du givre commitée par Margaux pour s'assurer que c'est bien physique et que ça ne perturbe pas le modèle actuel de Mars. Le cas échéant, activer l'option par défaut dans le modèle. - Un volontaire ? Commiter des fichiers "launch" (scripts bash) dans deftank pour lancer les runs sur CICLAD (commande "qsub launch") - François, Ehouarn et Antoine : faire grosse simulation sur OCCIGEN pour fournir des profils intéressants de sorties de modèle, et maintenir les bases de données à l'avenir. Commencer par répondre à Kevin.
14th November 2019
Présents:
========
Sur place: Christophe Francois, Lola
Antoine, Antony, Margaux, Joseph
En retard: Aymeric, Anni
Pas de visio.
Tour de table:
=============
- Antony: Probleme zrecast lié à la mémoire (?) ne voulais pas tourner sur baker, avec un fichier de 8 Go. Cela a marché sur penn... Message d'erreur : "processus arrêté..." (ce n'était pas résolu par l'istruction "limit -s unlimited"). Poursuit ces travaux sur Aerofreinage. doit apprendre a faire tourner le GCM...
- Antoine: Rapport de stage de fin d'étude rendu. Soutenance le lundi 2 decembre à 14h. Simulation avec ondes orographiques avec un run type MCD5.
- Margaux: rentre de 2 semaines de vacances. Soutenance confirmée le 29/11 à 14h, rapports reçus.
- Joseph: Probleme des gaz non-condensables résolus: les starts des "mauvaises" simulations étaient trop riches en CO2.
La variable "watercap" (budget de glace pérenne des glaciers) a été codée. En cours de test.
- Aymeric:
-Développement des outils de calcul de diagnostique dynamique. A nommer.
A modifier pour se comporter comme un utilitaire sans avoir à modifier le
code python.
- Doit passer planetoplot (actuellement en Python 2) en Python 3.
Propose de refonder planetoplot sur "xarray" voire de remplacer
planetoplot par xarray...
- Christophe: progrès sur la conservation de la masse de CO2 : c'est lié à la condensation atmosphérique. Discussion: cela pourrait etre du à la gestion des particules de glace de CO2 qui ne sont pas comptés comme il faut dans le bilan de pression à la surface :
- Francois: Information sur les nouvelles mesures de Ar, N2, O2, CO par Curiosity.
A faire:
=======
- Antoine: faire tourner le simulateur d'observable sur les stats et comparer: 1) MCS moyenné sur des mois de 30 degres LS 2) MCS simulé avec run MCD5.3 3) MCS simulé avec la nouvelle simulation avec nouvelle poussière MCD6alpha.
- Margaux: commiter la prise en compte de la chaleur latente pour les temperature de la glace d'eau en surface. => utilisation par Joseph pour les paleoclimats
- Joseph: changer la lecture de solperi dans "controle" des starts en Lsperi => Changer ecriture et lecture des starts ; test à la lecture pour convertir la valeur de solperi en Lsperi équivalente. Faire en sorte que le Lsperi actuel donne exactement le solperi de d'habitude. A discuter avec Ehouarn ?
- Francois: envoie des transparents de la réunion GCM Mars de Mardi 12/11
- Ticket a ajouter : modifier les valeurs de références de AR, N2, CO, O2 dans les start de références.
7th November 2019
Présents:
========
En visio:
Sur place: Loic, Franck L. Christophe, Francois, Ehouarn,
Antoine, Antony, Aymeric, Joseph
News:
====
- Projet Venus Climate Database (contrat ESA) va voir le jour.
- Embauche (au 01/01/2020) de Thomas Pierron ingé. CDD
- Jérémy Leconte a deux ingés CDD qui démarent sur le générique/Dynamico;
Lancer des POULDP Générique ?
- Projet (Francois & Ehouarn) de faire une réunion Q&A du GCM Mars +
bonnes pratiques...
- Félicie viendra au LMD le 14 novembre matin 10h pour réunion sur son projet
- Francisco sera à Paris du jeudi 28 novembre au 8 décembre
Commissions récentes:
====================
- r2173 : Ondes non-orographiques, lecture de RUWMAX (== EPMAX) dans .def
- r2174 : commissions modèle générique (.def pour cas de référence Saturne)
Tour de table:
=============
- Antoine: tests sur les ondes orographiques, l'effet est croissant jusqu'en
haut du modèle (120km) vers les poles. NB: le modèle est trop chaud par
rapport à MCS (même sans changer le forçage des ondes de gravité), donc il
y a par ailleur un problème poussière/températures.
Pas trop d'effet sur les ondes de marée thermique (par rapport à MCS)
- Joseph: a fait un ticket "watercap" pour référencer les dépots de glace
d'eau anciennes. A testé son "newstart", va le commettre.
- Franck a son login au CINES!
Souhaite utiliser l'albedo du modèle pour la chimie; repèrer les cas
bare ground / h20/CO2 ice sur la valeur de l'albedo seule
- Christophe: toujours des soucis de conservation du CO2
- Loic: toujours le pb de D/H qui varie (alors qu'il ne devrait pas),
pb reproduit en 1D, le pb semble lié aux nuages.
- Ehouarn: simulations pour Gabriella. MCDv6 alpha à faire
- Aymeric: a un équivalent python de l'outils streamfunction.F90
=> streamfunction ne marche plus avec les diagfis. => Ticket!
A aussi un outil qui calcule la PV!
- Antony: Compare les observations TGO (densité) avec les sorties modèles
(via zrecast).
- Francois: RAS. A quelques tickets sur le feu...
A faire:
=======
- Ajouter Antony à la liste pouldp (Ehouarn)
- Aymeric: mettre à jour le README des outils pour indiquer/pointer sur
ses outils.
31th October 2019
Présents:
========
En visio: Francisco
Sur place: Christophe, Vincent, Joseph, Antoine, Francois, Ehouarn, Loic
Guest star: Félicie
News:
====
- Antony Delavoie arrive la semaine prochaine (Ingé. LMD)
Commissions récentes:
====================
- aucune cette semaine
Tour de table:
=============
- Félicie: coucher de soleil sur Mars (installation itinérante);
projet d'appli ou fond d'écran ? Lié d'une façon ou autre
à position /situation réelle.
+ projets de bas-reliefs lumineux.
- Christophe: perte de masse dans nuages de CO2 avec lors de la condensation
- Joseph: - enrichissement des gas non condensables => regarder dans les
runs MCD que les diagnostiques sont bons (au site MSL vs global) ?
- va regarder la paramétrisation de l'albedo du givre
- Loic: - ajouter un traceur surfacique (épaisseur de "vielle glace") pour
pouvoir vérifier la conservation de l'eau.
- problèmes avec le rapport D/H. Chereche à vérifier si tout
marche bien sans fractionnement. Ca n'est pas toujours le cas.
pb dans les nuages et à la surface (tests/seuils de conservation
de la quantité de traceur).
- Antoine: - rédige son rapport (deadline début novembre)
- ondes de gravité orographiques : fait des simus en changeant
le coefficient de la trainée (entre 0.1 -référence-, 1. -ref.
de la littérature- et 10 => la simu plante au mois #7)
- Francisco : avec les ondes non oro de Gabriella: problèmes sur les
valeurs par défaut dans le code => trop instable.
En effet, Ehouarn et Gabriella sont en train de regarder
quels doivent être les bons paramètres
- Vincent : analyse des données MOLA => séparations en groupes de nuages.
- Francois : affaires en cours comme la semaine dernière...
- Ehouarn : simulations ondes de gravité non-orographiques avec Gabriella.
simulations MCDv6.0 alpha finies.
A faire:
=======
- Francois et Christophe regardent ensemble les histoires de CO2 au Bocal.
- Joseph: commettre le newstart => à voir au Bocal.
Tickets:
=======
- Rien de neuf.
24th October 2019
Présents:
========
Sur place: Aymeric, Christophe, Vincent, Joseph, Jean-Yves, Antoine,
Margaux, Francois
News:
====
- L'artiste Félicie d'Estienne d'Orves passera au LMD le 31 matin pour discuter du projet "Eyes on Mars". Rendez vous à 10h
Commissions svn récentes:
========================
- r2171: Mise en commun du z2sig.def haute resolution
- r2172 : modif de arch-centos
Tour de table:
=============
- Francois : Nouvel utilitaire "tracercolumn" compile mal sur son PC/gfortran. A suivre.
- Francois travaille sur les données CO PFS => NetCDF compliqué
- Idée Francois: simulateurs d'image vu par les Cameras type VMC, HRSC, EXI, etc... Ticket pour commencer : mise en place d'un choix d'altitude pour la visualisation au format "sphere" sur la MCD en ligne.
- Christophe: travaille sur perte de masse de CO2 liée aux nuages CO2
- Sans sédimentation: problème résolu en en remodifiant un test dans co2cloud (qui avait déjà été modifié récemment...)
- Avec sédimentation : il reste un problème de conservation de la masse. François: "Etes vous sur du diagnostic ?". A suivre.
- Margaux a testé le lessivage de la glace d'eau par les nuages de CO2 dans la nuit polaire. Prometteur !
- Joseph : a finalisé un newtart pour changer la concentration de CO2 de manière cohérente.
Tickets:
=======
=> Nouveau : retuning ondes orographique (Lott) pour comparaison MCS
(pour Antoine)
17th October 2019
Présents:
========
En visio: Francisco
Sur place: Christophe, Vincent, Joseph, Antoine, Lola, Franck,
Francois, Margaux, Anni
News:
====
- Félicie d'Estienne d'Orves passera au LMD le 30 ou 31 pour discuter
du projet "Eyes on Mars".
Commissions svn récentes: 2168 et +
========================
- r2168: Ne concerne pas Mars
- r2169: Mise du z2sig.def nouveau dans "deftank"
- r2170: Update de la Chimie (FL+FGG) pour tourner
sans ionosphère + regroupement d'options de la photochimie
Tour de table:
=============
- Ehouarn: simulations MCDv6alpha. Problèmes avec les simulations dust storm.
simulations ondes non orographiques de Gabriella
- Franck: couplages photolyse-GCM: il faudrait prendre en compte l'albedo variable
du sol.
- Antoine: simulations pour tester le lessivage de la poussière par le CO2
tests avec tau variable ou fixe (0.5). Ca semble marcher dans le bon sens.
Paramètre clé: le "scavenging ratio" dont la valeur est a optimiser (10,50?)
A venir: lessivage de l'eau par le CO2...
Le simulateur d'observables MCS est amélioré avec sorties des
temperatures de nuit (poussières à venir)
- Francois: Nouvel utilitaire (fait avec Mohammed) pour lire les traceurs et en
déduire la colonne "gencol" => "tracercolumn". A intégrer dans les outils.
Joseph va tester.
- Margaux: Rédaction du manuscript de thèse terminée!
- Joseph: (ticket #13) Quand les gas non-condenasbles ne sont pas présents, newstart
(option "composition") ne mets pas correctement (implicitement) leurs valeurs.
Va corriger.
- Vincent: analyse les données MOLA; pas (encore) sur le GCM.
- Christophe: Se lance en 1D sur la microphysique, recherche les pertes de masse
identifiéespar Luca. Sedimentation avec vlz_fi en cause?
Il y a d'autres processus qui
semblent aussi ne pas conserver la masse...
- Francisco: va implémenter la chimie de HDO et les ions de l'eau.
A faire:
=======
- Ajouter un z2sig.def haute résolution (~150) (Margaux)
- nettoyer dans physiq les référence à tauref /vs tauref_scenario (Ehouarn)
et vérifier le tau utilisé (visible?) dans le calcul des photolyses (Franck)
Tickets:
=======
#5 : Franck regarde cette histoire de conservation des traceurs
(test sur 10 sols, probablement trop court)
Ajout #14: écrire dans le diagfi la longueur d'onde de ref (TES ou MCS) de la dsod
10th October 2019
Présents:
========
En visio:
Sur place: Christophe, Aymeric, Anni, Franck L., Joseph, Margaux,
Antoine
News:
====
- Rappel: Formation LMDZ, inscrivez vous vite!
Commisions svn récentes: 2166 et +
=======================
- r2166: => correction pour tourner avec opacité constante et injection
1ere commission d'Antoine!
- r2167: => nettoyage du code: on vire l'ancienne dynamique LMDZ.MARS
=> Reste à mettre à jour la doc!
Tour de table:
=============
Franck: Francisco a envoyé un fix pour tourner sans la ionosphère.
Va tester puis commettre.
Antoine: Le supercycling de la chimie change t'il beaucoup les résultats?
Simulations avec lessivage en cours
A Faire:
=======
- mettre à dispo le "nouveau" z2sig.def pour tous (Ehouarn)
Tickets:
=======
#1: MCDv6.0alpha : des simulations sont en cours (sans lessivage)
=> envoyer l'info à tous (simu sur Occigen)
=> Penser à ouvrir des comptes pour tous sur Occigen (relancer Ehouarn)
#2: Franck est dessus
#6: cf rev 2167 , doc encore à faire. => Ehouarn ajoute le mail de Lola à
ce sujet pour rappel de points à éclaircir/améliorer.
3th October 2019
Présents:
========
En visio: Francisco
Sur Place: Christophe, Vincent, Margaux, Aymeric, Antoine, Francois, Joseph,
Antony (avec un mois d'avance), Lola, Loic
News:
====
- Rappel: il y a les formations Modélisation en novembre
et LMDZ en décembre 9-11.
- Info: la MCD est pluggée dans JMARS
Commisions svn récentes: (2165)
=======================
- r2165 => fichiers d'architecture pour Jean Zay
Tour de table:
=============
- Antoine:
1) Le lessivage de la poussière par la glace de CO2 dans la nuit polaire
version simple: on lessive une proportion de la poussière lorsque le
CO2 condense. => ajout via un flag "scavco2cond"
Test 1D sur 30 sol dans la nuit polaire prometteurs
2) Le supercycling de la chimie : la chimie représente ~40% du GCM
sans supercycling. NB: watercloud ~17% et radiatif lwmain ~5%
avec supercycling, la chimie prend 19% du temps total.
- Francisco : travil en cours si les ions sont dans traceur.def alors
on ne fait pas la ionosphère.
- Joseph: regarde l'évolution des gas non condensables (suite du travail
de stage de Mohammed). C'est complètement lié aux fichiers start!!
suivant qu'on ait 95% ou 97% de CO2 dans le système...
- Ehouarn: Simu MCDv6. Sur Occigen ~10h par mois.
(c'était ~6h par mois en configurations MCDv5.3)
Ajout de tickets:
=================
Info générale: les tickets tracs sont là:
http://web.lmd.jussieu.fr/trac-planeto/
chacun pourra ouvrir/éditer les tickets trac:
A priori les identifiants/mdp svn doivent fonctionner.
Sinon on va voir (avec Laurent) pour ouvrir des comptes aux nouveaux
Ajout ticket #8 : faire une doc sur vdifc
Ajout ticket #9 : documenter mesoinstall
Ajout ticket #10 : calcul de la constante des gas pour mélanges variables
Attachments (19)
- report_pouldp_06022020 (3.5 KB) - added by 5 years ago.
- report_pouldp_30012020 (5.1 KB) - added by 5 years ago.
- report_pouldp_23012020 (3.9 KB) - added by 5 years ago.
- report_pouldp_16012020 (4.2 KB) - added by 5 years ago.
- report_pouldp_09012020 (3.8 KB) - added by 5 years ago.
- report_pouldp_02012020 (1.8 KB) - added by 5 years ago.
- report_pouldp_19122019 (2.5 KB) - added by 5 years ago.
- report_pouldp_12122019 (1.9 KB) - added by 5 years ago.
- report_pouldp_05122019 (2.0 KB) - added by 5 years ago.
- report_pouldp_28112019 (2.1 KB) - added by 5 years ago.
- report_pouldp_21112019 (4.0 KB) - added by 5 years ago.
- report_pouldp_14112019 (2.5 KB) - added by 5 years ago.
- report_pouldp_07112019 (2.2 KB) - added by 5 years ago.
- report_pouldp_31102019 (2.1 KB) - added by 5 years ago.
- report_pouldp_24102019 (1.4 KB) - added by 5 years ago.
- report_pouldp_17102019 (2.4 KB) - added by 5 years ago.
- report_pouldp_10102019 (1.1 KB) - added by 5 years ago.
- report_pouldp_03102019 (1.7 KB) - added by 5 years ago.
- report_pouldp_13022020 (3.3 KB) - added by 5 years ago.
Download all attachments as: .zip
