Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Links to checkpoint62r: | (view) (download) (annotate) |
Sticky Tag: |
rename iicekey as ticekey to avoid conflict with pkg/seaice
Some small changes to avoid "potential conflict" warning. Works ok on weddell under MPI with nPx=6
Fix several store dirs.
More modifs toward adjoint.
Initial changes for adjoint of pkg/thsice
fix syntax error
vectorization leads to major modifications - main i,j-loops have been split into many smaller loop sections in order to move k-loop outside of i,j-loops - a lot of local scalar variables removed and replaced by 2D variables - inlined thsice_reshape_layers (which has become obsolete, but not yet removed) - all CADJ statements are commented out, because they do not mean anything anymore - lot's of commented lines remain, I am not sure how many of those we want to keep.
avoid unbalanced quote (single or double) in commented line
change comments
- Now hsMax limits snow heitght by turning it into ice following flooding scheme (and now conserving energy) - Slight rewriting of the formula in thsice_calc_thick.F - Parameter rhowi eliminated to be replaced by floodFac = (rhosw-rhoi)/rhos --> changes results of global_ocean.cs32x15.icedyn and global_ocean.cs32x15.thsice
use continuous power-law function for partition of energy between lateral melting/freezing and thinning/thickening ; old code using step function still available with #undef THSICE_FRACEN_POWERLAW
prevent (or reduce) lateral melting if sea-ice fraction is too close to iceMaskMin
ensure area > iceMaskMin : Pb appears after changing (May-04) to melt only laterally if thin ice.
fix few little problems: in thsice_calc_thickn.F: 1) growth vertically (and not laterally) if iceFrac == iceMaskMax 2) melt laterally only if hIce < hThinIce (as the comments say) and in thsice_extend.F: 3) allow to form ice even when iceFrac == iceMaskMax (by increasing thickness) 4) start to form ice as soon as the minimum ice-volume is reachable.
if already ice, thickness of new ice to form is not larger than hNewIceMax
add description of how the partitioning between lateral and vertical melting/growth works.
rename few parameters: himin -> hIceMin himin0 -> hThinIce hihig -> hThickIce i0 -> i0swFrac transCoef -> bMeltCoef frac_energy -> fracMelting and add: hNewIceMax, fracFreezing, dhSnowLin
First set of modifs for TAF-ing thsice.
- put i,j loops inside S/R: THSICE_ALBEDO, THSICE_SOLVE4TEMP, THSICE_EXTEND and THSICE_CALC_THICKN - split thsice_step_fwd.F in 2 S/R: thsice_step_temp.F & thsice_step_fwd.F
get rid of "goto" statement
include header file EEPARAMS.h
improve "debug" writing.
maintain snow precipitation unchanged (even if Tsf >=0).
major changes in pkg/thsice: allows atmospheric model (AIM) to use thsice. - split thsice_therm.F in 2 S/R: thsice_solve4temp.F & thsice_calc_thickn.F - move most of the ocean & bulk_force interface in thsice_main.F - add a "slab ocean" component to be used with atmospheric model
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 |