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 | # ==================== |
| 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 | Integr_GeoPot=1, |
| 23 | 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 | saltAdvScheme=77, |
| 35 | & |
| 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 | nTimeSteps=10, |
| 48 | #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 | monitorFreq=1.0, |
| 57 | & |
| 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 |