Parent Directory
| 
Revision Log
| 
Revision Graph
| Links to HEAD: | (view) (download) (annotate) |
| Links to checkpoint59a: | (view) (download) (annotate) |
| Sticky Tag: |
Updates for NLFS adjoint.
2nd-Order Moment Advection Scheme (Prather, 1986): first check-in - enable by setting #define GAD_ALLOW_SOM_ADVECT (in GAD_OPTIONS.h) - used without limiter (AdvScheme=80) or with Prather limiter (AdvScheme=81)
add GAD_OPTIONS.h (to have DISABLE_MULTIDIM_ADVECTION def/undef) and fix argument list of CALC_WSURF_TR call.
o add calc_wsurf_tr.F to diff.list o slight changes in store dirs in do_oceanic_phys (M. Mazloff)
include ifndef ALLOW_AUTODIFF around linear free surface correction
Second wave: correction of tracer source/sink due to Linear Free surface
start to implement deep-atmosphere and/or anelastic formulation
Hidden recomputation for GMREDI BOLUS
Fix line length problem for NLFS adjoint.
Remove a "hidden" recomputation in bottom_ctrl_5x5/ But does not reconcile G r a d accuracy w.r.t. reference.
make a local copy of velocity to pass (like u,v,r_Trans) to tracer advection S/R
Modifications for bottom topography control o replace hFacC by _hFacC at various places o replace ALLOW_HFACC_CONTROL by ALLOW_DEPTH_CONTROL o add non-self-adjoint cg2d_nsa o update autodiff support routines o re-initialise hfac after ctrl_depth_ini o works for 5x5 box, doesnt work for global_ocean.90x40x15
Bracket Kuz, Kvz stores
Fix ptracers adjoint o avoid extensive recomputatations o fix missing re-init. of gptr (missed by TAF)
* use flag useOffLine rather than CPP #ifndef ALLOW_OFFLINE * does not skip calls in initialization stage if using OffLine mode
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)
Adjust store directives to changes in gad_calc_rhs
add IF (usePtracers) that was missing.
First step for a NLFS adjoint o initially suppress rStar (new flag DISABLE_RSTAR_CODE) o new init. routines for calc_r_star, calc_surf_dr o still need to deal with ini_masks_etc o testreport seemed happy
Allow to apply AB on T,S rather than on AB(gT,gS): - implemented within #ifdef ALLOW_ADAMSBASHFORTH_3 - use the same arrays (gtNm,gsNm) to hold tracer field at previous time-steps (if AB(T,S)) and tendencies (if AB(gT,gS)).
remove MINIMAL_TAVE_OUTPUT option.
remove call to diagnostics_fill_surf_flux (badly placed, within bi,bj loop)
o fix Samar's sloppiness: all code needs to be <73 characters long
Added support for pkg/matrix
do not refer to gxNm1 if unnecessary (to work also with AB-3)
added storing of kappaRTr, kappaRk
o add surface flux diagnostics: TAUX, TAUY, TFLUX, TRELAX, TICE, SFLUX, SRELAX
change argument list of S/R impldiff.F: tracerIdentity replace deltaT - allow to implement deltaT function of level k - make diagnostics easier
implement Implicit Vertical advection for pTracers
Added CPP option MINIMAL_TAVE_OUTPUT for minimal time-averaged output: S, T, U, V, W, ETA, and phiHydLow.
depth convergence accelerator: replace deltaTtracer by dTtracerLev(k)
allow to compile without generic_advdiff pkg.
removed unused arrays
Change arrays for total vertical diffusivity (2.D explicit, 3.D implicit)
Explicit vertical diffusivity: use a 2.D vertical diffusivity array
use flags: fluidIsAir/Water, usingP/ZCoords instead of buoyancyRelation
o bringing adjoint up to date for sheduled c55
o packages PP81, MY82, and GGL90 have variable diffusivity - include them to the list/logical variableK
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
o adding offline package
max number of passive tracers is now defined in PTRACERS_SIZE.h
try to keep an accurate diagnostic of timeave_surf_flux
take out 1rst part of thermodynamics (KPP,GM,slope ..) into do_oceanic_phys.F
Keeping up with JMC's latest modifs.
prepare splitting of thermodynamics: store convective counter in common block
T & S: separate Vert.Advec.Scheme from horizontal Advec.Scheme
o Adjoint for EBM package (needs further testing) o bug fix for parallel EBM in ebm_zonalmean.F
Commiting new energy balance model to repository o package is pkg/ebm o verif. is verification/global_ocean_ebm o references are in ebm_driver.F Will need long integration testing.
Bringing flow and store directives up-to-date
o limit timeave output for hi-res integrations
rewrite (as in MultiDimAdvec) explicit tracer stepping (gad_calc_rhs.F) to work with implicit vertical advection and AB
add calls for implicit vertical direction (advection & diffusion)
but keep impldiff for implicit diffusion & viscosity only.
don't call external_forcing_surf in atmospheric simulation
use the new thermodynamic Sea-Ice pkg: thSIce
add PTRACERS_OPTIONS.h (needed in forward mode, for NUMBER_OF_PTRACERS)
additional changes for FREEZE: - new S/R FREEZE_SURFACE only apllied to surface level. - add run-time parameter "useOldFreezing" to use the old version "FREEZE"
o modifications to make FREEZE flux visible to pkg/kpp
- moved surfaceTendencyTice from pkg/seaice to main code
- FREEZE moved to FORWARD_STEP
- subroutine FREEZE now limits only surface temperature
this means new output.txt for global_ocean.90x40x15,
global_ocean.cs32x15, and global_with_exf, but note
that results for these three experiments remain
bit-identical to before if allowFreezing=.FALSE.
o added surface flux output variables to TIMEAVE_STATVARS
o time-averaged output for pkg/ptracers
o merging from ecco-branch o minor CPP options update
o cleanup: convert '#ifndef DISABLE_DEBUGMODE"' to '#ifdef ALLOW_DEBUG"'
o initialisation of rFlx extended to full array (required by TAF) and shifted to thermodynamics o removed PTRACERS.h in ptracers routine o added surfacetendencyPtr to S/R parameter list pracers_forcing
o undid all of the cp51 checkin pending some ongoing code cleanups and discussion
o modifications to make FREEZE flux visible to pkg/kpp
- moved surfaceTendencyTice from pkg/seaice to main code
- FREEZE & EXTERNAL_FORCING_SURF moved to FORWARD_STEP
- subroutine FREEZE now limits only surface temperature
(this means new output.txt for global_ocean.90x40x15,
global_ocean.cs32x15, and global_with_exf)
o added surface flux output variables to TIMEAVE_STATVARS
adjusted some flow directives
o first check-in for the "branch-genmake2" merge
o verification suite as run on shelley (gcc 3.2.2):
Wed Oct 8 23:42:29 EDT 2003
T S U V
G D M c m s m s m s m s
E p a R g m m e . m m e . m m e . m m e .
N n k u 2 i a a d i a a d i a a d i a a d
2 d e n d n x n . n x n . n x n . n x n .
OPTFILE=NONE
Y Y Y Y 13 16 16 16 0 16 16 16 16 16 16 16 16 13 12 0 0 pass adjustment.128x64x1
Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 0 0 16 16 0 0 pass adjustment.cs-32x32x1
Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 22 0 16 16 22 0 pass adjust_nlfs.cs-32x32x1
Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O advect_cs
Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O advect_xy
Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O advect_xz
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass aim.5l_cs
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass aim.5l_LatLon
Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass exp0
Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass exp1
Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass exp2
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass front_relax
Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass global_ocean.90x40x15
Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL global_ocean.cs32x15
Y Y Y Y 6 11 12 13 13 12 13 16 13 9 9 9 9 10 9 9 11 FAIL global_ocean_pressure
Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass global_with_exf
Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16 9 16 pass hs94.128x64x5
Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass hs94.cs-32x32x5
Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL ideal_2D_oce
Y Y Y Y 8 16 16 16 16 16 16 16 16 13 13 8 16 16 16 16 16 FAIL internal_wave
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass inverted_barometer
Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL lab_sea
Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL natl_box
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass plume_on_slope
Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass solid-body.cs-32x32x1
fix Problem with bulk_force & therm_seaice.
Bringing code up to date for AD o remove some IF-statements which cause excessive dependencies o provide interface for ADM*TLM
o Mods and bug fixes to pkg/cal, pkg/exf, etc., needed for computation of tracer Green's fucntions for ocean inversion project.
o Mods and bug fixes to pkg/cal, pkg/exf, etc. needed for computation of tracer Green's fucntions for ocean inversion project.
o Mods and bug fixes to pkg/cal and pkg/exf needed for computation of tracer Green's fucntions for ocean inversion project.
initialize tendency to zero : proposed by Martin to fix the bug when setting temp/salt advection=.false.
put thermSeaIce CALL(s) within if (useThermSeaIce) block
o introducing integer flag debugLevel o introducing pathname variable mdsioLocalDir for mdsio
o disentangled ALLOW_PTRACERS using new ALLOW_GCHEM
Preparing next differentiable checkpoint and sync of MAIN vs. ecco-branch (updating store after changes in checkpoint50b_post, plus still messing around with init. sequence).
Extended pkg/debug and instrumented main code to help track down fatal errors.
fix bug in convective adj. diagnostic (ConvectCount is a local array !) and add time-average diagnostic of U*V and diff_Kr*theta
Changes to restore differentiability of code w.r.t. previous tag (mostly adding new routines to make list and replacing pressure by totPhiHyd).
remove unused phiHyd
Dummy init. of visbeckK for TAF.
o few modif.'s to get latest version adjointed (mainly kick out code in ini_linear_phisurf) o modif's to run adjoint with exactConserv o case GM_BOLUS_ADVEC should be cleaned S/R gmredi_slope_psi should be cleaned o verification/carbon now has exactConserv=.TRUE.
Adjusting one storing for case GM_BOLUS_ADVECT
modification needed for the new AIM package (aim_v23)
Two packages: bulk_force (Bulk forcing) and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2 Earlier integration from Stephaine Dutkiewicz and Patrick Heimbach. Two ifdef statements for compile time, ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE Two switches in data.pkg to turn on at run-time: cat data.pkg # Packages &PACKAGES useBulkForce=.TRUE., useThermSeaIce=.TRUE., & WARNING: useSEAICE and useThermSEAICE are mutually exclusive. The bulk package requires an additional parameter file with two namelists, data.ice and data.blk. c ADAPTED FROM: c LANL CICE.v2.0.2 c----------------------------------------------------------------------- c.. thermodynamics (vertical physics) based on M. Winton 3-layer model c.. See Bitz, C. M. and W. H. Lipscomb, 1999: "An energy-conserving c.. thermodynamic sea ice model for climate study." J. Geophys. c.. Res., 104, 15669 - 15677. c.. Winton, M., 1999: "A reformulated three-layer sea ice model." c.. Submitted to J. Atmos. Ocean. Technol. c.. authors Elizabeth C. Hunke and William Lipscomb c.. Fluid Dynamics Group, Los Alamos National Laboratory c-----------------------------------------------------------------------
Two packages: bulk_forcing (Bulk forcing) and thermodynamic_seaice (adapted from LANL CICE.v2.0.2) Earlier integration from Stephaine Dutkiewicz and Patrick Heimbach. Two ifdef statements for compile time, ALLOW_TSEAICE and ALLOW_BULKFORMULA Two switches in data.pkg to turn on at run-time: cat data.pkg # Packages &PACKAGES useBulkf=.TRUE., useTSeaIce=.TRUE., & The bulk package requires an additional parameter file with two namelists. cat data.blk &BULKF_PARM01 RainFile= 'ncep_precip_m_cubed.bin', SolarFile= 'ncep_downsolr_cubed.bin', AirTempFile= 'ncep_tair_cubed.bin', AirhumidityFile='ncep_qair_g_cubed.bin', LongwaveFile= 'ncep_netlw_cubed.bin', UWindFile= 'ncep_uwind_cubed.bin', VWindFile= 'ncep_vwind_cubed.bin', WspeedFile= ' ', RunoffFile= ' ', QnetFile= ' ', EmPFile= 'ncep_emp_calc_cubed.bin', CloudFile= 'ncep_totalcloud_cubed.bin', & &BULKF_PARM02 qnet_off=0.0, empmr_off=0.0, conservcycle=311040000., & c ADAPTED FROM: c LANL CICE.v2.0.2 c----------------------------------------------------------------------- c.. thermodynamics (vertical physics) based on M. Winton 3-layer model c.. See Bitz, C. M. and W. H. Lipscomb, 1999: "An energy-conserving c.. thermodynamic sea ice model for climate study." J. Geophys. c.. Res., 104, 15669 - 15677. c.. Winton, M., 1999: "A reformulated three-layer sea ice model." c.. Submitted to J. Atmos. Ocean. Technol. c.. authors Elizabeth C. Hunke and William Lipscomb c.. Fluid Dynamics Group, Los Alamos National Laboratory c----------------------------------------------------------------------- WARNING: useSEAICE and useTSEAICE are mutually exclusive. todo: thermodynamic.F should be reviewed and cleaned up a bit.
differentiable version of checkpoint46n_post o external_fields_load now part of differentiation list o pressure needs multiple storing; would be nice to have store_pressure at beginning or end of forward_step, e.g. by having phiHyd global (5-dim.) (NB: pressure is needed for certain cases in find_rho, which is also invoked through convective_adjustment). o recomputations in find_rho for cases 'JMD95'/'UNESCO' or 'MDJWF' are OK. o #define ATMOSPHERIC_LOADING should be differentiable o ini_forcing shifted to begining of initialise_varia
Merging from release1_p8 branch: Adding package parameters and hooks for new seaice package.
* S/R aim_initialise.F replace S/R aim_init from file aim_do_inphys.F: - read AIM physics parameters from a file (data.aimphys) - set defaults values = F.Molteni paper (Clim.Dyn., 2002) - phiHyd removed from arg. list of S/R aim_do_atmos_physics
* split calc_exact_eta in 2 S/R : integr_continuity & update_etaH * move wVel computation at the end of the time step, in S/R integr_continuity * create specific S/R to exchange T,S before DYNAMICS (for stagger time step) * update timeave pkg for wVel diagnostic ; put convertEmP2rUnit in PARAMS.h
o Include a new diagnostic variable phiHydLow for the ocean model
- in z-coordinates, it is the bottom pressure anomaly
- in p-coordinates, it is the sea surface elevation
- in both cases, these variable have global drift, reflecting the mass
drift in z-coordinates and the volume drift in p-coordinates
- included time averaging for phiHydLow, be aware of the drift!
o depth-dependent computation of Bo_surf for pressure coordinates
in the ocean (buoyancyRelation='OCEANICP')
- requires a new routine (FIND_RHO_SCALAR) to compute density with only
Theta, Salinity, and Pressure in the parameter list. This routine is
presently contained in find_rho.F. This routine does not give the
correct density for 'POLY3', which would be a z-dependent reference
density.
o cleaned up find_rho
- removed obsolete 'eqn' from the parameter list.
o added two new verification experiments: gop and goz
(4x4 degree global ocean, 15 layers in pressure and height coordinates)
Merging from release1_p5 (cf. tag-index for checkpoint46).
Oly must be > 3 when using an advection scheme other than 2nd O and staggerTimeStep & multiDimAdvection & GM_AdvForm & not.GM_AdvSeparate Add an exchange of gT,gS in this case.
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
Included CPP option SINGLE_LAYER_MODE to configure barotropic setup (Martin Losch).
extend the computation domain of Rho to cover the whole domain. (necessary when using GM-Advect form + stagger time step with Olx=Oly=3)
o Modified initialisations to break adjoint dependencies o removed some store directives o added options files for KPP, GMREDI
add GM-bolus transport to Eulerian transport to advect Tracers (T,S, ...)
Bug fix: Missing #ifdef's for when PTRACERS are not enabled.
Added PTRACERS package
This allows an arbitrary number of passive tracers to be integrated
forward simultaneously with the dynamicaly model.
+ Implemented so far:
- basic forward algorithm (time-stepping, advection, diffusion, convection)
- I/O and checkpointing
- GM/Redi *but* using the GM/Redi coefficient of Salt
+ Not implemented so far:
- KPP
- OBCS
+ No specific example supplied (yet) but global_ocean.90x40x15 has the
necessary data.ptracer file. Simply use -enable=ptracers and uncomment
line in data.pkg. PTRACER01 then reproduces Salt exactly.
+ This package is disabled by default since it increases storage.
add argument myIter to S/R obcs_calc
put the k loop inside S/R gmredi_calc_tensor
fix diagnostic of convective adjustment (IVDC) (broken since thermo-/dynamics split)
Removed CPP connection between DISABLE_MULTIDIM_ADVECTION and ALLOW_AUTODIFF_TAMC. o ie. removed #ifdef ALLOW_AUTODIFF_TAMC from GAD_OPTIONS.h o updated comments in thermodynamics.F to reflect this change
Switched sense of CPP macro for controlling multi-dimensional advection: o DISABLE_MULTIDIM_ADVECTION is set in GAD_OPTIONS.h o automatically set if differentiating code (comments around call to gad_advection point out how to re-enable it) o this avoids needing to add the former macro to CPP_OPTIONS.h - reason for this is there's no point in any of the new schemes without it.
Fixed AD-related problems: o Store directives up-to-date with re-arranged Adams-Bashforth (mainly thermodynamics.F) o New store directives for multi-dim. advection schemes * new CPP flag ALLOW_MULTI_DIM_ADVECTION * new common block and key passkey (mainly gad_advection.F) o Modified store directives for split of dynamics/thermodynamics for the case ALLOW_KPP o Cleaned argument list for timestep_tracer.F
Deleted apostrophies in comments: "Isn't" breaks cpp...
Bringing comments up to data and formatting for document extraction.
Corrected parameter lists for case tr1stepping.
Re-arranged sequence of operations for Adams-Bashforth o this does not change numbers o this makes it very easy to extract forcing/diffusion out of ABII by changing calling sequence in calc_gt, calc_gs,... Key modifications: o new s/r: ADAMS_BASHFORTH2 gT=3/2*gT-1/2*gTnm1 o changed TIMESTEP_TRACER from gTnm1=t+dt*(3/2*gT-1/2*gTnm1) to gT=T+dt*gT o changed CALC_GT,CALC_GS & CALC_GTR1 to calcuate "gT" defined by new timestep_tracer (ie. including forcing, ABII, N-L F-S, etc...) now calls ADAMS_BASHFORTH2 and FREESURF_RESCALE_G o changed CYCLE_TRACER appropriately T=gT only Other associated mods: o new s/r: FREESURF_RESCALE_G applies non-linear free-surface term this used to be in TIMESTEP_TRACER o added myIter as argument to CALC_GS,CALC_GT,CALC_GTR1
Reinstate calls to advection routines after gratuitous commenting removed them in previous update! :^)
Include initialisation of gTracer fields to break dependency for TAF.
Added multi-dimensional form of advection o available only for single step schemes (ie. can't be used with ABII) o stable for max(cfl_u,cfl_v,cfl_w)<=1 (without cfl_u+cfl_v+cfl_w <=1) o selected using multiDimAdvection=.T. (default) o had to hack some existing routines to work on local arrays
Added run-time control of advection schemes.
- advection scheme determines method of forward integration.
- unfortunately, we have to use integers in "data" since ENUM_CENTERED_2ND
doesn't mean anything to fortran
- defaults are centered second
- output differs due to these mods! This is due to the g77 optimization.
I have tested that using -ffloat-store, these mods do not affect
the output so am confident about changes.
T S U V
C D M c m s m s m s m s
n p a R g m m e . m m e . m m e . m m e .
f n k u 2 i a a d i a a d i a a d i a a d
g d e n d n x n . n x n . n x n . n x n .
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 22 16 16 pass adjustment.128x64x1
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass adjustment.cs-32x32x1
Y Y N N -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- N/O aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass aim.5l_LatLon
Y Y N N -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- N/O aim.5l_zon-ave
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 22 16 pass exp0
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 22 16 pass exp1
Y Y Y Y 13 16 16 16 16 16 16 16 13 16 13 13 13 13 13 13 16 pass exp2
Y Y Y Y 12 16 16 13 16 16 16 16 16 13 16 16 16 16 13 13 16 FAIL exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass front_relax
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 11 16 pass hs94.128x64x5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 22 16 pass hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 pass hs94.cs-32x32x5
Y Y Y Y 16 16 16 22 16 16 16 16 16 16 16 22 16 16 16 16 16 pass internal_wave
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass natl_box
Y Y Y Y 16 16 16 16 16 13 16 13 16 16 16 16 16 16 16 16 16 pass solid-body.cs-32x32x1
Modifications related to split into thermodynamics.F, dynamics.F o missing initialisations in dynamics.F added o some fields no longer needed in dynamics/thermodynamics deleted o split of calc_diffusivity.F into calc_viscosity.F (plus split of kpp_calc_diff.F into kpp_calc_visc.F) o Modifications of some store directives for TAF
Split dynamics.F into dynamics.F and thermodynamics.F - idea is to make algorithm more transparent??? - probably less efficient - has exchanges at end of thermodynamics.F (which are needed if using staggered time-stepping with the cube OR using AIM.
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 |