--- MITgcm/doc/tag-index 2002/11/12 20:59:27 1.103 +++ MITgcm/doc/tag-index 2003/02/26 03:29:08 1.137 @@ -1,8 +1,264 @@ -$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.103 2002/11/12 20:59:27 heimbach Exp $ +$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.137 2003/02/26 03:29:08 jmc Exp $ Notes on tags used in MITgcmUV ============================== +untagged +o improve "ini_p_ground" when using finite difference form for calc_phi_hyd +o Added net flux to downward flux conversion to pkg/exf/exf_getffields.F +o Added SEAICE_initialHEFF to pkg/seaice + +checkpoint48g_post +--- EOS: +o compute locally the pressure for use in EOS: UNESCO, JMD95P, MDJWF + store total Potential in totPhyHyd for diagnostics & EOS. +o fix restart and overlap Pb when using Z-coord and EOS funct. of P +--- Pressure/geopotential gradient: +o use Ro_surf & R_low instead of hFac in cal_phi_hyd: +o NLFS + p-coord.: eta*Alpha' is not longer lagging 1.time-step behind + ( changes results of global_ocean_pressure ) +o NLFS + z-coord.: add eta*Rho' contribution. +o accurate phiHyd when using interface-W at the middle between 2 centers + ( changes results of ideal_2D_oce ) +o includes r* 2nd term in potential gradient. + +checkpoint48f_post +Merging from release1_p12: +o Modifications for using pkg/exf with pkg/seaice + - improved description of the various forcing configurations + - added basic radiation bulk formulae to pkg/exf + - units/sign fix for evap computation in exf_getffields.F + - updated verification/global_with_exf/results/output.txt +o Added pkg/sbo for computing IERS Special Bureau for the Oceans + (SBO) core products, including oceanic mass, center-of-mass, + angular, and bottom pressure (see pkg/sbo/README.sbo). +o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F + to avoid negative values in shallow regions. + - updated verification/natl_box/results/output.txt + - updated verification/lab_sea/results/output.txt +o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F +o Added useSingleCpuIO option (see PARAMS.h). +o Updated useSingleCpuIO option in mdsio_writefield.F to + work with multi-field files, e.g., for single-file pickup. +o pkg/seaice: + - bug fix in growth.F: QNET for no shortwave case + - added HeffFile for specifying initial sea-ice thickness + - changed SEAICE_EXTERNAL_FLUXES wind stress implementation +o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf, + kpp_transport_t.F, forward_step.F, and the_main_loop.F +o pkg/seaice: + - adjoint-friendly modifications + - added a SEAICE_WRITE_PICKUP at end of the_model_main.F + +checkpoint48e_post +o New S/R for diagnostic of bottom pressure (phi0surf was missing in c48d_post) +o dynamics: change definition of computational domain & adapt mom_fluxform + accordingly ; when viscA4=0, allows to run the dynamics with Olx=Oly=2. +o following exp works with Olx=Oly=2 (and pass through checkbounds option): + adjust_nlfs.cs-32x32x1 adjustment.128x64x1 adjustment.cs-32x32x1 + exp0 exp1 exp2 exp4 exp5 hs94.cs-32x32x5 hs94.128x64x5 hs94.1x64x5 + global_ocean.90x40x15 global_with_exf global_ocean_pressure + inverted_barometer solid-body.cs-32x32x1 + (high order AdvScheme and/or multiDimAd works with Olx=Oly=2 ? not sure) + +checkpoint48d_post +o in preparation for r*: + a) use pre-computed gradient of hydrostatic potential: + changes in timestep.F & mom_cdscheme.F affects results of ideal_2D_oce + b) move phi0surf from calc_phi_hyd to calc_grad_phi_hyd: + => affects results of glob_oce_pressure (different truncation error) + update output.txt of ideal_2D_oce & global_ocean_pressure + +checkpoint48d_pre +o in preparation for r*: + new S/R (calc_grad_phi_hyd.F) to compute Hydrostatic potential gradient. + pass the 2 comp. of the grad. as arguments to momentum S/R. + but for the moment, only used if it does not change the results. +o add finite volume form of calc_phi_hyd for buoyancyRelation OCEANIC & OCEANICP + +checkpoint48c_post +Modifications in treatment of initial flux forcing controls +to avoid extra storing of flux fields. + +checkpoint48c_pre +o Added tidal_basin_2d experiment (currently 1D!) +o can use r* Coordinate with OBC (but Orlanski Rad. BC not yet implemented) + +checkpoint48b_post +o r* coordinate added in #ifdef NONLIN_FRSURF block. + (modification to pressure gradient not yet implemented) +o mom_fluxform: new S/R mom_calc_rtrans.F computes vertical transports + (for U & V) that are used for vertical advection of momentum. + +checkpoint48a_post +o In preparation for r*, use temp. array to compute vertical transport + in flux-form momentum pkg. => change results of few test-exp. + Update output of: exp4, exp5, front_relax, global_ocean_pressure, + lab_sea, natl_box & plume_on_slope + +checkpoint48 +o fix small bug that was causing lot of warning (TIMER_STOP 'UPDATE_CG2D') + +checkpoint47j_post +o add new exp. (ideal_2D_oce) to test "exotic" parameters and options. +o moved adexch from forward_step to addummy_in_stepping +o added exch in ctrl_map_... +o aveFreq.=0. at end of cost_final to avoid overwrite in + reverse checkpointing loop +o Yet more changes in gmredi: + * adgmredi_calc_tensor + avoiding all recomputation of gmredi_slope_limit + * adgmredi_x/y/rtransport + added flag for excessive storing to avoid recomp. of + u/v/rtans, dTdx/y/z + -> this is not really necessary and very memory-consuming + * adgmredi_slope_psi: + consistency with gmredi_slope_limit in treatment of GM_slopeSqCutoff + * gmredi_slope_limit + re-activated full calculation of taperfct for case 'ac02' +o updating verification/ adjoints of carbon, global_ocean.90x40x15 + +checkpoint47i_post +o GM_EXCLUDE_() replace GM_TAPER_(): by default, all the code is compiled. +o add GM_Small_Number & GM_slopeSqCutoff as run-time parameters (data.gmredi) + default values = 1.D-12 & 1.D+48 (=> give same results as checkpoint47f_post) +o cg2dTargetResWunit: define tolerance using deltaTmom (and not deltaTfreesurf) +o gmredi: fix few bugs. + +checkpoint47h_post (duplicate tag for PH's c48 tag, in order to move c48) +This will be a good checkpoint to plug into the ecco-branch. +o few modif.'s to get latest version adjointed + (mainly kick out code in ini_linear_phisurf) +o modif's to run adjoint with exactConserv +o case GM_BOLUS_ADVEC should be cleaned + S/R gmredi_slope_psi should be cleaned +o verification/carbon now has exactConserv=.TRUE. + +checkpoint47g_post +o time-average AIM output: becomes independent from statvars time-average + (can be used with #undef ALLOW_TIMEAVE) ; add new variables (e.g., RH) +o After TAF bug fix: + Tested adjoint of gmredi package. + Tests were done with checkpoint47 rather than with latest + checkpoint47f_post since the latter is broken for the adjoint + verification/carbon/ contains both ADM and TLM gradient checks. +Here they are, before they get changed/lost/stolen. +Mostly modif.'s to fix numerical sensitivities. +Gradient checks OK for +- GM_taper_scheme: + * clipping + * ac02 + * linear + * glw91 + * dm95 + * ldd97 +- GMREDI_OPTIONS: + * GM_VISBECK_VARIABLE_K + * GM_NON_UNITY_DIAGONAL + * GM_EXTRA_DIAGONAL + * GM_BOLUS_ADVEC +in conjunction with data.gmredi parameters to be checked in +in a few minutes under verification/carbon/code/ + +checkpoint47f_post +Merging from release1_p10: +o modifications for using pkg/exf with pkg/seaice + - pkg/seaice CPP options SEAICE_EXTERNAL_FORCING + and SEAICE_EXTERNAL_FLUXES + - pkg/exf CPP options EXF_READ_EVAP and + EXF_NO_BULK_COMPUTATIONS + - usage examples are Experiments 8 and 9 in + verification/lab_sea/README + - verification/lab_sea default experiment now uses + pkg/gmredi, pkg/kpp, pkg/seaice, and pkg/exf + +checkpoint47e_post +o created a tag prior to creating a branch + +checkpoint47d_post +o fixed bug in the oceanic pressure coordinates code: vertical viscosity + at the bottom boundary had an erroneous half slip boundary condition +o added bottom drag and no slip boundary condition capabilities to + oceanic pressure coordinates code +o added new verification experiment global_ocean_pressure, this experiment is + also described in the tutorial section of the manual +checkpoint47d_pre + +checkpoint47c_post +o Added a verification experiment for a global cubed-sphere ocean run + - verification/global_ocean.cubed32x32x30 + - uses pkg/therm_seaice and pkg/bulk_forcing + +o allow to run AIM physics with SPEEDY input files (from Franco Molteni) +o allow a more accurate definition of Ro_Surf (selectFindRoSurf=1) + when using P-coordinate; only implemented for atmospheric config. +o OCEANICP & realFreshWater: include P-E direct effect on wVel ; + NOTES: requires option NONLIN_FRSURF to be "#define". +o update advect_xz/results/output.txt (left from checkpoint44g_post) + +checkpoint47b_post +Merging from release1_p9: +o pkg/seaice + - removed GOTO's and added taf directives + - double precision constants to reduce the g77 (Linux) + to F77 (SGI) differences reported in release1_p8 +o tools/genmake + - added SGI options +o verification/testscript + - updated to that of checkpoint47a_post +o verification/global_ocean.90x40x15/input/eedata + - modified for SGI f77 compatibility +o verification/lab_sea + - added description of sea-ice model + - added missing matlab routines + - added test of thermodynamics parallelization + +checkpoint47a_post +o new pkg: aim_v23 + = F. Molteni atmos.physics (SPEEDY, ver23) adapted to MITgcm. + for now, keep the same surface forcing as before. + - Part-Cell implemented into AIM; check that Heat & Water are conserved + - aim.5l_cs: use new aim pkg (run 10yr & get better results than before) + - aim.5l_LatLon & aim.5l_Equatorial_Channel are still using old aim pkg +o add diagnostic of surface correction term in monitor output + +o bulk_force and therm_seaice +- Two packages: bulk_force (Bulk forcing) + and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2 +- Earlier integration from Stephaine Dutkiewicz + and Patrick Heimbach. +- Two ifdef statements for compile time, + ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE + +checkpoint47 +differentiable version of checkpoint46n_post +o bug fix in quasihydrostaticterms +o * "clean" adjoint code (in terms of extensive recomputations) + can now be obtained for all GMREDI options (i.e. for + - GM_VISBECK_VARIABLE_K + - GM_NON_UNITY_DIAGONAL + - GM_EXTRA_DIAGONAL + - GM_BOLUS_ADVEC ) + * However, wrong gradient check problem remains unsolved. + * New CPP options have been introduced for different + tapering schemes +o external_fields_load now part of differentiation list +o pressure needs multiple storing; + would be nice to have store_pressure at beginning or + end of forward_step, e.g. by having phiHyd global (5-dim.) + (NB: pressure is needed for certain cases in find_rho, + which is also invoked through convective_adjustment). +o recomputations in find_rho for cases + 'JMD95'/'UNESCO' or 'MDJWF' are OK. +o #define ATMOSPHERIC_LOADING should be differentiable +o ini_forcing shifted to begining of initialise_varia +o Incorporating QNVS line search routines into MITgcm + (this is separate code, and therefore not under pkg) + lsopt/, optim/ +o Updated verification/carbon/ + - #define GM_VISBECK_VARIABLE_K + - GM_taper_scheme = 'ldd97' + checkpoint46n_post Merging from release1_p8: o verification/natl_box: