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

Diff of /MITgcm/doc/tag-index

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

revision 1.53.4.32.2.3 by dimitri, Sat Mar 29 17:35:57 2003 UTC revision 1.86 by mlosch, Wed Sep 18 16:38:01 2002 UTC
# Line 3  $Header$ Line 3  $Header$
3      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
4      ==============================      ==============================
5    
6  Branch icebear  checkpoint46g_post
7    o Include a new diagnostic variable phiHydLow for the ocean model
8  29 mar 03 (cvs tag icebear5)    - in z-coordinates, it is the bottom pressure anomaly
9   Some additional options for pkg/cost/cost_ice.F    - in p-coordinates, it is the sea surface elevation
10   Updates to verification/natl_box_adjoint    - in both cases, these variable have global drift, reflecting the mass
11        drift in z-coordinates and the volume drift in p-coordinates
12  7 mar 03    - included time averaging for phiHydLow, be aware of the drift!
13   Branched off ecco_c44_e27  o depth-dependent computation of Bo_surf for pressure coordinates
14   pkg/ctrl/ctrl_init.F    in the ocean (buoyancyRelation='OCEANICP')
15     Fix for determination of startrec and endrec.    - requires a new routine (FIND_RHO_SCALAR) to compute density with only
16   pkg/cal/cal_toseconds.F      Theta, Salinity, and Pressure in the parameter list. This routine is
17     Fix for negative intervals shorter than 1 day.      presently contained in find_rho.F. This routine does not give the
18   pkg/cost/cost_ice.F, cost.h, cost_readparms.F      correct density for 'POLY3', which would be a z-dependent reference
19     Added costIceStart and costIceEnd.      density.
20   verification/natl_box_adjoint.  o cleaned up find_rho
21     Updated test experiments    - removed obsolete 'eqn' from the parameter list.
22    o added two new verification experiments: gop and goz
23  ====================================================================    (4x4 degree global ocean, 15 layers in pressure and height coordinates)
24    checkpoint46g_pre
25  ecco_c44_e27  
26  o exf:  checkpoint46f_post
27    - Enable initialisation of forcing fields to constant  o Enable tangent linear (forward mode) gradient checks:
28      (runtime) values.    - extended active file handling to g_... files
29    - in exf_getffields.F    - added TANGENT_SIMULATION to theSimulationMode
30      Reduce i-/j-loop to interior domain, discarding overlaps.    - extended grdchk package accordingly
31      That also fixes wrong TAF-key computations for key_1, key_2  o added tangent linear setup in verification/carbon/code_ftl
32      with bulf formulae.  o added adjoint setup in verification/front_relax/code_ad
33    - exf_init.F modify #ifdef for exf_init_evap  
34    - exf_getffieldrec.F, ctrl_getrec.F  checkpoint46e_post
35      The following INT-usages are not safe:  o Added new equation of state -> MDJWF
36        fldsecs  = int(fldsecs/fldperiod)*fldperiod    - EOS of McDougall et al., 2002, JAOT, submitted
37        fldcount = int(fldsecs/fldperiod) + 1    - caveat: the equation of state is only valid for a smaller (more
38      and were modified.      realistic?) range of values than JMD95P/Z and UNESCO
39    - added repeatPeriod (from MAIN branch)    - added masks to the calculation of pressure in store_pressure
40  o grdchk: updated and merged    - added more check values for density in check_eos (ini_eos.F), some of
41    - has both ADM and TLM checks      the old check values are out of the range of the MDJWF-eos, so don't
42    - works for single- and multi-proc.      expect perfect matches for those
43    - output cleaned  checkpoint46e_pre
44    - worked successfully for parallel DIVA  
45    Needed one update in EEPARAMS.h and updates in ctrl package  checkpoint46d_post
46    More modif's may be necessary for ECCO code to support TLM.  o fixed store_pressure to work with both buoyancy relation = 'OCEANIC' and
47  o autodiff: checkpoint_lev...    'OCEANICP', also initialised field pressure correctly in ini_eos in the
48    - updated storing for bulk formulae (some arrays were missing)    case of pressure coordinates. eosType='JMD95Z' in
49    - added state for seaice package    combination with buoyancyRelation='OCEANICP' now causes an error.
50    - synchronized with main branch  o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
51  o DIVA:    in find_alpha/beta.
52    Modif.s in  checkpoint46d_pre
53    - cost_init.F  
54    - cost_averagesfields.F  o Changed p = pressure(i,j,k,bi,bj)  to  p = pressure(i,j,kRef,bi,bj)
55    - ecco_the_main_loop.F    so that JMD95Z and JMD95P give approptiate static stability in find_rho.
56    - ecco_the_model_main.F  
57    and added cost_final_restore.F cost_final_store.F  checkpoint46c_post
58    o With these changes, it is possible to fully separate  o Changes necessary for ocean in p-coordinates
59      the first forward run to compute the cost function (mdthe_main_loop)   -  Added new buoyancy relation = 'OCEANICP'
60      and adjoint runs over individual DIVA intervals (adthe_main_loop)   -  Added new parameters = gravitySign (this used to be contained inside
61      into separate model executions.      the factor dRdZ which I added when we first switched to R coordinates).
62      This enables better tuning of DIVA to given queue sizes.   X GM/Redi is not compatible (yet)
63    o Loss of cost function value in consecutive DIVA runs is   X bottom drag and no-slip need to be debugged.
64      fixed by storing to 'costfinal' file.  checkpoint46c_pre
65    o Appropriate call of ctrl_unpack/_pack,grdchk_main at  
66      beginning/end of adjoint integration is accounted for.  checkpoint46b_post
67    The current changes apply to TAF-1.4.26.  o Added new equation of state -> JMD95Z and JMD95P
68    TAF has one remaining bug, requiring one manual change    - EOS of Jackett and McDougall, 1995, JAOT
69    (see MITgcm manual).    - moved all EOS parameters into EOS.h
70  o bug fix plus modif for ubar, vbar in cost_averagesfields    - new routines ini_eos.F, store_pressure.F
71  o cal: synchronized between MAIN and ecco-branch  o Added UNESCO EOS, but not recommended because it requires
72      in-situ temperature (see JMD95)
73  icebear4  o Modified formatting for knudsen2.f in utils/knudsen2 and added
74  o Added cost_ice for sea-ice volume sensitivity experiment.    unesco.f to be used with POLY3
75  o ecco_the_model_main.F: SEAICE_WRITE_PICKUP must precede WRITE_CHECKPOINT  checkpoint46b_pre
76  o Reduced derivative singularities in pkg/seaice.  
77  o Updated verification/natl_box_adjoint  checkpoint46a_post
78    o Added (atmospheric) pressure loading
79  icebear3    - new field in FFIELDS.h, etc...
80  o Removed spurious iikey and idkey    - new cpp flag ATMOSPHERIC_LOADING
81  o pkg/exf  o Changed hFacC to _hFacC in calc_phi_hyd.F
82    - bug fix in treatment of downward radiation  o Added SHORTWAVE_HEATING to some files for consistency
83  o pkg/seaice  checkpoint46a_pre
84    - added initial sea-ice thickness parameter SEAICE_initialHEFF  
85  o updated verification/natl_box_adjoint  checkpoint46
86    further merges from release1_p5
87  icebear2  o new ctrl package
88  o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,    - adopted from ECCO environment to enable optimization
89    kpp_transport_t.F,and ecco_the_main_loop.F    - added Eliassen Palm fluxes to controls
90  o adjoint-friendly modifications to pkg/seaice  o added Eliassen Palm flux controls to gradient check package
91  o added verification/lab_sea as per release1_p12_pre  o cost package
92      - Compatible with new ctrl package
93  ecco_c44_e26    - added Eliassen Palm cost hooks
94  o Merging pkg/exf and pkg/seaice from cvs tag ecco_ice2    - modif's of existing cost functions
95  o exf: Modifications to correctly handle evap field for adjoint  o Adjoint-related bug fixes in kpp:
96  o grdchk:    - kpp_calc: sore of kpphbl avoids recomputation/call to S/R kppmix
97    Added code to perform gradient checks for bulk formulae/atmos. state    - kpp_routines: store of Rib avoids partial recomputation bug of TAF.
98    (atemp, aqh, uwind, vwind).  o autodiff package
99    NOTE:    - flow directives:
100    This package is out of synch with c48 package.      + new for S/R do_the_model_io
101    The latter also has the tangent linear gradient checks.      + enhanced for S/R checkpoint
102    Need to be merged.      + new gmredi.flow
103  o verification/natl_box_adjoint:    - Introduced CPP option ALLOW_AUTODIFF_MONITOR to
104    * added new setup: code_bulk/, input_bulk, results_bulk      disable adcommon.h/g_common.h by default.
105      (Experiment 2, described in README)    - fixed problem for adjoint of global_max
106      uses bulk formulae and atmos. state as control instead of  o modified verification/carbon/
107      air-sea fluxes.  o New setup for adjoint of global_ocean
108      Performs gradient checks for first 4 elements of atemp control.  o added hooks for Eliassen Palm fluxes (dynamics.F, ini_fields.F)
109    * changed standard verification to perform  o cleaned TAF keys iikey, idkey (dynamics.F, thermodynamics.F)
110      gradient checks on first 4 elements on hflux control.  o incorporated hooks for sponge layer code (external_forcing.F)
111    * gmredi switched off (will come after merge with c48).  
112    checkpoint45d_post
113  ecco_c44_e25  o import 2 fixes from release1_p5 (CG2D.h & packages_readparms.F)
114  o Merging some GM changes from main trunk (c47) to ecco-branch.  o use recip_dx*,recip_dy* instead of /dx*,/dy* in orlanski_E,W,N,S
115  o Updating verif. natl_box_adjoint with GMRedi.    (affects plume_on_slope "cg2d" output)
116  o Adding setup with explicit differentiaton of pressure solver.  o add new parameter deltaTfreesurf for free-surface time stepping
117  o Bug fix: missing OBCS_OPTIONS.h for sponge layer code.  o exchange gT,gS when using staggerTimeStep & GM & Oly<4
118  o (Re-?) Commiting old bug fixes in kpp to ecco-branch  o NONLIN_FRSURF: rescale gNm1 to get a better conservation with A-B
119    (kpp_calc.F, kpp_routines.F)  o change phiMin in exp1 (to agree with documentation)
120  o exf:  o add new flags:
121    updated external forcing package   - T,S forcing inside/outside Adams-Bashforth
122    - cleaned old exfa stuff   - turn off temp,salt advection and/or forcing
123    - bug fixes   - for each tracer, define internal flag for multiDimAdvection & A-B
124      (missing OBCS_OPTIONS.h in two routines)  o monitor: change definition of KE to check conservation of Energy
125    - enable easy to use "no forcing".  
126    - added exf I/O for atmospheric loading  checkpoint45c_post
127    - added exf I/O for runoff data  o Removed f90 comments ("!")
128    - transfered scaling between exf <-> MITgcm to exf namelist  o Modified initialisations to enable exf package for MITgcm
129      (without invoking ECCO_PACKAGE).
 ecco_c44_e24  
