/[MITgcm]/MITgcm/doc/tag-index

Diff of /MITgcm/doc/tag-index

Parent Directory | Revision Log | View Revision Graph Revision Graph | View Patch Patch

-revision 1.2063 by gforget,
Sun Sep 14 15:22:40 2014 UTC
+revision 1.2114 by jmc,
Wed Jan 14 20:49:27 2015 UTC
 Line 4 
 $Name$
      Notes on tags used in MITgcmUV
      ==============================
+ o model/src:
+   - move (from ini_forcing.F) initialisation of all arrays in header file
+     FFIELDS.h into new S/R ini_ffields.F ; this fix broken restart
+     (since Nov 2012) with option selectAddFluid=1.
+ o pkg/shelfice:
+   - change units (+ document) of newly added field "shelfIceMassDynTendency"
+   - upgrade from _RS to _RL shelficeMass & shelficeLoadAnomaly arrays.
+ o pkg/obcs: re-organize obcs_calc_stevens.F to prepare for passive tracers
+ o pkg/seaice: two new runtime parameters:
+   - SEAICE_cStar replaces the hard wired "20" in the strength formulation,
+     long overdue
+   - SEAICE_tensilFac: preparation for Koenig-Beatty+Holland (2012)
+     parameterization of tensil stress for fast ice (does not do anything yet)
+ o pkg/seaice/seaice_check.F: left over from restricted addtive Schwarz method
+     for LSR: move checking overlaps out of JFNK-if-block
+ o pkg/atm_phys:
+   - move copy of grid and dynamical fields out of atm_phys_driver.F into
+     new S/R ATM_PHYS_DYN2PHYS.
+   - fix level height for Linear Free-Surf case (to be consistent with
+     fixed-in-time pressure levels).
+ o pkg/mom_common (bottomdrag):
+   - average bottomDragFld (ALLOW_BOTTOMDRAG_CONTROL) from grid-cell center
+     to correct location (U & V point) in S/R MOM_U/V_BOTTOMDRAG;
+     update output (adm & tlm) of experiment global_ocean.90x40x15.bottomdrag.
+ o pkg/shelfice (u/v_drag):
+   - add one argument (the other velocity component) to S/R SHELFICE_U/V_DRAG;
+   - use SHELFICEDragLinear & SHELFICEDragQuadratic in shelfice_v_drag.F
+     (now similar to shelfice_u_drag.F)
+   - remove condition on bottomDragTerms when calling SHELFICE_U/V_DRAG
+     from mom_vecinv.F (similar to mom_fluxform calls).
+ o pkg/mom_common (bottomdrag):
+   - add one argument (the other velocity component) to S/R MOM_U/V_BOTTOMDRAG;
+   - add options for alternative discretisation of velocity norm in quadratic
+     bottom drag: averaging of grid-cell-center KE to U.point and V.point
+     (selectBotDragQuadr=0, original code) can produce noisy patterns; this
+     improves when averaging only the other component (e.g. vVel for U bottom
+     drag) using wet-point method (selectBotDragQuadr=2) or not (=1).
+   - add options (bottomVisc_pCell) to account for partial-cell in bottom
+     friction (no-slip BC) thus getting the true distance from the bottom;
+     off by default since this also reduces model stability (explicit bottom
+     friction). Changes produce different machine truncation with linear-drag
+     + no-slip BC: update output of experiments exp4.stevens (10 digits) and
+     front_relax (bvp & mxl, 10 & 11 digits).
+ o model/src & mom pkgs:
+   - add one more level to vertical viscosity local arrays (Nr+1, previously Nr)
+     since no-slip bottom BC uses viscosity @ k+1 to update velocity @ level k
+   - for now and until vertical mixing scheme are updated to fill up level Nr+1,
+     just copy level Nr value to level Nr+1.
+ o doc:
+   - add notes form Abhisek Chakraborty on how to close Heat & Salt budgets.
+ o pkg/streamice:
+   - shelfice_init_varia.F, shelfice_init_fixed.F, shelfice_readparms.F,
+     SHELFICE.h: move array shelficemass to init_varia, initialise
+     from array shelficemassinit, which is set in init_fixed and
+     has option to be set from file; introduce array
+     shelficeMassDynTend array, set thru file in init_fixed
+   - NEW FILE: shelfice_step_icemass.F: update shelficemass from either
+     shelfIceFreshWaterFlux and shelficeMassDynTend or H_streamice in
+     pkg/streamice
+   - shelfice_thermodynamics.F: call shelfice_step_icemass if
+     SHELFICEMassStepping = .true.
+   - shelfice_ad_diff.list: add new subroutine shelfice_step_icemass.f
+ checkpoint65h (2014/12/17)
+ o verification/global_ocean.90x40x15, global_ocean_ebm:
+   - update to generic controls.
+ o pkg/ctrl:
+   - ctrl_map_ini_genarr.F: test for igen>0
+   - ctrl_init_variables.F: add call to CTRL_MAP_INI_GENTIM2D.
+   - move gentim2d part of ctrl_map_ini_genarr.F to
+     new separate routine : ctrl_map_ini_gentim2d.F
+   - ctrl_cost_driver.F,ctrl_init.F ctrl_map_ini_genarr.F:
+     add time variable weights and cyclic controls (time mean,
+     seasonal cycle, etc) capabilities to gentim2d
+ o pkg/ecco:
+   - move reading of error, conversion to weights from
+     ecco_cost_init_fixed.F to ecco_toolbox.F (ecco_readwei)
+ o global_oce_llc90, cs32:
+   - activte sbo.
+   - CTRL_SIZE.h : add maxCtrlProc.
+   - activate the new genctrl capabilities.
+   - update following reorganization of global_oce_input_fields.tar.gz
+ o pkg/ctrl/ctrl_ad_diff.list:
+   - add missing ctrl_init_rec.f
+ o tools/OAD_support/keepOriginal:
+   - add ctrl_init_rec.F
+ o model/src/the_main_loop.F, pkg/ctrl/ctrl_map_ini_genarr.F:
+   - avoid recomputations of xx_gentim2d_dummy.
+ o momentum diagnostics:
+   - move filling of diagnostics Um_Diss & Vm_Diss from mom_fluxform.F and
+     mom_vecinv.F to timestep.F (to include Smag-3D contribution).
+ o pkg/seaice/seaice_lsr.F:
+   - fix computation of BU/BV, make sure they are non-zero, only necessary
+     for cubed sphere exchanges and SEAICE_OLx/y>0.
+     This makes calling s/r fill_cs_corner_rl obsolete (removed).
+ o verification/global_ocean.cs32x15/input.seaice
+   - add test for strong implicit coupling and restricted addtive Schwarz
+     methods for LSR
+   - update results/output.seaice.txt
+ o verification/tutorial_global_oce_biogeo:
+   - AD & OAD set-up: switch off ALLOW_OLD_VIRTUALFLUX (like in fwd set-up)
+     and set PTRACERS_EvPrRn(1,2) to zero (for DIC & ALK);
+     update output_adm,_oadm,_tlm.txt
+ o pkg/dic:
+   - change default to "#undef ALLOW_OLD_VIRTUALFLUX" in DIC_OPTIONS.h
+   - add a stop when ALLOW_OLD_VIRTUALFLUX is defined and PTRACERS_EvPrRn
+     for tracer 1 & 2 is set (to avoid double counting of dilution effect).
+ o verification/tutorial_global_oce_biogeo:
+   - AD & OAD set-up: do not set PTRACERS_EvPrRn for tracer 1 & 2 since
+     dilution is already added by pkg/dic ALLOW_OLD_VIRTUALFLUX ;
+     update output_adm,_oadm,_tlm.txt
+ o eesupp/src:
+   - add fill_cs_corner_rl.F (_RL-copy of fill_cs_corner_rs.F)
+     and call it in seaice_lsr when using nonzero SEAICE_OLx/y
+ o pkg/seaice: introduce strongly implicit coupling, intended for stabilizing
+   LSR (following Hutchings et al. 2004)
+   - introduce zetaZ as a global variable (requires adjustments in
+     seaice_evp.F and seaice_jfnk.F) and compute analoguously to etaZ
+     in seaice_calc_viscosities.F
+   - new flag SEAICEuseStrImpCpl (default off)
+   - add new terms zetaZ*du/dy and zetaZ*dv/dx on both sides of the
+     momentum equations
+   - JFNK-solver (preconditioner) adjusted to allow using this method
+     (where it does not have any positive effect)
+ o pkg/seaice
+   - add overlap for restrictive additive Schwartz method also to LSR, but
+     turn it off by default (because it only seems to work for non-cubed
+     sphere/llc topologies)
+ o pkg/autodiff
+   - add extra exchanges in addummy_in_steppging also for mnc output
+ o pkg/seaice:
+   - add extra exchanges before dumping adjoint fields in seaice_ad_dump.F
+     to avoid stripes in the adjoint output (similar to addummy_in_stepping.F)
+ o pkg/obcs:
+   - add option to specify domain-connected piece Id along OB though files
+     (for now, only used to diagnose PhiVel; default =1 -> global correction)
+   - modify obcs_diag_balance.F to process each connected-part separately;
+     exclude OB where Id = 0 but keep the full divergence there (-> constant
+     Psi along this portion of OB and all inflow derived from grad(Phi))
+   - add example of OB-conect files in verification/so_box_biogeo/input.
+ o pkg/diagnostics:
+   - improve treatment of Open-Boundaries in Velocity-Potential calculation:
+   a) use non zero matrix element across OB and solve for Phi 1 grid point
+      outside OB (i.e, where tracer OBCS are applied).
+   b) global imbalance of inFlow/outFlow across OB: instead of modifying
+      velocity across OB, spread it uniformly along the OB in the RHS.
+ checkpoint65g (2014/11/19)
+ o pkg/diagnostcis:
+   - add preconditioner off-diagonal factor (diagCG_pcOffDFac) as run-time param
+ o pkg/my82,pp81,kl10:
+   - fix ALLOW_3D_DIFFKR implementation (previously had no effect): use 3-D
+     field diffKr as background diffusivity to set minimum diffusivity value.
+ o pkg/seaice:
+   - add new parameter SEAICEuseEVPrev (for revisited EVP by
+     Bouillon et al 2013)
+   - add SEAICE_evpAlpha/Beta to namelist
+   - modify seaice_evp code to able to use SEAICEuseEVPrev. These changes
+     involve trivial re-ordering of divisions/multiplications that should not
+     have any effects, but ... they do change the results of lab_sea.hb87 so
+     that only 2 digits of agreement remain, in line with the general
+     instability/fragility of the evp-code. For now, I do not care and just
+     update the verification experiment
+ o pkg/ctrl : add generic pre-processing (xx_genarr2d_preproc etc)
+   - remove xx_genarr2d_numsmooth etc
+   - use xx_genarr2d_preproc etc instead
+   - carry out smooth_correl2D at initialization stage rather than during
+     time stepping.
+   - get smoothOpNb from xx_gentim2d_preproc
+ o global_oce_cs32 :
+   - add test of generic ecco capabilities.
+ o pkg/exf/exf_getsurfacefluxes.F :
+   - add to tmpUE/tmpVN rather than re-initializing them (conventional behavior).
+ o pkg/ctrl,ecco :
+   - omit un-necessary arguments to ctrl_cost_gen.F unless ECCO_CTRL_DEPRECATED
+     is defined.
+ o pkg/seaice: add diagnostics for Delta (deformation parameter)
+ o model/src:
+   - when ALLOW_3D_DIFFKR is defined, add a Warning if any of the vertical
+     diffusivity (background) is specified but ignored.
+ o verification:
+   - in testreport, convert hostname to lower case for default output dir name
+ o model, pkg : add useCTRL, useECCO parameters (on/off run time switch).
+ o pkg/dic: add modifications from Hajoon:
+   - add iron sedimentation flux function of the flux of PO4 out of bottom layer
+   - light attenuation: add self-shading effect by climatological Chlorophyll
+     (read from file).
+ o pkg/seaice/seaice_preconditioner.F : fix cpp flags, so that code also
+     compiles when SEAICE_ALLOW_DYNAMICS is undefined
+ o global_oce_llc90,cs32 : ecco_v4 verification exps now use separate
+     AUTODIFF_OPTIONS.h, CTRL_OPTIONS.h, ECCO_OPTIONS.h rather than
+     ECCO_CPPOPTIONS.h (i.e. replacing the deprecated, circumvoluted,
+     approach where CPP_OPTIONS.h includes an omnipotent
+     ECCO_CPPOPTIONS.h that superseeds pkgs '_OPTIONS.h').
+ o pkg/ctrl/ctrl_map_ini_genarr.F : include GMREDI_OPTIONS.h where
+     ALLOW_3D_KAPGM, ALLOW_3D_KAPREDI (to be added later) belong.
+ o pkg/seaice/SEAICE_SIZE.h : replace ALLOW_AUTODIFF_TAMC with ALLOW_AUTODIFF
+     to avoid needing AUTODIFF_OPTIONS.h anytime SEAICE_SIZE.h is included
+     (it seems that THSICE_SIZE.h, PTRACERS_SIZE.h have the same issue...)
+ o model/src, pkg/smooth, profiles, salt_plume, seaice, exf, ecco, ctrl :
+   add ECCO_OPTIONS.h, AUTODIFF_OPTIONS.h, CTRL_OPTIONS.h when needed
+   - ECCO_OPTIONS.h is needed when including ecco_cost.h, ecco.h
+   - AUTODIFF_OPTIONS.h is needed when including tamc.h, tamc_keys.h
+   - CTRL_OPTIONS.h is needed when including ctrl.h, optim.h, etc
+ o pkg/ctrl,ecco,exf :
+   - remove variables associated with NON-generic costs,
+     unless ECCO_CTRL_DEPRECATED is defined.
+   - remove NON-generic NON-obcs controls,
+     unless ECCO_CTRL_DEPRECATED is defined.
+ o pkg/ctrl : added ctrl_local_params.h
+ o pkg/ecco : added ecco_local_params.h
+ o verification : added ECCO_CTRL_DEPRECATED to most as and oad
+     experiments, which will progressively be removed once the
+     experiments are revised to use only generic controls/costs.
+ checkpoint65f (2014/10/13)
+ o pkg/ctrl, ecco, seaice : minor changes in headers
+ o global_oce_llc90,cs32 : ecco_v4 verification exps now test generic
+     cost function and control capabilities added since checkpoint65c.
+ o pkg/ctrl: regroup all obcs ctrl variables in CTRL_OBCS.h.
+ o pkg/ecco: increase modularity and improve logic
+   - ecco_toolbox.F (new) : basic operations that are often used,
+     including ecco_zero, ecco_cp, ecco_cprsrl, ecco_diffmsk,
+     ecco_addcost, ecco_add, ecco_div, ecco_readbar.
+   - cost_gencal.F (new) : determine filnames and records.
+   - cost_genread.F, cost_generic.F : condense code using ecco_toolbox.
+   - cost_generic.F : use spzeroloc rather than hard-coded -9999,
+     and exploit the above functions leading to more compact code.
+   - cost_genread.F : rm un-necessary arguments, re-arrange args order.
+   - ecco_check.F, ecco_readparms.F : bring logic up to date (using_tpj,
+     using_ers, using_gfo, using_mdt is now only defined within gencost)
+ o pkg/ecco: extend generic cost capabilities
+   - pkg/ecco : added 3D gencost (ALLOW_GENCOST3D), added gencost
+     capabilities via gencost_preproc (+ _c,_i,_r related params),
+     gencost_posproc (+ _c,_i,_r related params), gencost_nrecperiod,
+     gencost_outputlevel arguments to cost_generic
+   - pkg/ecco/cost_genread.F (new) : interface between cost_generic
+     and active_read, making room for gencost_preproc to take effect.
+     Added preproc code ('climmon') to replace GENERIC_BAR_MONTH.
+   - pkg/ecco/cost_generic.F :
+     compute and store misfit maks*(model-data) in localdif that is
+     then used in cost and for output, output misfit if outlev.GT.0
+   - pkg/ecco/cost_averagesfields.F : depend on gencost_barskip
+     rather gencost_barfile for deciding whether to write barfile
+   - pkg/ecco/ecco_cost_driver.F : extract call to cost_profiles from cost_hyd
+   - pkg/ecco/ecco_cost_init_fixed.F : extract gencost weights from
+     ecco_cost_weights, re-activate ecco_summary, determine gencost_barskip
+   - pkg/ecco/ecco_summary.F : report gencost variables
+   - pkg/ecco/cost_gencost_customize.F : treat case of scatterometer costs
+ o pkg/sea/SEAICE_OPTIONS.h : add ALLOW_COST_ICE, ALLOW_SEAICE_COST_SMR_AREA
+ o pkg/ecco: internal re-organization
+   - resolve un-necessary dependancies within pkg/ecco and to pkg/ctrl
+     (rm includes of optim.h and ctrl.h, replace optimcycle with eccoiter,
+     mv relevant variables definition from ecco_cost.h to ecco.h,
+     make sure gencost routines only include ecco.h)
+   - throughout : add missing CCP switches, rm unused variables,
+     and rm un-necessary CPP switches
+   - throughout gencost routines : replace modelstartdate with
+     gencost_startdate when needed
+   - cost_bp_read.F : use arguments, instead of common blocs
+     so that it can be used also by (new) cost_gencost_bpv4.F
+   - ecco_init_varia.F : initialize RHOsumGlob_0, VOLsumGlob_0
+   - ecco_phys.F : compute etanFull
+   - ecco_readparms.F : rename using_topex as using_tpj, add
+     bpv4-grace to gencost list, initialize added parameters
+   - ecco_summary.F : include SIZE.h needed for augmented ecco.h
+ checkpoint65e (2014/09/28)
+ o pkg/obcs:
+   - replaced the just-added CPP options with run-time variables
+     (see pkg/obcs/OBCS_PARAMS.h for details)
+ o pkg/obcs:
+   - add additional CPP options for sponge boundary layers
+     (see pkg/obcs/OBCS_OPTIONS.h for details)
+ o pkg/atm_phys:
+   - implement simple (obliquity only, assuming circular planet orbit)
+     seasonal cycle for incoming insolation ; used when run-time parameter
+     "select_incSW" is set to 1 ; by default (=0), uses old unchanged formula.
+ o pkg/exf, pkg/ecco, pkg/ctrl:
+   - add ctrlUseGen and .NOT.ctrlUseGen exclusive brackets
+   - activate gentim2d for xx_atemp,xx_aqh,xx_swdown,xx_lwdown,xx_precip
+   - implement gentim2d for xx_tauu, xx_tauv
+ o pkg/ctrl, pkg/ecco :
+   - fix and complement the generic controls implementation
+ o pkg/exf/exf_getclim.F, exf_getffields.F, exf_getsurfacefluxes.F,
+     pkg/shelfice/shelfice_thermodynamics.F,
+     verification/hs94.1x64x5/code_ad/ctrl_map_gentim2d.F,
+     verification/hs94.1x64x5/code_oad/ctrl_map_gentim2d.F :
+   - add weight in CTRL_GET_GEN calls
+ o pkg/ctrl :
+   - ctrl_get_gen.F : add weight as input parameter, and apply
+     it directly rather than using smooth_correl2Dw (to be retired).
+   - ctrl_map_gentim2d.F : add wgentim2d in CTRL_GET_GEN call
+ o pkg/ctrl :
+   - introduce run time switch ctrlUseGen (in ctrl.h, ctrl_check.F,
+     ctrl_init_variables.F, ctrl_readparms.F, ctrl_summary.F)
+   - re-activate ctrl_summary (call in ctrl_init.F)
+ o pkg/ctrl/ctrl_cost_gen.F :
+   - fix CCP brackets, remove un-necessary myiter, mytime
+ o pkg/ecco, pkg/ctrl :
+   - remove cost_genctrl.F (now pkg/ctrl/ctrl_cost_driver.F)
+   - remove cost_forcing_gen.F (now pkg/ctrl/ctrl_cost_gen2d.F)
+   - add ctrl_cost_driver.F (replaces pkg/ecco/cost_genctrl.F)
+   - add ctrl_cost_gen2d.F (replaces pkg/ecco/cost_forcing_gen.F)
+ checkpoint65d (2014/09/15)
+ o pkg/ctrl/CTRL_GENARR.h, ctrl_readparms.F :
+   - add xx_genarr2d_bounds,xx_genarr3d_bounds,
+     xx_genarr2d_numsmooth, xx_genarr3d_numsmooth
+ o pkg/ctrl/ctrl_map_ini_genarr.F :
+   - remove compilation failure
+   - replace hard coded example
+   - use ctrl_map_ini_gen2D, 3D
+ o pkg/ctrl/ctrl_pack.F, ctrl_unpack.F :
+   - on/off switch based upon xx_genarr2d_weight,3d
  o src/packages_init_variables.F :
    - regroup pkg/estimation calls
    - remove CTRL_MAP_INI, CTRL_MAP_INI_ECCO
-Line 267 
 o pkg/seaice: more clean up
+Line 582 
 o pkg/seaice: more clean up
      nearly a factor of 10 (but keep in mind there is mostly seaice code in
      this experiment):
      > size mitgcmuv_ad
      text    data        bss     dec        hex
      5100659 1216  220458952  225560827  d71c8fb
      5339425 1200 1868783272 1874123897 6fb4dc79 (SEAICE_LSR_ADJOINT_ITER)
      but only a moderate increase of 36% for global_ocean.cs32x15:
      text    data        bss     dec        hex

 Legend:



Removed from v.1.2063
 


changed lines


 
Added in v.1.2114
 Legend:



Removed from v.1.2063
 


changed lines


 
Added in v.1.2114
-Removed from v.1.2063
+Added in v.1.2114

	ViewVC Help
Powered by ViewVC 1.1.22