1 
Change log for branch "pre38". 
2 

3 
o cosine(latitude) scaling of viscosity and diffusion 
4 
 moved declaration of "cosFac" from calc_mom_rhs.F to GRID.h 
5 
 now calculated in ini_spherical_grid.F (set to 1 in ini_cartesian_grid.F) 
6 
 "Method III" (preferred) available for biharmonic term 
7 
 anisotropic (preferred) scaling available 
8 
 added scaling for diffusion (same features as for viscosity above) 
9 
x thought for the future: lateral GM/Redi fluxes are currently scaled. 
10 
If we scale them we would have to scale the other tensor components 
11 
in order to retain the adiabatic property of scheme. 
12 
o combined initializations of G's and other timestepping arrays 
13 
 created ini_garrays.F to initialize G's to zero 
14 
 combined uVel + vVel initialization into ini_vels.F 
15 
 deleted ini_uvel.F + ini_vvel.F 
16 
 removed ini_gw.F since no longer necessary 
17 
o rearranged order of arrays in /GRID_R/ in GRID.h 
18 
g77 suddenly started complaining about padding. Odd but SGI concurred! 
19 
o fixed CPP_EEOPTIONS.h in hs94.128x64x5 
20 
o cg2d: added CPP flag for using "regular" overlaps instead of Olx=Oly=1 
21 
o added inheritence of disable option to package subdirectories 
22 
o testscript has skip option and puts full output into output.txt 
23 
o added pkg/mom_fluxform: conservative form of momentum equations 
24 
 *should* produce identical results to original calc_mom_rhs.F 
25 
o added pkg/mom_vecinv: vectgor invariant form of momentum equations 
26 
 noslip side conditions not implemented. Don't try to use them! 
27 
o define separatetly surface and bottom position  set maskC as 3D global array 
28 
o modify Surface correction term to deal with topography in Atmosphere 
29 
o prepare implementation of nonLinear Free Surface: 
30 
 add parameters and CPP_options for nonLinear Free Surface 
31 
 move etaNm1 to main common block 
32 
o allow to recompute eta to conserve exactly the total volume. 
33 
o added fCoriG to PARAMS.h, ini_cori.F, FCORI_MACROS.h 
34 
 needed for VI momentum equations which are broken without. 
35 
o CG2D: allow to switchoff RHS normalisation setting Residual_Target in W_unit 