130  o modifications to gradient check package (Martin Losch)  o modifications to gradient check package (Martin Losch)
131    - enable centered differences    - enable centered differences
132    - modified format of standard output    - modified format of standard output
133  o bug fix for nwettile contents  
134    so far assumed global sum  in first entry (1,1,k),  checkpoint45b_post
135    needed for global control vector  o cleaned exchange calls in cg2d
136    -> introduced new variable nwet[c/w/s]global(nr)  o Included CPP option ALLOW_AUTODIFF_MONITOR in forward_step
137  o shifts of calls to cal/exf package to be accessible  o included CPP option SINGLE_LAYER_MODE in thermodynamics
138    in MITgcm stand-alone    to configure barotropic setup (Martin Losch)
139  o shifted packages_init_variables before convective_adjustment_ini  o moved some initialisations within bi,bj loop in dynamics
140    (corresponding update of TAF tape initialisation).    (TAF flow dependency)
141  o removed auxiliary field in impldiff  o in initialise_varia, moved call packages_init_variables
142    (caused numerical differences for GMRedi cases)    before convective_adjustment_init
143  o updated obcs control part (G. Gebbie)    (corresponding adjustment of TAF tape initialisation).
144  o removed some f90 comments ("!")  o in calc_diffusivity/calc_viscosity extended domain to full overlap.
145  o Bug fixes and updates in cost package  
146   * cost_final: correct multiplier mult_tauv, mult_sflux  checkpoint45a_post
147   * cost_hyd: added call cost_ssh  o compute Rho over the whole domain (extend iMin,jMin indices)
148   * cost_init_barfiles.F: replaced arrays by tmpfld's  o specific exchange for X-slice domain (case Ny=1)
149   * cost_readtopexmean.F: ersmask -> tpmeanmask  o argument futureIter of S/R obcs_calc missing in exp4 & internal_wave /code
150   * cost_weights.F: corrected wtheta[2], wsalt[2]  o solve_for_pressure : set cg2d_x to zero where OBCS are applied
151  o cleaned exchanges in cg2d.F  
152    checkpoint44h_post = checkpoint45
153  ecco_c44_e23  o Fix to restore differentiability:
154  merged patches between release1_p1 and release1_p2:    place at which packages_init_variables is called in
155  o fixed array boundary problem for adjoint    initialise_varia is crucial for reverse control flow.
156    (related to partial recomputation within subroutine  o modifications in GMRedi to enable stable adjoint
157    related to dynamics -> calc_viscosity -> mom_...    (corresponding modif's in thermodynamics, dynamics)
158  o fixed recomputation in thermodynamics  o added missing hook for sponge layer code in external_forcing
159  o prepare (but not complete) recomputation problem in impldiff  o modified test cost functions
160  o fixed problem for adjoint of global_max  o storing in gad_advection (plus removal of write(0,*)!
161  o switched off sponge layer code in external_forcing  o missing headers TR1.h added in convective_adjustment routines
162  Bug fix in ctrl_mask_set_yz (spotted by G. Gebbie)  o errorMessageUnit set non-equal zero in eeboot
163      (conflict on some platforms, e.g. IBM SP3)
164  ecco_c44_e22  o modified carbon verif.
165  o Finally, the merges with the Scripps cost function terms:  
166    Various additions, modifications, quality checks.  checkpoint44h_pre
167  o Plus: Included switch ALLOW_AUTODIFF_MONITOR to handle  o vertical grid option: allow to put the Interface at the middle between
168    monitoring of adjoint state.    2 cell-centers ; replace delZ in pkg/kpp by drF.
169  o Plus: removed comm_stats, ini_communication_patterns from diff. list  o GM Advective form: Tracers are advected using the residual transport (=
170    (no more need to "keep" MPI headers for adjoint generation)    Euler+GM-bolus); set param GM_AdvSeparate=T to return to previous form.
171  o Plus: exf changes:  
172    o added run-off routine  checkpoint44g_post
173    o additional storing needed for ustress,vstess  o fix surface correction term in multiDimAdvection (affects plume_on_slope
174      because of state-dependent quality check    results) ; ifdef missing in thermodynamics.F
175    o added missing initialisations  
176  o Plus: genmake with MAKEDEPEND variable was still missing in ecco-branch  checkpoint44f_pre,checkpoint44f_post
177  NB: - verification/natl_box_adjoint/ 240 steps run  o added PTRACERS package
178        * looks ok under Linux (twain.lcs.mit.edu)    This allows an arbitrary number of passive tracers to be integrated
179        * looks ok on MIT SUN's (fjord.mit.edu)    forward simultaneously with the dynamicaly model.
180        * NaN in reverse mode on Scripps SUN's at iter=214    + Implemented so far:
181        * NaN in reverse mode on IBM SP3 (BlueHorizon) at iter=204      - basic forward algorithm (time-stepping, advection, diffusion, convection)
182      - 1 year forward int. global 2x2 deg. produces similar cost functions      - I/O and checkpointing
183        * all cost terms of new code are slightly smaller (0.5 to 3%).      - GM/Redi  *but*  using the GM/Redi coefficient of Salt
184      + Not implemented so far:
185  ecco_c44_e21      - KPP
186  o reversed order of ctrl/cost initialisation      - OBCS
187  o gmredi_slope_psi.F: Inconsistent types in MIN/MAX expressions    + No specific example supplied (yet) but global_ocean.90x40x15 has the
188  o modified to enable ctrl field swapping      necessary data.ptracer file. Simply use -enable=ptracers and uncomment
189  o bug fix for global sum of ncvarlength in ctrl_init      line in data.pkg.
190  o no more call to ctrl_unpack from ctrl_init    + This package is disabled by default since it increases storage.
191  o no more write of ctrl vector in ctrl_pack  
192    however, init. of control vector from ctrl_unpack  
193  o various bug fixes & consistent nobcs handling in  o testscript: added -cleanup option
194    ctrl_pack/unpack routines (by G. Gebbie)  o Shap_filt: a) new shap_filt S/R to use no-slip BC with S2 filter
195  o optim package: Added obcs control part for lsopt I/O               b) enable to filter 2D fields.
196    (by G. Gebbie)  
197  o missing call to obcs_diags  checkpoint44e_pre,post
198  o moved calls in packages_init_fixed to packages_readparms  o re-wrote convective adjustment to use pre-calculated weights
199    because some quantities are needed earlier.    Affects: model/src/convective_adjustment.F, convective_adjustment_ini.F
200      Adds:    model/src/convective_weights.F, convectively_mixtracer.F
201  ecco_c44_e20    Original calls to convect() still in place.
202  o before major merges from the Scripps setup  
203    here are two small updates:  chkpt44d_pre,post
204    - field swapping for surface flux controls enabled  o added missing EXCLUDE_MONITOR flags
205      (remains t.b.d. for atmos. state plus obcs)  o changed "e" to "_d" in gmredi_slope_limit, gmredi_slope_psi
206    - bug fix in active_file.F for xz-/yz- routines    (incompatible typ in MIN/MAX expressions caused problems
207    - added missing xz-/yz- routines to active_file_ad.F    on IBM SP3)
208    o in genmake added variable MAKEDEPEND
209  ecco_c44_e19    plus resetting for case SunOS
210  o merged necessary patches:  o added timer_stats.c routine for IBM SP3
211    - AD-related changes for GMRedi  o removed variables in dynamics
212    - fixes i KPP (delZ -> drF)  !!! I forgot to create tag chkpt44d_pre !!!
213    - hook to OBCS songe layer code in external_forcing  !!! But chkpt44d_pre = chkpt44c_post    !!!
214    - errorMessageUnit non-zero in eeboot.F  
215    - included GMREDI, EXACT_CONSERV to natl_box_adjoint verif.  o real fresh water flux implemented with non-linear free-surface.
216    
217  ecco_c44_e18  chkpt44c_pre,post
218  o sponge layer code for obcs was missing  o few fix (mask in shap_s2, EmPmR in external_field_load,
219    - added hook in external_forcing.F    USE_NATURAL_BCS in solve_for_P);
220    - added storing for TAF  o add arguments myIter & myTime to S/R obcs_calc & solve_for_P
221  o new routine cost_init_barfiles to initialise  
222    adjoint fields adtbar, adsbar, adpsbar  checkpoint44b_pre/_post
223    (needed in conjunction with corrected ladinit for active I/O)  o merge of relevant stuff from the ecco-branch:
 o cleaned up and corrected order in packages_readparms, ecco_readparms  
 o corrected and simplified weights for obcs controls  
   - weight fields simplified  
   - bug fix in ctrl_pack  
   - no more write of control vector in ctrl_pack  
 o i/j index bug fixes in cost_ssh_mean  
   
 ecco_c44_e17  
 o bug fixes in cost package  
   * wrong identifier for active I/O of monthly mean psbar  
     (cost_averagesfields, cost_ssh) xx_psbar_dummy  
   * cleaned all identifier for active I/O of monthly fields  
   * wrong ladinit in various cost routines  
 o merge with interval (chkpt44c_post - chkpt44d_post)  
   * added missing EXCLUDE_MONITOR flags  
   * changed "e" to "_d" in gmredi_slope_limit, gmredi_slope_psi  
     (incompatible typ in MIN/MAX expressions caused problems  
     on IBM SP3)  
   * in genmake added variable MAKEDEPEND  
     plus resetting for case SunOS  
     added timer_stats.c routine for IBM SP3  
   * removed variables in dynamics  
   
 ecco_c44_e16  
 o new ECCO checkpoint, the first since ecco_c43_e15  
   which is branched from MITgcm.  
   All changes were merged in between (see ecco-branch-mod? tags).  
 o last changes:  
   * assuming path for TAF and TAMC are known (e.g. in $path)  
   * genmake:  
     - added tools/xmakedepend to circumvene shortcoming of  
       SUN's makedepend (size too small)  
     - invoking ../tools/xmakedepend instead of makedepend  
     - changed compiler flags for case SunOS  
   
 ecco-branch-mod5  
 o mod's to reproduce open boundary run  
 o enable field swapping for OB's and climatology  
 o update store directives for OB control  
 o added option ALLOW_OBCS_PRESCRIBE  
 o call of getobcs now out of obcs_calc (param. list modified)  
   
 ecco-branch-mod4  
 o this is the second part of changes necessary to reconcile  
   the natl_box_adjoint verification (with and without grdchk).  
   >>>!!! CHANGES BETWEEN THE TAGS ecco-branch-mod3 !!!<<<  
   >>>!!!                      AND ecco-branch-mod4 !!!<<<  
   >>>!!! SHOULD BE MERGED INTO THE MAIN TRUNK      !!!<<<  
224    - genmake: removed $S64 overwrite for case SunOS    - genmake: removed $S64 overwrite for case SunOS
225    - pkg/exf: update and corrections for field swapping and obcs    - pkg/exf: update and corrections for field swapping and obcs
226    - pkg/ecco: parameter lists for the_model_main, the_main_loop    - pkg/ecco: parameter lists for the_model_main, the_main_loop
# Line 277  o this is the second part of changes nec Line 228  o this is the second part of changes nec
228    - pkg/autodiff: added flow directives for obcs, mdsio_gl_slice    - pkg/autodiff: added flow directives for obcs, mdsio_gl_slice
229                    updated checkpointing_lev... lists for obcs                    updated checkpointing_lev... lists for obcs
230    - model/src: minor changes in forward_step, plot_field    - model/src: minor changes in forward_step, plot_field
231                   added directive for divided adjoint in the_main_loop
232    - pkg/mdsio: added mdsio_gl_slice    - pkg/mdsio: added mdsio_gl_slice
233    
234  ecco-branch-mod3  o check parameters & config (chkpt44a_pre,post)
235  o this is the first part of changes necessary to reconcile  o OBC and NonLin_FrSurf.
   natl_box_adjoint verification (with and without grdchk).  
   This set of changes is relevant for ECCO only.  
   The next set/tag of changes (ecco-branch-mod4)  
   should be merged into the main trunk.  
   - update some options  
   - parameter lists of ecco_the_model_main and ecco_the_main_loop  
     differed from those of the MITgcm counterparts.  
   
 ecco-branch-mod2  
 o update ecco-branch-mod1 to checkpoint44  
 o exclude from updates:  
   - verification/  
   - adjoint/  
   - pkg/cost  
   - pkg/ctrl  
   - pkg/ecco  
   - pkg/obcs  
   - model/src/the_model_main.F  
   - model/src/the_main_loop.F  
236    
237  checkpoint44  checkpoint44
 o apply Shapiro filter on T&S before calc_PhiHyd (stagger time stepping)  
238  o fix bug in mom_vi_del2uv  o fix bug in mom_vi_del2uv
239    o select when filters are applied ; add options to zonal_filter (data.zonfilt)
240    o gmredi: fix Pb in the adiabatic form ; add options (.e.g. Bolus advection)
241    o update AIM experiments (NCEP input files)
242    o improve and extend diagnostics (Monitor, TimeAve with NonLin-FrSurf)
243  o added some stuff for AD  o added some stuff for AD
244    
 ecco-branch-mod1  
 o based on c43  
 o replaced packages cost, ctrl, ecco, obcs  
 o removed MITgcm verif.'s, added ECCO verif.'s  
 o added lsopt, optim for line search optimization.  
   
245  checkpoint43  checkpoint43
246  o added two new advection tests  o added two new advection tests
247  o added global 4x4 degree ocean test  o added global 4x4 degree ocean test

Legend:
Removed from v.1.53.4.32.2.3  
changed lines
  Added in v.1.86

  ViewVC Help
Powered by ViewVC 1.1.22