--- MITgcm/doc/tag-index 2005/12/31 01:21:41 1.596 +++ MITgcm/doc/tag-index 2006/03/20 14:49:24 1.635 @@ -1,9 +1,162 @@ -$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.596 2005/12/31 01:21:41 jmc Exp $ +$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.635 2006/03/20 14:49:24 jmc Exp $ $Name: $ Notes on tags used in MITgcmUV ============================== +o add S/R to compute NH-metrics terms & Coriolis f' in vertical momentum + (were missing in NH mode, and wrong in quasihydrostatic mode). + waiting for Chris updated version of calc_gw to add those CALLs. +o fix bug in quasihydrostatic S/R ; + add a test for this option: global_ocean.90x40x15.qshyd +o new S/R "SET_PARMS": set parameters that might depend on the use of some pkgs. + (called from INITIALISE_FIXED, after INI_PARMS & PACKAGES_READPARAMS) +o add on/off flag useOffLine for package OFFLINE +o seaice: + - modify advection so that it works also on the C-grid + - extend loop ranges in seaice_lsr and seaice_dynsolver; fixes + problem in multi-tile runs with C-grid + - add flags for clipping ice velocities and/or masking solver forcing + - fix drag (ocean on ice) formulation and a bug that caused noise in the + C-grid solver + - added a few comments to seaice_diffusion.F + - separate C-grid and B-grid versions more cleanly (UVM and seaiceMaskU/V + are now exclusive) + - add a new ocean-ice stress coupling method (Hibler and Bryan, 1987) which + required re-ordering the code a little: + + make DAIRN global variable (defined in SEAICE) + + move computation of the viscosities eta and zeta into a new separate + routine (because I need to recompute them in seaice_ocean_stress.F) + +checkpoint58b_post +o AIM+thSIce: fix sign of Short-Wave through the Ice ; + update output.txt of exp. aim.5l_cs.thSI +o various mnc cleanups and improvements: + + shrink lookup tables by factor of ~4 + + better error reporting when running out of lookup space + + able to handle longer path/file names (up to 500 chars) +o adjoint: + - Another overhaul of store dirs. for NLFS to eliminate "hidden" recomputations. + - TBD: "hidden" mom_vecinv recomp. in dynamics +o change forcing_In_AB to affects both T,S forcing and Momentum forcing + (allow to differentiate between forcing components using new integer flags: + momForcingOutAB=1/0 & tracForcingOutAB=1/0) + and add new flag to put Dissipation tendency out of Adams-Bashforth. +o ecco + - enable trend removal from SSH anomaly obs. +o adjoint + - a fix for PmEpR when #undef EXACT_CONSERV + - fix of "hidden" recomp. in thermodynamics (no TAF warnings) + - try to fix "hidden" recomp. in forward_step for NLFS + (NLFS seems broken sometimes after I checked it in) +o seaice: + - add c-grid version of the dynamic solver: + + three new routines that are the c-grid substitute lsr, dynsolver + and ostres + + put a few fields that were local to dynsolver into global common + blocks, so that I can move the computation of stresses etc into + seaice_lsr (saves coding but may break the adjoint). + - replace more hardwired parameters by runtime parameters + - add ice masks that mask the rhs of the implicit solvers where there + is no ice (commented out in seaice_dynsolver, because i am not sure + if this works properly), eventually this should replace the clipping + of ice velocities in seaice_dynsolver to +/-40cm/s. +o remove the BUILD_INFO.h include from CPP_EEOPTIONS.h and explicitly + include it only where its needed +o move balancing of fluxes before the call of external_forcing_surf +o seaice: + - correct units of turning angles (now they are really in degrees) + - move computation of GWATX/Y out of main routine into dynsolver + - initialize GWATX/Y in seaice_init +o Apply mask on Fresh-Water flux (needed for SSH forcing). + update output.txt of exp. global_ocean_ebm & cpl_aim+ocn (affects cg2d residual) +o seaice: add runtime flags for turning angles (default to zero) +o Adding runtime parameter to exf/, ctrl/, to remove on-the-fly global mean + and trend from exf forcing and ctrl control adjustments +o fix variable aliasing in gad_calc_rhs that cannot be handled by TAF + (since checkpoint57x_post) +o change the default rotation period to 86164.s = the Earth rotation +o 1rst implementation of Implicit IGW using the 3-D solver (use3Dsolver=T) + and based on a reference stratification. Tested within hs94.cs-32x32x5. +o seaice advection: add code for choosing between the default (2nd order + central differences) and multiDim advection schemes provided by + pkg/generic_advdiff +o seaice: + - removed all references to SEAICE_GRID.h fields from + dynsolver, advect, diffus, lsr, without changing the results of lab_sea + (cubed-sphere configuration will change a little) + - grid computations in seaice_init are obsolete now. + - cleaned up code a little so that is easier to read (hopefully). + - added new (better?) discretizations of default advection and diffusion + (turn on with SEAICEuseFluxForm) + - replace 1 in 3d global fields by kSurface in growth.F + +checkpoint58a_post +o seaice: + - bug fix in advect.F +o shelfice: + - separate shelfice load anomaly from pload (breaks with time + dependent forcing), => introduce constant field shelficeLoadAnomaly. + Its default is 0., but it may be computed more cleverly from (unknown) + t- and s-profiles (tRef, sRef) and actual EOS. For now this has to be + done offline. A good approximation of the pressure load anomaly is + necessary to avoid large initial adjustment processes underneath + deep-reaching shelfice. + - small fix in ISOMIP shelfice_thermodynamics + - add more sophisticated thermodynamics of Hellmer (1989), this is + now the default + - give a few constants more intuitive names +o add code to balance EmPmP and Qnet at the end of do_ocean_physics. Useful + if bulk formulae are used in long integration (especially EmPmR). Turn + on with balanceEmPmR = .true. or balanceQnet = .true. in data, PARM01 + if balancePrintMean, the imbalance that is substracted is printed to + STDOUT. +o fix bug (typo) in obcs_apply_eta.F +o add a new package shelfice which enables us to model the circulation + underneath shelf-ice with draughts of hundreds of meters. + - works + - to do: + include improved thermodynamics + + check whether kSurfC (from model/inc/SURFACE.h) can be used + instead of new field kTopC + + can we use topoFile, Ro_surf and existing code in ini_depth.F + instead of new field R_shelfIce and new code in ini_depth.F + and ini_mask_etc.F? + + adjust vertical mixing schemes? +o seaice pkg: + - add a cpp-flag to include a simple formula for a salinity dependent + freezing point of sea water at the surface + - change the contents of SEAICE_TEST_ICE_STRESS_1 so that zonal and + meridional stresses are now different (and hopefully correct) +o AIM pkg: change threshold on large-scale condensation to conserve Heat. + update results of aim.5l_LatLon test ; fix title & units of AIM diagnostics + and change 3 diagnostics names (OSR,OSRCLR,KM -> TSR,TSRCLR,DRAG) +o pkg/diagnostics: add options and parameters for regional-statistics + diagnostics ; define regions by reading a region-mask from a file + (might add other options later) +o pkg/bulk_force: + - add option and parameters to use AIM surface-flux formulae. + - Change loading part: S/R BULKF_FIELDS_LOAD only take care of bulkf_fields ; + others forcing fields (defined in PARM05, parameter file "data") are loaded + from S/R EXTERNAL_FIELDS_LOAD, as usually (with useBulkforce = T or F). + - initialise all bulkf_fields in bulkf_init.F ; do in-lining of exf_bulkcdn.F + - use the right EXCH call for uwind,vwind (to work on CS-grid). + - re-arrange header files (move parameters from BULKF.h to BULKF_PARAMS.h) + and parameters (note: calcWindStress replaces .NOT.readwindstress). + - cleaner interface with pkg/thSIce (new S/R THSICE_GET_PRECIP). + +checkpoint58 +o add "--embed-source" option to genmake2 which, when the supporting + tools can be compiled (as determined by genmake2) will turn on the + embed_files package which then embeds the entire MITgcm source code + (*.[fFhc] + Makefile) used for the build within the executable + - requested by CNH but off by default + - adds a paltry <9MB to the mitgcmuv executable in most cases + - only writes the output when useEMBEDSRC in eedata is true +o clean up 3-dim. relaxation codes: + * Fix I/O inconsistency in pkg/rbcs: replace precFloat32 by readBinaryPrec + * Remove 3-dim. relaxation code from pkg/exf (now use only pkg/rbcs) + * Thanks to Tom Haine for testing! + checkpoint57z_post o make CG3D solver compatible with Free-surface at k > 1 (p-coordinate): affects output (@ truncation level) of exp. plume_on_slope ; updated.