28 | 28 | |
29 | 29 | ''' (I.1)''' standardisation of the name of quantities derived from horizontal scale factors (many modules of NEMO)[[BR]]The chosen names of quantities derived from horizontal scale factors are:[[BR]]- inverse of scale factors : r1_e1t, r1_e2t, r1_e1u, r1_e2u, r1_e1v, r1_e2v, r1_e1f, r1_e2f ;[[BR]]- surfaces and their inverse: e1e2t, e1e2u , e1e2v, e1e2f and, r1_e1e2t , r1_e1e2u , r1_e1e2v, r1_e1e2f ;[[BR]]- ratio (used in traldf...) : e2_e1u , e1_e2v(I.1) generalize of the use of surface scale factors ;[[BR]]These quantities are defined in memory in dom_oce.F90, and set to their proper value in domhgr.F90 |
96 | | '''II.2 —''' '''Advective trends''' :TRA/traadv.... ===>>> OK see '''revision [https://forge.ipsl.jussieu.fr/nemo/log/branches/2015/dev_r5721_CNRS9_NOC3_LDF?rev=5737 5770] & [https://forge.ipsl.jussieu.fr/nemo/log/branches/2015/dev_r5721_CNRS9_NOC3_LDF?rev=5737 5771]''' |
97 | | |
98 | | - traadv_muscl2.F90 removed[[BR]]- traadv_muscl.F90 renamed traadv_mus.F90[[BR]]- traadv_cen2.F90 replaced by traadv_cen which provides 2nd and 4th order centered scheme on both horizontal and vertical directions.[[BR]]- name TVD has been changed to the more commonly used FCT (Flux Correction Transport)[[BR]]- FCT offers now 2nd and 4th order centered scheme on both horizontal and vertical directions.[[BR]]- UBS offers now in option either a 2nd order FCT or a 4th order COMPACT scheme in the vertical direction[[BR]]- NB in case of 4th order in the vertical, CEN, UBS and FCT use a compact 4th order scheme[[BR]]- key_esopa removed from traadv.F90 |
| 96 | '''II.2 —''' '''Advective trends''' :TRA/traadv.... ===>>> OK see '''revision [https://forge.ipsl.jussieu.fr/nemo/log/branches/2015/dev_r5721_CNRS9_NOC3_LDF?rev=5737 5770] & [https://forge.ipsl.jussieu.fr/nemo/log/branches/2015/dev_r5721_CNRS9_NOC3_LDF?rev=5737 5771] ===>> SETTE OK with revision [https://forge.ipsl.jussieu.fr/nemo/log/branches/2015/dev_r5721_CNRS9_NOC3_LDF?rev=5737 5774] & [https://forge.ipsl.jussieu.fr/nemo/log/branches/2015/dev_r5721_CNRS9_NOC3_LDF?rev=5737 5775]''' |
| 97 | |
| 98 | * traadv_muscl2.F90 removed[[BR]]- traadv_muscl.F90 renamed traadv_mus.F90[[BR]]- traadv_cen2.F90 replaced by traadv_cen which provides 2nd and 4th order centered scheme on both horizontal and vertical directions.[[BR]]- name TVD has been changed to the more commonly used FCT (Flux Correction Transport)[[BR]]- FCT offers now 2nd and 4th order centered scheme on both horizontal and vertical directions.[[BR]]- UBS offers now in option either a 2nd order FCT or a 4th order COMPACT scheme in the vertical direction[[BR]]- NB in case of 4th order in the vertical, CEN, UBS and FCT use a compact 4th order scheme[[BR]]- key_esopa removed from traadv.F90 |
| 127 | ''' III. Phasing of viscous trends''' |
| 128 | |
| 129 | '''III.1 —''' '''LDF on dynamics''' - viscous trends : ===>>> OK see '''revision .... and ..... with SETTE OK''' |
| 130 | |
| 131 | ''operator'': [[BR]]- laplacian operator: re-entrant laplacian for iso-level operators (dynldf_lap);[[BR]]- bilaplacian operator: dynldf_blp.F90 module: use of the re-entrant iso-level laplacian [[BR]]- suppression of traadv_eiv.F90 ; traldf_bilap.F90 ; traldf_bilapg.F90 ;[[BR]]- rotb and rotn arrays suppressed (even from the restart file) [[BR]]- hdivn and hdivb removed from restart too [[BR]]- divcur.F90 becomes divhor and only computes the horizontal divergence |
| 132 | |
| 133 | [[BR]]''lateral mixing coefficient'':[[BR]]- new user interface (namelist) : |
| 134 | |
| 135 | {{{ |
| 136 | !----------------------------------------------------------------------- |
| 137 | &namdyn_ldf ! lateral diffusion on momentum |
| 138 | !----------------------------------------------------------------------- |
| 139 | ! ! Type of the operator : |
| 140 | ! ! no diffusion: set ln_dynldf_lap=..._blp=F |
| 141 | ln_dynldf_lap = .false. ! laplacian operator |
| 142 | ln_dynldf_blp = .false. ! bilaplacian operator |
| 143 | ! ! Direction of action : |
| 144 | ln_dynldf_lev = .false. ! iso-level |
| 145 | ln_dynldf_hor = .false. ! horizontal (geopotential) |
| 146 | ln_dynldf_iso = .false. ! iso-neutral |
| 147 | ! ! Coefficient |
| 148 | nn_ahm_ijk_t = 0 ! space/time variation of eddy coef |
| 149 | ! ! =-30 read in eddy_viscosity_3D.nc file |
| 150 | ! ! =-20 read in eddy_viscosity_3D.nc file |
| 151 | ! ! = 0 constant |
| 152 | ! ! = 10 F(k)=c1d |
| 153 | ! ! = 20 F(i,j)=F(grid spacing)=c2d |
| 154 | ! ! = 30 F(i,j,k)=c2d*c1d |
| 155 | ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) |
| 156 | rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] |
| 157 | rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] |
| 158 | rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] |
| 159 | ! |
| 160 | ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) |
| 161 | / |
| 162 | }}} |
| 163 | [[BR]]- default value of the namelist: ''no operator chosen'': ln_dynldf_lap=dyndf_blp=0 No lateral diffusion applied on momentum[[BR]]- create the 2D and 3D ahtu and ahtv file for ORCA_R2 '''===>>> TO BE DONE for ORCA_R1''' . nn_ahm_ijk_t must be set to -20 or -30 for ORCA2 and ORCA1 (read in a file) [[BR]]- both laplacian and bi-laplacian coef. are always 3D arrays defined at T- and F- points. Bi-laplacian eddy coefficient defined as SQRT(bhm) '''====>>> TO BE DONE: explore the possibility of u- and v points'''[[BR]]- redesign of space variation of the eddy viscosity: suppression of: ldfdyn_substitute.h90 ; ldfdyn_c1d.h90 , ldfdyn_c2d.h90 and ldfdyn_c3d.h90 (replaced by ldfc1d_c2d.F90, a module share with ldftra.F90) ; ldfdyn_oce.F90 (now included in ldfdyn.F90)[[BR]]- remove ldfdyn_smag.F90 (suppression of smagorinsky option) '''==>>> TO BE re-add in the new framework by someone...'''[[BR]]- suppression of key_esopa in traldf.F90 |
| 164 | |
| 165 | [[BR]]'''IV. Solve the pending issues''' |
| 166 | |
| 167 | ''' (IV.1) ''' specification of eddy viscosity/diffusivity coefficients in case of grid-size dependency.[[BR]] Issue: define the grid-size dependency in such a way that does not require the Equator being inside the domain.[[BR]] Solution: add a scalar in the coordinate file which is equal to the grid-size at the equator even if Equator is not in the domain[[BR]] and use this scalar to compute the grid-size dependency of the coefficients. [[BR]]'''(IV.2) ''' specification of lateral friction at the coast.[[BR]] Issue: the different values of fmask at the coast have been removed so that this lateral boundary condition doesn't affect the advective term (''' double check that''' !). fmask is still used for laplacian and bilaplacian viscosities. Nevertheless, the magnitude of lateral friction depends on the magnitude of ahm (bhm, reap.). It seems better to define this magnitude [[BR]] Solution: add a scalar in the coordinate file which is equal to the grid-size at the equator even if Equator is not in the domain[[BR]] and use this scalar to compute the grid-size dependency of the coefficients. [[BR]]'''(IV.3)''' Add a proper coding of Smagorinsky eddy viscosity and diffusivity |
| 168 | |
| 169 | ... |
| 170 | |
| 171 | '''remaining to do :''' |
| 172 | |
| 173 | * '''IV.1 '''In case of variation of the 2D coeff. as a function of the grid size, use a reference to a 1 degree grid size for ah on both tracers and momentum |
| 174 | * '''IV.2''' management of no-slip boundary condition as fmask no more used in dynvor.F90 .... |
| 175 | * '''IV.3''' re-introduce in a consistent way with the re-design LDF if decided... |
128 | | |
129 | | : |
130 | | |
131 | | ''' III. Phasing of viscous trends''' |
132 | | |
133 | | ''' (III.1) LDF on dynamics''' |
134 | | |
135 | | * rotb and rotn arrays suppressed (even from the restart file) [[BR]]- hdivn and hdivb removed from restart too [[BR]]- divcur.F90 becomes divhor and only computes the horizontal divergence |
136 | | |
137 | | ... |
138 | | |
139 | | ''' IV. Solve the pending issues''' |
140 | | |
141 | | ''' (IV.1) ''' specification of eddy viscosity/diffusivity coefficients in case of grid-size dependency.[[BR]] Issue: define the grid-size dependency in such a way that does not require the Equator being inside the domain.[[BR]] Solution: add a scalar in the coordinate file which is equal to the grid-size at the equator even if Equator is not in the domain[[BR]] and use this scalar to compute the grid-size dependency of the coefficients. ''' (IV.2) ''' specification of lateral friction at the coast.[[BR]] Issue: the different values of fmask at the coast have been removed so that this lateral boundary condition doesn't affect the advective term (''' double check that''' !). fmask is still used for laplacian and bilaplacian viscosities. Nevertheless, the magnitude of lateral friction depends on the magnitude of ahm (bhm, reap.). It seems better to define this magnitude [[BR]] Solution: add a scalar in the coordinate file which is equal to the grid-size at the equator even if Equator is not in the domain[[BR]] and use this scalar to compute the grid-size dependency of the coefficients. |
142 | | |
143 | | ... |
144 | | |
145 | | '''remaining to do :''' |
146 | | |
147 | | * management of no-slip boundary condition as fmask no more used in dynvor.F90 .... |
148 | | * create the 2D ahtu and ahtv file for ORCA_R2 (and also for ORCA_R1) |
149 | | * In case of variation of the 2D coeff. as a function of the grid size, use a reference to a 1 degree grid size for ah on both tracers and momentum |
150 | | * suppress Smagorinsky option which does not work any more. Can be re-introduce in a consistent way with the re-design LDF if decided... |
153 | | ''' V. Finalize the simplification''' |
154 | | |
155 | | ''' (V.1) add''' in NEMOGCM/TOOLS a directory a program that read a coordinates.nc file apply a reduction of scale factors in some straits, and create the associated extended coordinates file (coordinates_e1e2u_v.nc).[[BR]] create the extended coordinates files form ORCA 2°, 1°, and 0.5°. ===>>> Done for ORCA 2° and 1°, ''' 1/2° and 1/4° missing''' ... |
| 179 | [[BR]] ''' V. Finalize the simplification''' |
| 180 | |
| 181 | ''' (V.1) add''' in NEMOGCM/TOOLS directory a program that read a coordinates.nc file apply a reduction of scale factors in some straits, and create the associated extended coordinates file (coordinates_e1e2u_v.nc).[[BR]]'''(V.2) create''' the extended coordinates files form ORCA 2°, 1°, and 0.5°. ===>>> Done for ORCA 2° and 1°, ''' 1/2° and 1/4° missing''' ...[[BR]]'''(V.3) create''' in NEMOGCM/TOOLS directory a program that create the ORCA2 and ORCA1 eddy viscosity''' ????''' |