| 1 | Instructions for setting up a cube-sphere integration with sea-ice. | Instructions for setting up a cube-sphere integration with sea-ice. | 
| 2 |  |  | 
|  | Forcing fields and other input files needed for the CS510 |  | 
|  | configuration can be obtained from |  | 
|  | http://ecco2.org/products/model_setups/cs510.html |  | 
|  |  |  | 
|  | ================= |  | 
|  |  |  | 
|  | UPDATES |  | 
|  |  |  | 
|  | July 24, 2007. Cube72 |  | 
|  | July 8, 2007. Cube66: SEAICE_MULTICATEGORY, SEAICE_VARIABLE_FREEZING_POINT, |  | 
|  | rhonil=1027.5, [temp/salt]AdvScheme=7, convertFW2Salt=-1, LSR_ERROR=1e-6 |  | 
|  | useRealFreshWaterFlux=.TRUE., runoff_lat_inc=179*1, SEAICE_salinity=10 |  | 
|  | bathyFile=GEBCO_510x6x510_ver06_dig.bin, HeffFile=heff_file8 |  | 
|  | SEAICE_availHeatFrac=.5, useGMRedi=.TRUE. |  | 
|  | April 24, 2007. Cube65: Bryan and Lewis 1979 diffusivity |  | 
|  | (cube67, cube68, and cube69 were carried out with cube65 code) |  | 
|  | April 7 , 2007. Cube64: era40_fullres forcing |  | 
|  |  |  | 
|  |  |  | 
| 3 | ########################################################################## | ########################################################################## | 
| 4 | # Get forcing files, etc. | # Get forcing files, etc. | 
| 5 | # Example below is for cube65 on /nobackup2a/menemenl | # Example below is for cube88 on pleiades | 
| 6 |  |  | 
| 7 | cd /nobackup2a/menemenl | mkdir cube94 | 
| 8 | mkdir cube72 | cd cube94 | 
| 9 | cd cube72 | ln -sf /nobackupp10/hzhang1/cs510/run_template . | 
| 10 | ln -sf ../run_template . | ln -sf /nobackupp10/hzhang1/forcing/jra25 . | 
| 11 | ln -sf ../ncep_rgau . | ln -sf /nobackupp10/hzhang1/forcing/cube78_forcing . | 
|  | ln -sf ../era40_ecmwf_blend . |  | 
|  | ln -sf ../cube59_GPCP . |  | 
| 12 |  |  | 
| 13 |  |  | 
| 14 | ########################################################################## | ########################################################################## | 
| 15 | # getting the code from anonymous CVS server | # getting the code from anonymous CVS server | 
| 16 |  |  | 
|  | bash or sh shell: |  | 
|  | $ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' |  | 
|  | $ 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" ) |  | 
|  |  |  | 
| 17 | cvs co MITgcm_contrib/high_res_cube/README.cs510 | cvs co MITgcm_contrib/high_res_cube/README.cs510 | 
| 18 | cvs co MITgcm_contrib/high_res_cube/code-mods | cvs co MITgcm_contrib/high_res_cube/code-mods | 
| 19 | cvs co MITgcm_contrib/high_res_cube/input | cvs co MITgcm_contrib/high_res_cube/input | 
| 21 |  |  | 
| 22 |  |  | 
| 23 | ########################################################################## | ########################################################################## | 
| 24 | # 510*510*6*50, 216-cpu cube sphere on altix with ifort | # 510*510*6*50, 450-cpu cube sphere on pleiades with ifort | 
| 25 |  |  | 
| 26 | cd MITgcm | cd MITgcm | 
| 27 | mkdir build run | mkdir build run | 
| 28 | cd build | cd build | 
| 29 | \rm * | \rm * | 
| 30 | \cp ../../MITgcm_contrib/high_res_cube/code-mods/* . | \cp ../../MITgcm_contrib/high_res_cube/code-mods/* . | 
| 31 | \cp ../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* . | \cp ../../MITgcm_contrib/high_res_cube/code-mods/s450t_34x102/SIZE.h . | 
| 32 | \cp SIZE.h_216 SIZE.h | \rm README SIZE.h_mpi | 
| 33 | ../tools/genmake2 -of ../../MITgcm_contrib/high_res_cube/code-mods/linux_ia64_ifort+mpi_altix_nas | ../tools/genmake2 -of ../tools/build_options/linux_amd64_ifort+mpi_ice_nas | 
| 34 | make depend | make depend | 
| 35 | make -j | make -j 16 | 
| 36 | cd ../run | cd ../run | 
| 37 | \rm * | \rm * | 
| 38 | \cp ../../MITgcm_contrib/high_res_cube/input/* . | \cp ../../MITgcm_contrib/high_res_cube/input/*data* . | 
| 39 | ln -sf ../../run_template/* . | \cp ../../MITgcm_contrib/high_res_cube/input/modpickup . | 
| 40 | \rm pickup* *192_94* *.m D*1 E* G*all.bin L* W* h*[6-7] PHCT* | \cp ../../MITgcm_contrib/high_res_cube/input/job* . | 
| 41 |  | \mv data.seaice_hr data.seaice | 
| 42 | \mv data.exf_hr data.exf | \mv data.exf_hr data.exf | 
| 43 | \mv data_hr data | \mv data_hr data | 
| 44 | ln -sf ../../run_template/pickup.0000000216.cube71 pickup.0000000216 | \ln -sf ../../run_template/tile* . | 
| 45 | ln -sf ../../run_template/pickup_seaice.0000000216.cube66 pickup_seaice.0000000216 | \ln -sf ../../run_template/G*dig* . | 
| 46 |  | \ln -sf ../../run_template/*Init.cube78 . | 
| 47 | \cp ../build/mitgcmuv . | \cp ../build/mitgcmuv . | 
| 48 | \mv job_altix cube72 | qsub job_cs510_pleiades450 | 
| 49 | qsub cube72 |  | 
| 50 |  | ================= | 
| 51 |  |  | 
| 52 |  | UPDATES | 
| 53 |  |  | 
| 54 |  | 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 | 
| 99 |  | and era40_ecmwf_blend/cube59_GPCP forcing | 
| 100 |  | November 3, 2007. Cube 76 - same as cube77 but 1979-2006 integration | 
| 101 |  | August 31, 2007.  Cube75 - the second blend. | 
| 102 |  | August 9, 2007. Cube74 | 
| 103 |  | August 6, 2007. Cube73 | 
| 104 |  | July 24, 2007. Cube72 | 
| 105 |  | July 8, 2007. Cube66: SEAICE_MULTICATEGORY, SEAICE_VARIABLE_FREEZING_POINT, | 
| 106 |  | rhonil=1027.5, [temp/salt]AdvScheme=7, convertFW2Salt=-1, LSR_ERROR=1e-6 | 
| 107 |  | useRealFreshWaterFlux=.TRUE., runoff_lat_inc=179*1, SEAICE_salinity=10 | 
| 108 |  | bathyFile=GEBCO_510x6x510_ver06_dig.bin, HeffFile=heff_file8 | 
| 109 |  | SEAICE_availHeatFrac=.5, useGMRedi=.TRUE. | 
| 110 |  | April 24, 2007. Cube65: Bryan and Lewis 1979 diffusivity | 
| 111 |  | (cube67, cube68, and cube69 were carried out with cube65 code) | 
| 112 |  | 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 |  | removed UVEL and VVEL from diagnostics | 
| 118 |  |  | 
| 119 |  | January 5, 2006: cube51, like cube50 but with bottomDragQuadratic = 0.0012 | 
| 120 |  |  | 
| 121 |  | January 5, 2006: cube50, like cube48 but with snow advection | 
| 122 |  |  | 
| 123 |  | December 30, 2006: updates for cube49, the first | 
| 124 |  | "blend" integration. | 
| 125 |  |  | 
| 126 |  | December 11, 2006: preparing for cube47, same as cube43 | 
| 127 |  | but with Arctic monthly river runoff climatology. | 
| 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 |  | Bug fixes in growth.F routine and new set of diagnostics | 
| 136 |  |  | 
| 137 |  | October 15, 2006: 216-CPU code for cube40 | 
| 138 |  |  | 
| 139 |  | September 28, 2006: 216-CPU code for cube38 | 
| 140 |  |  | 
| 141 |  | September 4, 2006 | 
| 142 |  | 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 |  | August 8, 2006: 216-CPU code and input used for cube37 | 
| 147 |  | using "June 10, 2006" MITgcm code. | 
| 148 |  |  | 
| 149 |  | for online T/S profile diagnostics all that is needed is | 
| 150 |  | profiles_init_fixed | 
| 151 |  | profiles_inloop | 
| 152 |  | see pkg/ecco/ecco_cost_init_fixed.F | 
| 153 |  | pkg/ecco/cost_averagesfields.F |