/[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.1576 by jmc, Fri Jan 20 01:30:45 2012 UTC revision 1.1596 by gforget, Mon Feb 13 23:45:41 2012 UTC
# Line 4  $Name$ Line 4  $Name$
4      Notes on tags used in MITgcmUV      Notes on tags used in MITgcmUV
5      ==============================      ==============================
6    
7    o pkg/seaice :
8      - simplify McPhee tapering, and allow to use it with SEAICEturbFluxFormula.EQ.
9      1 or 2, using newly added run time parameter SEAICE_availHeatTaper.
10      - fix d_AREAbyATM, d_AREAbyOCN, d_AREAbyICE diags I broke in r1.148.
11      - add permanent SItflux diag that corresponds to TFLUX but includes
12      ice+snow. Hence SItflux-TFLUX should match the ice+snow heat budget.
13      - allow activation/testing of a fix for suspected missing term in
14      ocn-ice heat budget (to be confirmed). To test this, you want to
15      undef SEAICE_DISABLE_HEATCONSFIX, and then set the run time param
16      SEAICEheatConsFix to .TRUE. that also allows the 'SIaaflux' diagnostic.
17      - no change in results, fwd or ad.
18    o global_ocean.cs32x15/results: three changes following from the pkg/seaice commit below.
19      - switch to variable freezing point.
20           Results change : cg2d -> 1 digits; ad -> 3-6-6 digits.
21      - bug fix : the sublimation term that was missing is now activated.
22           Results change : cg2d -> 2 digits; ad -> 8-6-6 digits.
23      - bug fix : flooding sign term for SIsal0.NE.0.
24           Result change : cg2d -> 8 digits; ad -> 11-12-9 digits.
25    o pkg/seaice:
26      - unified freezing point treatment : the old SEAICE_VARIABLE_FREEZING_POINT
27     is now the default, and the old default constant freezing point is recovered with
28     SEAICE_tempFrz0    = -1.96,
29     SEAICE_dTempFrz_dS = 0.,
30      - retire SEAICE_freeze that was the old way of specifycing the constant freezing point.
31      - remove ALLOW_SEAICE_FLOODING brackets; run time switch is already there.
32      - bug fix (thanks to M. Losch) : the sublimation term that was missing is now
33     activated. To reproduce old results that had this bug define SEAICE_DISABLE_SUBLIM.
34      - bug fix (silly me) : flooding sign term (seaice_growth.F r 1.149line 1211)
35     was wrong for SIsal0.NE.0. Changes cs32x15 results (see above item for detail).
36      - bug fix (thanks to O. Jahn) : area loss for melting  got messed up in
37     seaice_growth r149 for legacy branch (no results change).
38      - introduce SEAICE_CAP_SUBLIM : caps sublimation heat flux in solve4temp (code from I. Fenty).
39      - results did not change except for global_ocean.cs32x15, mostly due to the switch
40      to variable freezing point (see above item for detail).
41    o pkg/thsice:
42      - fix and simplify interface to salt_plume pkg ; diagnose directly amount of
43        seawater freezing in thsice_calc_thickn (not used).
44    o pkg/seaice: turbulent flux term specification via run time params.
45      - replace CPP switches with run time switches to facilitate testing
46      - retired CPPs : MCPHEE_OCEAN_ICE_HEAT_FLUX, GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR
47      - added runtime parameter : SEAICEturbFluxFormula
48      - move Mcphee constants to SEAICE_PARAMS.h
49      - seaice_readparms.F : read turbulent flux term params (SEAICEturbFluxFormula;
50        SEAICE_availHeatFrac/Frz, or SEAICE_gamma_t/_frz if apply); remove
51        seaice_summary call (moved to seaice_init_fixed)
52      - seaice_init_fixed.F : set SEAICE_availHeatFrac/Frz depending on
53        SEAICEturbFluxFormula, and seaice_summary call (moved from seaice_readparms)
54      - seaice_summary.F : add missing report of read turbulent flux term set up (and SITRACER).
55      - seaice_check.F : check range; inform and stop if retired CPPs are used (incl. SEAICE_AGE).
56      - seaice_growth.F : unified use of SEAICE_availHeatFrac/Frz, without
57            CPP switch, same computation as before for the different options.
58      - change results at the truncation level.
59    o pkg/seaice:
60      - remove un-used parameters LAD & MAX_TICE ; add new run-time parameters:
61        useMaykutSatVapPoly (default=F), postSolvTempIter (default=2) and
62        SEAICE_wetAlbTemp (default=-1e-3). Remove option SEAICE_SOLVE4TEMP_LEGACY.
63      - update verification/*/input*/data.seaice:
64        to reproduce previous results with former #defined SEAICE_SOLVE4TEMP_LEGACY
65        code, set useMaykutSatVapPoly=T & postSolvTempIter=0 ;
66        and to reproduce previous results with former non-SOLVE4TEMP_LEGACY code,
67        set SEAICE_wetAlbTemp=0 & SEAICE_snowThick=0.
68    o pkg/seaice/seaice_solve4temp.F (SOLVE4TEMP_LEGACY):
69      - remove LEGACY code for solving for tsurf (A1,A2,A3) but maintain the same
70        algorithm (same choice: useMaykutPolySatVap=T, postSolvTempIter=0);
71        difference in results only due to machine truncation.
72      - remove MAX_TICE (tsurf is always =< TMELT anyway); keep MIN_TICE if using
73        MaykutPolySatVap; keep MIN_LWDOWN.
74      - adapt SEAICE_MODIFY_GROWTH_ADJ code (untested) to non-legacy formulation.
75      - update output of verification exp. lab_sea (4 fwd: 11 cg2d digits except hb87,
76        + 3 AD: 6 digits except evp) and offline_exf_seaice.seaice (12 digits).
77    o pkg/salt_plume:
78      - isolate saltPlumeFlux array exch into a specific S/R called from
79        do_oceanic_phys.F
80      - add interface to pkg/salt_plume (filling of saltPlumeFlux) in pkg/thsice.
81    o verification:
82      - move untested set-up to MITgcm_contrib/verification_other.
83    o pkg/seaice:
84      - retire old way of SEAICE_AGE. Now replaced with ALLOW_SITRACER.
85    o pkg/seaice:
86      - seaice_solve4temp.F (keeping results unchanged):
87        * add 2 future run-time params to select which saturation vapor pressure
88          formula is used and how post-iteration fluxes are set.
89          Temporary set these 2 params according to option SOLVE4TEMP_LEGACY;
90        * New piece of code (not used) to update surf. atmos. heat flux F_ia using
91          the linearized expression to stay consistent with tsurf finding.
92      - seaice_freedrift.F avoid floating exception in sqrt(y) since y=sqrt(x*x)-x
93        can be < 0.
94    o pkg/seaice/seaice_solve4temp.F (SOLVE4TEMP_LEGACY undef):
95      - remove from F_ia (and from it's derivative dFiDTs1) the contribution of
96        conductive heat flux F_c , and add it explicitly when updating tsurf as
97        solution of: Fc = Fia + d/dT(Fia - Fc) *delta.tsurf
98        so that now F_ia has a consistent meaning through the entire routine.
99        This causes differences at truncation level: update output of fwd exp.
100        global_ocean.cs32x15.seaice.
101    o pkg/seaice:
102      - lots of editing/cleaning in seaice_solve4temp.F + better documentation
103        of arrays, sign/direction of fluxes; no effect on results.
104      - replace hard coded 273.16 (=previous value of celsius2K) by celsius2K
105        in seaice_budget_ocean.F & seaice_solve4temp.F (SOLVE4TEMP_LEGACY part);
106        note: this should have been done before changing celsius2K default.
107      - to keep results unchanged, add back old value of celsius2K in verification
108        exp. input*/data : lab_sea (AD), offline_exf_seaice (fwd+AD) & seaice_obcs
109      - replace hard coded seawater freezing temp in seaice_solve4temp.F (part
110        SOLVE4TEMP_LEGACY) of 271.20 (=previous value of celsius2K + default
111        value of SEAICE_freeze) by celsius2K+SEAICE_freeze so that seawater
112        freezing temp is consistent with value in seaice_growth.F ; this causes
113        differences at machine truncation level: update lab_sea AD output (3/4).
114    o pkg/seaice:
115      - mask (in exf_init_fixed.F) atemp,aqh,evap,precip,runoff,lwflx,areamask
116        and climsst/sss even when using seaice (to prevent calculation that might
117        use SST on land); does not affect any verification exp.
118      - apply mask to taux,tauy (in seaice_get_dynforcing.F); update output
119        of lab_sea.fd & lab_sea.hb87 and global_ocean.cs32x15.seaice_dynmix.
120    o pkg/seaice/seaice_growth.F: add mask to avoid non-zero ice thickness
121      over land (non-legacy code)
122    o vectorize salt_plume package:
123      - add extra code to salt_plume_frac and salt_plume_tendency_apply_s
124      - this code computes plumebk everywhere (also over land), thus requires
125        code to avoid divisions by zero and catch the case of kLev=Nr.
126      - efficiency of salt_plume_frac is still sub-optimal
127    o tools (genmake2):
128      - fix makefile NAME specification (option: -makefile NAME) by adding
129        "-f $MAKEFILE" in few places where it was missing including makedepend
130        command (contribution from Paul).
131    
132    checkpoint63i (2012/01/24)
133    o model/src:
134      - add factorized versions of some EOS code to find_rho.F,
135        define USE_FACTORIZED_EOS to turn it on (default for TARGET_NEC_SX)
136    o model/src:
137      - add celsius2K to namelist PARM01; remove unused recip_rhoNil from PARAMS.h;
138      - change default celsius2K from 273.16 to 273.15
139      - to keep results unchanged, add back old value of celsius2K in verification
140        exp. input*/data : aim.5l_cs, 1D_ocean_ice_column, lab_sea (fwd),
141         global_ocean.cs32x15(seaice) & cpl_aim+ocn(input_atm).
142    o model/src: add a little initialisation trick to find_rho.F and find_alpha.F
143      to make TAF generate vectorizable code
144    o pkg/seaice/seaice_evp.F: add code (copied from CICE) that avoids underflows.
145       Turn on by defining SEAICE_EVP_ELIMINATE_UNDERFLOWS. Most compilers
146       have flags that do this more efficiently.
147    o pkg/diagnostics/diagstats_calc.F
148      - for TARGET_NEC_SX fix the treatment of the scaling factor tmpFac=scaleFact
149  o pkg/obcs (with exch2):  o pkg/obcs (with exch2):
150    - derive OB indices in overlap region from OB-InsideMask (in obcs_init_fixed.F)    - derive OB indices in overlap region from OB-InsideMask (in obcs_init_fixed.F)
151      and remove code (in obcs_readparms.F) which was based on EXCH of indices      and remove code (in obcs_readparms.F) which was based on EXCH of indices
Legend:
Removed from v.1.1576  
changed lines
  Added in v.1.1596
  ViewVC Help
Powered by ViewVC 1.1.22