Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Links to checkpoint62s: | (view) (download) (annotate) |
Sticky Tag: |
undoing 1.50 to 1.51 modification
call OBCS_COPY_TRACER for advected seaice variable (AREA, HEFF, HSNOW, HSALT and ICEAGE) (should fix Pb when using high-order scheme with OB near tile edge)
For OBCS, close southern and western boundaries for sea ice solver based on myYGlobalLo and myXGlobalLo instead of OB_Js and OB_Iw.
- rename: SEAICE_DIAGS.h --> SEAICE_TAVE.h - use simpler (no level index) cumulative-time counter: SEAICE_timeAve(bi,bj) - cumulate fields only if doing time-ave output (SEAICE_taveFreq > 0)
o Adding I. Fenty seaice cost extensions (after minor cleanup) - smrsst, smrsss - based on modified seaice_cost_driver (not yet checked into CVS) - new runtime "clamp" parameters SEAICE_clamp_salt, SEAICE_clamp_theta o Adjust a few seaice STORE directives (following some recent modifs).
make maximum zeta a runtime parameter
add missing "_d 0" in HSNOW initialisation
add bj in exch2 arrays and S/R
Changed initialisation and re-initialisations of seaiceMaskU/V
cosmetic change: replace uVel(i,j,1,bi,bj) with uVel(i,j,kSrf,bi,bj) etc.
third and step of replacing 3D versions of UICE,VICE,HEFF,AREA by 2D versions. Keep your fingers crossed.
first steps in cleaning up the 3-time levels of UICE,VICE,HEFF,AREA: 1. introduce fields for time step (N-1 = *Nm1) 2. do not use levels 2 and 3 anymore (but they are still around) these changes probably break the adjoint, but the pickup is still working, so I am stopping now to give P. a chance to fix the adjoint if necessary and then I will get rid of the 3D fields altogether carefully taking into account the pickup (not clear how).
o adding JM's doc/diags_changes.txt: a very useful summary that I refer to all the time o adding pkg/seaice diagnostics for atmospheric and oceanic heat flux and for atmospheric and oceanic ice growth
- move initialization of grid variables from init_varia and init_fixed where possible
Adding stressDivergenceX/Y to store.
new header file "W2_EXCH2_SIZE.h" coming with new W2-Exch2 topology code
change macros (EXCH & GLOBAL_SUM/MAX) sufix _R4/_R8 to _RS/_RL when applied to _RS/_RL variable
forgot a few instances of 330., when replace this number by a runtime parameter SEAICE_rhoSnow. There is still one instance in cost_ice_test.F, which I am not touching
change a few defaults: B-grid flags (not tested in verification, so no effect) turn on flooding by default turn on advection of snow by default
add 4 coeffient-2D-fields for the metric terms. For now they are defined in SEAICE.h and computed in seaice_init_varia.F, eventually they should be moved to newly created SEAICE_GRID.h (together with other grid parameter fields) and computed in seaice_init_fixed.F. That will be a different step.
added sea ice age tracer; #define SEAICE_AGE in SEAICE_PARAMS.h
Make frWtrAtm global for FW balancing in presence of sea-ice, and I need it outside of seaice code. Also move DIAGNOSTICS_FILL(frWtrAtm ... (hope I'm not making a betise).
More _IF fixes.
next fix: put all previous changes within #ifdef SEAICE_ALLOW_TD_IF
More fixes.
Add Ian Fenty sea-ice thermodynamics modifications. Currently this is in separate routines to enable alternative testing. (relevant routines are ..._if.F) CPP switch is via #define SEAICE_ALLOW_TD_IF I.F. code doesnt contain latest SEAICE_SALINITY code
fixed a comment
complete the functionality of starting from pickup even if nIter0=0 (did not work so far because of to generous if statement)
removed uneeded arrays KGEO, GWATX, and GWATY from SEAICE_CGRID code
o pkg/seaice: removed SEAICE_FFIELDS.h and seaice_get_forcing.F seaice forcing fields can now be read only through pkg/exf
always initialise all variables in common block.
o Added open boundary conditions capability for seaice HEFF and AREA For time being this capability requires pkg/exf, pkg/obcs, and pkg/seaice UICE and VICE are masked at the southern and western open boundaries, if any, of the domain. Eventually needs to be extended to SALT and SNOW.
define SEAICE_salinity as a fraction of the model's surface level salinity
initial HSALT set to ICE2WATR*rhoConstFresh*SEAICE_salinity
replace hardcoded lower limit for zeta (zMin) by runtime parameter
adding ice salinity HSALT as a prognostic variable
changes for using pkg/seaice with pkg/obcs
Close southern and western boundaries for UICE/VICE when useOBCS=.TRUE.
added missing exchange for pkg/seaice OBCS mask
fix typo in latest check-in
added a missing #inclulde "OBCS_OPTIONS.h" (bug fix by An Nguyen)
oops
applied OBCS T/S mask to HEFFM sea ice mask This fixes pkg/seaice OBCS problem but in my opinion this information should eventually be transferred to the GRID.h arrays: hFacC, maskC, etc.
Adding AREA, HEFF, HSNOW as control variables.
- fixed problem with HSNOW initialization to avoid thick snow when AREA is initilized with small values - does not change lab_sea results. - added HsnowFile initialization capability Modified Files: doc/tag-index pkg/seaice/SEAICE_PARAMS.h pkg/seaice/seaice_init_varia.F pkg/seaice/seaice_readparms.F
initialise PRESS0, zmin, zmax as a function of AREA and HEFF
- make it possible to read AREA from a file, too.
make pLoad & sIceLoad always available (but still only used if ATMOSPHERIC_LOADING is #define).
- rearrange fields for evp a little to avoid having to call a: too many exchanges and b: exchanges for stress trensor components which we do not have.
o add new interface for atmospheric wind forcing (C-grid only) - new routine - two new field (taux, tauy for stress over ice) - removed one global field (DAIRN) for C-grid case only - small clean up - even remembered to update seaice_ad_diff.list
Initially set HSNOW to zero intead of 20 cm when HEFF is zero. This does not change lab_sea experiments, which are initialized with HEFF=1 everywhere, but it "will" change results for any experiment that are initialized with partial open water domains.
o fix multi-category seaice: - change cpp flag SEAICE_MULTILEVEL to more meaningful name: SEAICE_MULTICATEGORY - fix short wave heat flux o replace field areaLoc by scalar variable
o split seaice_init into _fixed, _varia o add first version of diagnostics interface (may need improvements, additions, etc)
This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, select a symbolic revision name using the selection box, or choose 'Use Text Field' and enter a numeric revision.
ViewVC Help | |
Powered by ViewVC 1.1.22 |