Instructions for setting up a quarter degree integration forced 1. by NCEP fields and 2. by SIO fields. ########################################################################## # getting the code from anonymous CVS server cvs co MITgcm_contrib/quarter_degree_global cvs co MITgcm_code cvs co MITgcm/verification/global_with_exf cvs co MITgcm/verification/testreport ############################################################################# # 1-CPU verification/global_with_exf including pkg/diagnostics cd MITgcm/verification/global_with_exf/build \rm * \cp ../code/* . \cp ../../../../MITgcm_contrib/quarter_degree_global/code/packages.conf . \cp ../../../../MITgcm_contrib/quarter_degree_global/code/ECCO_CPPOPTIONS.h . \cp ../../../../MITgcm_contrib/quarter_degree_global/code/CPP_OPTIONS.h . \cp ../../../../MITgcm_contrib/quarter_degree_global/code/DIAGNOSTICS_SIZE.h . \cp ../../../../MITgcm_contrib/quarter_degree_global/code/SIZE.h_test1 . \mv SIZE.h_test1 SIZE.h ../../../tools/genmake2 make depend make -j cd .. mkdir run cd run \rm * \cp ../build/mitgcmuv . \ln -sf ../input/* . \rm data \cp ../../../../MITgcm_contrib/quarter_degree_global/input/data_test data \cp ../../../../MITgcm_contrib/quarter_degree_global/input/data.pkg . \cp ../../../../MITgcm_contrib/quarter_degree_global/input/data.kpp . \cp ../../../../MITgcm_contrib/quarter_degree_global/input/data.diagnostics . mitgcmuv >& output.txt there is comparison output in ../../../../MITgcm_contrib/quarter_degree_global/results/output.txt ############################################################################# # mpi test verification/global_with_exf including pkg/diagnostics on columbia cd MITgcm/verification/global_with_exf/build \rm * \cp ../code/* . \cp ../../../../MITgcm_contrib/quarter_degree_global/code/packages.conf . \cp ../../../../MITgcm_contrib/quarter_degree_global/code/ECCO_CPPOPTIONS.h . \cp ../../../../MITgcm_contrib/quarter_degree_global/code/CPP_OPTIONS.h . \cp ../../../../MITgcm_contrib/quarter_degree_global/code/SIZE.h_test . \mv SIZE.h_test SIZE.h ../../../tools/genmake2 -of ../../../tools/build_options/linux_ia64_ifort+mpi_altix_nas make depend make -j cd .. mkdir run cd run \rm * \cp ../build/mitgcmuv . \ln -sf ../input/* . \rm data \cp ../../../../MITgcm_contrib/quarter_degree_global/input/data_test data mpirun -np 2 mitgcmuv there is comparison output in ../../../../MITgcm_contrib/quarter_degree_global/results/STDOUT.0000 ########################################################################## # 1440x1088 204-cpu (120x12)x(64x17) 1/4-deg integration on altix with efc # requires input files in subdirectory iter69 to be at same level as MITgcm # and MITgcm_contrib. A copy of iter69 resides on columbia10:/nobackup10 # The integration length is 4018 days, 1992-2002. # First three days need to be integrated using deltat=200. cd MITgcm mkdir bin exe cd bin \rm * \cp ../../MITgcm_contrib/quarter_degree_global/code/* . ../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_nas make depend make -j cd ../exe \rm * \cp ../bin/mitgcmuv . \cp ../../MITgcm_contrib/quarter_degree_global/input/* . \ln -sf ../../iter69/* . qsub job_altix ########################################################################## # 2880x2176 480-cpu (96x30)x(136x16) 1/4-deg integration on altix with efc # requires input files in sudirectory iter69 to be at same level as MITgcm # and MITgcm_contrib. A copy of iter69 with 1/8-deg files resides on # columbia18:/nobackup18/menemenl/eight. # The integration length is 3 days. # First three days need to be integrated using deltat=100. cd MITgcm mkdir bin exe cd bin \rm * \cp ../../MITgcm_contrib/quarter_degree_global/code/* . \cp SIZE.h_8_480 SIZE.h ../tools/genmake2 -of ../tools/build_options/linux_ia64_efc+mpi_altix make depend make -j cd ../exe \rm * \cp ../bin/mitgcmuv . \cp ../../MITgcm_contrib/quarter_degree_global/input/* . \ln -sf ../../iter69/* . \cp data.exf_clim_8 data.exf_clim \cp data_SIO_8 data qsub job_altix ########################################################################## # 2880x2176 960-cpu (96x30)x(68x32) 1/4-deg integration on altix with efc # requires input files in sudirectory iter69 to be at same level as MITgcm # and MITgcm_contrib. A copy of iter69 with 1/8-deg files resides on # columbia18:/nobackup18/menemenl/eight. # The integration length is 3 days. # First three days need to be integrated using deltat=100. module purge module load modules scsl.1.5.0.0 intel-comp.8.0.070 mpt.1.11-85 setenv MPI_DSM_DISTRIBUTE cd MITgcm mkdir bin exe cd bin \rm * \cp ../../MITgcm_contrib/quarter_degree_global/code/* . \cp SIZE.h_8_960 SIZE.h ../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_nas make depend make -j cd ../exe \cp ../bin/mitgcmuv . \cp ../../MITgcm_contrib/quarter_degree_global/input/* . \ln -sf ../../iter69/* . \cp data.exf_clim_8 data.exf_clim \cp data_SIO_8 data modpickup_8 mpirun -v -a a2k c18 480 runon 11-490 mitgcmuv : c20 480 runon 11-490 mitgcmuv mpirun -v -a a2k c18 508 runon 11-490 mitgcmuv : c20 508 runon 11-490 mitgcmuv on columbia18: ln -sf /nobackup18/menemenl/eighth /tmp/eighth on columbia20: ln -sf /nobackup20/menemenl/eighth /tmp/eighth cd ~/exe \cp /tmp/eighth/MITgcm/exe/mitgcmuv960 . \cp /tmp/eighth/MITgcm_contrib/quarter_degree_global/input/* . \ln -sf /tmp/eighth/iter69/* . \cp data.exf_clim_8 data.exf_clim \cp data_SIO_8 data \ln -sf /tmp/eighth/MITgcm/exe/pick* . mpirun -v -a a2k -d /tmp/eighth/MITgcm/exe c18 480 runon 11-490 mitgcmuv960 : \ c20 480 runon 11-490 mitgcmuv960 >&! output.txt & \cp data.test data scp data:columbia17:/nobackup/menemenl/eighth/MITgcm/exe scp data:columbia19:/nobackup/menemenl/eighth/MITgcm/exe scp data:columbia20:/nobackup/menemenl/eighth/MITgcm/exe mpirun -v -a a2k -d /nobackup/menemenl/eighth/MITgcm/exe \ c18 480 runon 11-490 mitgcmuv1440 : \ c19 480 runon 11-490 mitgcmuv1440 : \ c20 480 runon 11-490 mitgcmuv1440 \ >&! output.test1440 mkdir tests mv STDOUT.0000 output.test1440 tests/stdout1440 \rm STD* mpirun -v -a a2k -d /nobackup/menemenl/eighth/MITgcm/exe \ c17 480 runon 11-490 mitgcmuv1920 : \ c18 480 runon 11-490 mitgcmuv1920 : \ c19 480 runon 11-490 mitgcmuv1920 : \ c20 480 runon 11-490 mitgcmuv1920 \ >&! output.test1920 mv STDOUT.0000 output.test1920 tests/stdout1920 \rm STD* \cp data.run data scp data:columbia17:/nobackup/menemenl/eighth/MITgcm/exe scp data:columbia19:/nobackup/menemenl/eighth/MITgcm/exe scp data:columbia20:/nobackup/menemenl/eighth/MITgcm/exe /opt/sgi/mpt/1.11-85/bin/mpirun.real -v -a a2k \ -d /nobackup/menemenl/eighth/MITgcm/exe \ c17 480 runon 11-490 mitgcmuv1920 : \ c18 480 runon 11-490 mitgcmuv1920 : \ c19 480 runon 11-490 mitgcmuv1920 : \ c20 480 runon 11-490 mitgcmuv1920 \ >&! output.run \cp data.test data scp data columbia19:/nobackup/menemenl/eighth/MITgcm/exe scp data columbia20:/nobackup/menemenl/eighth/MITgcm/exe mpirun -v -a a2k -d /nobackup/menemenl/eighth/MITgcm/exe \ c18 480 runon 9-511 mitgcmuv1440 : \ c19 480 runon 9-511 mitgcmuv1440 : \ c20 480 runon 9-511 mitgcmuv1440 \ >&! output.test1440 mkdir tests mv STDOUT.0000 output.test1440 tests/stdout1440 \rm STD* \cp data.run data scp data columbia17:/nobackup/menemenl/eighth/MITgcm/exe scp data columbia19:/nobackup/menemenl/eighth/MITgcm/exe mpirun -v -a a2k -d /nobackup/menemenl/eighth/MITgcm/exe \ c17 480 runon 9-511 mitgcmuv1440 : \ c18 480 runon 9-511 mitgcmuv1440 : \ c19 480 runon 9-511 mitgcmuv1440 \ >&! output.run fails for 1/8 > viscA4GridMin=3.e-4, > viscA4GridMax=4.e-3, > viscC4Leith=0.05, > deltat=300 fails: < viscA4GridMin=6.e-4, < viscA4GridMax=8.e-3, < viscC4Leith=0.1, fails: < viscA4GridMin=6.e-4, < viscA4GridMax=8.e-3, < viscC4Leith=0.2, fails: < viscA4GridMin=1.e-3, < viscA4GridMax=1.e-2, < viscC4Leith=0.2, fails: < viscA4GridMin=1.e-3, < viscA4GridMax=1.e-2, < viscC4Leith=0.3, fails: < viscA4GridMin=1.2e-3, < viscA4GridMax=1.6e-2, < viscC4Leith=0.2, fails: < viscA4GridMin=1.2e-3, < viscA4GridMax=1.6e-2, < viscC4Leith=0.4, fails < viscA4GridMin=1.2e-3, < viscA4GridMax=1.2e-2, < viscC4Leith=1.2, fails: < viscA4GridMin=2e-3, < viscA4GridMax=1e-2, < viscC4Leith=0.1, fails: < viscA4GridMin=2e-3, < viscA4GridMax=1e-2, < viscC4Leith=0.2, fails: < viscA4GridMin=4e-3, < viscA4GridMax=8e-3, < viscC4Leith=0.2, fails: viscA4GridMin=5e-3, viscA4GridMax=1e-2, viscC4Leith=.2, starting with time step 50112, try: (this mean day 174 becomes day 116) viscA4GridMin=6.e-4, viscA4GridMax=8.e-3, viscC4Leith=0.1, deltat=200 %%%%%%%%%%%%%%%%%%%%%%%%%%%%% 1/8 timing tests: columbia20 scsl.1.5.0.0 c18 intel-comp.8.0.070 240-CPU default mpt.1.11-85 ------------------------------------------------ MAIN LOOP 7115 5066 EXF_GETFORCING 372 38 DO_OCEANIC_PHYS 1221 1272 STATEKPP 509 532 KPPMIX 370 374 DYNAMICS 1527 1429 SOLVE_FOR_PRESSURE 1804 206 BLOCKING_EXCHANGES 117 115 THERMODYNAMICS 1166 1100 DO_THE_MODEL_IO 173 177 WRITE_CHECKPOINT 341 333 columbia20 scsl.1.5.0.0 c18 intel-comp.8.0.070 480-CPU default mpt.1.11-85 c19 default c17 default, other processes c19 new ------------------------------------------------------------------------- MAIN LOOP 4024 2917 5379 8895 3571 EXF_GETFORCING 218 33 880 261 227 DO_OCEANIC_PHYS 622 619 621 913 620 STATEKPP 269 265 268 310 265 KPPMIX 190 188 189 291 189 DYNAMICS 747 718 745 1456 714 SOLVE_FOR_PRESSURE 966 119 1703 2748 900 BLOCKING_EXCHANGES 64 70 392 375 70 THERMODYNAMICS 553 582 566 1510 584 DO_THE_MODEL_IO 178 156 149 184 152 WRITE_CHECKPOINT 364 332 864 420 333 CPU 240 480 960 1440 1920 ------------------------------------------------------------- MAIN LOOP 5066 2917 2227 1551 1276 EXF_GETFORCING 38 33 126 60 83 DO_OCEANIC_PHYS 1272 619 320 390 186 STATEKPP 532 265 138 181 77 KPPMIX 374 188 98 112 53 DYNAMICS 1429 718 313 319 173 SOLVE_FOR_PRESSURE 206 119 483 169 310 BLOCKING_EXCHANGES 115 70 53 29 28 THERMODYNAMICS 1100 582 294 279 172 DO_THE_MODEL_IO 177 156 131 202 155 WRITE_CHECKPOINT 333 332 302 524 352 4 clusters timeave Dec 16 08:17 Ttave.0000051840.data Dec 16 08:47 Ttave.0000053136.data Dec 16 09:17 Ttave.0000054432.data Dec 16 09:53 Ttave.0000055728.data Dec 17 14:08 Ttave.0000075168.data Dec 17 14:33 Ttave.0000076464.data Dec 17 14:58 Ttave.0000077760.data Dec 17 15:29 Ttave.0000079056.data 1 cluster timeave Dec 16 22:14 Ttave.0000059616.data Dec 16 23:16 Ttave.0000060912.data Dec 17 00:19 Ttave.0000062208.data Dec 17 01:27 Ttave.0000063504.data Dec 17 02:29 Ttave.0000064800.data Dec 17 03:31 Ttave.0000066096.data Dec 17 04:39 Ttave.0000067392.data Dec 17 05:42 Ttave.0000068688.data Dec 17 06:44 Ttave.0000069984.data 2 cluster timeave Dec 17 21:19 Ttave.0000086832.data Dec 17 21:55 Ttave.0000088128.data 36 Dec 17 22:30 Ttave.0000089424.data 35 Dec 17 23:12 Ttave.0000090720.data 42 Dec 17 23:47 Ttave.0000092016.data 35 Dec 18 00:23 Ttave.0000093312.data 36 Dec 18 01:04 Ttave.0000094608.data 41 Dec 18 01:39 Ttave.0000095904.data 35 2 cluster diagnostics Dec 18 02:26 THETA.0000086832.data Dec 18 03:05 THETA.0000088128.data 39 Dec 18 03:43 THETA.0000089424.data 38 Dec 18 04:27 THETA.0000090720.data 44 Dec 18 05:06 THETA.0000092016.data 39 Dec 18 05:45 THETA.0000093312.data 39 Dec 18 06:29 THETA.0000094608.data 44 Dec 18 07:08 THETA.0000095904.data 39 %%%%%%%%%%%%%%%%% e1=readbin('test_oldstyle/ETAtave.0000000018.data',[90 40]); e2=readbin('ETAN.0000000018.data',[90 40]); e16=readbin('test_oldstyle/Eta.0000000016.data',[90 40]); e17=readbin('test_oldstyle/Eta.0000000017.data',[90 40]); e18=readbin('test_oldstyle/Eta.0000000018.data',[90 40]); minmax(e1-e2) minmax(e2-(e16+e17+e18)/3) e1=readbin('test_oldstyle/sFluxtave.0000000018.data',[90 40]); e2=readbin('SFLUX.0000000018.data',[90 40]); minmax(e1-e2) e1=readbin('test_oldstyle/tFluxtave.0000000018.data',[90 40]); e2=readbin('TFLUX.0000000018.data',[90 40]); minmax(e1-e2) e1=readbin('TRELAX.0000000018.data',[90 40]); subplot(211),mypcolor(e1');colorbar e1=readbin('SRELAX.0000000018.data',[90 40]); subplot(211),mypcolor(e1');colorbar e1=readbin('TICE.0000000018.data',[90 40]); subplot(211),mypcolor(e1');colorbar e1=readbin('test_oldstyle/wVeltave.0000000018.data',[90 40 15]); e2=readbin('WVEL.0000000018.data',[90 40 15]); e16=readbin('test_oldstyle/W.0000000016.data',[90 40 15]); e17=readbin('test_oldstyle/W.0000000017.data',[90 40 15]); e18=readbin('test_oldstyle/W.0000000018.data',[90 40 15]); minmax(e1-e2) minmax(e2-(e16+e17+e18)/3)