| 1 | dimitri | 1.1 | # ==================== | 
| 2 |  |  | # | Model parameters | | 
| 3 |  |  | # ==================== | 
| 4 |  |  | # | 
| 5 |  |  | # Continuous equation parameters | 
| 6 |  |  | &PARM01 | 
| 7 |  |  | tRef=15*20., | 
| 8 |  |  | sRef=15*35., | 
| 9 |  |  | viscAh =3.E5, | 
| 10 |  |  | viscAr =1.E-3, | 
| 11 |  |  | diffKhT=0., | 
| 12 |  |  | diffK4T=0., | 
| 13 |  |  | diffKrT=3.E-5, | 
| 14 |  |  | diffKhS=0., | 
| 15 |  |  | diffK4S=0., | 
| 16 |  |  | diffKrS=3.E-5, | 
| 17 |  |  | ivdc_kappa=10., | 
| 18 |  |  | implicitDiffusion=.TRUE., | 
| 19 |  |  | gravity=9.81, | 
| 20 |  |  | rhonil=1035., | 
| 21 |  |  | eosType='JMD95Z', | 
| 22 |  |  | #staggerTimeStep=.TRUE., | 
| 23 |  |  | vectorInvariantMomentum=.TRUE., | 
| 24 |  |  | implicitFreeSurface=.TRUE., | 
| 25 |  |  | exactConserv=.TRUE., | 
| 26 |  |  | select_rStar=2, | 
| 27 |  |  | nonlinFreeSurf=4, | 
| 28 |  |  | hFacInf=0.2, | 
| 29 |  |  | hFacSup=2.0, | 
| 30 |  |  | useRealFreshWaterFlux=.TRUE., | 
| 31 |  |  | allowFreezing=.TRUE., | 
| 32 |  |  | hFacMin=.1, | 
| 33 |  |  | hFacMinDr=20., | 
| 34 |  |  | readBinaryPrec=64, | 
| 35 |  |  | globalfiles=.true. | 
| 36 |  |  | tempAdvScheme=33 | 
| 37 |  |  | saltAdvScheme=33 | 
| 38 |  |  | & | 
| 39 |  |  |  | 
| 40 |  |  | # Elliptic solver parameters | 
| 41 |  |  | &PARM02 | 
| 42 |  |  | cg2dMaxIters=200, | 
| 43 |  |  | #cg2dTargetResidual=1.E-9, | 
| 44 |  |  | cg2dTargetResWunit=1.E-14, | 
| 45 |  |  | & | 
| 46 |  |  |  | 
| 47 |  |  | # Time stepping parameters | 
| 48 |  |  | &PARM03 | 
| 49 |  |  | nIter0=72000, | 
| 50 |  |  | nTimeSteps=3, | 
| 51 |  |  | deltaT   =1200., | 
| 52 |  |  | abEps = 0.1, | 
| 53 |  |  | forcing_In_AB=.FALSE., | 
| 54 |  |  | pChkptFreq  =1728000., | 
| 55 |  |  | #chkptFreq   =311040000., | 
| 56 |  |  | taveFreq    =1244160000., | 
| 57 |  |  | dumpFreq    =622080000., | 
| 58 |  |  | monitorFreq =31104000., | 
| 59 |  |  | periodicExternalForcing=.TRUE., | 
| 60 |  |  | externForcingPeriod=2592000., | 
| 61 |  |  | externForcingCycle=31104000., | 
| 62 |  |  | # 2 months restoring timescale for temperature | 
| 63 |  |  | tauThetaClimRelax = 5184000., | 
| 64 |  |  | # 2yrs restoring timescale for salinity | 
| 65 |  |  | tauSaltClimRelax = 62208000., | 
| 66 |  |  | monitorFreq =1., | 
| 67 |  |  | & | 
| 68 |  |  |  | 
| 69 |  |  | # Gridding parameters | 
| 70 |  |  | &PARM04 | 
| 71 |  |  | usingCartesianGrid=.FALSE., | 
| 72 |  |  | usingSphericalPolarGrid=.FALSE., | 
| 73 |  |  | usingCurvilinearGrid=.TRUE., | 
| 74 |  |  | delR= 50., 70., 100., 140., 190., | 
| 75 |  |  | 240., 290., 340., 390., 440., | 
| 76 |  |  | 490., 540., 590., 640., 690., | 
| 77 |  |  | & | 
| 78 |  |  |  | 
| 79 |  |  | # Input datasets | 
| 80 |  |  | &PARM05 | 
| 81 |  |  | bathyFile      ='bathy_Hmin50.bin', | 
| 82 |  |  | hydrogThetaFile='lev_T_cs_15k.bin', | 
| 83 |  |  | hydrogSaltFile ='lev_S_cs_15k.bin', | 
| 84 |  |  | zonalWindFile  ='trenberth_taux.bin', | 
| 85 |  |  | meridWindFile  ='trenberth_tauy.bin', | 
| 86 |  |  | thetaClimFile  ='lev_surfT_cs_12m.bin', | 
| 87 |  |  | saltClimFile   ='lev_surfS_cs_12m.bin', | 
| 88 |  |  | surfQFile      ='shiQnet_cs32.bin', | 
| 89 |  |  | EmPmRFile      ='shiEmPR_cs32.bin', | 
| 90 |  |  | & |