| 1 | Instructions for setting up a cube-sphere integration with sea-ice. | Instructions for setting up a cube-sphere integration with sea-ice. | 
| 2 |  |  | 
| 3 | Forcing fields and other input files needed for the CS510 | ########################################################################## | 
| 4 | configuration can be obtained from | # Get forcing files, etc. | 
| 5 | http://ecco2.org/products/model_setups/cs510.html | # Example below is for cube88 on pleiades | 
| 6 |  |  | 
| 7 |  | mkdir cube94 | 
| 8 |  | cd cube94 | 
| 9 |  | ln -sf /nobackupp10/hzhang1/cs510/run_template . | 
| 10 |  | ln -sf /nobackupp10/hzhang1/forcing/jra25 . | 
| 11 |  | ln -sf /nobackupp10/hzhang1/forcing/cube78_forcing . | 
| 12 |  |  | 
| 13 |  |  | 
| 14 |  | ########################################################################## | 
| 15 |  | # getting the code from anonymous CVS server | 
| 16 |  |  | 
| 17 |  | cvs co MITgcm_contrib/high_res_cube/README.cs510 | 
| 18 |  | cvs co MITgcm_contrib/high_res_cube/code-mods | 
| 19 |  | cvs co MITgcm_contrib/high_res_cube/input | 
| 20 |  | cvs co MITgcm_code | 
| 21 |  |  | 
| 22 |  |  | 
| 23 |  | ########################################################################## | 
| 24 |  | # 510*510*6*50, 450-cpu cube sphere on pleiades with ifort | 
| 25 |  |  | 
| 26 |  | cd MITgcm | 
| 27 |  | mkdir build run | 
| 28 |  | cd build | 
| 29 |  | \rm * | 
| 30 |  | \cp ../../MITgcm_contrib/high_res_cube/code-mods/* . | 
| 31 |  | \cp ../../MITgcm_contrib/high_res_cube/code-mods/s450t_34x102/SIZE.h . | 
| 32 |  | \rm README SIZE.h_mpi | 
| 33 |  | ../tools/genmake2 -of ../tools/build_options/linux_amd64_ifort+mpi_ice_nas | 
| 34 |  | make depend | 
| 35 |  | make -j 16 | 
| 36 |  | cd ../run | 
| 37 |  | \rm * | 
| 38 |  | \cp ../../MITgcm_contrib/high_res_cube/input/*data* . | 
| 39 |  | \cp ../../MITgcm_contrib/high_res_cube/input/modpickup . | 
| 40 |  | \cp ../../MITgcm_contrib/high_res_cube/input/job* . | 
| 41 |  | \mv data.seaice_hr data.seaice | 
| 42 |  | \mv data.exf_hr data.exf | 
| 43 |  | \mv data_hr data | 
| 44 |  | \ln -sf ../../run_template/tile* . | 
| 45 |  | \ln -sf ../../run_template/G*dig* . | 
| 46 |  | \ln -sf ../../run_template/*Init.cube78 . | 
| 47 |  | \cp ../build/mitgcmuv . | 
| 48 |  | qsub job_cs510_pleiades450 | 
| 49 |  |  | 
| 50 | ================= | ================= | 
| 51 |  |  | 
| 52 | UPDATES | UPDATES | 
| 53 |  |  | 
| 54 | November 22, 2007. Cube 78 - blend 3 test | March 19, 2012. Cube94 on pleiades | 
| 55 |  | This integration is driven by JRA25 (as was cube92) | 
| 56 |  | but it includes shelfice (as did cube86), | 
| 57 |  | it is based on MITgcm > checkpoint63k (that is post ICE workshop code), | 
| 58 |  | it uses optimized Arctic sea ice parameters from | 
| 59 |  | http://ecco2.jpl.nasa.gov/data1/arctic/output/newcode_050109/JRA25_OPT3/input/data.seaice | 
| 60 |  | http://ecco2.org/manuscripts/2011/NguyenJGR2011.pdf | 
| 61 |  | it uses Southern Ocean ice and ice shelf parameters from JSO24 | 
| 62 |  | /data14/scratch/hzhang/CP/JSO24/input | 
| 63 |  | (bathymetry and initial conditions tbd), | 
| 64 |  | and it uses pkg/bbl and frazil instead of down_slope | 
| 65 |  |  | 
| 66 |  | July 27, 2010.  Cube88 on pleiades - same as cube84 but with: | 
| 67 |  | 1. JRA25 instead of cube84 blend forcing | 
| 68 |  | 2. useSALT_PLUME  = .TRUE., (off in Southern Ocean) | 
| 69 |  | 3. useDOWN_SLOPE  = .TRUE., | 
| 70 |  | 4. split north/south data.sec17saice parameters from regional optimizations | 
| 71 |  | 5. initialize from T/S/U/V/ETAN and sea ice from cube78 Jan 1992 pickup | 
| 72 |  | 6. 3-day 3D diagnostics | 
| 73 |  | 7. SEAICE_no_slip = .TRUE. to keep code stable | 
| 74 |  |  | 
| 75 |  | June 30, 2009.  Cube86 -same as cube85 but with: | 
| 76 |  | 1. updated Southern OCean Bathymetry | 
| 77 |  | 2. ice shelf cavities | 
| 78 |  |  | 
| 79 |  | June 16, 2009. Cube85 -same as cube84 but with: | 
| 80 |  | 1. split north/south data.seaice parameters from regional optimizations | 
| 81 |  | 2. includes pkg/salt_plume in the Arctic and pkg/down_slope globally | 
| 82 |  | 3. initialize from T/S/U/V/ETAN and sea ice files | 
| 83 |  | 4. 3-day 3D diagnostics | 
| 84 |  | 5. SEAICE_no_slip=.TRUE. to keep code stable | 
| 85 |  |  | 
| 86 |  | June 13, 2008. Cube 81 - same as cube79 but with | 
| 87 |  | 1. new values for data.seaice and data.exf | 
| 88 |  | 2. DIFFKR_2_20_1_lat6070_cube81, which has diffKrS/T=3.55e-7 in Arctic | 
| 89 |  | 3. GMRedi turned off | 
| 90 |  |  | 
| 91 |  | December 22, 2007. Cube80 - same as cube79 but with | 
| 92 |  | SEAICE_BREAKS_VERIFICATION (SST relaxation depends on SIarea) | 
| 93 |  | ALLOW_CLIMSST_RELAXATION to Reynold's daily 1/4-deg AVHRR analysis | 
| 94 |  | ALLOW_CLIMSST_RELAXATION to monthly 1-deg WOA05 analysis | 
| 95 |  | December 2, 2007. Cube79 - same as cube78 (blend3) but with | 
| 96 |  | useRelativeWind=.TRUE. and SEAICE_gamma_t=259200 | 
| 97 |  | November 22, 2007. Cube 78 - blend 3 | 
| 98 | November 6, 2007. Cube 77 - same as cube75 but with SEAICE_SALINITY | November 6, 2007. Cube 77 - same as cube75 but with SEAICE_SALINITY | 
| 99 | and era40_ecmwf_blend/cube59_GPCP forcing | and era40_ecmwf_blend/cube59_GPCP forcing | 
| 100 | November 3, 2007. Cube 76 - same as cube77 but 1979-2006 integration | November 3, 2007. Cube 76 - same as cube77 but 1979-2006 integration | 
| 111 | (cube67, cube68, and cube69 were carried out with cube65 code) | (cube67, cube68, and cube69 were carried out with cube65 code) | 
| 112 | April 7 , 2007. Cube64: era40_fullres forcing | April 7 , 2007. Cube64: era40_fullres forcing | 
| 113 |  |  | 
| 114 |  | February 6, 2006: cube56, same as cube 50 with SEAICEuseFlooding | 
| 115 |  |  | 
| 116 | ########################################################################## | February 1, 2006: back to cube50 default (bottomDragQuadratic = 0.002) | 
| 117 | # Get forcing files, etc. | removed UVEL and VVEL from diagnostics | 
|  | # Example below is for cube78 on /nobackup2a/menemenl |  | 
| 118 |  |  | 
| 119 | cd /nobackup2a/menemenl | January 5, 2006: cube51, like cube50 but with bottomDragQuadratic = 0.0012 | 
|  | mkdir cube78 |  | 
|  | cd cube78 |  | 
|  | ln -sf ../run_template . |  | 
|  | ln -sf ../cube78_forcing . |  | 
| 120 |  |  | 
| 121 | ########################################################################## | January 5, 2006: cube50, like cube48 but with snow advection | 
|  | # getting the code from anonymous CVS server |  | 
| 122 |  |  | 
| 123 | bash or sh shell: | December 30, 2006: updates for cube49, the first | 
| 124 | $ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' | "blend" integration. | 
|  | $ cvs login |  | 
|  | ( enter the CVS password: "cvsanon" ) |  | 
|  |  |  | 
|  | tcsh or csh shell: |  | 
|  | $ setenv CVSROOT ':pserver:cvsanon@mitgcm.org:/u/gcmpack' |  | 
|  | $ cvs login |  | 
|  | ( enter the CVS password: "cvsanon" ) |  | 
| 125 |  |  | 
| 126 | cvs co MITgcm_contrib/high_res_cube/README.cs510 | December 11, 2006: preparing for cube47, same as cube43 | 
| 127 | cvs co MITgcm_contrib/high_res_cube/code-mods | but with Arctic monthly river runoff climatology. | 
|  | cvs co MITgcm_contrib/high_res_cube/input |  | 
|  | cvs co MITgcm_code |  | 
| 128 |  |  | 
| 129 |  | December 5, 2006: 216-CPU code for cube46 | 
| 130 |  | Changed diagnostics to accomodate JMCs latest and some | 
| 131 |  | changes to code-mods/CPP_EEOPTIONS.h and W2_OPTIONS.h | 
| 132 |  | to speed up code. | 
| 133 |  |  | 
| 134 | ########################################################################## | November 7, 2006: 216-CPU code for cube43 | 
| 135 | # 510*510*6*50, 216-cpu cube sphere on altix with ifort | Bug fixes in growth.F routine and new set of diagnostics | 
| 136 |  |  | 
| 137 | cd MITgcm | October 15, 2006: 216-CPU code for cube40 | 
|  | mkdir build run |  | 
|  | cd build |  | 
|  | \rm * |  | 
|  | \cp ../../MITgcm_contrib/high_res_cube/code-mods/* . |  | 
|  | \cp ../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* . |  | 
|  | \cp SIZE.h_216 SIZE.h |  | 
|  | ../tools/genmake2 -of ../../MITgcm_contrib/high_res_cube/code-mods/linux_ia64_ifort+mpi_altix_nas |  | 
|  | make depend |  | 
|  | make -j |  | 
|  | cd ../run |  | 
|  | \rm * |  | 
|  | \cp ../../MITgcm_contrib/high_res_cube/input/* . |  | 
|  | ln -sf ../../run_template/* . |  | 
|  | \rm pickup* *192_94* *.m D*1 E* G*all.bin L* W*10.bin h*[6-7] PHC* |  | 
|  | \mv data.exf_hr data.exf |  | 
|  | \mv data_hr data |  | 
|  | ln -sf ../../run_template/pickup.0000000216.cube78 pickup.0000000216 |  | 
|  | ln -sf ../../run_template/pickup_seaice.0000000216.cube78 pickup_seaice.0000000216 |  | 
|  | ln -sf ../../run_template/pickup.0000000216.meta . |  | 
|  | ln -sf ../../run_template/pickup_seaice.0000000216.meta . |  | 
|  | \cp ../build/mitgcmuv . |  | 
|  | \mv job_altix cube78 |  | 
|  | qsub cube78 |  | 
| 138 |  |  | 
| 139 |  | September 28, 2006: 216-CPU code for cube38 | 
| 140 |  |  | 
| 141 | ########################################################################## | September 4, 2006 | 
| 142 | # generating initial pickup files for cube74 | Modified 32*32*6*15 experiments to use grid_cs32.* files from | 
| 143 |  | verification/tutorial_held_suarez_cs/input instead of | 
| 144 |  | tile* files from verification/global_ocean.cs32x15/input | 
| 145 |  |  | 
| 146 | data.cal | August 8, 2006: 216-CPU code and input used for cube37 | 
| 147 | # start on December 30 because of bug that does | using "June 10, 2006" MITgcm code. | 
|  | # not allow a January 1, hour 0 start |  | 
|  | startDate_1=19881230, |  | 
|  |  |  | 
|  | data |  | 
|  | # run for 4 days with 2-min time step |  | 
|  | niter0=0, |  | 
|  | endtime=345600., |  | 
|  | deltaT=120., |  | 
|  | pChkptFreq=345600., |  | 
|  | monitorFreq=1., |  | 
|  |  |  | 
|  | data.seaice |  | 
|  | # set SEAICE_salinity to 33 to avoid initialization schock if |  | 
|  | # ice is created or melts during first time step |  | 
|  | HeffFile           = 'heff_file8', |  | 
|  | SEAICE_salinity    = 33., |  | 
| 148 |  |  | 
| 149 | \mv pickup.0000002880.data ../../run_template/pickup.0000000216.cube74 | for online T/S profile diagnostics all that is needed is | 
| 150 | \mv pickup_seaice.0000002880.data ../../run_template/pickup_seaice.0000000216.cube74 | profiles_init_fixed | 
| 151 |  | profiles_inloop | 
| 152 |  | see pkg/ecco/ecco_cost_init_fixed.F | 
| 153 |  | pkg/ecco/cost_averagesfields.F |