Parent Directory
| 
Revision Log
| 
Revision Graph
| Links to HEAD: | (view) (download) (annotate) |
| Links to checkpoint65k: | (view) (download) (annotate) |
| Sticky Tag: |
Fix a few comments + small changes in the code (does not change results of global_ocean.gm_k3d)
For the k3d module in gmredi: 1/ Fixed bug in velocity rotation & conversion to A grid. 2/ Fixed bug in barotropic velocity calculation. 3/ Changed a couple of logicals.
- rename uMean,vMean --> uEulerMean,vEulerMean - fix arguments of few DIAGNOSTICS_FILL calls in gmredi_k3d.F
no needs to include CPP_OPTIONS.h (already there from GMREDI_OPTIONS.h)
for K3D parameterisation: fixed up minor bugs relating to finding the eigenvalues/vectors
for the K3D parameterisation: made sure that all variables are initialised properly (i.e. assigned some value everywhere
K3D: Added upper and lower bounds to renormalisation factor as namelist variables. Fixed small bug with interpolation to U & V grid. Also added extra documentation and a diagnostic.
fixed problem with floating point exceptions when calculating the renormalisation factor
K3d: Added option to keep diffusivity constant in the mixed layer. Also added some diagnostics.
1) cleaned up and renamed namelist parameters to better reflect their purpose, 2) added some new diagnostics 3) imposed a maximum length for urms length scale and 4) imposed a minimum length for eddy length scale
Fixed bugs in 1) Eady growth rate calculation for k3d 2) diagnostic description for residual model
1/Introduced namelist logical to turn the PV sheet on and off. 2/ Fixed bug when GM_K3D_likeGM=.TRUE. which gave the wrong diffusivity to the isoneutral diffusion tensor + added some helpful comments.
In K3d (for gmredi) 1/facility for a minimum bound on urms length scale introduced, 2/ extra diagnostics, 3/ changes to the Eady growth rate calculation
Changes for the residual model.
added an extra diagnostic to GM_K3D
avoid unbalanced single or double quote (CPP warnings)
1/ Fixed units in diagnostics of Xi. 2/ Made calculation of beta more robust 3/ Fixed bug in calculation Eady growth rate. The bug was related to the depth of integration.
Namelist option to ignore beta in the calculation of grad(q) for GM_K3D
Changes associated with the PV eddy closure (GM_K3D) include: o To improve efficiency the call to solve for the eigenvectors only happen once every GM_K3D_vecFreq seconds. This required the following changes: - read and write pickup files for the eigenvectors and deformation radius (gmredi_read_pickup.F and gmredi_write_pickup.F) - making the number of modes (GM_K3D_NModes) a parameter which must be specified at compile time in GMREDI.h - A new namelist variable, GM_K3D_vecFreq - Added modesC, modesW, modesS and Rdef to the common block o If the CPP option use_lapack is undefined, then a WKB approximation to the eigenvectors and deformation radius is now used (although, it seems unstable; so for the moment an error is raised in gmredi_check if GM_K3D is defined but use_lapack is not). o Changed gmredi_calc_eigs returns the deformation radius rather than the deformation wavenumber o Fixed bug in calculation of tfluxX and tfluxY for the instance where the surface layer is the depth of the water column. o Added warning messages if there are problems with calculating eigenmodes and eigenvectors o Cleaned up code o Improved documentation o Rationalised diagnostics o Added some extra startup checks (gmredi_check)
1/ use calc_oce_mxlayer to calculate the mixed layer depth for the PV eddy closure in gmredi_k3d. 2/ added smoothing option for constant K when using the PV closure (gmredi_k3d). 3/ added some extra checks for when GM_useK3D 4/ cleaned up gmredi_k3d a bit
added a debugging option to run the eddy PV closure in a GM limit (GM_K3D_likeGM)
Added a new eddy diffusivity parameterisation pkg/gmredi. More detailed description in tag-index.
This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, select a symbolic revision name using the selection box, or choose 'Use Text Field' and enter a numeric revision.
| ViewVC Help | |
| Powered by ViewVC 1.1.22 |