New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
WorkingGroups/AGRIF (diff) – NEMO

Changes between Version 8 and Version 9 of WorkingGroups/AGRIF


Ignore:
Timestamp:
2014-03-27T09:43:54+01:00 (10 years ago)
Author:
clevy
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • WorkingGroups/AGRIF

    v8 v9  
    1 [[TOC(heading=AGRIF,AGRIF/*, depth=1)]] 
     1[[PageOutline]] 
     2 
     3[[TOC(heading=AGRIF,AGRIF/*, depth=2)]] 
    24 
    35= '''AGRIF''' = 
     
    5860 
    5961== REPORT OF THE WORKING GROUP MEETING (January 30, 2014) == 
    60 Participants : Rachid Benshila, Daley Calvert, Julien Calandreau, Jérome Chanut, Laurent Debreu, Christian Ethe, Rachel Furner, Claire Levy, Gurvan Madec, Sébastien Masson 
    61  
    62   
     62Participants : Rachid Benshila, Daley Calvert, Julien Calandreau, Jérome Chanut, Laurent Debreu, Christian Ethe, Rachel Furner, Claire Levy, Gurvan Madec, Sébastien Masson 
    6363 
    64641) Discuss the update of the 3.6 version 
    6565 
    66   
    67  
    68 History : 
    69  
    70 -       AGRIF version (and corresponding NEMO subroutines) has been updated by Julien Calandreau (INRIA) and Laurent 
    71  
    72 -       Jérome has made a lot of improvements and bug fixes in his Gilbraltar Strait configuration 
    73  
    74 -       Same for Sébastien in the Pulsation ANR project 
    75  
    76 -       Other developments in LIM, PISCES 
    77  
    78   
    79  
    80 Bugs that have been identified and corrected since version 3.4 : 
    81  
    82 -       zwork array (J. Chanut): a global array  « zwork » was defined in the libmpp module and (due to a bug in the AGRIF converter) this led to several problems and in particular the update step was not performed … 
    83  
    84 -       Diffusion coefficients in the corners of the sponge layern were wrong (J. Chanut) 
    85  
    86 -       Scale factors were not correctly computed in the ghost cells (where it is assumed that we are in full steps) (J. Chanut) 
     66History : 
     67 
     68 * AGRIF version (and corresponding NEMO subroutines) has been updated by Julien Calandreau (INRIA) and Laurent 
     69 
     70 * Jérome has made a lot of improvements and bug fixes in his Gilbraltar Strait configuration 
     71 
     72 * Same for Sébastien in the Pulsation ANR project 
     73 
     74 * Other developments in LIM, PISCES 
     75 
     76Bugs that have been identified and corrected since version 3.4 : 
     77 
     78 * zwork array (J. Chanut): a global array  « zwork » was defined in the libmpp module and (due to a bug in the AGRIF converter) this led to several problems and in particular the update step was not performed … 
     79 
     80 * Diffusion coefficients in the corners of the sponge layern were wrong (J. Chanut) 
     81 
     82 * Scale factors were not correctly computed in the ghost cells (where it is assumed that we are in full steps) (J. Chanut) 
    8783 
    8884è same problem for OBC 
    8985 
    90 -       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. 
     86 * 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. 
    9187 
    9288è Same problem for GLS and for OBC 
    9389 
    94 -       Asselin Filter (J. Chanut) : when updating the « after » (n+1) fields, the corresponding corrections of the « now » fields linked with the Asselin filter were not done. 
    95  
    96 -       arrays out of bound in procnames in MPI (wrong declarations of AGRIF profiles in agrif_declare_variables: jpi, jpj have been replaced by ncli,nlcj) 
    97  
    98   
    99  
    100 Bugs  / Problems that have been identified 
    101  
    102 ·      The bathymetry connection check is not complete (add a check that the coarse grid bathy have been correctly updated) 
    103  
    104 ·      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). 
    105  
    106 ·      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. 
    107  
    108   
    109  
    110 - Use of AGRIF for degradation online (problème de north pole, coarsening) 
    111  
    112 - Average operators with update values 
    113  
    114 - High order operators with coastlines 
    115  
    116   
    117  
    118 - BDY problem have been solved 
    119  
    120   
     90 * Asselin Filter (J. Chanut) : when updating the « after » (n+1) fields, the corresponding corrections of the « now » fields linked with the Asselin filter were not done. 
     91 
     92 * arrays out of bound in procnames in MPI (wrong declarations of AGRIF profiles in agrif_declare_variables: jpi, jpj have been replaced by ncli,nlcj) 
     93 
     94Bugs  / Problems that have been identified 
     95 
     96·      The bathymetry connection check is not complete (add a check that the coarse grid bathy have been correctly updated) 
     97 
     98·      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). 
     99 
     100·      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. 
     101 
     102 * Use of AGRIF for degradation online (problème de north pole, coarsening) 
     103 
     104 * Average operators with update values 
     105 
     106 * High order operators with coastlines 
     107 
     108 * BDY problem have been solved 
    121109 
    122110Rachid is ready do the merge (as soon the beta version is fixed) 
    123111 
    124   
    125  
    1261122) What are the known problems / bugs? 
    127113 
    128   
    129  
    130 -       Number of ghost cells 
    131  
    132 o   Antartic problem (periodicity): require only one ghost cell instead of 2 
    133  
    134 o   Number of required ghost cells for high order advection schemes may be larger than 2 
    135  
    136 o   -> Remove hardcoding of the number of ghost cells, do the required modifications to the AGRIF “profiles”  (Laurent) 
    137  
    138 -       North fold problem ? 
    139  
    140 -       Timings for each grid 
    141  
    142   
     114 * Number of ghost cells 
     115 
     116o   Antartic problem (periodicity): require only one ghost cell instead of 2 
     117 
     118o   Number of required ghost cells for high order advection schemes may be larger than 2 
     119 
     120o   -> Remove hardcoding of the number of ghost cells, do the required modifications to the AGRIF “profiles”  (Laurent) 
     121 
     122 * North fold problem ? 
     123 
     124 * Timings for each grid 
    143125 
    1441263) What are the main next developments? 
    145127 
    146   
    147  
    148 -       Time splitting (OK for version 3.6, end of the year, Jérome) 
    149  
    150 o   only consider forward integration of the forcing mode permits to ensure volume conservation 
    151  
    152 o   (refluxing or conservative interpolation at the boundary). 
    153  
    154 o   Tells the users that we have only AGRIF working in time splitting version in the forward version of the code. 
    155  
    156 -       Parallel runs (Laurent) 
    157  
    158 -       Update of the parent grid rhs before barotropic step(s) 
    159  
    160   
     128 * Time splitting (OK for version 3.6, end of the year, Jérome) 
     129 
     130o   only consider forward integration of the forcing mode permits to ensure volume conservation 
     131 
     132o   (refluxing or conservative interpolation at the boundary). 
     133 
     134o   Tells the users that we have only AGRIF working in time splitting version in the forward version of the code. 
     135 
     136 * Parallel runs (Laurent) 
     137 
     138 * Update of the parent grid rhs before barotropic step(s) 
    161139 
    1621404) Version 5.0 (time and space schemes in the system). 
    163141 
    164   
    165  
    166   
    167  
    168 -       Vertical grid refinement (Laurent) 
    169  
    170 o   Probably before version 5 
    171  
    172 -       Online degradation using AGRIF 
    173  
    174 -       Agrif with variable volume 
    175  
    176 o   (update of the volume, and modification of boundary conditions) (Jérome). 
    177  
    178   
     142 * Vertical grid refinement (Laurent) 
     143 
     144o   Probably before version 5 
     145 
     146 * Online degradation using AGRIF 
     147 
     148 * Agrif with variable volume 
     149 
     150o   (update of the volume, and modification of boundary conditions) (Jérome). 
    179151 
    1801525) Nesting tools: 
    181153 
    182   
    183  
    184 Two stages:  
    185  
    186 a) for 3.6 version   
    187  
    188 -       agree on all the points: bathymetry smoothing, update of coarse grids (take the maximum as default) 
    189  
    190 -       Incorporate corrections done by Jean Marc 
    191  
    192 -       The nesting tools are only essential for the creation of bathymetry: 
    193  
    194 o   Do a list available tools for 2D/3D interpolation 
    195  
    196 o   For 3.6 keep the 3D fields interpolation in nesting tools 
    197  
    198   
     154Two stages: 
     155 
     156a) for 3.6 version 
     157 
     158 * agree on all the points: bathymetry smoothing, update of coarse grids (take the maximum as default) 
     159 
     160 * Incorporate corrections done by Jean Marc 
     161 
     162 * The nesting tools are only essential for the creation of bathymetry: 
     163 
     164o   Do a list available tools for 2D/3D interpolation 
     165 
     166o   For 3.6 keep the 3D fields interpolation in nesting tools 
    199167 
    200168b) next version 
    201169 
    202   
    203  
    204 -       The definition of vertical coordinates should be done outside NEMO 
    205  
    206 -       Nesting tools for the s coordinates 
    207  
    208   
     170 * The definition of vertical coordinates should be done outside NEMO 
     171 
     172 * Nesting tools for the s coordinates 
    209173 
    2101745) Documentation: 
    211175 
    212 -       Update of the nesting tools 3.6 documentation 
    213  
    214 -       Update of the users guide: describes what is possible or not with AGRIF. 
    215  
    216 -       Documentation for nesting in NEMO/AGRIF: starts from the notes of Sebastien, then circulate the powerpoint 
    217  
    218 -       GMD paper 
    219  
    220   
     176 * Update of the nesting tools 3.6 documentation 
     177 
     178 * Update of the users guide: describes what is possible or not with AGRIF. 
     179 
     180 * Documentation for nesting in NEMO/AGRIF: starts from the notes of Sebastien, then circulate the powerpoint 
     181 
     182 * GMD paper 
    221183 
    2221846) Standard configuration and sanity tests 
    223185 
    224   
    225  
    226186The standard configuration should: 
    227187 
    228 -       Include more than one level of embedding (test of the recursive update starting from the finer resolution grid(s)) 
    229  
    230 -       Uses time splitting 
    231  
    232 -       Include a zoom in a sea ice area (Labrador Sea, LIM3) 
    233  
    234 -       Include a zoom with PISCES 
    235  
    236 -       Zoom in Perou Chili 
    237  
    238   
     188 * Include more than one level of embedding (test of the recursive update starting from the finer resolution grid(s)) 
     189 
     190 * Uses time splitting 
     191 
     192 * Include a zoom in a sea ice area (Labrador Sea, LIM3) 
     193 
     194 * Include a zoom with PISCES 
     195 
     196 * Zoom in Perou Chili 
    239197 
    240198Sanitiy tests 
    241199 
    242 -       Zoom with a mesh refinement of 1 
    243  
    244 -       Test constancy preservation for tracers 
    245  
    246   
     200 * Zoom with a mesh refinement of 1 
     201 
     202 * Test constancy preservation for tracers 
    247203 
    248204Longer visions of having AGRIF zoom in AMM