| 4 |
Notes on tags used in MITgcmUV |
Notes on tags used in MITgcmUV |
| 5 |
============================== |
============================== |
| 6 |
|
|
| 7 |
|
o pkg/mdsio: use MAX_LEN_FNAM (instead of hard coded 80) for file-names |
| 8 |
|
fix the singleCpuIO RS version of MDSWRITEFIELD |
| 9 |
|
and add a WARNING if using globalFiles=T with multi-processors. |
| 10 |
|
o write grid-file(s) in only one format (MDS or MNC): default=MNC (if |
| 11 |
|
using MNC); can force MDS output by setting writegrid_mnc=F in data.mnc |
| 12 |
|
o tracer advection scheme: |
| 13 |
|
add DST-2 (=Lax-Wendroff) & 1rst order upwind schemes (the way to get |
| 14 |
|
those 2 schemes before, changing GAD_FLUX_LIMITER.h + to recompile, |
| 15 |
|
was not very convenient). |
| 16 |
|
add implicit version for DST-3 & DST3_Flux_Limit vertical advection; |
| 17 |
|
(further tests needed for the implicit DST & Flux-Limiter schemes). |
| 18 |
|
o PTRACERS: added namelist parameter PTRACERS_ref in analogy to tRef |
| 19 |
|
and sRef for convenience, handy, when you want to initialize |
| 20 |
|
ptracers in a large domain with a constant non-zero value. Also, |
| 21 |
|
untangle ptracers_init.F a little. |
| 22 |
|
|
| 23 |
|
checkpoint57v_post |
| 24 |
|
o change dome exp. to use variable viscosity in side-drag terms. |
| 25 |
|
o remove all the files from the cfc_offline test since they are now |
| 26 |
|
incorporated into tutorial_cfc_offline which has been producing |
| 27 |
|
identical results for a few days of testing |
| 28 |
|
o add new parameter: sideDragFactor to enable half-slip-side BC |
| 29 |
|
and change side-drag S/R to use viscosity from mom_calc_visc |
| 30 |
|
(can switch back by setting sideDragFactor=0.); affects results at |
| 31 |
|
truncation level; regenerate output for exp.: global_ocean.cs32x10(x3), |
| 32 |
|
global_ocean_ebm, global_ocean_pressure, matrix_example & MLAdjust(x5). |
| 33 |
|
o vorticity & strain are now computed without masking, to account for |
| 34 |
|
no-slip-side BC in variable viscosity calculation ; mask is applied |
| 35 |
|
anyway later on, before computing advection/coriolis terms. |
| 36 |
|
o modify exp4 to test new obcs support for reading from external files |
| 37 |
|
and for passive tracers |
| 38 |
|
o fix behavior of ptracers_monitor when monitor output is written to a |
| 39 |
|
netcdf-file; still not independent of monitor.F (writes in the same |
| 40 |
|
file). |
| 41 |
|
o change default value of tRef (put a uniform, meaningful value for ocean |
| 42 |
|
& atmos) and initialize sRef (was missing, as Martin mentioned). |
| 43 |
|
o Redefine viscC?Leith*. Change the coefficients of Leith viscosity |
| 44 |
|
to be proportional to enstrophy dissipation length. Thus, |
| 45 |
|
viscC?leith?_new=pi*(viscC?leith?_old)^(1/3) |
| 46 |
|
This is now consistent with Smagorinsky and with the manual. |
| 47 |
|
|
| 48 |
|
checkpoint57u_post |
| 49 |
|
o OBCS and PTRACERS: add open boundary support for passive tracers |
| 50 |
|
- either use homogenous (pseudo) v.Neumann conditions or prescribe |
| 51 |
|
OB-values from file; this is not different from the way theta and salinity |
| 52 |
|
are treated |
| 53 |
|
- however, Orlanski-radiation conditions are not supported, and the model |
| 54 |
|
will stop if you use pTracers and Orlanski at the same time. |
| 55 |
|
- beefed up the rountine obcs_external_fields_load: now only those open |
| 56 |
|
boundary values are overwritten with values from files for which there |
| 57 |
|
are really files, otherwise the OB-fields remain untouched. This makes |
| 58 |
|
it possible to use different OBs at different ends of the domain (as |
| 59 |
|
with EXF) |
| 60 |
|
o reduce number of exch calls in shap_filt (s2). |
| 61 |
|
o CS-grid: fix divergence contribution to Leith viscosity. Affects (@ |
| 62 |
|
truncation level) 4/5 MLAdjust test results (11&12 digits match). |
| 63 |
|
|
| 64 |
|
checkpoint57t_post |
| 65 |
|
o test biharmonic viscosity + no-slip with exp: global_ocean.cs32x15.viscA4 |
| 66 |
|
o finish adding brackets (coriolis S/R, vecinv) in multi-terms sum to get |
| 67 |
|
same truncation, independent of the face number: affects affects several |
| 68 |
|
test. exp.: Update global_ocean.cs32x15 (+.thsice) & MLAdjust (all 5). |
| 69 |
|
o fix bug in mom_v_sidedrag (biharmonic viscosity was counted 2 times); |
| 70 |
|
note that biharmonic viscosity + no-slip was not tested in any experiment. |
| 71 |
|
o fix bug in non-linear free-surf side-drag ; generate new output for |
| 72 |
|
global_ocean.cs32x15(+.thsice) & global_ocean_pressure. |
| 73 |
|
o make mnc honor the writeBinaryPrec flag for most kinds of output |
| 74 |
|
o eesupp: split EEPARAMS into EESIZE/EEPARAMS and add a vector global sum |
| 75 |
|
o minor change in averaging W (mom_vi_?_vertshear.F) that affects truncation |
| 76 |
|
error; update output from global_ocean.cs32x15(+.thsice) & MLAdjust (all 5). |
| 77 |
|
o add brackets in multi-terms sum to get same truncation, independent of |
| 78 |
|
the face number (manage to get exactly same value of gU,gV at the cs-edge) |
| 79 |
|
affects several test. exp. Update global_ocean.cs32x15 (+.thsice). |
| 80 |
|
o Clean, new version of mom_calc_visc. Adds a new feature, viscAhRemax |
| 81 |
|
and viscA4Remax, which allow one to limit the minimum viscosity by specifying |
| 82 |
|
a maximum grid-scale Reynolds number. Also, changes all length calculations |
| 83 |
|
in the case of grid anisotropy to geometric mean of dx**2 and dy**2. |
| 84 |
|
Suggested values of viscosity parameters found in mom_calc_visc. |
| 85 |
|
Many new viscosity diagnostics available. |
| 86 |
|
|
| 87 |
|
checkpoint57s_post |
| 88 |
|
o Add clearing-house for viscosity calculations, mom_calc_visc and |
| 89 |
|
correct numerous bugs in LeithD and Smagorinsky discretizations. |
| 90 |
|
mom_calc_visc now allows for Smagorinsky, Leith and LeithD to be |
| 91 |
|
used simultaneously, and to be used in mom_vecinv and mom_fluxform. |
| 92 |
|
Mom_vecinv can be used with standard viscous calculation (mom_vi_hdissip.F) |
| 93 |
|
and with Strain-Tension viscous term calculation (mom_hdissip.F). |
| 94 |
|
The Strain-Tension form is now called with useStrainTensonVisc=.true. |
| 95 |
|
ViscAhStrain and viscAhTension are now retired. |
| 96 |
|
o diagnostics "VISCAHD" & "VISCA4D" replace "VISCAH" & "VISCA4". |
| 97 |
o In response to JMC's bug report, genmake2 now ignores any soft-links |
o In response to JMC's bug report, genmake2 now ignores any soft-links |
| 98 |
that it finds in the *current* directory when searching for files. |
that it finds in the *current* directory when searching for files. |
| 99 |
Soft-links found in any other source directories will still work as |
Soft-links found in any other source directories will still work as |
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch