Parent Directory | Revision Log | Revision Graph
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
1 | adcroft | 1.2 | # ==================== |
2 | # | Model parameters | | ||
3 | # ==================== | ||
4 | # | ||
5 | # Continuous equation parameters | ||
6 | &PARM01 | ||
7 | tRef=281.089, 289.321, 307.191, 331.377, 432.413, | ||
8 | sRef=5*0., | ||
9 | viscAr=0.E3, | ||
10 | viscAh=0.E6, | ||
11 | viscA4=0.E17, | ||
12 | no_slip_sides=.FALSE., | ||
13 | no_slip_bottom=.FALSE., | ||
14 | diffKhT=0.E3, | ||
15 | diffKrT=0., | ||
16 | diffK4T=0.E17, | ||
17 | diffKrS=0.E2, | ||
18 | diffKhS=0.E3, | ||
19 | diffK4S=0.E17, | ||
20 | buoyancyRelation='ATMOSPHERIC', | ||
21 | eosType='IDEALGAS', | ||
22 | jmc | 1.5 | Integr_GeoPot=1, |
23 | adcroft | 1.2 | gravity=9.81, |
24 | rhonil=1.0, | ||
25 | rhoConst=1.0, | ||
26 | rigidLid=.FALSE., | ||
27 | implicitFreeSurface=.TRUE., | ||
28 | staggerTimeStep=.TRUE., | ||
29 | nonHydrostatic=.FALSE., | ||
30 | readBinaryPrec=64, | ||
31 | writeBinaryPrec=64, | ||
32 | hFacMin=1.0, | ||
33 | zonal_filt_lat=45.0, | ||
34 | adcroft | 1.6 | saltAdvScheme=77, |
35 | adcroft | 1.2 | & |
36 | |||
37 | # Elliptic solver parameters | ||
38 | &PARM02 | ||
39 | cg2dMaxIters=200, | ||
40 | #cg2dTargetResidual=1.E-12, | ||
41 | cg2dTargetResWunit=1.E-17, | ||
42 | & | ||
43 | |||
44 | # Time stepping parameters | ||
45 | &PARM03 | ||
46 | niter0=25920, | ||
47 | adcroft | 1.4 | nTimeSteps=10, |
48 | adcroft | 1.2 | #endTime=31104000., |
49 | deltaT=300.0, | ||
50 | abEps=0.1, | ||
51 | cAdjFreq=0., | ||
52 | pChkptFreq=2592000.0, | ||
53 | chkptFreq=0.0, | ||
54 | dumpFreq=2592000.0, | ||
55 | taveFreq=2592000.0, | ||
56 | adcroft | 1.4 | monitorFreq=1.0, |
57 | adcroft | 1.2 | & |
58 | |||
59 | # Gridding parameters | ||
60 | &PARM04 | ||
61 | usingCartesianGrid=.FALSE., | ||
62 | usingSphericalPolarGrid=.TRUE., | ||
63 | delX=128*2.8125, | ||
64 | delY=64*2.8125, | ||
65 | delR=1.E4, 2.5E4, 3.E4, 2.E4, 1.5E4, | ||
66 | Ro_SeaLevel=1.E5, | ||
67 | groundAtK1=.TRUE., | ||
68 | phiMin=-90., | ||
69 | rSphere=6370.E3, | ||
70 | & | ||
71 | |||
72 | # Input datasets | ||
73 | &PARM05 | ||
74 | bathyFile='orography_fullcell.bin', | ||
75 | & |
ViewVC Help | |
Powered by ViewVC 1.1.22 |