Warning, /verification/so_box_biogeo/README is written in an unsupported language. File is not indexed.
view on githubraw file Latest commit 698b072c on 2014-08-27 21:38:16 UTC
698b072c19 Jean*0001 Southern-Ocean box with Biochemistry, using Open-Boundary Conditions
0002 (pkg/obcs) at Northern, Eastern and Western edges of the domain.
0003 ======================================================================
0004
0005 This experiment illustrates and tests the use of package DIC with OBCS.
0006
0007 The configuration (e.g., resolution), model parameters and forcing are
0008 almost identical to tutorial_global_oce_biogeo expect that the horizontal
0009 domain is limited to a sub-domain around Drake passage with open-boundary
0010 conditions coming from the last year of a 2 yrs global simulation and
0011 initial conditions taken from t=1.yr of this same simulation.
0012 This enable to compare directly the results of the first year simulation
0013 of this regional set-up with the 2nd year results of the global set-up run.
0014
0015 Configure and compile the code:
0016 cd build
0017 ../../../tools/genmake2 -mods ../code [-of my_platform_optionFile]
0018 make depend
0019 make
0020 cd ..
0021
0022 To run:
0023 cd run
0024 ln -s ../input/* .
0025 ln -s ../build/mitgcmuv .
0026 ./mitgcmuv > output.txt
0027 cd ..
0028
0029 There is comparison output in the directory:
0030 results/output.txt
0031
0032 ----------------------------------------------------------------------
0033 To generate inital and open boundary conditions:
0034 a) The global set-up (using executable from: tutorial_global_oce_biogeo/build)
0035 was run for 2 years using model-parameter from inp_global/.
0036 The only differences with the ones used in tutorial_global_oce_biogeo are:
0037 - CD-Scheme is turned off (useCDscheme=F), since it is not implemented for OBCS ;
0038 - as a consequence, horizontal viscosity (viscAh) is increased from 2.E5 to 3.E5;
0039 - convective-adjustment diffusivity (ivdc_kappa) is reduced from 100. to 10.m^2/s;
0040 - implicit vertical viscosity is turned off (not needed);
0041 - free-surface exact conservation is turned on (exactConserv=T);
0042 - use the Avective-Form of GM (in data.gmredi: GM_AdvForm=T) instead of Skew-Flux;
0043
0044 to repeat this global simutlation, from so_box_biogeo/
0045 mkdir run_glob ; cd run_glob
0046 ln -s ../inp_global/* .
0047 ../inp_global/prepare_run
0048 ln -s ../../tutorial_global_oce_biogeo/build/mitgcmuv .
0049 ./mitgcmuv > output.glob
0050 cd ..
0051
0052 b) Then regional model forcing files, initial conditions and open-boundary
0053 conditions files are generated using the matlab script: "mk_box_input.m"
0054 cd inp_global/
0055 matlab
0056 >> mk_box_input