Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Sticky Tag: |
new options 1) to account for true vertical distance (including hFac) in vertical viscous flux and diffusive flux ; 2) to increase vertical mixing near surface and/or bottom where partial cell is too thin ;
add run-time param to allow to turn off tidal forcing.
new run-time param (useMin4hFacEdges) to select method for setting hFacC
- add specific run-time param to select level of printed plot-field-maps, set by default to debugLevel. Remove internal param "printDomain".
comment out writing of "pickup_write_immed" (never used + not in namelist)
new input parameter "rwSuffixType" (read from PARM03 namelist) to select file suffix format; by default (=0), use iteration number;
cleaning
add new parameter "selectImplicitDrag" for implicit treatment of bottom drag
- new parameters "top_Pres" & "seaLev_Z" (replacing Ro_SeaLevel and recently added phi0Ref) to set vertical axis origin and phiRef origin according to coordinate and fluid type.
rename internal parameter "useDynP_inEos_Zc" to "storePhiHyd4Phys" (since it's also used with pkg/atm_phys) and add few more parameters (but still inactive, src code yet to come)
forgot to add new run-time flag (applyExchUV_early) in namelist (PARM03)
adjust POLY3 coeffs summary
add run-time parameters to select bottom-drag options
remove "recip_Cp" from common block ; add thetaConst (for pStar)
fix reported value of hFacMinDr
useSingleCpuInput separates single-CPU input from single-CPU output
add run-time parameters for isotropic 3-D Smagorinsky (Smag3D)
change include order (for OpenAD): MOM_VISC.h before SET_GRID.h
print value of viscAhW,viscA4W and (new) logical parameters: useHarmonicVisc, useBiharmonicVisc & useVariableVisc.
extend length of drC & recip_drC from Nr to Nr+1 (forgot to check-in this one)
add run-time parameter to select CG2D solver minimum-residual solution (in case of a poor convergence)
- allow to read-in time-constant addMass field from file (addMassFile); - add support to account for frictional heating (run-time switch & CPP option).
report value of params balanceQnet/EmPmR, doTheta/SaltClimRelax and balanceTheta/SaltClimRelax
add logical flag, true if using Cubed-Sphere Exch with CS-corners inside the domain
add coeffs (1-D array) for vertical biharmonic diffusion on T & S.
add an internal flag, true when free-surface is always and everywhere at level k=1 ; for now, just set it to "usingZCoords"
- only report Linear EOS parameters if using Linear EOS - add celsius2K ; improve description of rhoConst,rhoNil,rhoConstFresh.
move delX,delY to new header file (SET_GRID.h) and adjust length to 1rst face dimensions.
new run-time switch "doResetHFactors" (default=false) to reset (or not) thickness factors (hFac*) at the beginning of forward_step.F
new internal flag "calc_wVelocity" (= momStepping .OR. exactConserv)
report printResidualFreq and improve debLev[*] description
print out useEnergyConservingCoriolis flag value; add "V.I" to description of vector-invariant specific parameters.
fix typo in description
add new param to scale curvilinear grid distance and area when using a different rSphere.
print out "nEndIter" value (was missing).
- remove useConstantF, useBetaPlaneF, useSphereF and replace them with integer parameter "selectCoriMap" (=0,1,2). - add parameter "fPrime" for constant second Coriolis coeff (default=0).
first check-in of sigma (and hybrid-sigma) coordinate code
Adding parameters to specify temperature and salinity of addMass array separately from temp_EvPrRn and salt_EvPrRn
remove unused variables
add parameter for Non-Hydrostatic free-surface option (in development)
forgot to include output of new parameter useSRCGSolver and a test
epsAB_CD (CD-code AB parameter): default is set to abEps
allow to specify a vertical profile of vertical viscosity
changed to pass when compiling with strick checking of arguments across S/R
call WRITE_1D_R* with the right type of argument
- 1 argument added to S/R PRINT_LIST_* - call WRITE_0D_RL (instead of WRITE_0D_R8) to print "RL" parameters
rename thetaMin,phiMin -> xgOrigin,ygOrigin (temporary backward compatibility in ini_parms.F, until next checkpoint)
Add smoothing parameter for mixed-layer diagnostics
dRhoSmall (for mixed-layer diagnostics) now a run-time paramater
update bottom-drag units
initialisation and option setting for mass source/sink of fluid
add description for selectVortScheme=3
update description of "select_rStar" (but still not very explicit)
set flags for Dynamically Active Tracers
new parameter "selectVortScheme" to select which scheme to use for vorticty-term (in Vector Inv.); remove "SadournyCoriolis" flag.
unexpected values (mass2rUniti x 3) for the 3 "Euler angle" ; fixed.
add a flag writePickupAtEnd (default=true) to be able to suppress writing a pickup at the last timestep
add parameter rotateGrid and Euler angles EulerPhi/Theta/Psi. For usingSphericalGrid, this allows to define the rotated grid coordinates via phiMin/thetaMin/dxSpacing/dySpacing, etc., but then re-compute the geographical coordinates according to the inverse of the rotation defined by the Euler angles, so that Coriolis parameter and online-interpolation by exf works with the geographical coordinates on XC/YC/XG/YG. The three Euler angles PhiEuler, ThetaEuler, and PsiEuler define the rotation about the original z-axis (of an sphere centered cartesian grid), the new x-axis, and the new z-axis, respectively. Naturally, this feature does not work with all packages, so the some combinations are prohibited in config_summary (flt, flt_zonal, ecco, profiles), because there the coordinates are assumed to be regular spherical grid coordinates.
- fix units of bottom drag Coeffs - distinguishe between 'm' or 'Pa' for vertical viscosity/diffusivity coeffs
add "monitorSelect" to reduce number of monitored fields (default unchanged)
add new run-time parameter, "pickupStrictlyMatch": when reading pickup, check and stop if pickup-files do not stricly match; default=T.
use mass2rUnit & rUnit2mass as units conversion factor (for surface forcing) and remove horiVertRatio (& recip_horiVertRatio);
add temperaure criteria for Mixed-layer diagnostic
use WRITE_0D_C to print character parameter value.
define and use reference profile for w <-> omega conversion in NH-code.
initial check in: correction of tracer source/sink due to Linear Free surface
print out "the_run_name".
start to implement deep-atmosphere and/or anelastic formulation
put the right units for dx?,dy? and rA? ; and write globalArea.
add useSphereF, useBetaPlaneF and use3dCoriolis to namelist PARAM01 ; change default of useSphereF to T (if set to F => read Coriolis from files);
fix typo.
stop when mom_fluxform/mom_vecinv are not compiled but monStepping=T
print value of quasiHydrostatic, metricTerms, useNHMTerms & momDissip_In_AB
change forcing_In_AB to affects both T,S forcing and Momentum forcing (allow to differentiate between forcing components using new integer flags: momForcingOutAB=1/0 & tracForcingOutAB=1/0) and add new flag to put Dissipation tendency out of Adams-Bashforth.
add internal parameter "use3Dsolver"
add new parameter (implicitIntGravWave) even if the code is not yet ready; (hard to maintain a different version of theses files during testing phase)
add write/readBinaryPrec, globalFiles & useSingleCpuIO to the summary.
add 3 lines at the end (easier to search in STDOUT)
new parameter: doAB_onGtGs to (always) apply Adams.Bashforth on tracer tendencies (rather than on tracers <- coming soon).
new parameter: sideDragFactor to enable half-slip-side BC.
add KEscheme selector.
Allow use of old scheme for viscous lengths with useAreaViscLength=.TRUE. Default: .FALSE.
Parameter and config_summary changes for mom_calc_visc.F
"dumpInitAndLast" replaces "nodumps" (& partly MINIMAL_TAVE_OUTPUT)
echo debug-levels.
rkSign & usingPCoords replace -rkFac & groundAtK1 (<- removed)
add flag for momentum vertical advection (upwindShear)
add parameters for AB-3 ; remove readPickupWithTracer & writePickupWithTracer
add baseTime parameter = model base time (time origin)
make the comments shorter.
Oops, fixed line widths...
Added Harmonic Smagorinsky Viscosity. May not work with strongly anisotropic grids.
Expanded Leith viscosity to follow Ahgridmin and Ahgridmax, and select between 'full' and approximate form. Also added leith-like viscosity that is proportional to grad(div.v_h).
values of mnc flags (pickup_write_mnc, pickup_read_mnc, snapshot_mnc, monitor_mnc) were not printed; fix it.
depth convergence accelerator: replace deltaTtracer by dTtracerLev(k)
remove some obsolete params.
fix typo error
add separated viscosity for Divergence and Vorticity
set flags: fluidIsAir, fluidIsWater, usingPCoords, usingZCoords according to buoyancyRelation
allow to set a vertical profile of vertical diffusivity for T & S
o I am *sick* of moving these variables around. But Jean-Michel has all but threatened a hissy fit if they aren't removed from PARAMS.h. So now here they are *back* in MNC_PARAMS.h where they were just a few days ago.
o move useMNC and related runtime switches to PARAMS.h
o fix name error: "monitor_mdsio" --> "monitor_stdio" - add "monitor_grid" output for ini_cori.F,ini_parms.F
o overhaul of IO so that we now have flags for MDSIO and/or MNC - all verification tests compile and run with linux_ia32_g77 - defaults are compatible with current input files--nothing should change if you were not previously using MNC - MNC output has been added in numerous places (eg. timeave) but there are still a few writes not yet do-able with MNC (this is in progress) - flags now allow for either/or/both use of MDSIO and MNC and documentation will soon follow - numerous small formatting cleanups for ProTeX
o make the mon_iotype setup in monitor the poster child for MNC/MDSIO flag handling
print out no_slip_bottom flag
Only update cg2d preconditioner every "cg2dPreCondFreq" iter. (default=1)
needs valid U,V in halo region for multi-Dim-Advec; do Exch(U,V) if staggerTimeStep
merged cylindrical coord configuration and rotating_tank exp
Added new flag to fix aim.5l_cs - useJamartMomAdv distingiushes the non-linear terms from linear Coriolis when using Jamart wet points
Added vertical diffusivity profile (T/S) due to Bryan and Lewis, 1979. New parameters: diffKrBL79surf - surface diffusivity diffKrBL79deep - deep diffusivity diffKrBL79Ho - turning depth for arctan function diffKrBL79scl - depth scale for arctan function This diffusivity is added to all other diffusivities. The defaults are set so as to give zero diffusivity.
Added Sadourny discretization of Coriolis in V.I. mode - moved some PARAMETERS from mom_*_coriolis.F to PARAMS.h - re-enabled use of omega3 in mom_vecinv.F
Added variable viscosity for the vector invariant equations based on Leith, 1968, Phys. Fluids (10) 1409-1416 - the use of the variable viscosty in the no-slip boundary conditions has not been implemented (but should be) - new parameters viscC2leith and viscC4leith are non-dimensional - I decided to modulate the variable viscosuty with the same viscAhMax and viscA4max; ideally we should have another maximum based on dx^2/dt etc.
Added new grid-dependent viscosity and biharmonic dissipation terms: o new parameters = viscAhGrid, viscAhMax, viscA4Grid and viscA4Max o new stuff in mom_vi_hdissip.F which is bound to break the adjoint!!!! This has been "verified" against old experiments **BUT** there is a known issue involving exch2 - results are dependent on using exch1 or exch2 and we don't know which is right (if either).
add run-time parameters for implicit vertical advection.
fix typo error.
delp,delz viscAz,viscAp ... are only local variables (and no longer in a common bloc).
* reverse gravitySign back to "logical" value: +1 when R increases downward * add a parameter to apply surface relaxation only in a latitude band
specific volume (atmospheric geopotential) includes water vapor effect
Merging from c51
new flag "useCDscheme" (default=F);
set the Earth rotation either directly (omega) or from the rotation period
merging c49 and e27
minor cleaning.
first step in r* implementation: add param "select_rStar"
* allows a more accurate definition of Ro_Surf (selectFindRoSurf=1) when using P-coordinate; only implemented for atmospheric config.
Fixed inconsistent arguments that show up using run-time checking with IFC. - fool_the_compiler() called with an argument more often than not! - config_summary() was use vector I/O for scalars. To make model work with run-time checking in general you will need to not enable it for mdsio_byteswap_r?.F also.
o cleaned up the use of rhoNil and rhoConst. - rhoNil should only appear in the LINEAR equation of state, everywhere else rhoNil is replaced by rhoConst, e.g. find_rho computes rho-rhoConst and the dynamical equations are all divided by rhoConst o introduced new parameter rhoConstFresh, a reference density of fresh water, to remove the fresh water flux's dependence on rhoNil. The default value is 999.8 kg/m^3 o cleanup up external_forcing.F and external_forcing_surf.F - can now be used by both OCEANIC and OCEANICP
o Added new equation of state -> JMD95Z and JMD95P - EOS of Jackett and McDougall, 1995, JPO - moved all EOS parameters into EOS.h - new routines ini_eos.F, store_pressure.F o Added UNESCO EOS, but not recommended because it requires in-situ temperature (see JMD95) o Modified formatting for knudsen2.f in utils/knudsen2 and added unesco.f to be used with POLY3
Add new flags: * T,S forcing outside Adams-Bashforh * temp,salt Advection and Forcing (turn on/off) * for each tracer: internal flag for multiDimAdvection & A-B
Changes encapsulated by checkpoint43a-release1mods and chkpt44d_post on the main trunk. These are: o added missing EXCLUDE_MONITOR flags o changed "e" to "_d" in gmredi_slope_limit, gmredi_slope_psi (incompatible typ in MIN/MAX expressions caused problems on IBM SP3) o in genmake added variable MAKEDEPEND plus resetting for case SunOS o added timer_stats.c routine for IBM SP3 o removed variables in dynamics o real fresh water flux implemented with non-linear free-surface. o few fix (mask in shap_s2, EmPmR in external_field_load, USE_NATURAL_BCS in solve_for_P); o add arguments myIter & myTime to S/R obcs_calc & solve_for_P o merge of relevant stuff from the ecco-branch: - genmake: removed $S64 overwrite for case SunOS - pkg/exf: update and corrections for field swapping and obcs - pkg/ecco: parameter lists for the_model_main, the_main_loop harmonized between ECCO and MITgcm - pkg/autodiff: added flow directives for obcs, mdsio_gl_slice updated checkpointing_lev... lists for obcs - model/src: minor changes in forward_step, plot_field added directive for divided adjoint in the_main_loop - pkg/mdsio: added mdsio_gl_slice o check parameters & config (chkpt44a_pre,post) o OBC and NonLin_FrSurf. o fix bug in mom_vi_del2uv o select when filters are applied ; add options to zonal_filter (data.zonfilt) o gmredi: fix Pb in the adiabatic form ; add options (.e.g. Bolus advection) o update AIM experiments (NCEP input files) o improve and extend diagnostics (Monitor, TimeAve with NonLin-FrSurf) o added some stuff for AD These were merged with cvs co -r release1 -P MITgcm cd MITgcm cvs update -kk cvs update -j checkpoint43a-release1mods -j chkpt44d_post -d -P -kk
o define cell Center vertical spacing and then put Interface at the middle
add parameters to use real fresh water flux (replace USE_NATURALBCS option) with non-linear free-surface.
Bringing comments up to data and formatting for document extraction.
use gBaro instead of gravity to allow External mode reduced gravity test
remove printing of gBaro (no longer used)
add 2 parameters to define a Crank-Nickelson barotropic time-stepping
o Added printing of key grid variables in config_summary.F and removed write(0,*) output of these variables in ini_spherical_polar_grid.F o Added two new routines to do consistently formatted output of lines of constant X or Y for an XY variable. New routines are in read_write.F
Made sure each .F and .h file had the CVS keywords Header and Name at its start. Most had header but very few currently have Name, so lots of changes!
Merged changes from branch "branch-atmos-merge" into MAIN (checkpoint34) - substantial modifications to algorithm sequence (dynamics.F) - packaged OBCS, Shapiro filter, Zonal filter, Atmospheric Physics
Packaged OBCS (open boundary conditions) o source code now in pkg/obcs o controlled by namelist in data.obcs o includes Orlanski radiation condition (in same directory) o affects exp4 and internal_wave
vertical discretization for atmospheric run, with k=1 at the surface: add parameter: Ro_SeaLevel = Total pressure at Sea-Level ; add flag: groundAtK1 = True if surface(k=1) at lower boundary (ground) set horiVertRatio to Gravity*rhoConst ;
Removed GM/Redi parameters as part of packaging process.
Modifications to include TAMC directives, tape key computations and initialisations to make code TAMC compatible. Routines the_model_main.F and initialise_fixed.F are left unchanged for the moment. (P.H.)
Changed some erroneous declarations of myCurrentTime (sometimes myTime) from REAL to _RL. This probably was caught before because most people are using -r8 compile option in conjuction with REAL4_IS_SLOW which converts everything to real*8. Spotted by C.E.
Modifications for non-hydrostatic ability + updates for open-boundaries.
Added IMPLICIT NONE in a lot of subroutines. Also corrected the recip_Rhonil bug: we didn't set it in ini_parms.F
Changes to allow for atmospheric integration builds of the code
Changes to support - g77 compilation under Linux - LR(1) form of 64-bit is D or E for constants - Modified adjoint of exch with adjoint variables acuumulated.
Consistent isomorphism changes
Consistent isomorphism changes
Consistent isomorphism changes
Isomorphism consistency changes
Added print of deltaTclock - the model clock
Minor changes to correct bugs with multi-process mode of operation
Added flags for turning off momentum equation and/or temperature equation
Added extra IO features - runtime on/off flags - more reporting of configuration
Added changes to support implicit free-surface. - included gBaro a "barotropic" gravity that can be set differently to the g.rhoprime gravity. - discovered and fixed coding error in dynamics loop. Per tile temporaries that needed correct initial values were not being reset for each tile.
Added support for binary IO of model fields for restart and/or postprocessing
Further changes to reporting of model configuration on standard output
Improvements to I/O and feedback info.
*** empty log message ***
*** empty log message ***
Further $Id to $Header conversions
Imported working release
Initial revision
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 |