--- MITgcm/verification/lab_sea/README 2003/05/22 20:49:14 1.1.2.10 +++ MITgcm/verification/lab_sea/README 2002/12/28 10:11:11 1.4 @@ -24,8 +24,8 @@ vwindFile = 'v10m.labsea79' # 10-m meridional wind atempFile = 'tair.labsea1979' # 2-m air temperature aqhFile = 'qa.labsea1979' # 2-m specific humidity - lwdownFile = 'flo.labsea1979' # downward longwave radiation - swdownFile = 'fsh.labsea1979' # downward shortwave radiation + lwfluxFile = 'flo.labsea1979' # longwave radiation + swfluxFile = 'fsh.labsea1979' # shortwave radiation precipFile = 'prate.labsea1979' # precipitation evapFile = 'evap.labsea1979' # evaporation @@ -84,21 +84,6 @@ make mv mitgcmuv mitgcmuv_exf_forcing -1-CPU executable with SEAICE_EXTERNAL_FORCING, without READ_EVAP - cd ../../../verification/lab_sea/input - ln -sf ../code/SIZE.h . - ln -sf ../code/CPP_EEOPTIONS.h . - ln -sf ../code/CPP_OPTIONS.h . - ln -sf ../code/ECCO_CPPOPTIONS.BULK.h ECCO_CPPOPTIONS.h - ln -sf ../code/SEAICE_OPTIONS.h . - ../../../tools/genmake -makefile - ==> on alhena use: - ../../../tools/genmake -platform=o2k_noopt -makefile - make clean - make depend - make - mv mitgcmuv mitgcmuv_compute_evap - Generating 1-CPU executable with SEAICE_EXTERNAL_FLUXES: cd ../../../verification/lab_sea/input ln -sf ../code/SIZE.h . @@ -142,13 +127,13 @@ To run Experiment 1: cd ../../../verification/lab_sea/input ln -sf DATA.10hours data - ln -sf DATA.SEAICE.lsr data.seaice + ln -sf DATA.SEAICE.adi data.seaice mitgcmuv_1 >&! output.txt mkdir exp1 mv *tave.0000000010.data exp1 Use matlab script lookat_exp1.m to compare the output -of exp1 with that from release1_p16 sea-ice code: +of exp1 with that from release1_p10 sea-ice code: cd ../../../verification/lab_sea/matlab matlab lookat_exp1 @@ -157,13 +142,13 @@ Instructions for running Experiment 2 ===================================== -This is a 1-cpu test of the ADI solver. The solution is -compared to that of experiment 1, which used the LSR solver. +This is a 1-cpu test of the LSR solver. The solution is +compared to that of experiment 1, which used the ADI solver. To run Experiment 2: cd ../../../verification/lab_sea/input ln -sf DATA.10hours data - ln -sf DATA.SEAICE.adi data.seaice + ln -sf DATA.SEAICE.lsr data.seaice mitgcmuv_1 >&! output.txt mkdir exp2 mv *tave.0000000010.data exp2 @@ -195,17 +180,16 @@ mkdir exp3b mv *tave.0000000001.data exp3b -Use matlab script lookat_exp3.m to look at -results from exp3a and exp3b: +Use matlab script lookat_exp3.m to compare +the output of exp2 to that of exp1: cd ../../../verification/lab_sea/matlab matlab lookat_exp3 -Both solvers should and do converge to solutions -that are constant at each latitude. But note that -the two solutions are not exactly the same because -the two solvers represent slightly different -approximations of the viscous/plastic equations. +If the dynamic ice solvers are correct, they should converge to a +constant solution. The tests above show that both the ADI and the LSR +solver are unable to handle periodic domains and therefore that they +cannot be "correctly" parallelized. Instructions for running Experiment 4 @@ -217,7 +201,7 @@ To run Experiment 4: cd ../../../verification/lab_sea/input ln -sf DATA.10hours data - ln -sf DATA.SEAICE.lsr data.seaice + ln -sf DATA.SEAICE.adi data.seaice mpirun -np 2 mitgcmuv_2x1 mkdir exp4 mv *tave.0000000010.data exp4 @@ -228,12 +212,12 @@ matlab lookat_exp4 -The accuracy of the dynamic solvers is increased by -reducing parameter LSR_ERROR for the LSR solver and -by increasing parameter NPSEUDO for the ADI solver. -The computational cost/accuracy trade-off is quite -high for the ADI solver but relatively small for -the LSR solver, which is the default solver. +The inaccuracy at the tile boundary can be decreased by increasing +NPSEUDO in data.seaice at the expense of computation time. +Also as sea-ice ages (thickens) the difference at the tiles becomes +increasingly small. For solver pkg/seaice/adi.F a value of NPSEUDO=10 +appears adequate for forward integrations but cannot be used for sea-ice +adjoint model. Work is underway to fix this problem. Instructions for running Experiment 5 @@ -245,7 +229,7 @@ To run Experiment 5: cd ../../../verification/lab_sea/input ln -sf DATA.2years data - ln -sf DATA.SEAICE.lsr data.seaice + ln -sf DATA.SEAICE.adi data.seaice ln -sf evap.labsea1979 evap.labsea1980 ln -sf flo.labsea1979 flo.labsea1980 ln -sf fsh.labsea1979 fsh.labsea1980 @@ -262,10 +246,10 @@ matlab lookat_exp5 -The comparison here is just for fun. -The configuration is too coarse and artificial -to reproduce realistic sea-ice conditions. - +Disclaimer: +The comparison here is just for fun, not really supposed +to look anything like the data. +Otherwise it would put a lot of people out of business :-) Instructions for running Experiment 6 ===================================== @@ -287,7 +271,6 @@ matlab lookat_exp6 - Instructions for running Experiment 7 ===================================== @@ -322,7 +305,7 @@ To run Experiment 8: cd ../../../verification/lab_sea/input ln -sf DATA.10hours_exf data - ln -sf DATA.SEAICE.exf data.seaice + ln -sf DATA.SEAICE.adi_exf data.seaice mitgcmuv_exf_forcing >&! output.txt mkdir exp8 mv *tave.0000000010.data exp8 @@ -344,7 +327,7 @@ To run Experiment 9: cd ../../../verification/lab_sea/input ln -sf DATA.10hours_exf data - ln -sf DATA.SEAICE.exf data.seaice + ln -sf DATA.SEAICE.adi_exf data.seaice mitgcmuv_exf_fluxes >&! output.txt mkdir exp9 mv *tave.0000000010.data exp9 @@ -354,24 +337,3 @@ cd ../../../verification/lab_sea/matlab matlab lookat_exp9 - - -Instructions for running Experiment 10 -====================================== - -This is a 1-cpu, 10-hour integration similar to exp8, but -with bulk formula computation of evaporation fields. - -To run Experiment 10: - cd ../../../verification/lab_sea/input - ln -sf DATA.10hours_exf data - ln -sf DATA.SEAICE.exf data.seaice - mitgcmuv_compute_evap >&! output.txt - mkdir exp10 - mv *tave.0000000010.data exp10 - -Use matlab script lookat_exp10.m to compare -the output of exp10 to that of exp8: - cd ../../../verification/lab_sea/matlab - matlab - lookat_exp10