/[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.525 by edhill, Sat Jul 30 00:58:40 2005 UTC revision 1.628 by jmc, Mon Mar 13 04:44:25 2006 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    checkpoint58b_post
8    o AIM+thSIce: fix sign of Short-Wave through the Ice ;
9      update output.txt of exp. aim.5l_cs.thSI
10    o various mnc cleanups and improvements:
11      + shrink lookup tables by factor of ~4
12      + better error reporting when running out of lookup space
13      + able to handle longer path/file names (up to 500 chars)
14    o adjoint:
15      - Another overhaul of store dirs. for NLFS to eliminate  "hidden" recomputations.
16      - TBD: "hidden" mom_vecinv recomp. in dynamics
17    o change forcing_In_AB to affects both T,S forcing and Momentum forcing
18     (allow to differentiate between forcing components using new integer flags:
19      momForcingOutAB=1/0 & tracForcingOutAB=1/0)
20     and add new flag to put Dissipation tendency out of Adams-Bashforth.
21    o ecco
22      - enable trend removal from SSH anomaly obs.
23    o adjoint
24      - a fix for PmEpR when #undef EXACT_CONSERV
25      - fix of "hidden" recomp. in thermodynamics (no TAF warnings)
26      - try to fix "hidden" recomp. in forward_step for NLFS
27        (NLFS seems broken sometimes after I checked it in)
28    o seaice:
29     - add c-grid version of the dynamic solver:
30       + three new routines that are the c-grid substitute lsr, dynsolver
31         and ostres
32       + put a few fields that were local to dynsolver into global common
33         blocks, so that I can move the computation of stresses etc into
34         seaice_lsr (saves coding but may break the adjoint).
35     - replace more hardwired parameters by runtime parameters
36     - add ice masks that mask the rhs of the implicit solvers where there
37       is no ice (commented out in seaice_dynsolver, because i am not sure
38       if this works properly), eventually this should replace the clipping
39       of ice velocities in seaice_dynsolver to +/-40cm/s.
40    o remove the BUILD_INFO.h include from CPP_EEOPTIONS.h and explicitly
41      include it only where its needed
42    o move balancing of fluxes before the call of external_forcing_surf
43    o seaice:
44      - correct units of turning angles (now they are really in degrees)
45      - move computation of GWATX/Y out of main routine into dynsolver
46      - initialize GWATX/Y in seaice_init
47    o Apply mask on Fresh-Water flux (needed for SSH forcing).
48      update output.txt of exp. global_ocean_ebm & cpl_aim+ocn (affects cg2d residual)
49    o seaice: add runtime flags for turning angles (default to zero)
50    o Adding runtime parameter to exf/, ctrl/, to remove on-the-fly global mean
51      and trend from exf forcing and ctrl control adjustments
52    o fix variable aliasing in gad_calc_rhs that cannot be handled by TAF
53      (since checkpoint57x_post)
54    o change the default rotation period to 86164.s = the Earth rotation
55    o 1rst implementation of Implicit IGW using the 3-D solver (use3Dsolver=T)
56      and based on a reference stratification. Tested within hs94.cs-32x32x5.
57    o seaice advection: add code for choosing between the default (2nd order
58      central differences) and multiDim advection schemes provided by
59      pkg/generic_advdiff
60    o seaice:
61      - removed all references to SEAICE_GRID.h fields from
62        dynsolver, advect, diffus, lsr, without changing the results of lab_sea
63        (cubed-sphere configuration will change a little)
64      - grid computations in seaice_init are obsolete now.
65      - cleaned up code a little so that is easier to read (hopefully).
66      - added new (better?) discretizations of default advection and diffusion
67        (turn on with SEAICEuseFluxForm)
68      - replace 1 in 3d global fields by kSurface in growth.F
69      
70    checkpoint58a_post
71    o seaice:
72      - bug fix in advect.F
73    o shelfice:
74      - separate shelfice load anomaly from pload (breaks with time
75        dependent forcing), => introduce constant field shelficeLoadAnomaly.
76        Its default is 0., but it may be computed more cleverly from (unknown)
77        t- and s-profiles (tRef, sRef) and actual EOS. For now this has to be
78        done offline. A good approximation of the pressure load anomaly is
79        necessary to avoid large initial adjustment processes underneath
80        deep-reaching shelfice.
81      - small fix in ISOMIP shelfice_thermodynamics
82      - add more sophisticated thermodynamics of Hellmer (1989), this is
83        now the default
84      - give a few constants more intuitive names
85    o add code to balance EmPmP and Qnet at the end of do_ocean_physics. Useful
86      if bulk formulae are used in long integration (especially EmPmR). Turn
87      on with balanceEmPmR = .true. or balanceQnet = .true. in data, PARM01
88      if balancePrintMean, the imbalance that is substracted is printed to
89      STDOUT.
90    o fix bug (typo) in obcs_apply_eta.F
91    o add a new package shelfice which enables us to model the circulation
92      underneath shelf-ice with draughts of hundreds of meters.
93     - works
94     - to do: + include improved thermodynamics
95              + check whether kSurfC (from model/inc/SURFACE.h) can be used
96                instead of new field kTopC
97              + can we use topoFile, Ro_surf and existing code in ini_depth.F
98                instead of new field R_shelfIce and new code in ini_depth.F
99                and ini_mask_etc.F?
100              + adjust vertical mixing schemes?
101    o seaice pkg:
102     - add a cpp-flag to include a simple formula for a salinity dependent
103       freezing point of sea water at the surface
104     - change the contents of SEAICE_TEST_ICE_STRESS_1 so that zonal and
105       meridional stresses are now different (and hopefully correct)
106    o AIM pkg: change threshold on large-scale condensation to conserve Heat.
107       update results of aim.5l_LatLon test ; fix title & units of AIM diagnostics
108       and  change 3 diagnostics names (OSR,OSRCLR,KM -> TSR,TSRCLR,DRAG)
109    o pkg/diagnostics: add options and parameters for regional-statistics
110       diagnostics ; define regions by reading a region-mask from a file
111       (might add other options later)
112    o pkg/bulk_force:
113     - add option and parameters to use AIM surface-flux formulae.
114     - Change loading part: S/R BULKF_FIELDS_LOAD only take care of bulkf_fields ;
115       others forcing fields (defined in PARM05, parameter file "data") are loaded
116       from S/R EXTERNAL_FIELDS_LOAD, as usually (with useBulkforce = T or F).
117     - initialise all bulkf_fields in bulkf_init.F ; do in-lining of exf_bulkcdn.F
118     - use the right EXCH call for uwind,vwind (to work on CS-grid).
119     - re-arrange header files (move parameters from BULKF.h to BULKF_PARAMS.h)
120       and parameters (note: calcWindStress replaces .NOT.readwindstress).
121     - cleaner interface with pkg/thSIce (new S/R THSICE_GET_PRECIP).
122    
123    checkpoint58
124    o add "--embed-source" option to genmake2 which, when the supporting
125      tools can be compiled (as determined by genmake2) will turn on the
126      embed_files package which then embeds the entire MITgcm source code
127      (*.[fFhc] + Makefile) used for the build within the executable
128      - requested by CNH but off by default
129      - adds a paltry <9MB to the mitgcmuv executable in most cases
130      - only writes the output when useEMBEDSRC in eedata is true
131    o clean up 3-dim. relaxation codes:
132      * Fix I/O inconsistency in pkg/rbcs: replace precFloat32 by readBinaryPrec
133      * Remove 3-dim. relaxation code from pkg/exf (now use only pkg/rbcs)
134      * Thanks to Tom Haine for testing!
135    
136    checkpoint57z_post
137    o make CG3D solver compatible with Free-surface at k > 1 (p-coordinate):
138      affects output (@ truncation level) of exp. plume_on_slope ; updated.
139    o fix the pickup behavior for obcs in combination with ptracers
140    o fix a bug in ptracers_external_fields_load
141    o no AB-2 at the first iteration for Gw (consistent with Gu,Gv,Gt,Gs ...)
142      affects results of NH experiments (but all our NH tests have Gw[it=1]==0 !!)
143    o transfered surface NH pressure to eta field (if exactConserv).
144    o gmredi: fix bugs in slope tapering ldd97 (depthZ had the wrong sign)
145         and update lab_sea output.txt file.
146    o completing 3-dim. relaxation code: additional term
147         1/tau*( theta(i,j,k) - thetaStar(i,j,k) )
148      in S/R EXTERNAL_FORCING_T (file external_forcing.F)
149      * At compile time, in "EXF_OPTIONS.h" set
150          #define  ALLOW_CLIMTEMP_RELAXATION
151          #define  ALLOW_CLIMSALT_RELAXATION
152      * At runtime in "data", under &PARM03 set (e.g., units in seconds)
153          tauThetaClimRelax3Dim = 15552000.0,
154          tauSaltClimRelax3Dim  = 15552000.0,
155      * 3 cases:
156          1) Use a constant-in-time 3-dim. field
157          2) Provide 12 3-dim. fields corresponding to Jan., ... , Dec.
158          3) Real-time relax. to time-evolving prescribed field
159    o Adding Laplacian-type smoothness constraint on
160      initial and boundary (forcing) controls (ic, bc)
161      * enable at compile-time via
162        ALLOW_SMOOTH_IC_COST_CONTRIBUTION, ALLOW_SMOOTH_BC_COST_CONTRIBUTION
163      * enable at run-time via mult_smooth_ic = 1, mult_smooth_bc = 1
164      * Absolute scale somewhat unclear so far
165      * no area weighting was included since no other cost term does so far
166    
167    checkpoint57y_post
168    First step for a NLFS adjoint
169    o initially suppress rStar (new flag DISABLE_RSTAR_CODE)
170    o new init. routines for calc_r_star, calc_surf_dr
171    o still need to deal with ini_masks_etc
172    o testreport seemed happy
173    checkpoint57y_pre
174    
175    o obcs: always reset wVel & eta (to zero or prescribe it) at the OBs.
176    o add new package rbcs (for relaxing ptracers at boundaries/depth/surface etc)
177    o add new offset feature to offline pkg (plus minor cleanup)
178    o cfc fixes for final year of integration
179    o clean up GCHEM_OPTIONS.h/DIC_OPTIONS.h (all options should
180      be in GCHEM_OPTIONS).
181    
182    checkpoint57x_post
183    o fix testreport Pb in multiple test experiment.
184    o cleanup our namelist syntax in verification:
185      - make it more F95 standards compliant (some work still needed)
186      - allows the vast majority of our tests to work with gfortran
187    o Non-Hydrostatic:
188     - put all NH variables (formerly in DYNVARS.h & GW.h) in NH_VARS.h
189     - fix exchange calls for CS-grid using the new EXCH_S3D_RL.
190    o Allow to apply AB on T,S rather than on AB(gT,gS):
191     - implemented within #ifdef ALLOW_ADAMSBASHFORTH_3
192     - use the same arrays (gtNm,gsNm) to hold tracer field at previous
193       time-steps (if AB(T,S)) and tendencies (if AB(gT,gS)).
194     - need a "clever" read_pickup to restart from the other AB formulation
195    o move local commom bloc /GlobalLo/ (used for singleCpuIO) in EESUPPORT.h
196    o remove some unused variables (reduces number of compiler warnings)
197    o adjoint
198      * Disable useVariableViscosity for adjoint for time being
199        (there are a number of extensive recomputations)
200      * Initialize certain variables needed for adjoint
201    
202    checkpoint57w_post
203    o pkg/mdsio: use MAX_LEN_FNAM (instead of hard coded 80) for file-names
204                 fix the singleCpuIO RS version of MDSWRITEFIELD
205       and add a WARNING if using globalFiles=T with multi-processors.
206    o write grid-file(s) in only one format (MDS or MNC): default=MNC (if
207      using MNC); can force MDS output by setting writegrid_mnc=F in data.mnc
208    o tracer advection scheme:
209      add DST-2 (=Lax-Wendroff) & 1rst order upwind schemes (the way to get
210       those 2 schemes before, changing GAD_FLUX_LIMITER.h + to recompile,
211       was not very convenient).
212      add implicit version for DST-3 & DST3_Flux_Limit vertical advection;  
213       (further tests needed for the implicit DST & Flux-Limiter schemes).
214    o PTRACERS: added namelist parameter PTRACERS_ref in analogy to tRef
215      and sRef for convenience, handy, when you want to initialize
216      ptracers in a large domain with a constant non-zero value. Also,
217      untangle ptracers_init.F a little.
218    
219    checkpoint57v_post
220    o change dome exp. to use variable viscosity in side-drag terms.
221    o remove all the files from the cfc_offline test since they are now
222      incorporated into tutorial_cfc_offline which has been producing
223      identical results for a few days of testing
224    o add new parameter: sideDragFactor to enable half-slip-side BC
225      and change side-drag S/R to use viscosity from mom_calc_visc
226      (can switch back by setting sideDragFactor=0.); affects results at
227      truncation level; regenerate output for exp.: global_ocean.cs32x10(x3),
228       global_ocean_ebm, global_ocean_pressure, matrix_example & MLAdjust(x5).
229    o vorticity & strain are now computed without masking, to account for
230      no-slip-side BC in variable viscosity calculation ; mask is applied
231      anyway later on, before computing advection/coriolis terms.
232    o modify exp4 to test new obcs support for reading from external files
233      and for passive tracers
234    o fix behavior of ptracers_monitor when monitor output is written to a
235      netcdf-file; still not independent of monitor.F (writes in the same
236      file).
237    o change default value of tRef (put a uniform, meaningful value for ocean
238      & atmos) and initialize sRef (was missing, as Martin mentioned).
239    o Redefine viscC?Leith*.  Change the coefficients of Leith viscosity
240      to be proportional to enstrophy dissipation length.  Thus,
241       viscC?leith?_new=pi*(viscC?leith?_old)^(1/3)
242      This is now consistent with Smagorinsky and with the manual.
243    
244    checkpoint57u_post
245    o OBCS and PTRACERS: add open boundary support for passive tracers
246      - either use homogenous (pseudo) v.Neumann conditions or prescribe
247        OB-values from file; this is not different from the way theta and salinity
248        are treated
249      - however, Orlanski-radiation conditions are not supported, and the model
250        will stop if you use pTracers and Orlanski at the same time.
251      - beefed up the rountine obcs_external_fields_load: now only those open
252        boundary values are overwritten with values from files for which there
253        are really files, otherwise the OB-fields remain untouched. This makes
254        it possible to use different OBs at different ends of the domain (as
255        with EXF)
256    o reduce number of exch calls in shap_filt (s2).
257    o CS-grid: fix divergence contribution to Leith viscosity. Affects (@
258      truncation level) 4/5 MLAdjust test results (11&12 digits match).
259    
260    checkpoint57t_post
261    o test biharmonic viscosity + no-slip with exp: global_ocean.cs32x15.viscA4
262    o finish adding brackets (coriolis S/R, vecinv) in multi-terms sum to get
263      same truncation, independent of the face number: affects affects several
264      test. exp.: Update global_ocean.cs32x15 (+.thsice) & MLAdjust (all 5).
265    o fix bug in mom_v_sidedrag (biharmonic viscosity was counted 2 times);
266      note that biharmonic viscosity + no-slip was not tested in any experiment.
267    o fix bug in non-linear free-surf side-drag ; generate new output for
268      global_ocean.cs32x15(+.thsice) & global_ocean_pressure.
269    o make mnc honor the writeBinaryPrec flag for most kinds of output
270    o eesupp: split EEPARAMS into EESIZE/EEPARAMS and add a vector global sum
271    o minor change in averaging W (mom_vi_?_vertshear.F) that affects truncation
272      error; update output from global_ocean.cs32x15(+.thsice) & MLAdjust (all 5).
273    o add brackets in multi-terms sum to get same truncation, independent of
274      the face number (manage to get exactly same value of gU,gV at the cs-edge)
275      affects several test. exp. Update global_ocean.cs32x15 (+.thsice).
276    o Clean, new version of mom_calc_visc.  Adds a new feature, viscAhRemax
277      and viscA4Remax, which allow one to limit the minimum viscosity by specifying
278      a maximum grid-scale Reynolds number.  Also, changes all length calculations
279      in the case of grid anisotropy to geometric mean of dx**2 and dy**2.
280      Suggested values of viscosity parameters found in mom_calc_visc.
281      Many new viscosity diagnostics available.
282    
283    checkpoint57s_post
284    o Add clearing-house for viscosity calculations, mom_calc_visc and
285      correct numerous bugs in LeithD and Smagorinsky discretizations.  
286      mom_calc_visc now allows for Smagorinsky, Leith and LeithD to be
287      used simultaneously, and to be used in mom_vecinv and mom_fluxform.
288      Mom_vecinv can be used with standard viscous calculation (mom_vi_hdissip.F)
289      and with Strain-Tension viscous term calculation (mom_hdissip.F).  
290      The Strain-Tension form is now called with useStrainTensonVisc=.true.
291      ViscAhStrain and viscAhTension are now retired.
292    o diagnostics "VISCAHD" & "VISCA4D" replace "VISCAH" & "VISCA4".
293    o In response to JMC's bug report, genmake2 now ignores any soft-links
294      that it finds in the *current* directory when searching for files.
295      Soft-links found in any other source directories will still work as
296      they did before.
297    o add HAVE_SETRLSTK define and useSETRLSTK flag (eedata) which calls
298      a C routine to unlimit the stack size
299    o mnc-ify pkg/aim snapshot vars
300    o mnc-ify pkg/land as requested by Daniel
301    o various changes to mnc including:
302      - all files use the new "BASENAME[[.ITER].{t|f}NUM].nc" format
303      - output can now be grouped so that all files within a group
304          change the ITER portion of their names in lock-step together
305      - can now read ("global") PER-FACE (in addition to PER-TILE) files
306          and works with both EXCH1 and EXCH2 (but needs more testing)
307      - writing works for all verification test cases w/ g77 on Linux
308    
309    checkpoint57r_post
310    o make MNC variable initialization happen earlier
311    o add diagnostics for (almost) each term of the momentum Eq.
312    o remove SST,SSS,UVEL_k2,VVEL_k2 diagnostics.
313    o move pkg timeave out of pkg-groups "oceanic" & "atmospheric" (pkg/pkg_groups)
314    o remove MINIMAL_TAVE_OUTPUT option ;
315      add "dumpInitAndLast" to replace "nodumps" (and MINIMAL_TAVE_OUTPUT).
316    o move some basic S/R (used elsewhere in the code) from pkg/mdsio to eesupp/src
317      allows to compile without pkg/mdsio (& test successfully hs94.1x64x5)
318    o add MNC output for adjoint of EXF (adexf), SEAICE (adseaice)
319    o add time-dependent SST, SSS control
320    o cost terms in ecco/ modularized (new cost_generic.F)
321    o cost_averages routines are modularized (new cost_averagesgeneric.F)
322    o new infrastructure for seaice model vs. obs misfits in seaice/
323      * driver from the_main_loop.F is seaice_cost_driver.F
324      * first term is AREA vs. SMR observation misfit (ALLOW_SEAICE_COST_SMR_AREA)
325      * gradient check in Exp. 6a of natl_box_adjoint/
326    o natl_box_adjoint input_seaice for Lab. Sea now available online via wget
327      (see README file)
328    
329    checkpoint57q_post
330    o fix initialisation Pb in calc_gw (get NANs in exp5 with g77)
331    o fix annoying months-old bug in mnc/mdsio output for seaice pkg
332      - bug was reported by Jean-Michel and Ian Fenty
333    o fix initialization Pb in AIM-Ocean coupled model (aim_sice2aim.F)
334      update exp. cpl_aim+ocn output.
335    o Implement 4-level checkpointing: AUTODIFF_4_LEVEL_CHECKPOINT
336      Needed for high-res. runs with large global files
337      to be able to limit tapelev size
338    o Changed some tests for adjoint:
339      * carbon/ now tests AUTODIFF_2_LEVEL_CHECKPOINT
340      * natl_box_adjoint/code_bulk/ tests AUTODIFF_4_LEVEL_CHECKPOINT
341      * global_ocean.90x40x15/ back to default 3-level checkpointing
342    o Modifications to enable scripted COMMON to MODULE conversion
343      * verification/OpenAD/code/ has necessary scripts and description;
344      * converts, compiles and runs successfully;
345    
346    checkpoint57p_post
347    o reading of delXfile & delYfile moved from ini_parms.F to ini_grid.F
348    o move ${PKG}_MNC_INIT from ${PKG}_READ_PARAMS to ${PKG}_INIT_FIXED
349      (already the case for some pkgs, including recent MNC_init, e.g. thsice)
350    o call CALC_GW from DYNAMICS (instead of from FORWARD_STEP)
351    o Dissipation & phiHyd gradient are always added to gU,gV in timestep.F
352      (was already the case for dissipation with mom_vecinv,
353       and also the case for grad.PhiHyd when staggered-timeStep)
354      This will allow to put dissipation out-off the AB time-stepping.
355     Unfortunately, affects truncation error. Update output of exp. that "fails"
356      (cg2d matching number) and leave the others for later: front_relax(10),
357      global_ocean_pressure(10),ideal_2D_oce(10),internal_wave(9),lab_sea(12),
358      matrix_example(7!),natl_box(12),plume_on_slope(11),vermix(12x3,11x2).
359    
360    checkpoint57o_post
361  o fix MNC/MDSIO flags in PTRACERS so that (1) it works, and  o fix MNC/MDSIO flags in PTRACERS so that (1) it works, and
362    (2) its ready to be documented during the 20050801 DocFest    (2) its ready to be documented during the 20050801 DocFest
363  o MNC-ify the gmredi package -- tested and works w/ and w/o mnc  o MNC-ify the gmredi package -- tested and works w/ and w/o mnc
364    for ideal_2D_oce and lab_sea on linux_ia32_g77    for ideal_2D_oce and lab_sea on linux_ia32_g77
365    o set to #undef CTRL_SET_OLD_MAXCVARS_30
366      in global1x1_tot/code_seaice_adj to access indices 32, 33
367      (worth another tag)
368    
369  checkpoint57n_post  checkpoint57n_post
370  o adding precip control xx_precip 32  o adding precip control xx_precip 32

Legend:
Removed from v.1.525  
changed lines
  Added in v.1.628

  ViewVC Help
Powered by ViewVC 1.1.22