--- MITgcm/doc/tag-index 2002/09/05 20:49:32 1.83 +++ MITgcm/doc/tag-index 2002/12/27 14:26:14 1.112 @@ -1,9 +1,252 @@ -$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.83 2002/09/05 20:49:32 mlosch Exp $ +$Header: /home/ubuntu/mnt/e9_copy/MITgcm/doc/tag-index,v 1.112 2002/12/27 14:26:14 cheisey Exp $ Notes on tags used in MITgcmUV ============================== +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: + updating new external_fields_load routine +o New package: pkg/seaice + Sea ice model by D. Menemenlis (JPL) and Jinlun Zhang (Seattle). + The sea-ice code is based on Hibler (1979-1980). + Two sea-ice dynamic solvers, ADI and LSR, are included. + In addition to computing prognostic sea-ice variables and diagnosing + the forcing/external data fields that drive the ocean model, + SEAICE_MODEL also sets theta to the freezing point under sea-ice. + The implied surface heat flux is then stored in variable + surfaceTendencyTice, which is needed by KPP package (kpp_calc.F and + kpp_transport_t.F) to diagnose surface buoyancy fluxes and for the + non-local transport term. Because this call precedes model + thermodynamics, temperature under sea-ice may not be "exactly" at + the freezing point by the time theta is dumped or time-averaged. + N.B.: fairly large differences in accuracy occur across + different platforms/compilers; comparison between + g77 (Linux) and F77 (SGI Origin 2000/MIPSpro) + generated output gives: + T S U V + C D M c m s m s m s m s + n p a R g m m e . m m e . m m e . m m e . + f n k u 2 i a a d i a a d i a a d i a a d + g d e n d n x n . n x n . n x n . n x n . + + Y Y Y Y 5 5 7 7 7 8 10 9 6 6 6 6 7 5 7 5 7 FAIL lab_sea +o GAD: + - generated new common blocks to account for call of + same gad routines with differing traceridentities + (needed to modify tracerIdentity indices in GAD.h) + - generated separate common blocks for case useCubedSphereExchange + (Department of Futurology) + - parameter lists to gmredi_?transport: added tracerIdentity + - added new key indices to tamc.h +o external_fields_load: + - added this routine to TAF list + - needed to make some common blocks global and additional storing + along the same lines as exf package (checkpoint_lev?_directives.h) +o exf: + updated external forcing package + - bug fixes carried over from ecco-branch + (missing OBCS_OPTIONS.h in two routines) + - enable easy to use "no forcing". + - added exf I/O for atmospheric loading + - added exf I/O for runoff data + - transfered scaling between exf <-> MITgcm to exf namelist +o Adding new verification experiment global_with_exf. + This verification experiment is almost identical to the + global_ocean.90x40x15 experiment, except that it uses + the exf (and cal) package instead of the + external_fields_load routine (cf README file). + +checkpoint46m_post +o modified testscript to return -- when there are no lines of "%MON" output +o added new flag "quasiHydrostatic" for QH mode + and new routine quasihydrostaticterms() +o added 2*Omega*Cos(phi)*W to u equations (non-hydrostatic Coriolis term) + New routine: mom_u_coriolis_nh() +o fixed sign of non-hydrostatic metric terms in u and v equations. + Depend on correct sign of gravitySign. +o added work.pc* to .cvsignore in aim directories using ifc +o introduced flag for controllin non-hydrostatic metric terms +o removed gravitySign from PARM04 +o switched to g77/ifc from pgf77 in aim experiments + +checkpoint46l_post +o made convective adjustment work with pressure coordinates: + - changed the direction of k-loop in convective_adjustment.F for the + case of pressure coordinates (OCEANICP,ATMOSPHERIC buoyancyRelation) + - adjusted the reference pressure k-index in convective_adjustment.F + - adjusted the convection condition in convect.F (in analogy to + calc_ivdc.F) + - convective_adjustment no longer computes anything on the halos + - removed the warnings about negative salinity from find_rho.F and + find_alpha.F; instead the new routine look_for_neg_salinity, called + at the beginning of find_rho, find_alpha, and find_beta, does a + check of the entire slice, if CPP-option + CHECK_SALINITY_FOR_NEGATIVE_VALUES is defined +checkpoint46l_pre + +checkpoint46k_post +o fixed the verification/global_ocean.90x40x15 experiment: + - new bathymetry (the world according to A., JMC, and M.) + - new initial fields and forcing fields (*.bin files) + - new POLY3.COEFFS (for the next release one should switch to a full + equation of state: JMD95P or MDJWF) + - fixed several errors and redundancies in the data file + - experiment uses looped cells + - added matlab directory with diagnostic scripts for plotting of output + +o S/R aim_initialise.F replace S/R aim_init from file aim_do_inphys.F: + - read AIM physics parameters from a file (data.aimphys) + - set defaults values = F.Molteni paper (Clim.Dyn., 2002) + +checkpoint46j_post +o split calc_exact_eta in 2 S/R : integr_continuity & update_etaH +o move wVel computation at the end of the time step, in S/R integr_continuity +o create specific S/R to exchange T,S before DYNAMICS (for stagger time step) +o update timeave pkg for wVel diagnostic. +checkpoint46j_pre + +checkpoint46i_post +o Clean up AIM package (and keep the results unchanged): + include CPP_OPTION and use IMPLICT NONE in all routines ; + declare all the variables _RL ; use _d 0 for all numerical constants. + use ifdef ALLOW_AIM everywhere. And now AIM can be used with g77 ! + +checkpoint46h_post +o cleaned up the use of rhoNil and rhoConst. + - rhoNil should only appear in the LINEAR equation of state, everywhere + else rhoNil is replaced by rhoConst, e.g. find_rho computes rho-rhoConst + and the dynamical equations are all divided by rhoConst +o introduced new parameter rhoConstFresh, a reference density of fresh + water, to remove the fresh water flux's dependence on rhoNil. The default + value is 999.8 kg/m^3 +o cleanup up external_forcing.F and external_forcing_surf.F + - can now be used by both OCEANIC and OCEANICP +checkpoint46h_pre + + Added code to convert surface volume flux (fresh water) into + a mass flux when using P coordinates in the ocean (OCEANICP). + Note: It assumes you have set rho0=rhoConst=density of fresh water. + +checkpoint46g_post +o Include a new diagnostic variable phiHydLow for the ocean model + - in z-coordinates, it is the bottom pressure anomaly + - in p-coordinates, it is the sea surface elevation + - in both cases, these variable have global drift, reflecting the mass + drift in z-coordinates and the volume drift in p-coordinates + - included time averaging for phiHydLow, be aware of the drift! +o depth-dependent computation of Bo_surf for pressure coordinates + in the ocean (buoyancyRelation='OCEANICP') + - requires a new routine (FIND_RHO_SCALAR) to compute density with only + Theta, Salinity, and Pressure in the parameter list. This routine is + presently contained in find_rho.F. This routine does not give the + correct density for 'POLY3', which would be a z-dependent reference + density. +o cleaned up find_rho + - removed obsolete 'eqn' from the parameter list. +o added two new verification experiments: gop and goz + (4x4 degree global ocean, 15 layers in pressure and height coordinates) +checkpoint46g_pre + +checkpoint46f_post +o Enable tangent linear (forward mode) gradient checks: + - extended active file handling to g_... files + - added TANGENT_SIMULATION to theSimulationMode + - extended grdchk package accordingly +o added tangent linear setup in verification/carbon/code_ftl +o added adjoint setup in verification/front_relax/code_ad + +checkpoint46e_post o Added new equation of state -> MDJWF - EOS of McDougall et al., 2002, JAOT, submitted - caveat: the equation of state is only valid for a smaller (more @@ -12,7 +255,7 @@ - added more check values for density in check_eos (ini_eos.F), some of the old check values are out of the range of the MDJWF-eos, so don't expect perfect matches for those -checkpoint47d_pre +checkpoint46e_pre checkpoint46d_post o fixed store_pressure to work with both buoyancy relation = 'OCEANIC' and