| 1 | == REPORT OF THE WORKING GROUP MEETING (January 30, 2014) == #REPORTOFTHEWORKINGGROUPMEETINGJanuary302014 |
| 2 | Participants : Rachid Benshila, Daley Calvert, Julien Calandreau, Jérome Chanut, Laurent Debreu, Christian Ethe, Rachel Furner, Claire Levy, Gurvan Madec, Sébastien Masson |
| 3 | |
| 4 | === 1) Discuss the update of the 3.6 version === #a1Discusstheupdateofthe3.6version |
| 5 | ==== History : ==== |
| 6 | * [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF] version (and corresponding NEMO subroutines) has been updated by Julien Calandreau (INRIA) and Laurent Debreu |
| 7 | |
| 8 | * Jérome Chanut has made a lot of improvements and bug fixes in his Gilbraltar Strait configuration |
| 9 | |
| 10 | * Same for Sébastien Masson in the Pulsation ANR project |
| 11 | |
| 12 | * Other developments in LIM, PISCES |
| 13 | |
| 14 | ==== Bugs that have been identified and corrected since version 3.4 : ==== |
| 15 | * zwork array (J. Chanut): a global array « zwork » was defined in the libmpp module and (due to a bug in the [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF] converter) this led to several problems and in particular the update step was not performed … |
| 16 | |
| 17 | * Diffusion coefficients in the corners of the sponge layern were wrong (J. Chanut) |
| 18 | |
| 19 | * Scale factors were not correctly computed in the ghost cells (where it is assumed that we are in full steps) (J. Chanut) |
| 20 | |
| 21 | -> same problem for OBC |
| 22 | |
| 23 | * TKE (S. Masson): near the boundaries, spatial averaging operations led to wrong Kz in the first cell inside the high resolution domain. Extrapolation is possible but does not fulfill the requirement that with, a mesh refinement factor of 1, the obtained solution of the HR grid is identical to a solution computed without a zoom. S. Masson has implemented exchanges of TKE between coarse and fine grids. |
| 24 | |
| 25 | -> Same problem for GLS and for OBC |
| 26 | |
| 27 | * Asselin Filter (J. Chanut) : when updating the « after » (n+1) fields, the corresponding corrections of the « now » fields were not done. |
| 28 | |
| 29 | * arrays out of bound in procnames in MPI (wrong declarations of [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF] profiles in agrif_declare_variables: jpi, jpj have been replaced by ncli,nlcj) |
| 30 | |
| 31 | ==== Bugs / Problems that have been identified ==== #BugsProblemsthathavebeenidentified |
| 32 | · The bathymetry connection check is not complete (add a check that the coarse grid bathy have been correctly updated) |
| 33 | |
| 34 | · Reproducibility is not guaranteed (due to MPI problem) when ln_spc_dyn is set to true in namelist (which means that we take 0 as special value for dynamics in interpolations/updates). |
| 35 | |
| 36 | · Recursive updates. When more than 1 zoom level are set, the update should start from the highest resolution grid and go recursively to the root grid. Laurent implemented changes that makes this possible. |
| 37 | |
| 38 | * Use of [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF] for degradation online (problème de north pole, coarsening) |
| 39 | |
| 40 | * Average operators with update values |
| 41 | |
| 42 | * High order operators with coastlines |
| 43 | |
| 44 | * BDY problem have been solved |
| 45 | |
| 46 | ==== Improvements ==== #Improvements |
| 47 | * [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF] and time splitting (OK for version 3.6, end of the year, Jérome) |
| 48 | * only consider forward integration of the forcing mode permits to ensure volume conservation |
| 49 | * without update of the parent grid right hand side |
| 50 | * Performance of parallel runs |
| 51 | * allows the distribution of the processors among the fine grids |
| 52 | |
| 53 | '''''Rachid Benshila''' is ready be in charge of the merge'' |
| 54 | |
| 55 | === 2) What are the known open problems ? (not to be solved for 3.6) === #a2Whataretheknownopenproblemsnottobesolvedfor3.6 |
| 56 | * Number of ghost cells ('''''L. Debreu''''') |
| 57 | * Antartic problem (periodicity): require only one ghost cell instead of 2 |
| 58 | * Number of required ghost cells for high order advection schemes may be larger than 2 |
| 59 | * -> Remove hardcoding of the number of ghost cells, do the required modifications to the [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF] “profiles” (Laurent) |
| 60 | |
| 61 | * North fold problem |
| 62 | |
| 63 | * Timings for each grid |
| 64 | |
| 65 | === 3) What are the main next developments? === #a3Whatarethemainnextdevelopments |
| 66 | * Coupling between grids at the level of the barotropic mode (in time splitting version) (Jérome) |
| 67 | * requires some modification in the organization of the step routine of NEMO |
| 68 | |
| 69 | * Vertical grid refinement (Laurent) |
| 70 | * Online coarsening for tracers using [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF] (L. Debreu, G. Madec) |
| 71 | * [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF] with variable volume |
| 72 | |
| 73 | === 4) Nesting tools: === |
| 74 | Two stages: |
| 75 | |
| 76 | a) for 3.6 version |
| 77 | |
| 78 | * agree on all the points: bathymetry smoothing, update of coarse grids (take the maximum as default) |
| 79 | |
| 80 | * Incorporate corrections done by Jean Marc |
| 81 | |
| 82 | * The nesting tools are only essential for the creation of bathymetry: |
| 83 | * forcings can be interpolated on the fly |
| 84 | * -> Do a list available tools for 3D interpolation (initialization / restart fields) |
| 85 | * For 3.6 keep the 3D fields interpolation in nesting tools |
| 86 | |
| 87 | b) next version |
| 88 | |
| 89 | * The definition of vertical coordinates should be done outside NEMO |
| 90 | |
| 91 | * Nesting tools for the s coordinates |
| 92 | |
| 93 | ==== 5) Documentation: ==== |
| 94 | * Update of the nesting tools 3.6 documentation |
| 95 | |
| 96 | * Update of the users guide: describes what is possible or not with [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF]. |
| 97 | |
| 98 | * Documentation for nesting in NEMO/[http://www-ljk.imag.fr/MOISE/AGRIF AGRIF]: starts from the notes of Sebastien, then circulate the powerpoint |
| 99 | |
| 100 | * GMD paper |
| 101 | |
| 102 | === 6) Standard configuration and sanity tests === #a6Standardconfigurationandsanitytests |
| 103 | The standard configuration should: |
| 104 | |
| 105 | * Include more than one level of embedding (test of the recursive update starting from the finer resolution grid(s)) |
| 106 | |
| 107 | * Uses time splitting |
| 108 | |
| 109 | * Include a zoom in a sea ice area (Labrador Sea, LIM3) |
| 110 | |
| 111 | * Include a zoom with PISCES |
| 112 | |
| 113 | * Zoom in Perou Chili |
| 114 | |
| 115 | Sanitiy tests |
| 116 | |
| 117 | * Zoom with a mesh refinement of 1 |
| 118 | |
| 119 | * Test constancy preservation for tracers |
| 120 | |
| 121 | Longer visions of having [http://www-ljk.imag.fr/MOISE/AGRIF AGRIF] zoom in AMM |
| 122 | |
| 123 | S coordinates zoom inside a z coordinate parent grid. |