Changeset 26545
- Timestamp:
- 11/04/21 17:30:01 (3 years ago)
- Location:
- issm/trunk-jpl/test
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/MITgcm/code_4003/SIZE.h
r26544 r26545 48 48 & OLy = 3, 49 49 & nSx = 1, 50 & nSy = 4,50 & nSy = 1, 51 51 & nPx = 1, 52 & nPy = 1,52 & nPy = 4, 53 53 & Nx = sNx*nSx*nPx, 54 54 & Ny = sNy*nSy*nPy, … … 63 63 & MAX_OLY = OLy ) 64 64 65 INTEGER nobcs66 PARAMETER ( nobcs = 4 )67 -
issm/trunk-jpl/test/MITgcm/input_4003/data
r24125 r26545 1 # :::====================1 # ==================== 2 2 # | Model parameters | 3 3 # ==================== … … 5 5 # Continuous equation parameters 6 6 &PARM01 7 # momstepping=.FALSE., 8 # tempstepping=.FALSE., 9 # saltstepping=.FALSE., 10 Tref = 90*-1.9, 11 Sref = 90*34.4, 7 tRef = 90*-1.9, 8 sRef = 90*34.4, 9 eosType='JMD95Z', 10 HeatCapacity_Cp = 3974.0, 11 rhoConst=1030., 12 gravity=9.81, 12 13 viscAr=1.E-3, 13 #viscAh=600.0, 14 viscAhGrid=0.2, 15 viscA4Grid=0.02, 14 viscAh= 300., 16 15 no_slip_sides=.FALSE., 17 16 no_slip_bottom=.FALSE., 18 diffKhT= 100.0, 19 #diffKhT= 600.0, 17 diffKhT= 10., 20 18 diffKrT=5.E-5, 21 diffKhS= 100.0, 22 #diffKhS= 600.0, 19 diffKhS= 10., 23 20 diffKrS=5.E-5, 24 21 bottomDragQuadratic=2.5E-3, 25 staggerTimestep = .true., 26 tempAdvScheme=30, 27 saltAdvScheme=30, 28 vectorInvariantMomentum = .true., 29 ### momImplVertAdv=.true., 30 tempImplVertAdv=.true., 31 saltImplVertAdv=.true., 32 eosType='JMD95Z', 33 #tAlpha = 2.E-4, 34 #sBeta = 7.4E-4, 35 HeatCapacity_cp = 3974.0, 36 rhoConst=1000., 37 rhoNil=1000., 38 gravity=9.81, 39 gBaro=9.81, 40 implicitDiffusion = .true., 41 implicitViscosity = .true., 42 rigidLid=.FALSE., 22 staggerTimestep = .TRUE., 23 tempAdvScheme=77, 24 saltAdvScheme=77, 25 vectorInvariantMomentum = .TRUE., 26 # momImplVertAdv=.TRUE., 27 #tempImplVertAdv=.TRUE., 28 #saltImplVertAdv=.TRUE., 29 implicitDiffusion = .TRUE., 30 implicitViscosity = .TRUE., 31 selectImplicitDrag = 2, 43 32 implicitFreeSurface=.TRUE., 44 useRealFreshWaterFlux = . true.,33 useRealFreshWaterFlux = .TRUE., 45 34 exactConserv=.TRUE., 46 # start nlfs with rstar 47 #select_rStar=2, 35 # start NLFS (without rstar) 48 36 nonlinFreeSurf=4, 49 hFacInf=0. 2,37 hFacInf=0.05, 50 38 hFacSup=2.0, 51 # end nlfs with rstar52 hFacMin=0. 2,39 # end NLFS 40 hFacMin=0.10, 53 41 #integr_GeoPot = 1, 54 #ivdc_kappa = 0.005., 55 #implicitDiffusion=.true., 56 #implicitViscosity=.true., 57 #useCDscheme=.true., 58 #nonHydrostatic=.true., 59 #allowFreezing = .true., 42 ivdc_kappa = 1., 43 useJamartWetPoints=.TRUE., 44 selectBotDragQuadr=1, 45 #bottomVisc_pCell=.TRUE., 46 #- 60 47 readBinaryPrec=64, 61 #writeBinaryPrec=64, 62 debuglevel = 5, 63 selectCoriMap = 0, 64 f0 = 0.0, 48 writeBinaryPrec=64, 65 49 useSingleCpuIO=.TRUE., 66 globalFiles=.TRUE., 67 useCoriolis=.TRUE., 50 #globalFiles=.TRUE., 51 debuglevel = 2, 52 plotLevel = 0, 68 53 & 69 54 … … 72 57 cg2dMaxIters=300, 73 58 cg2dTargetResidual=1.E-11, 74 # cg2dMinColumnEps = 1. 75 # pReleaseVisc = 1. 76 #cg3dMaxIters=40, 77 #cg3dTargetResidual=1.E-13, 59 #printResidualFreq=100, 78 60 & 79 61 80 62 #Time stepping parameters 81 63 &PARM03 82 niter0=11395, 83 nTimeSteps=1440, 84 # nTimeSteps=20, 85 # startTime=0., 86 # startTime = 0., 87 #startTime=62208000., 88 #startTime=15552000., 89 #endTime=31104000., 90 #endTime=86400., 91 # deltaT=300.0, 92 deltaT=60.0, 93 ### deltaT=1200.0, 94 forcing_In_AB = .false., 95 #rCD = 0.9896, 64 #nIter0=2880, 65 #nTimeSteps=288, 66 #startTime=0., 67 #endTime=2592000., 68 deltaT=300.0, 69 forcing_In_AB = .FALSE., 96 70 abEps=0.1, 97 cAdjFreq = -1., 98 pChkptFreq=3000, 99 chkptFreq=3000.0, 100 # pChkptFreq=600.0, 101 # chkptFreq=600.0, 102 dumpFreq=3000, 103 # monitorFreq=15552000., 104 # monitorSelect=0, 105 monitorFreq=60., 71 pChkptFreq= 2592000., 72 chkptFreq = 2592000., 73 dumpFreq = 86400., 74 monitorFreq= 7200., 75 monitorSelect=1, 76 monitorFreq=1800., 77 dumpFreq = 10800., 78 #- for testing purpose: 79 nIter0=2898, 80 nTimeSteps=20, 81 #dumpFreq = 1200., 82 monitorFreq=1., 106 83 & 107 84 108 85 # Gridding parameters 109 86 &PARM04 110 useMin4hFacEdges=.TRUE.,111 87 usingCartesianGrid=.FALSE., 112 88 usingSphericalPolarGrid=.TRUE., 113 #delR = 10., 10., 10., 10., 10., 10., 10., 10.01, 10.03, 10.11,114 #10.32, 10.8, 11.76, 13.42, 16.04, 19.82, 24.85, 31.1, 38.42, 46.5,115 #55., 63.5, 71.58, 78.9, 85.15, 90.18, 93.96, 96.58, 98.25, 99.25,116 #100.01, 101.33, 104.56, 111.33, 122.83, 139.09, 158.94, 180.83, 203.55, 226.5,117 #249.5, 272.5, 295.5, 318.5, 341.5, 364.5, 387.5, 410.5, 433.5, 456.5,118 89 delR=90*10., 119 delX= 3*.125,90 delX=1*.125, 120 91 delY=200*.0078125, 121 # delYFile='dlat.bin', 122 # delY=30*10e3, 123 xgOrigin = -105.5, 124 ygOrigin = -75.4457, 92 xgOrigin = 0., 93 ygOrigin = -75.5, 94 pCellMix_select=20, 95 pcellMix_delR = 10., 96 #interViscAr_pCell = .TRUE., 97 #interDiffKr_pCell = .TRUE., 98 pCellMix_viscAr = 90*1.E-3, 99 pCellMix_diffKr = 90*1.E-4, 100 # current default: diffKzT=5.E-5, viscAz=1.E-3, 125 101 & 126 102 127 103 # Input datasets 128 104 &PARM05 129 bathyFile='bathy .box',130 hydrogThetaFile='t heta.init',131 hydrogSaltFile ='salt.init',132 pSurfInitFile='etainit.round.bin', 105 bathyFile='bathy_flat.bin', 106 hydrogThetaFile='temp_ini.bin', 107 hydrogSaltFile ='salt_ini.bin', 108 # pSurfInitFile ='etai_ini.bin' 133 109 & -
issm/trunk-jpl/test/MITgcm/input_4003/data.shelfice
r26452 r26545 16 16 SHELFICEmergeThreshold = 0.10, 17 17 #-- 18 SHELFICEtopoFile='icetopo.bin', 19 SHELFICEwriteState = .TRUE., 18 SHELFICEtopoFile='shelficeTopo.Lin.bin', 19 SHELFICEmassFile='shelficeMass.Lin.bin', 20 SHELFICEMassDynTendFile='shelfice_dMdt.r02.bin', 21 #SHELFICEwriteState = .TRUE., 20 22 & -
issm/trunk-jpl/test/MITgcm/input_4003/data_uncoupled
r24123 r26545 1 # :::====================1 # ==================== 2 2 # | Model parameters | 3 3 # ==================== … … 5 5 # Continuous equation parameters 6 6 &PARM01 7 # momstepping=.FALSE., 8 # tempstepping=.FALSE., 9 # saltstepping=.FALSE., 10 Tref = 90*-1.9, 11 Sref = 90*34.4, 7 tRef = 90*-1.9, 8 sRef = 90*34.4, 9 eosType='JMD95Z', 10 HeatCapacity_Cp = 3974.0, 11 rhoConst=1030., 12 gravity=9.81, 12 13 viscAr=1.E-3, 13 #viscAh=600.0, 14 viscAhGrid=0.2, 15 viscA4Grid=0.02, 14 viscAh= 300., 16 15 no_slip_sides=.FALSE., 17 16 no_slip_bottom=.FALSE., 18 diffKhT= 100.0, 19 #diffKhT= 600.0, 17 diffKhT= 10., 20 18 diffKrT=5.E-5, 21 diffKhS= 100.0, 22 #diffKhS= 600.0, 19 diffKhS= 10., 23 20 diffKrS=5.E-5, 24 21 bottomDragQuadratic=2.5E-3, 25 staggerTimestep = .true., 26 tempAdvScheme=30, 27 saltAdvScheme=30, 28 vectorInvariantMomentum = .true., 29 ### momImplVertAdv=.true., 30 tempImplVertAdv=.true., 31 saltImplVertAdv=.true., 32 eosType='JMD95Z', 33 #tAlpha = 2.E-4, 34 #sBeta = 7.4E-4, 35 HeatCapacity_cp = 3974.0, 36 rhoConst=1000., 37 rhoNil=1000., 38 gravity=9.81, 39 gBaro=9.81, 40 implicitDiffusion = .true., 41 implicitViscosity = .true., 42 rigidLid=.FALSE., 22 staggerTimestep = .TRUE., 23 tempAdvScheme=77, 24 saltAdvScheme=77, 25 vectorInvariantMomentum = .TRUE., 26 # momImplVertAdv=.TRUE., 27 #tempImplVertAdv=.TRUE., 28 #saltImplVertAdv=.TRUE., 29 implicitDiffusion = .TRUE., 30 implicitViscosity = .TRUE., 31 selectImplicitDrag = 2, 43 32 implicitFreeSurface=.TRUE., 44 useRealFreshWaterFlux = . true.,33 useRealFreshWaterFlux = .TRUE., 45 34 exactConserv=.TRUE., 46 # start nlfs with rstar 47 #select_rStar=2, 35 # start NLFS (without rstar) 48 36 nonlinFreeSurf=4, 49 hFacInf=0. 2,37 hFacInf=0.05, 50 38 hFacSup=2.0, 51 # end nlfs with rstar52 hFacMin=0. 2,39 # end NLFS 40 hFacMin=0.10, 53 41 #integr_GeoPot = 1, 54 #ivdc_kappa = 0.005., 55 #implicitDiffusion=.true., 56 #implicitViscosity=.true., 57 #useCDscheme=.true., 58 #nonHydrostatic=.true., 59 #allowFreezing = .true., 42 ivdc_kappa = 1., 43 useJamartWetPoints=.TRUE., 44 selectBotDragQuadr=1, 45 #bottomVisc_pCell=.TRUE., 46 #- 60 47 readBinaryPrec=64, 61 #writeBinaryPrec=64, 62 debuglevel = 5, 63 selectCoriMap = 0, 64 f0 = 0.0, 48 writeBinaryPrec=64, 65 49 useSingleCpuIO=.TRUE., 66 globalFiles=.TRUE., 67 useCoriolis=.TRUE., 50 #globalFiles=.TRUE., 51 debuglevel = 2, 52 plotLevel = 0, 68 53 & 69 54 … … 72 57 cg2dMaxIters=300, 73 58 cg2dTargetResidual=1.E-11, 74 # cg2dMinColumnEps = 1. 75 # pReleaseVisc = 1. 76 #cg3dMaxIters=40, 77 #cg3dTargetResidual=1.E-13, 59 #printResidualFreq=100, 78 60 & 79 61 80 62 #Time stepping parameters 81 63 &PARM03 82 niter0=11395, 83 # nTimeSteps=25920000, 64 #nIter0=2880, 65 #nTimeSteps=288, 66 #startTime=0., 67 #endTime=2592000., 68 deltaT=300.0, 69 forcing_In_AB = .FALSE., 70 abEps=0.1, 71 pChkptFreq= 2592000., 72 chkptFreq = 2592000., 73 dumpFreq = 86400., 74 monitorFreq= 7200., 75 monitorSelect=1, 76 monitorFreq=1800., 77 dumpFreq = 10800., 78 #- for testing purpose: 79 nIter0=2898, 84 80 nTimeSteps=20, 85 # startTime=0., 86 # startTime = 0., 87 #startTime=62208000., 88 #startTime=15552000., 89 #endTime=31104000., 90 deltaT=300.0, 91 ### deltaT=1200.0, 92 forcing_In_AB = .false., 93 #rCD = 0.9896, 94 abEps=0.1, 95 cAdjFreq = -1., 96 pChkptFreq=3000, 97 chkptFreq=3000.0, 98 # pChkptFreq=600.0, 99 # chkptFreq=600.0, 100 dumpFreq=3000, 101 # monitorFreq=15552000., 102 # monitorSelect=0, 103 # monitorFreq=31104000., 81 #dumpFreq = 1200., 82 monitorFreq=1., 104 83 & 105 84 106 85 # Gridding parameters 107 86 &PARM04 108 useMin4hFacEdges=.TRUE.,109 87 usingCartesianGrid=.FALSE., 110 88 usingSphericalPolarGrid=.TRUE., 111 #delR = 10., 10., 10., 10., 10., 10., 10., 10.01, 10.03, 10.11,112 #10.32, 10.8, 11.76, 13.42, 16.04, 19.82, 24.85, 31.1, 38.42, 46.5,113 #55., 63.5, 71.58, 78.9, 85.15, 90.18, 93.96, 96.58, 98.25, 99.25,114 #100.01, 101.33, 104.56, 111.33, 122.83, 139.09, 158.94, 180.83, 203.55, 226.5,115 #249.5, 272.5, 295.5, 318.5, 341.5, 364.5, 387.5, 410.5, 433.5, 456.5,116 89 delR=90*10., 117 delX= 3*.125,90 delX=1*.125, 118 91 delY=200*.0078125, 119 # delYFile='dlat.bin', 120 # delY=30*10e3, 121 xgOrigin = -105.5, 122 ygOrigin = -75.4457, 92 xgOrigin = 0., 93 ygOrigin = -75.5, 94 pCellMix_select=20, 95 pcellMix_delR = 10., 96 #interViscAr_pCell = .TRUE., 97 #interDiffKr_pCell = .TRUE., 98 pCellMix_viscAr = 90*1.E-3, 99 pCellMix_diffKr = 90*1.E-4, 100 # current default: diffKzT=5.E-5, viscAz=1.E-3, 123 101 & 124 102 125 103 # Input datasets 126 104 &PARM05 127 bathyFile='bathy .box',128 hydrogThetaFile='t heta.init',129 hydrogSaltFile ='salt.init',130 pSurfInitFile='etainit.round.bin', 105 bathyFile='bathy_flat.bin', 106 hydrogThetaFile='temp_ini.bin', 107 hydrogSaltFile ='salt_ini.bin', 108 # pSurfInitFile ='etai_ini.bin' 131 109 & -
issm/trunk-jpl/test/MITgcm/input_4003/eedata
r24123 r26545 1 1 # Example "eedata" file 2 2 # Lines beginning "#" are comments 3 # nTx - No. threads per process in X 4 # nTy - No. threads per process in Y 3 # nTx :: No. threads per process in X 4 # nTy :: No. threads per process in Y 5 # debugMode :: print debug msg (sequence of S/R calls) 5 6 &EEPARMS 6 7 useCoupler=.TRUE., 7 8 debugMode=.TRUE., 8 9 & 9 # Note: Some systems use & as the 10 # namelist terminator. Other systems 11 # use a / character (as shown here). 10 # Note: Some systems use & as the namelist terminator (as shown here). 11 # Other systems use a / character. -
issm/trunk-jpl/test/MITgcm/input_4003/eedata_uncoupled
r24121 r26545 1 1 # Example "eedata" file 2 2 # Lines beginning "#" are comments 3 # nTx - No. threads per process in X 4 # nTy - No. threads per process in Y 3 # nTx :: No. threads per process in X 4 # nTy :: No. threads per process in Y 5 # debugMode :: print debug msg (sequence of S/R calls) 5 6 &EEPARMS 6 7 useCoupler=.FALSE., 7 # debugMode=.TRUE.,8 8 & 9 # Note: Some systems use & as the 10 # namelist terminator. Other systems 11 # use a / character (as shown here). 9 # Note: Some systems use & as the namelist terminator (as shown here). 10 # Other systems use a / character. -
issm/trunk-jpl/test/NightlyRun/test4003.m
r26543 r26545 3 3 % 4 4 %Script control parameters 5 steps=1:11; 6 steps=1:5; 7 steps=6; 5 steps=1:12; 6 steps=1:6; 8 7 final_time=1/365; 9 8 10 9 %To download and recompile MITgcm from scratch: 11 %!rm -rf ${ISSM_DIR}/test/MITgcm/install_dngoldberg12 %!rm -rf ${ISSM_DIR}/test/MITgcm/build/*13 %!rm -rf Models10 !rm -rf ${ISSM_DIR}/test/MITgcm/install_dngoldberg 11 !rm -rf ${ISSM_DIR}/test/MITgcm/build/* 12 !rm -rf Models 14 13 15 14 %Organizer … … 21 20 % {{{ Parameters: 22 21 if perform(org,'Parameters'), 23 Nx= 3; %number of longitude cells22 Nx=1; %number of longitude cells 24 23 Ny=200; %number of latitude cells 25 24 Nz=90; %number of MITgcm vertical cells 26 25 nPx=1; %number of MITgcm processes to use in x direction 27 nPy= 1; %number of MITgcm processes to use in y direction26 nPy=4; %number of MITgcm processes to use in y direction 28 27 xgOrigin=1; %origin of longitude 29 28 ygOrigin=1; %origin of latitude
Note:
See TracChangeset
for help on using the changeset viewer.