Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Links to checkpoint61s: | (view) (download) (annotate) |
Sticky Tag: |
Adding stressDivergenceX/Y to store.
One temp. store.
Add TAF option "kind" (or adjust "byte") to enable real*4 common blocks
Add timers for KPP_CALC and GGL90_CALC
Add some headers and (re-)init.
AD: back to version 1.72 ; + compute rhoInSitu if fluidIswater, otherwise set it to zero.
try to put back "if fluidIsWater" (needed for hs94 AD test)
A few changes following intro of new variable rhoInSitu
store 3-D (in-situ) density in commom block (this save 1 rho computation); re-arange rho diagnostics.
always call FREEZE_SURFACE if "allowFreezing" is set ; (+ stop in config_check if trying to use seaice pkgs with allowFreezing=T)
new pkg "down_slope".
replace calls to "FIND_RHO" with recent version "FIND_RHO_2D"
More flexibility in zeroadj.
GMREDI_WITH_STABLE_ADJOINT CPP option: special setup of gmredi for which the adjoint is stable enough for parameter optimization.
Arctic is different from global (prbbl bc of obcs)
More hidden recomputations, this time more thoroughly: avoid seaice_model and load_fields_driver
Remove this silly ifdef (what was I sinking...)
Some stuff needed for latest v3 production.
added check, readparms, and more diagnostice to pkg/salt_plume also changed package flag from useSaltPlume to useSALT_PLUME
added pkg/salt_plume by gathering code, which were previously spread around various files in model/inc and model/src results remain numerically identical to before, as a first step towards adding more options, etc., to this package
Fix a bug reported by An Nguyen: calc_salt_plume_depth requires that sigmaR be available. Bug fixed by adding flag "usePLUME", which is also a first step for addition of pkg/plume
Put back IF(useSEAICE) for adjoint (seems benign).
Adding code for seaice export (and others in future) sensitivity calculation.
Include SEAICE_SALINITY in adjoint.
Update NLFS adjoint.
Update NLFS adjoint.
Comment all relevant #ifndef ALLOW_AUTODIFF_TAMC that used to hide exch2 or cubed-sphere specific code (commented via 'cph-exch2')
o starting implementation of Duffy et al. (GRL 1999) salt plume scheme (#ifdef ALLOW_SALT_PLUME)
forgot to check-in this one (pkg/gmredi modification from yesterday).
One more store change (revert).
Re-arrange store directives.
forgot the #ifndef ALLOW_AUTODIFF_TAMC in the previous check-in
- move GMRedi after KPP and all vertical mixing scheme. - add call for computing mixed-layer depth (but hidden from TAF for now)
Seaice_model suddenly needs much more storing. Don't have time to check what has changed (again).
add call to KPP_DO_EXCH (after calling KPP_CALC). - fix conservation problem. - might be a temporary fix until we check carefully the computation domain of KPPviscAz.
add myIter to KPP_CALC call.
o Remove recomputations after latest modifs o EVP TLM and ADM gradient checks completely wrong (1.E34)
Fix KPP adjoint after introducing new fields KPPdiffKzS, KPPdiffKzT
Removing exf_clim code.
move EXF header files from lower_case.h to UPPER_CASE.h ;
First set of modifs for TAF-ing thsice.
First set of modifs for TAF-ing thsice.
o add calc_wsurf_tr.F to diff.list o slight changes in store dirs in do_oceanic_phys (M. Mazloff)
if fluid is not water, by-pass all the calls, except OBCS_CALC which was moved down to the end of bi,bj loops ;
o Small modifs (mostly porting for non-ecco setups)
Updating seaice adjoint, step 1 (everything, except SEAICE_EVP).
Adding store.
o ecco_readparms namelist was incomplete for new control terms o seaice load added to stored state o slightly more storing in do_oceanic_phys
Split seaice_init into _fixed, _varia
add hooks for atm2d package
move call to SEAICE_MODEL from FORWARD_STEP to DO_OCEANIC_PHYS (and Patrick fixed the adjoint)
use the right CPP option: "#ifdef ALLOW_OCN_COMPON_INTERF"
o Another overhaul of store dirs. for NLFS to eliminate "hidden" recomputations. o TBD: "hidden" mom_vecinv recomp. in dynamics
o Fix PmEpR store for #undef EXACT_CONSERV o Fix "hidden" recomp. in thermodynamics (no TAF warnings) o Try to fix NLFS adjoint (hidden recomp in forward_step)
move balanceEmPmR & balanceQnet out of the bi,bj loops
o move balancing of fluxes before the call of external_forcing_surf (thanks to Jean-Michel for catching this bug)
Change of EmPmR in external_forcing_surf necessitates re-arranging adjoint STORE (what a pain).
o add code to balance EmPmP and Qnet at the end of do_ocean_physics. Useful if bulk formulae are used in long integration (especially EmPmR). Turn on with balanceEmPmR = .true. or balanceQnet = .true. in data, PARM01 if balancePrintMean, the imbalance that is substracted is printed to STDOUT.
o add hooks for new package shelfice, painless
time-ave surf_flux & convect_adjust moved from thermodynamics to do_oceanic_phys
add diagnostics for Convective Adjusment.
- add call to diags_oceanic_surf_flux (ocean only) - remove unused variables.
add buoyancy (rho) flux diagnostics
use new S/R diagnostics_fill to fill-in diagnostics.
This set of changes restores TAMC(!) compatibility.
use flags: fluidIsAir/Water, usingP/ZCoords instead of buoyancyRelation
o Fixing directives for non-ECCO forward_step, after ECCO directives have been fixed to avoid recomp. of do_oceanic_phys o Works well for global_ocean.90x40x15, but recomp. issue remains for carbon
o some delicate re-shuffle of store directives to avoid one extra call of do_oceanic_physics o NB: this may break global_ocean adjoint temporarily, but it is clear how to fix it. Will do later, need this now.
modified argument list of CALL OBCS_CALC to myTime+deltaTclock
o bringing adjoint up to date for sheduled c55
o prepare addition of new packages GGL90 and OPPS - new parameters useGGL90 and useOPPS - include call to the new routines at the (hopefully) right places
add diagnostics od d.Rho/dr (=SigmaR).
o remove single quotes (eg. "don't"-->"do not") so that the on-line code browser does not get confused
o add calls for two new packages - pp81 (Packanowski and Philander, 1981), Richardson number and stratification dependent mixing - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
move ocean specific S/R calls (FREEZE_SURFACE, THSICE_MAIN, OCN_APPLY_IMPORT) from forward_step.F to do_oceanic_phys.F
replace surfaceTendency U,V,S,T,Tice by surfaceForcing U,V,S,T,Tice
max number of passive tracers is now defined in PTRACERS_SIZE.h
put atmospheric physics & state-vars diagnostics calls in 2 dedicated S/R.
take out 1rst part of thermodynamics (KPP,GM,slope ..) into do_oceanic_phys.F
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 |