--- MITgcm/doc/tag-index 2004/12/22 07:27:07 1.427 +++ MITgcm/doc/tag-index 2006/03/21 13:00:33 1.637 @@ -1,14 +1,631 @@ -$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.427 2004/12/22 07:27:07 dimitri Exp $ +$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.637 2006/03/21 13:00:33 mlosch Exp $ +$Name: $ Notes on tags used in MITgcmUV ============================== +o seaice: add an EVP solver following Hunke and Dukowicz (1997) and the + documentation of CICE. Turn on by defining SEAICE_ALLOW_EVP in + SEAICE_OPTIONS.h and SEAICEuseEVP=.true. in data.seaice. Works only + with SEAICE_CGRID defined. +o move calendarDumps from "data" to "data.cal" and clean-up the code + with a simple call to pkg/cal S/R: CAL_TIME2WRITE +o add useEXF & useCAL flags (for now, set in hard-coded way) +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. +o fix the pickup behavior for obcs in combination with ptracers +o fix a bug in ptracers_external_fields_load +o no AB-2 at the first iteration for Gw (consistent with Gu,Gv,Gt,Gs ...) + affects results of NH experiments (but all our NH tests have Gw[it=1]==0 !!) +o transfered surface NH pressure to eta field (if exactConserv). +o gmredi: fix bugs in slope tapering ldd97 (depthZ had the wrong sign) + and update lab_sea output.txt file. +o completing 3-dim. relaxation code: additional term + 1/tau*( theta(i,j,k) - thetaStar(i,j,k) ) + in S/R EXTERNAL_FORCING_T (file external_forcing.F) + * At compile time, in "EXF_OPTIONS.h" set + #define ALLOW_CLIMTEMP_RELAXATION + #define ALLOW_CLIMSALT_RELAXATION + * At runtime in "data", under &PARM03 set (e.g., units in seconds) + tauThetaClimRelax3Dim = 15552000.0, + tauSaltClimRelax3Dim = 15552000.0, + * 3 cases: + 1) Use a constant-in-time 3-dim. field + 2) Provide 12 3-dim. fields corresponding to Jan., ... , Dec. + 3) Real-time relax. to time-evolving prescribed field +o Adding Laplacian-type smoothness constraint on + initial and boundary (forcing) controls (ic, bc) + * enable at compile-time via + ALLOW_SMOOTH_IC_COST_CONTRIBUTION, ALLOW_SMOOTH_BC_COST_CONTRIBUTION + * enable at run-time via mult_smooth_ic = 1, mult_smooth_bc = 1 + * Absolute scale somewhat unclear so far + * no area weighting was included since no other cost term does so far + +checkpoint57y_post +First step for a NLFS adjoint +o initially suppress rStar (new flag DISABLE_RSTAR_CODE) +o new init. routines for calc_r_star, calc_surf_dr +o still need to deal with ini_masks_etc +o testreport seemed happy +checkpoint57y_pre + +o obcs: always reset wVel & eta (to zero or prescribe it) at the OBs. +o add new package rbcs (for relaxing ptracers at boundaries/depth/surface etc) +o add new offset feature to offline pkg (plus minor cleanup) +o cfc fixes for final year of integration +o clean up GCHEM_OPTIONS.h/DIC_OPTIONS.h (all options should + be in GCHEM_OPTIONS). + +checkpoint57x_post +o fix testreport Pb in multiple test experiment. +o cleanup our namelist syntax in verification: + - make it more F95 standards compliant (some work still needed) + - allows the vast majority of our tests to work with gfortran +o Non-Hydrostatic: + - put all NH variables (formerly in DYNVARS.h & GW.h) in NH_VARS.h + - fix exchange calls for CS-grid using the new EXCH_S3D_RL. +o Allow to apply AB on T,S rather than on AB(gT,gS): + - implemented within #ifdef ALLOW_ADAMSBASHFORTH_3 + - use the same arrays (gtNm,gsNm) to hold tracer field at previous + time-steps (if AB(T,S)) and tendencies (if AB(gT,gS)). + - need a "clever" read_pickup to restart from the other AB formulation +o move local commom bloc /GlobalLo/ (used for singleCpuIO) in EESUPPORT.h +o remove some unused variables (reduces number of compiler warnings) +o adjoint + * Disable useVariableViscosity for adjoint for time being + (there are a number of extensive recomputations) + * Initialize certain variables needed for adjoint + +checkpoint57w_post +o pkg/mdsio: use MAX_LEN_FNAM (instead of hard coded 80) for file-names + fix the singleCpuIO RS version of MDSWRITEFIELD + and add a WARNING if using globalFiles=T with multi-processors. +o write grid-file(s) in only one format (MDS or MNC): default=MNC (if + using MNC); can force MDS output by setting writegrid_mnc=F in data.mnc +o tracer advection scheme: + add DST-2 (=Lax-Wendroff) & 1rst order upwind schemes (the way to get + those 2 schemes before, changing GAD_FLUX_LIMITER.h + to recompile, + was not very convenient). + add implicit version for DST-3 & DST3_Flux_Limit vertical advection; + (further tests needed for the implicit DST & Flux-Limiter schemes). +o PTRACERS: added namelist parameter PTRACERS_ref in analogy to tRef + and sRef for convenience, handy, when you want to initialize + ptracers in a large domain with a constant non-zero value. Also, + untangle ptracers_init.F a little. + +checkpoint57v_post +o change dome exp. to use variable viscosity in side-drag terms. +o remove all the files from the cfc_offline test since they are now + incorporated into tutorial_cfc_offline which has been producing + identical results for a few days of testing +o add new parameter: sideDragFactor to enable half-slip-side BC + and change side-drag S/R to use viscosity from mom_calc_visc + (can switch back by setting sideDragFactor=0.); affects results at + truncation level; regenerate output for exp.: global_ocean.cs32x10(x3), + global_ocean_ebm, global_ocean_pressure, matrix_example & MLAdjust(x5). +o vorticity & strain are now computed without masking, to account for + no-slip-side BC in variable viscosity calculation ; mask is applied + anyway later on, before computing advection/coriolis terms. +o modify exp4 to test new obcs support for reading from external files + and for passive tracers +o fix behavior of ptracers_monitor when monitor output is written to a + netcdf-file; still not independent of monitor.F (writes in the same + file). +o change default value of tRef (put a uniform, meaningful value for ocean + & atmos) and initialize sRef (was missing, as Martin mentioned). +o Redefine viscC?Leith*. Change the coefficients of Leith viscosity + to be proportional to enstrophy dissipation length. Thus, + viscC?leith?_new=pi*(viscC?leith?_old)^(1/3) + This is now consistent with Smagorinsky and with the manual. + +checkpoint57u_post +o OBCS and PTRACERS: add open boundary support for passive tracers + - either use homogenous (pseudo) v.Neumann conditions or prescribe + OB-values from file; this is not different from the way theta and salinity + are treated + - however, Orlanski-radiation conditions are not supported, and the model + will stop if you use pTracers and Orlanski at the same time. + - beefed up the rountine obcs_external_fields_load: now only those open + boundary values are overwritten with values from files for which there + are really files, otherwise the OB-fields remain untouched. This makes + it possible to use different OBs at different ends of the domain (as + with EXF) +o reduce number of exch calls in shap_filt (s2). +o CS-grid: fix divergence contribution to Leith viscosity. Affects (@ + truncation level) 4/5 MLAdjust test results (11&12 digits match). + +checkpoint57t_post +o test biharmonic viscosity + no-slip with exp: global_ocean.cs32x15.viscA4 +o finish adding brackets (coriolis S/R, vecinv) in multi-terms sum to get + same truncation, independent of the face number: affects affects several + test. exp.: Update global_ocean.cs32x15 (+.thsice) & MLAdjust (all 5). +o fix bug in mom_v_sidedrag (biharmonic viscosity was counted 2 times); + note that biharmonic viscosity + no-slip was not tested in any experiment. +o fix bug in non-linear free-surf side-drag ; generate new output for + global_ocean.cs32x15(+.thsice) & global_ocean_pressure. +o make mnc honor the writeBinaryPrec flag for most kinds of output +o eesupp: split EEPARAMS into EESIZE/EEPARAMS and add a vector global sum +o minor change in averaging W (mom_vi_?_vertshear.F) that affects truncation + error; update output from global_ocean.cs32x15(+.thsice) & MLAdjust (all 5). +o add brackets in multi-terms sum to get same truncation, independent of + the face number (manage to get exactly same value of gU,gV at the cs-edge) + affects several test. exp. Update global_ocean.cs32x15 (+.thsice). +o Clean, new version of mom_calc_visc. Adds a new feature, viscAhRemax + and viscA4Remax, which allow one to limit the minimum viscosity by specifying + a maximum grid-scale Reynolds number. Also, changes all length calculations + in the case of grid anisotropy to geometric mean of dx**2 and dy**2. + Suggested values of viscosity parameters found in mom_calc_visc. + Many new viscosity diagnostics available. + +checkpoint57s_post +o Add clearing-house for viscosity calculations, mom_calc_visc and + correct numerous bugs in LeithD and Smagorinsky discretizations. + mom_calc_visc now allows for Smagorinsky, Leith and LeithD to be + used simultaneously, and to be used in mom_vecinv and mom_fluxform. + Mom_vecinv can be used with standard viscous calculation (mom_vi_hdissip.F) + and with Strain-Tension viscous term calculation (mom_hdissip.F). + The Strain-Tension form is now called with useStrainTensonVisc=.true. + ViscAhStrain and viscAhTension are now retired. +o diagnostics "VISCAHD" & "VISCA4D" replace "VISCAH" & "VISCA4". +o In response to JMC's bug report, genmake2 now ignores any soft-links + that it finds in the *current* directory when searching for files. + Soft-links found in any other source directories will still work as + they did before. +o add HAVE_SETRLSTK define and useSETRLSTK flag (eedata) which calls + a C routine to unlimit the stack size +o mnc-ify pkg/aim snapshot vars +o mnc-ify pkg/land as requested by Daniel +o various changes to mnc including: + - all files use the new "BASENAME[[.ITER].{t|f}NUM].nc" format + - output can now be grouped so that all files within a group + change the ITER portion of their names in lock-step together + - can now read ("global") PER-FACE (in addition to PER-TILE) files + and works with both EXCH1 and EXCH2 (but needs more testing) + - writing works for all verification test cases w/ g77 on Linux + +checkpoint57r_post +o make MNC variable initialization happen earlier +o add diagnostics for (almost) each term of the momentum Eq. +o remove SST,SSS,UVEL_k2,VVEL_k2 diagnostics. +o move pkg timeave out of pkg-groups "oceanic" & "atmospheric" (pkg/pkg_groups) +o remove MINIMAL_TAVE_OUTPUT option ; + add "dumpInitAndLast" to replace "nodumps" (and MINIMAL_TAVE_OUTPUT). +o move some basic S/R (used elsewhere in the code) from pkg/mdsio to eesupp/src + allows to compile without pkg/mdsio (& test successfully hs94.1x64x5) +o add MNC output for adjoint of EXF (adexf), SEAICE (adseaice) +o add time-dependent SST, SSS control +o cost terms in ecco/ modularized (new cost_generic.F) +o cost_averages routines are modularized (new cost_averagesgeneric.F) +o new infrastructure for seaice model vs. obs misfits in seaice/ + * driver from the_main_loop.F is seaice_cost_driver.F + * first term is AREA vs. SMR observation misfit (ALLOW_SEAICE_COST_SMR_AREA) + * gradient check in Exp. 6a of natl_box_adjoint/ +o natl_box_adjoint input_seaice for Lab. Sea now available online via wget + (see README file) + +checkpoint57q_post +o fix initialisation Pb in calc_gw (get NANs in exp5 with g77) +o fix annoying months-old bug in mnc/mdsio output for seaice pkg + - bug was reported by Jean-Michel and Ian Fenty +o fix initialization Pb in AIM-Ocean coupled model (aim_sice2aim.F) + update exp. cpl_aim+ocn output. +o Implement 4-level checkpointing: AUTODIFF_4_LEVEL_CHECKPOINT + Needed for high-res. runs with large global files + to be able to limit tapelev size +o Changed some tests for adjoint: + * carbon/ now tests AUTODIFF_2_LEVEL_CHECKPOINT + * natl_box_adjoint/code_bulk/ tests AUTODIFF_4_LEVEL_CHECKPOINT + * global_ocean.90x40x15/ back to default 3-level checkpointing +o Modifications to enable scripted COMMON to MODULE conversion + * verification/OpenAD/code/ has necessary scripts and description; + * converts, compiles and runs successfully; + +checkpoint57p_post +o reading of delXfile & delYfile moved from ini_parms.F to ini_grid.F +o move ${PKG}_MNC_INIT from ${PKG}_READ_PARAMS to ${PKG}_INIT_FIXED + (already the case for some pkgs, including recent MNC_init, e.g. thsice) +o call CALC_GW from DYNAMICS (instead of from FORWARD_STEP) +o Dissipation & phiHyd gradient are always added to gU,gV in timestep.F + (was already the case for dissipation with mom_vecinv, + and also the case for grad.PhiHyd when staggered-timeStep) + This will allow to put dissipation out-off the AB time-stepping. + Unfortunately, affects truncation error. Update output of exp. that "fails" + (cg2d matching number) and leave the others for later: front_relax(10), + global_ocean_pressure(10),ideal_2D_oce(10),internal_wave(9),lab_sea(12), + matrix_example(7!),natl_box(12),plume_on_slope(11),vermix(12x3,11x2). + +checkpoint57o_post +o fix MNC/MDSIO flags in PTRACERS so that (1) it works, and + (2) its ready to be documented during the 20050801 DocFest +o MNC-ify the gmredi package -- tested and works w/ and w/o mnc + for ideal_2D_oce and lab_sea on linux_ia32_g77 +o set to #undef CTRL_SET_OLD_MAXCVARS_30 + in global1x1_tot/code_seaice_adj to access indices 32, 33 + (worth another tag) + +checkpoint57n_post +o adding precip control xx_precip 32 +o adding swflux control xx_swflux 33 +o Replacing 10 forcing penalty cost routines by 1 generic. +o Implementing partitioning in time-mean and variable adjustment + for atmos. state. +o Time-mean constant error is now read in data.ecco (instead of hard-coded) +o in natl_box_adjoint: + * updated exp 1 to 10 + (but does not have time-mean vs. variable split results yet) +o updated global1x1/*_seaice_adj/ + to incorporate latest changes +o enable useSingleCpuIO for mdsio vector routines + (applies in particular to tapelev I/O) +o removed inefficient K-loop in all active_file*control*.F + +checkpoint57m_post +o re-arrange exp. adjustment.cs-32x32x1 (adjust_nlfs.cs-32x32x1 merge into) + and advect_cs (use EXCH-2 and standard cs-32 grid files). +o pkg/exch2: filling of face-corner halo regions is optional + (ifdef W2_FILL_NULL_REGIONS) and is turned off by default. + +checkpoint57l_post + tag before updating tutorials (pre-release-2) +o add cos & sin of grid-direction orientation angle (curvilinear grid only, + for output and interpolation of input fields onto model grid) +o fix radius in mom_v_metric_cylinder.F (wrong by dy/2). + Affects results of exp. rotating_tank ; update output.txt. +o pkg/diagnostics: 2 new S/R to fill a diagnostics using a scaling factor + and taking the square quantity and/or using a fraction-weight field. +o fix aim_diagnostics_init.F mate numbers. +o add mnc output capability to diagnostics/diagstat and update + our cvsignore files for the ACSII output generated + +checkpoint57k_post +o lab_sea now uses diagnostics package and provides a suite of + new EXF diagnostics output. + +checkpoint57j_post +o exf + * bug fix in exf_bulkformulae for psimh (spotted by Eliyahu Biton): + Bracket missing in eqn. for psimh (eqn. 7 of Large & Pond, 1982); + terms -2*ATAN(x) + pi/2 were added to both stable and unstable case. + * extended diagnostics output + - now includes all EXF variables + - hs, hl moved to common block + - separate S/R for diagnostics, called after ctrl vector is added + - name changed for EXFsw -> EXFswnet, EXFlw -> EXFlwnet + - EXFqnet contains true Qnet (i.e. includes qsw) + * updated verif. lab_sea, global_with_exf + * Rearranging hflux (exf_mapfields -> exf_getforcing) + to ensure that diagnostics sees a true EXFqnet +o add matlab scripts for CS-grid vectors (in utils/cs_grid) +o diagnostics pkg: change pointers so that 1 diagnostic can be used + several times with different frequencies. +o mnc-ify the thsice package as requested by Daniel Enderton + - the monitor--MNC output needs work + - many attributes need to be added (most are currently blank) +o mnc-ify theaim_v23 package as requested by Daniel Enderton + - many attributes need to be added (most are currently blank) +o remove "groundAtK1" & "rkFac" (replaced by usingPCoords & -rkSign) +o add Sadourny like discretization to the 4th.O vorticity advection scheme +o fix bug in rStar code (calc_r_star.F) when used with topo/bathy. + update output.txt from: aim.5l_cs (& .thSI), cpl_aim+ocn, + fizhi-gridalt-hs and global_ocean.cs32x15 (& .thsice). +o fix few little things to avoid floating error with -ftrapuv (ifort v9) + (ini_spherical_polar_grid.F, calc_r_star.F, monitor/mon_surfcor.F) + +checkpoint57i_post +o Diagnostics package was changed: diagnostic allocation is now deactivated + by setting an empty file name in the namelist (used to be by setting 0 frequency) +o momentum advection (VI form) options: activate 4th.O vorticity Coriolis ; + add 1rst Order upwind vertical Shear. +o add a quick fix for A/B checkpointing in OFFLINE mode +o add IGNORE_TIME define and -ignore_time option to genmake2 as an + end-run around the C/Fortran linking horror in eesup/src/timers.F +o added pkg/runclock: allows user to terminate run based on wall clock +o added non-hydrostatic scaling parameter: nh_am2 +o pkg/my82: futile attempt to make results compiler-independent by adding + some "_d 0", etc. no effect. +o for mnc output, fill the 'T' coordinate var with myTime and create a + separate 'iter' variable for iteration count +o pkg/diagnostics: + add new capability: output of Global/Regional & per-level statistics ; + tested against monitor output (including land & thSIce) on various exp. +o make seaice-pkg I/O type flags (MNC/MDS) more standard. + +checkpoint57h_done + tag after changing user customized files (DIAGNOSTICS_SIZE.h, + data.diagnostics, data.gchem, data.ptracers) is done. + +o pkg/diagnostics: prepare for Global/Regional statistics. + add 2nd namelist in data.diagnostics and update DIAGNOSTICS_SIZE.h +o remove "baseTime" (no used) from arg. list of DIFF_BASE_MULTIPLE + and rename it: DIFFERENT_MULTIPLE + +checkpoint57h_post +o move tIter0 from data.gchem to PTRACERS_Iter0 in data.ptracers +o pkg/diagnostics: change output frequency, now in seconds (instead + of Nb of iterations). This requires to change data.diagnostics. +checkpoint57h_pre + +o in mnc, let the user create the output dir if both mnc_outdir_date + and mnc_outdir_num are false +o add fflags() to pkg/diagnostics which allows one to specify, on a + per-file basis, the file precision [default behavior is unchanged] +o added diagnostics to pkg/kpp, including computation of mixed layer + depth based on a temperature/density criterion +o updated verification/natl_box to test the new pkg/kpp diagnostics +o bugs fixed (with useAbsVorticity=T & with useJamartMomAdv=T) +o test AB-3 code with 4th.Order advection scheme (advect_xy/input.ab3_c4) +o fix mnc flags in pkg/seaice so that they are "owned" by seaice, + not mnc -- this mistake was found by JMC and is now corrected + +checkpoint57g_post +o add Adams-Bashforth 3rd Order time stepping: + - default remains AB-2 (#undef ALLOW_ADAMSBASHFORTH_3) + - not yet implemented with NH, CD_Scheme, pTracers, MNC_pickup ... +checkpoint57g_pre + +o make lambda[Theta,Salt]ClimRelax spatially varying + (default remains spatially uniform constant based on tau) +o added new control variables + * init. uVel, vVel, etanN + * lambda[Theta,Salt]ClimRelax +o separate masks used for ctrl_pack/unpack 'from write_grid' output + (suggested by G. Forget) + +checkpoint57f_post +o add parameter "baseTime" = model base time = time at iteration zero + new function DIFF_BASE_MULTIPLE (replace DIFFERENT_MULTIPLE) to + account for baseTime shift +o move Coriolis arrays from PARAMS.h to GRID.h +checkpoint57f_pre + +o mnc-ifying the seaice package +o introduced an alternative length scale calculation for use in the limiting + maximum viscosity (off by default, see mom_vi_hdissip.F for details) +o S/R ptracers_calc_diff: argument list bug for CALL KPP_CALC_DIFF_S +o added modified Leith scheme with viscosity proportional to grad(Hdiv) + as well as grad(vort3). Made it compatible with both approximate form + for the cubed sphere and full form of these gradients. Added flag + useFullLeith to select between approximate and full form. +o add "d0" in getcon.F to reduce differences between compilers/platforms + update output.txt of exp. fizhi-gridalt-hs +o change port_rand funtion to allow to set the seed. +o fix obscure MNC bug involving conversion between reals and integers + - it caused no problems with certain (eg. gcc, ifc) compilers but the + new version is now more standards-compliant +o adjoint-related modifications + * gad_dst3_adv_?: to be able to use advscheme=30 in adjoint + * gad_calc_rhs: fool adjoint, revert to advscheme=30 in reverse when using 33 + * mom_vi_hdissip: to be able to use Leith scheme in adjoint + +checkpoint57e_post +o Adding eddy stress controls a la Ferreira et al. +checkpoint57e_pre (actually was tagged: eckpoint57e_pre) + +o fizhi pkg: fix AtoC and CtoA when called with same input & output array + update output.txt of exp. fizhi-gridalt-hs +o add more "Z" dimensions and associated coordinate vars to MNC +o pkg/exf: added twoDigitYear capability for option useExfYearlyFields +o parameter calendarDumps: when set, approximate months (30-31 days) and years + (360-372 days) for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq, + KPP_taveFreq, and freq in pkg/diagnostics are converted to exact calendar + months and years. Requires pkg/cal. +o diagnostics: add pickup functionality + - off by default for backwards compatibility + - current version only uses MDSIO -- MNC will be added +o allow to use pkg diagnostics for LAND pkg output & Shapiro Tendencies. +o pkg/shap_filt: fix for multi-tiles per face (CS-grid with EXCH2) + +checkpoint57d_post +o pkg/thsice used in the ocean: fix tile Pb related to sIceLoad forcing +o modify verification/ideal_2D_oce so that it tests the creation of mnc + output from the diagnostics package +o temporary fix in EXCH2_UV: add a 3rd call to fix the 24.tiles CS-grid set-up +o make genmake2 more portable by using ".$FS" instead of ".f" for more + of the test compilations (eg. IBM xlf on Mac OSX w/HFS filesys) +o add coordinate variables to the MNC output from the diagnostics package +o change MNC-flags default value: true when useMNC=TRUE +o add a "readgrid_mnc" flag to MNC and give model/src/ini_curvilinear_grid.F + the ability (off by default) to read the grid information using MNC +o lab_sea: use LSR_ERROR=1.E-12 for testing ; update output.txt +o diagnostics: fix a bug (sub-set of levels); keep full precision when averaging +o do snap-shot diagnostics at the middle of the time-interval (instead of + at the end). +o update exp. aim.5l_cs (use default LAND_OPTIONS.h file) and add + new test for AIM with thSIce + slab-ocean +o cpl_aim+ocn test: update output (optfile & default GMREDI_OPTIONS.h) +o move state variable diagnostics to the beginning of the time step ; + (note: still needs some work specially for snap-shot diagnostics) + duplicate & replace UVELVVEL diagnostics with UV_VEL_C & UV_VEL_Z + +checkpoint57c_post +o Change testing option (-ieee) to get closer agreement between g77 & ifort +o generate new output.txt files on faulks with g77 & the new -ieee option +checkpoint57c_pre + +o add "_d 0" in freeze ; change some "FLOAT" to "DFLOAT" and went back (TAF); + update output.txt of exp. global_ocean.90x40x15 +o fix in ini_eos.F (one "_d 0" was missing): now g77 & ifc give much closer + results; update global_ocean_pressure & global_ocean.cs32x15 output.txt +o add saltFlux(+File) forcing capability with periodicExternalForcing +o allow to use pkg diagnostics for GM-Redi & ThSIce output +o diagnostics_fill: fill in up to sNx+1 or sNy+1 (if U,V or Z pt) for mnc. +o improvements for the MNC output from DIAGNOSTICS: + - correct handling of the unlimitied dimension + - all dims now have associated coordinate variables + - variables located at U, V, and Corner ("Z") points are now written + with the correct horzontal dimensions (eg. 'Xp1', 'Yp1') so that + *all* data points are output +o added seaice_summary.F and removed obsolete ALLOW_SEAICE's from pkg/seaice +o added pkg/seaice/seaice_monitor.F o mdsio_readfield.F can now read files visible only to master MPI process. o pkg/seaice: some test code for ice-modified ocean stress o exf_getffields interpolation defaults to bilinear for all scalar forcing fields (remains bicubic for wind velocity and stress). This avoids, e.g., spurious negative numbers for precipitation and humidity. Will cause some small numerical differences for integrations using pkg/exf/exf_interp.F. +o pkg/ctrl/ + * Add ctrlvec diagnostics in pack/unpack for nondimensional I/O + May be enabled via doPackDiag + * introduce z-scaling of + + gradient (doZscalePack) and + + control (doZscaleUnpack) +o pkg/mdsio/ + * enable useSingleCpuIO for state estimation (ctrl_pack/unpack) mdsio_gl.F + * N.B.: leave tiled I/O for tapelev (i.e. mds...vector) checkpoint57b_post o simplify options in diagnostics_fill arguments (remove weird combinations)