/[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.1440 by jmc, Wed May 25 00:51:04 2011 UTC revision 1.1450 by jmc, Wed Jun 8 15:36:56 2011 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o model debug options:
8      - Switch reading of debugMode from S/R INI_PARMS (parameter file "data", stored
9        in PARAMS.h) to S/R EESET_PARMS (parameter file "eedata", stored in EEPARAMS.h).
10      - Refine debug-msg printing selection (3 new levels):
11        debugMode  : controls printing of debug messages (sequence of S/R calls).
12        debugLevel : controls printing of algorithm-relevant auxiliary/debug messages
13         debLevZero =  0 ; /* level of disabled aux. msg printing */
14         debLevA =  1 ; /* level of minimum  aux. msg printing */
15         debLevB =  2 ; /* level of low aux. print (report read-file opening)*/
16         debLevC =  3 ; /* level of moderate debug prt (most pkgs debug msg) */
17         debLevD =  4 ; /* level of enhanced debug prt (add DEBUG_STATS prt) */
18         debLevE =  5 ; /* level of extensive debug printing */
19        New debugLevel default: with default debugMode value (=F),
20          = 2 (debLevB) or = 1 (debLevA) in case ALLOW_AUTODIFF is defined.;
21          and if debugMode=T, default debugLevel = 4 (debLevD).
22        Relation with previous settings:
23         debLevB (=2) corresponds to former debLevA (=1)
24         debLevD (=4) corresponds to former debLevB (=2)
25        Note:
26         a) if debugLevel=0, skip WRITE_GRID call when nIter0 > 0 (same as before);
27         b) default debugLevel=1 is convenient in AD runs (avoid many file opening msgs);
28            to recover previous default, set debugLevel=2 (e.g. AD-test hs94.1x64x5)
29      - Add new parameter (printResidualFreq, read from PARM02 namelist) to print
30        detailed CG solvers convergence; off (=0) by default, except if debugLevel=5.
31    o pkg/obcs:
32      - change OBCSfixTopo default from True to False (should not change the results);
33        keep OBCSfixTopo=T in verification-exp. seaice_obcs.
34    o pkg/thsice:
35      - fix bug in thsice_calc_thickn.F (was introduced in version 1.18) in
36        2nd layer seaice sublimation (unlikely to happen)
37      - use Interior mask for OBCS in thsice_advection.F (similar to what is
38        done in gad_advection & seaice_advection).
39    o adjoint:
40      bug fix spotted by Gael Forget:
41      seaice_advection and thsice_advection still had
42      legacy CPP options to avoid CS corner updates for adjoint
43    o pkg/seaice:
44      => Change some of the defaults (2011/05/29):
45         - change SEAICEuseFluxForm default (from False to True)
46         - change DIFF1 default (from 0.004 m/s to 0. m/s)
47         And print warnings if DIFF1 is not set in data.seaice and centered
48          avection scheme is used without any diffusion.
49      also:
50      - fix seaice_diffusion calls in default AdvScheme part of seaice_advdiff.F;
51      - set SEAICEdiffKh[] default in the same way as SEAICEadvSch[]
52      - stop if DIFF1<>0 with SEAICEadvScheme <> 2, if DIFF1<>0 with any
53        other non-zero SEAICEdiffKh[] diffusivity, if SEAICEadvSalt=T with
54        #undef SEAICE_VARIABLE_SALINITY, if SEAICEadvAge=T with #undef SEAICE_AGE.
55      update results of test exp.:
56         lab_sea: output.fd (switch to AdvScheme=33 with diffKh=20)
57             output.hb87, output_adm.evp (SEAICEuseFluxForm=T and diffKh=200)
58             output_adm.noseaicedyn (SEAICEuseFluxForm=T)
59         offline_exf_seaice: output.seaicetd, output_adm (SEAICEuseFluxForm=T)
60         seaice_obcs: output (SEAICEuseFluxForm=T)
61      and prescribe old default for SEAICEuseFluxForm (=F) and DIFF1 (=0.004)
62          for standard fwd & AD lab_sea tests (to keep the same output).
63    o pkg/seaice/seaice_growth.F:
64      - use Ian Fenty's capping formulas for actual ice
65      thickness thoughout EVOLUTION branch.
66      - put treatment pathological case #2) in CPP brackets
67      (ALLOW_AVOID_INFINITESIMAL_AREA, undef by default).
68      - update global_ocean.cs32x15 and 1D_ocean_ice_colum results accordingly.
69    o pkg/seaice:
70      - seaice_diffusion.F : remove the min(DX,DY) factor permanently;
71      pass the Laplacian diffusivity (in m^2/s) diffKh as a parameter;
72      also pass fac as a param, which can be 1. (to update the tendency
73      field) or the time step (to update the field itself).
74      - seaice_advdiff.F: for multidim, call seaice_diffusion with
75      e.g. SEAICEdiffKhHeff and with fac=1; for single-dim, call
76      seaice_diffusion with e.g. SEAICEdiffKhHeff and with fac=time step.
77      - SEAICE_PARAMS.h, seaice_readparms.F, seaice_check.F :
78      hooks to handle SEAICEdiffKhHeff, SEAICEdiffKhArea, etc.
79    o pkg/seaice:
80      diffusion with centered advection scheme (diffus.F & advect.F):
81      - simplify S/R DIFFUS (put deltaT & DIFF1 factor outside, move EXCH outside)
82        to just compute the laplacian of input-field.
83      - unroll k loop (harmonic,biharmonic diff) in advect.F to save 1 exch call.
84      - fix flux-form version in diffus.F to conserve tracer.
85      update output of lab_sea.fd & lab_sea.hb87 (1rst set of modif affects
86       machine truncation) and lab_sea.salt_plume (conservative flux-form modif)
87    
88    checkpoint62y (2011/05/26)
89    o pkg/seaice:
90     - seaice_growth.F : in EVOLUTION branch, always include
91       a_QSWbyATM_cover in qnet (as it should be).
92     - verification/global_ocean.cs32x15 : update results accordingly
93       (output.seaice.txt output_adm.seaice.txt).
94    o pkg/seaice + obcs:
95     - use maskInW & maskInS in multi-dim advection and biharmonic diffusion:
96       equivalent of assuming no tracer gradient outside the OB region.
97       update output of seaice_obcs (uses DIFF1 <>0 => biharmonic diffusion)
98     - remove calls to OBCS_COPY_TRACER after EXCH calls (no longer needed
99       after changes above);
100     - obcs+seaice cleaning: apply seaice-tracer OB values to tracer field
101       only at OB location (instead of over a band of OLx width).
102    o pkg/seaice:
103     - fix wrong indices in old (non Flux-Form) diffusion default S/R diffus.F
104       (dxG(I+1, ...) changed to dxG(I, ...), no impact on Lat-Lon grid)
105  o pkg/obcs:  o pkg/obcs:
106    - split header file "OBCS.h" into 4 separated files:    - split header file "OBCS.h" into 4 separated files:
107      OBCS_PARAMS.h, OBCS_GRID.h, OBCS_FIELDS.h & OBCS_SEAICE.h      OBCS_PARAMS.h, OBCS_GRID.h, OBCS_FIELDS.h & OBCS_SEAICE.h
Legend:
Removed from v.1.1440  
changed lines
  Added in v.1.1450
  ViewVC Help
Powered by ViewVC 1.1.22