/[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.785 by dimitri, Wed Apr 11 00:00:47 2007 UTC revision 1.869 by jmc, Fri Sep 7 01:16:11 2007 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o testreport:
8      Allow to choose (per experiment) the list of tested output variables
9       and the main variable (used to issue a pass/FAIL ) ;
10      - for now, done by adding the specific list in file: [exp]/input/tr_checklist
11      - default is unchange: >cg2d< , (T,S,U,V,5 tracers)_[min,max,mean,StD] ;
12      Also modify printing of summary:
13       a) ">" NN "<" indicates the "main" variable (used for pass/FAIL)
14       b) echo the list of tested variables in the individual summary.txt of
15          each experiment.
16    o new S/R GLOBAL_SUM_TILE to compute global sum independently of number of
17      threads and number of processes (with #define GLOBAL_SUM_SEND_RECV)
18      (but still depend on the tile decompostion of the domain).
19    o pkg/seaice: change a discretization detail to make the B-grid code
20      "more symmetric". This change is likely to remove all "ridges/leads" that
21      where quasi-aligned with the grid in B-grid solutions. It may also be
22      the solution to the stress coupling issues with the B-grid.
23    o each tested exp: add a "run" directory with a standard ".cvsignore" file ;
24      remove .cvsignore from input dir. (since testreport no longer run here).
25    o pkg/thsice:
26     - cleaning: move common block out of "THSICE_2DYN.h" (removed) to "THSICE_VARS.h"
27     - fix reshaping of sea-ice just after advection:
28       melt all the ice only if ice-volume is too small.
29     - re-visit lateral/thinning repartition of melting (modifs from May 04, 2007):
30       prevent lateral melting if sea-ice fraction is too close to iceMaskMin ;
31     and update experiment global_ocean.cs32x15.icedyn (+start from pickup files).
32    o external_forcing_surf : extend the use of salt_EvPrRn & temp_EvPrRn (fresh-water
33      tracer content) to linFS or no RealFreshWaterFlux ; default remains unchanged.
34    o define mass2rUnit & rUnit2mass as units conversion factor, used for surface
35      forcing. mass2rUnit: from mass per unit area [kg/m2] to r-coordinate unit
36      (z-coord: =1/rhoConst ; p-coord: =gravity) ; rUnit2mass is just the reverse;
37      and remove horiVertRatio (& recip_horiVertRatio). Affects truncation error
38      of exp. global_oce_in_p: 10 16 16 14 16 13 16 16 16 13 13 13 16 13 13 12 16
39      -> update output.txt
40    o pkg/cd_code: add brackets in multiple terms sum (S/R CD_CODE_SCHEME):
41     this affects truncation error: -> update output.txt of exp:
42      global_ocean.90x40x15 (fw&adm), lab_sea (fw&adm) & lab_sea.lsr,
43      ideal_2D_oce, isomip.htd, natl_box, tutorial_global_oce_latlon.
44    
45    checkpoint59f
46    o pkg/autodiff: zero_adj resets an ajoint variable to zero (in adjoint mode)
47    o pkg/ctrl: ctrl_bound  forces parameter corrections to stay in
48      a bounded range (in forward mode)
49    o pkg/ecco, ctrl, autodiff and mom_common:
50     internal parameter controls, completements of implementations
51    o add few EXCH S/R (UV_BGRID, UV_DGRID, SM) for B & D-Grid velocity & 2nd Moment Sxy.
52    o use UV_BGRID exch for dxV,dyU: for now, only with pkg/exch2 & forward ;
53      (not seem any difference in testing results).
54    o pkg/cd_code: use the right EXCH calls (for CS-grid); but problem remains
55       at face corners.
56    o pkg/generic_advdiff:
57     - in GAD.h, define as parameters (iMinAdvR,iMaxAdvR,jMinAdvR,jMaxAdvR)
58      loop range indices for computing vertical advection tendency
59      and do the calculation in the interior only (no overlap) ;
60      for now, only implemented for SOM advection (to save CPU time).
61     - add argument "withSigns" to S/R FILL_CS_CORNER_TR_RL (needed for SOM_xy moment)
62    o exch2 adjoint
63      First set of changes:
64      - add exch2 hand-written adjoint templates
65      - exch2 hand-written no longer use TAF-flag 'nownew_arg'
66        This will require cleanup in pkg/autodiff/*.F
67        (changed routines currently kept separate in verif. code/)
68      - Comment all relevant #ifndef ALLOW_AUTODIFF_TAMC
69        that used to hide exch2 or cubed-sphere specific code
70        (commented via 'cph-exch2')
71    o pkg/seaice
72      - Close southern and western boundaries for UICE/VICE when useOBCS=.TRUE.
73    o pkg/bulk_force: fix loading into Qsw of Short-Wave flux; and update
74       output of global_ocean.cs32x15.thsice .
75    
76    checkpoint59e
77    o exch2: new S/R that only use the cube-tracer (center-position) exchange:
78      for B-grid vector (not yet tested) and C-grid vector (<- only used
79      when option W2_USE_R1_ONLY in W2_OPTIONS.h is defined).
80    o starting implementation of Duffy et al. (GRL 1999) salt plume scheme
81      (#ifdef ALLOW_SALT_PLUME)
82    o gather_2d/scatter_2d: removed halo regions
83    o gmredi: add option (GM_taper_scheme='fm07') for Ferrari & McWilliams 2007
84      scheme and the corresponding test case (under front_relax dir). For now,
85      only available with Skew-flux form and not yet tested in realistic set-up.
86      adjoint: store directives updated, but gradient checks deteriorated,
87      so maybe GM/Redi adjoint is now broken.
88    o seaice adjoint
89      * add AREA, HEFF, HSNOW as control variable
90        (should not change lab_sea results)
91      * CTRL_OPTIONS.h needs new flag: CTRL_SET_OLD_MAXCVARS_40
92    o pkg/seaice:
93      * added computation of saltFlux in seaice_growth;
94        For time being, sea ice salinity is assumed constant
95        contribution of snow flooding to freshening of sea ice is neglected.
96        SEAICE_salinity is set in SEAICE_PARM01 of data.seaice and it
97        defaults to 0.0 in order to maintain backward compatibility.
98        Typical values for SEAICE_salinity are 4 to 10 g/kg.
99      * replaced dRf(1) with dRf(kSurface)*hFacC(kSurface) in seaice_growth.F
100      * added SEAICE_availHeatFrac, which specifies fraction of surface level
101        heat content used to melt or grow ice
102      * commented out cap of HICE: it helps keep ice thinner
103      * reintroduced cap of HEFF via MAX_HEFF and SEAICE_CAP_HEFF option
104      * applied OBCS T/S mask to HEFFM sea ice mask
105    
106    checkpoint59d
107    o standard name: {PKG}_DIAGNOSTICS_STATE to fill diagnostics array
108     for {PKG} state variables, and called from DO_STATEVARS_DIAGS.
109     (keep {PKG}_DIAGNOSTICS_FILL for non-state variables diagnostics)
110    o pkg/gmredi (gmredi_readparms.F): add some missing "_d 0"
111     and change default value of parameter "GM_Small_Number";
112     update output of exp: lab_sea (x2), cfc_example, global_ocean.90x40x15,
113       global_ocean.cs32x15 (x4), global_ocean_ebm, global_with_exf,  
114       tutorial_global_oce_biogeo, tutorial_global_oce_latlon
115    o pkg/ecco: cost function when applying pkg/smooth to controls
116    o pkg/ctrl: add forcingPrecond run time parameter
117    o pkg/seaice:
118     - replace ZMIN=4.e8 with ZMIN=0. in dynsolver.F as in C-grid code
119     - fix bug in metric terms of u-equation in LSR solver (implicit part),
120       changes verification/lab_sea
121     - rearrange implicit part of u-equation in LSR solver to make it easier to
122       read, understand, and debug. This causes changes at the truncation level.
123     - fix (serious) bug in the boundary conditions (sign error in the implicit
124       part) of the LSR solver. Changes all LSR results.
125     - fixed problem with HSNOW initialization to avoid thick snow when AREA
126       is initilized with small values - does not change lab_sea results.
127     - added HsnowFile initialization capability
128    o pkg/smooth: addition to model/src and pkg/ctrl so allow
129     initialization of pkg/smooth and application to control vector
130    o pkg/seaice: fix another bug in the LSR no-slip boundary conditons
131    o dependencies: pkg/profiles requires pkg/cal
132    o pkg/profiles:
133     - added option ALLOW_PROFILES_GENERICGRID
134         -> pkg/profiles can thus be used with any grid (e.g. cs32x15) assuming
135         that interpolation information (grid points and coefficients) is provided
136         within the netcdf input files.
137         -> sample matlab scripts can be found in MITgcm_contrib/gael
138         that allow you to prepare netcdf input files.
139     - also modified: cleaner include statements, cleaner error print statements,
140       cleaner stops when error, and more tests of the netcdf input files consistency.
141    o additional test (in global_ocean.cs32x15) using seaice+thsice pkgs.
142    o pkg/seaice:
143     - move seaice-diagnostics from seaice_do_diags to do_statevars_diags; this
144       requires a new routine: seaice_diagnostics_fill.F
145     - make sure that seaice does not move if not forced (pressure replacement
146       method, Hibler and Ip, 1995), changes results
147     - and new damping criterion for EVP solver according to Hunke, JCP, 2001
148     - SEAICE_evpDampC = > 0. turns it on (615. is the recommended value).
149       It is turned off by default (=-1.), because it changes the results.
150     - make it possible to read AREA from a file, too.
151    
152    checkpoint59c
153    o add S/R and diagnostics for oceanic mixed-layer depth
154    o move GMRedi call after all vertical mixing schemes (including KPP)
155    o pkg/thsice:
156     - ensure iceMask > iceMaskMin : this problem appeared after May-04 modifs
157      (melt only laterally if thin ice): update output of global_ocean.cs32x15.thsice
158    o pkg/obcs:
159     - bug fix for balance code
160    o pkg/seaice:
161      - allow values of SEAICE_elasticParm > 1
162      - add new parameter SEAICE_evpTauRelax as an alternative
163        to SEAICE_elasticParm
164      - make lsr boundary conditions free slip (by default), changes result
165      - add code for no slip boundary conditions for lsr
166      - set minimum of zeta to zero (and not 4e8) (changes results)
167      - add a cap on AREA after advecting AREA if seaice_growth is not called
168    o change calculation of Visbeck K:
169      - no longer depend on tapering scheme; instead, compute slope (limited
170        by GM_maxSlope). this is equivalent to what was comuted before with
171        clipping or linear tappering.
172      - fix vertical averaging (wrong by 1/2 level)
173      - leave the old version (but will not be maintained) within
174        #ifdef OLD_VISBECK_CALC / #endif
175      - update output of exp. ideal_2D_oce
176    o change error catch in calc_r_star to improve vectorization
177    o exf runoff and seaice: remove ALLOW_SEAICE from exf pkg files and
178      add #define ALLOW_RUNOFF in SEAICE_OPTIONS.h (when not using exf).
179    o adjoint:
180      Re-organized adjoint checkpointing according to Matt Mazloff
181      (but with modifs/completions/cleanups)
182      Storing is now bundled into large arrays to optimized I/O
183      (and second step in cleanup of the_main_loop)
184    
185    checkpoint59b
186    o coupler (pkg/atm_ocn_coupler, pkg/ocn_compon_interf, pkg/atm2d):
187      - export/import (between ocean & coupler) fields needed for ocean DIC.
188    o pkg/seaice:
189      - get rid of three more exchanges seaice_calc_viscosities (required
190        different loop boundaries in seaice_lsr and seaice_ocean_stress)
191      - bug fix in seaice_ocean_stress (update of lab_sea required)
192      - reduce loop boundaries to necessary minimum also in seaice_evp
193      - go back to the old way of computing deltaC/Z as the adjoint likes that
194        better
195    o Refresh of tutorial_global_oce_latlon (includes changing salinity-like
196      passive tracer to an Age Tracer) + update output.txt
197    o pkg/exf:
198      - implement A-grid / C-grid selection for wind-stress input files.
199      - clean-up exf_wind.F (remove us, use wspeed instead ; add wStress)
200      for now, in exf_bulk_largeyeager04.F only :
201      - change implementation for case where windspeed + 2. components wind-stress
202        are read-in. AND use rhoConstFresh to convert evap.
203    o pkg/seaice:
204      - use correct EXCH_UV for wind-forcing (on A-grid).
205      - change SEAICE_OPTIONS.h to limit the number of option combinations (regarding
206        SEAICE_EXTERNAL_FORCING & SEAICE_EXTERNAL_FLUXES); add corresponding stops.
207    o pkg/exf:
208      - default EXF options are now set in EXF_OPTIONS.h (ifndef ALLOW_AUTODIFF_TAMC)
209      EXF_INTERPOLATION:
210      - merge global_with_exf/code version of exf_interp_read.F
211        to the main code ;
212      - Add new option: EXF_INTERP_USE_DYNALLOC to use Dynamical Allocation when
213        reading the file ; Important: This options needs to be defined to recover
214        previous code.
215        otherwise, use fixed size array to read in (size=exf_interp_bufferSize).
216      - no longer loose real*8 precision along the way when using exf_iprec=64.
217    o pkg/seaice: rearrangement of common blocks (really needs more cleaning up)
218      - change discretisation of evp-solver in order to avoid another exchange
219        (changes lab_sea) and clean up a little
220      - replace SQRT(MAX(deltaC,SEAICE_EPS)) by MAX(SQRT(deltaC),SEAICE_EPS)
221        help the adjoint (this also changes the results slightly)
222      - add a few "_d 0" where they had been missing in seaice_evp.F
223      - cleaning up: remove unused seaice_calc_rhs.F
224    o pkg/exf: some cleaning-up:
225      - change various "constants" into runtime parameters
226      - cleaned-up Large&Yeager04 routine which should eventually
227        replace exf_bulkformulae.F (changed names of S/R and CPP)
228        and merged various ALLOW_ATM_WIND options
229     only used with thsice pkg (in thsice_get_exf.F):
230      - implement fixed turbulent transfert Coeff over ice (as in Large&Yeager04)
231        which are used when useStabilityFct_overIce=TRUE (= the default)
232      - change default values of ice & snow emissivity and scale downward LW
233        by ice & snow emissivity.
234    o pkg/thsice:
235     fix few little problems, in thsice_calc_thickn.F:
236      - growth vertically (and not laterally) if iceFrac == iceMaskMax
237      - melt only laterally if hIce < hThinIce (as the comments say)
238      and in thsice_extend.F:
239      - allow to form ice even when iceFrac == iceMaskMax (by increasing thickness)
240      - start to form ice as soon as the minimum ice-volume is reachable.
241     update output of test-exp: aim.5l_cs(thSI), global_ocean.cs32x15(thsice),
242      and cpl_aim+ocn
243    o pkg/kpp
244     - move computation of surface related input fields to KPP into a new
245       subroutine kpp_forcing_surf.F
246     - remove wrong (within bi,bj loops) EXCH calls.
247     - little cleaning (add argument myThid in last position, cvs header ...)
248     - for now, fix conservation problem by exchanging KPPviscAz (may not
249       be the best fix)
250    o pkg/seaice
251     - check for inconsistent evp parameter/cpp-flag combination
252     - remove leftover errIO flat for obsolete name list error handling
253    o OpenAD
254     - changed verification/OpenAD to full baroclinic setup (including GM/Redi)
255    
256    checkpoint59a
257    o pkg/seaice: some parameter clarifications and a little bit of clean up
258    o pkg/thsice: rename run-time parameters:
259       himin       -> hIceMin
260       himin0      -> hThinIce
261       hihig       -> hThickIce
262       i0          -> i0swFrac
263       transCoef   -> bMeltCoef
264       frac_energy -> fracMelting
265      add new parameters: hNewIceMax, fracFreezing, dhSnowLin
266       with default values (for now) that do not change the results.
267      remove old version of thsice vars advection in pkg/seaice.
268    o in FFIELDS.h:
269      - replace #ifndef ALLOW_EXF  by #ifndef EXCLUDE_FFIELDS_LOAD
270      - make pLoad & sIceLoad always available (but still only used if
271         ATMOSPHERIC_LOADING is #define).
272    o pkg/kpp: replaced _KPP_RL with _RL
273    o bug fixes for vertical diffusivity computations when both KPP and
274        3D diffusivity arrays are used.
275    o pkg/seaice: fix bug in evp solver, while doing that
276      - change (improve?) discretization of stress computations for both
277        cgrid-lsr and evp solver
278      - add a new routine to compute strain rates
279      - the above changes lab_sea results
280      - add no slip conditions for evp solver (off by default), lsr-version
281        may follow
282      - fix exchange for shear at Z-points (does not change lab_sea)
283      - retire SEAICEuseEVP, instead turn on EVP by setting SEAICE_deltaTevp
284        (to force the user to pick a time step since there is not "safe" way
285        of choosing it)
286      - add principal stress diagnostics
287      - rearrange fields for evp a little to avoid having to call a: too many
288        exchanges and b: exchanges for stress trensor components which we do
289        not have.
290      - modify Hibler+Bryan (1987) stress coupling so that it reuses information
291        of evp, complete rewrite for lsr case (latter is does not yet work
292        in all cases)
293    o pkg/kpp: add code to have mixing in shelf ice caverns
294      - add myThid to all kpp routines (long overdue)
295    
296    checkpoint59
297    o pkg/exf:
298      - fix wind-stress location: defined on A-grid (if ALLOW_ATM_WIND &
299        ALLOW_BULKFORMULAE or USE_EXF_INTERPOLATION), otherwise, defined on C-grid
300      - move exf header files from lower_case.h to UPPER_CASE.h ;
301        add missing cvs Header & Name
302      - add adjoint routines and flow dirs for newly used exch_uv_agrid routines
303      - add namelist param. selectStressGridPosition, exf_monFreq
304      - retire many _CLIM routines, and remove data.exf_clim
305        + params moved to data.exf
306        + climsst, climsss now handled by exf_set_gen, i.e.
307          climss?period=0 refers to constant field,
308          climss?period=-12 to monthly forcing
309      - split up namelists in data.exf
310    o pkg/seaice:
311      - remove seaice_exch and seaice_exch_uv
312      - remove superfluous exchange of theta
313      - add myThid to seaice_budget_ocean/ice
314      - add new interface for atmospheric wind forcing (C-grid only)
315        + new routine
316        + two new field (taux, tauy for stress over ice)
317        + removed one global field (DAIRN) for C-grid case only
318    o pkg/thsise
319      - 2nd set of modifs for thsice adjoint in conjunction with exf
320    o enable runtime flag useEXF (useCAL is set to .T. if useEXF=.T.)
321    
322    checkpoint58y_post
323  o Added capability for latitudinal dependence of Bryan and Lewis, 1979  o Added capability for latitudinal dependence of Bryan and Lewis, 1979
324    vertical diffusivity, similar to that in MOM4; vertical diffusivity is    vertical diffusivity, similar to that in MOM4; turned on using
325    specified using diffKrBL79* diffKrBLEQ* and KbryanLewisLatTransition    #define ALLOW_BL79_LAT_VARY in CPP_OPTIONS.h and specified using
326      diffKrBL79*, diffKrBLEQ*, and BL79LatVary in runtime data file.
327  o Added capability for 3D specification of background vertical diffusion  o Added capability for 3D specification of background vertical diffusion
328    coefficient.  This is done by (i) #define ALLOW_3D_DIFFKR    coefficient.  This is done by (i) #define ALLOW_3D_DIFFKR
329    in CPP_OPTIONS.h and (ii) specifying file diffKrFile in runtime    in CPP_OPTIONS.h and (ii) specifying file diffKrFile in runtime

Legend:
Removed from v.1.785  
changed lines
  Added in v.1.869

  ViewVC Help
Powered by ViewVC 1.1.22