Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Links to checkpoint63k: | (view) (download) (annotate) |
Sticky Tag: |
rename parameters: SIsalFrac to SEAICE_saltFrac & SIsal0 to SEAICE_salt0
- update CPPs and run time params (SEAICE_tempFrz0 = -1.96, SEAICE_dTempFrz_dS = 0., to use constant freezing point). - no result chnages.
set useMaykutSatVapPoly=T & postSolvTempIter=0 to reproduce old results with former #defined SEAICE_SOLVE4TEMP_LEGACY code
- remove SEAICE_AGE and add two SITRACER in the foward tests.
- use pkg/exf default emissivities - use default SEAICEuseFluxForm with simple harmonic diffusivity (SEAICEdiffKhArea=200.) instead of DIFF1.
sort out emissivity and long wave radiation step 2: add emissivity for snow and ice as new run-time parameters that default to the corresponding exf values; this changes results in many verification experiments so that for backward compatibility these parameter are now set explicitly to inconsistent values (ugly)
use old default value for SEAICEuseFluxForm (=F) and DIFF1 (=0.004).
Updating of verification experiments now that SEAICE_SALINITY has been changed to SEAICE_VARIABLE_SALINITY and SEAICE_salinity has been changed to SIsalFRAC. Noticed that ICE2WATR was still being used to initialize the HSALT in the lab_sea verification experiments. Removed it from the code but that broke testreport. Hence, some lab_sea/results/*txt files had to be replaced. Now pretty sure that ICE2WATR has been complete neutralized. : Modified Files: : lab_sea/results/output_adm.evp.txt : lab_sea/results/output_adm.noseaicedyn.txt : lab_sea/results/output_adm.txt : lab_sea/input_ad/data.seaice lab_sea/input_ad.evp/data.seaice : lab_sea/input_ad.noseaicedyn/data.seaice : lab_sea/input.salt_plume/data.seaice : lab_sea/code_ad/SEAICE_OPTIONS.h lab_sea/code/SEAICE_OPTIONS.h : 1D_ocean_ice_column/input_ad/data.seaice : 1D_ocean_ice_column/input/data.seaice : 1D_ocean_ice_column/code_ad/SEAICE_OPTIONS.h : 1D_ocean_ice_column/code/SEAICE_OPTIONS.h : seaice_obcs/input/data.seaice : seaice_obcs/code/SEAICE_OPTIONS.h : global_ocean.cs32x15/code_ad/SEAICE_OPTIONS.h : global_ocean.cs32x15/code/SEAICE_OPTIONS.h : : offline_exf_seaice/code_ad/SEAICE_OPTIONS.h : offline_exf_seaice/code/SEAICE_OPTIONS.h : ----------------------------------------------------------------------
Verification experiment update 1) Remove "backward compatability only" lines from lab_sea, offline_exf_seaice, and seaice_obcs data.seaice files < # for backward compatibility only: < SEAICE_cpAir = 1.0039560439560439e+03, < SEAICE_lhSublim = 2.8340219780219775e+06, < SEAICE_rhoAir = 1.3E0, 2) Added rhoConstFresh = 999.8 in input*/data where no rhoConstFresh was previously defined 3) Update output*.txt where necessary Modified Files: 1D_ocean_ice_column/input/data.seaice 1D_ocean_ice_column/results/output.txt lab_sea/input_ad/data lab_sea/input_ad/data.seaice lab_sea/input_ad.evp/data lab_sea/input_ad.evp/data.seaice lab_sea/input_ad.noseaicedyn/data.seaice lab_sea/input_ad_seaice_ocean/data lab_sea/input_ad_seaice_ocean/data.seaice lab_sea/input_ad_seaice_only/data lab_sea/input_ad_seaice_only/data.seaice lab_sea/input/data lab_sea/input/data.seaice lab_sea/input.fd/data.seaice lab_sea/input.hb87/data lab_sea/input.hb87/data.seaice lab_sea/input.salt_plume/data lab_sea/input.salt_plume/data.seaice lab_sea/results/output_adm.evp.txt lab_sea/results/output_adm.noseaicedyn.txt lab_sea/results/output_adm.noseaice.txt lab_sea/results/output_adm.txt lab_sea/results/output.fd.txt lab_sea/results/output.hb87.txt lab_sea/results/output.salt_plume.txt lab_sea/results/output.txt offline_exf_seaice/input_ad/data offline_exf_seaice/input_ad/data.seaice offline_exf_seaice/input/data offline_exf_seaice/input.seaicetd/data offline_exf_seaice/input.seaicetd/data.seaice offline_exf_seaice/results/output_adm.txt offline_exf_seaice/results/output.seaicetd.txt seaice_obcs/input/data seaice_obcs/input/data.seaice seaice_obcs/results/output.txt
Update lab_sea and seaice_obcs verification experiments #1) Removed ICE2WATR from input/ and input_ad/data.seaice #2) Updated results/output.txt and output_adm.txt in those verif. experiments that used the incorrect ICE2WATR values : Modified Files: : lab_sea/input_ad/data.seaice lab_sea/input_ad.evp/data.seaice : lab_sea/input_ad.noseaicedyn/data.seaice : lab_sea/input_ad_seaice_ocean/data.seaice : lab_sea/input_ad_seaice_only/data.seaice : lab_sea/input/data.seaice lab_sea/input.fd/data.seaice : lab_sea/input.hb87/data.seaice : lab_sea/input.salt_plume/data.seaice : lab_sea/results/output_adm.evp.txt : lab_sea/results/output_adm.noseaicedyn.txt : lab_sea/results/output_adm.noseaice.txt : lab_sea/results/output_adm.txt : lab_sea/results/output.salt_plume.txt : seaice_obcs/input/data.seaice : seaice_obcs/results/output.txt : ----------------------------------------------------------------------
in order to reproduce the old results, many data.seaice now contain funny values for cpAir and latent heat of sublimation
set flags to their previous default value (SEAICEadvSalt,Age=F)
Changed ICE2WATR to SEAICE_rhoIce/rhoConst for consistency between dynamics and thermodynamics of pkg/seaice. Set ICE2WATR=0.92 to reproduce old output.txt in verification experiments. http://forge.csail.mit.edu/pipermail/mitgcm-devel/2007-October/003053.html
Change reference output for adjoint yet again (after checking that recent fixes reproduce old reference): o SEAICE_EXCLUDE_FOR_EXACT_AD_TESTING o SEAICE_salinity = 0.3 o SEAICE_availHeatFrac = 0.8
o Something seems wrong with EVP adjoint so revert to LSR adjoint o Add SEAICE_SALINITY adjoint
o Add SEAICEadvSnow, SEAICEuseFlooding to adjoint calculation (changes results) o Add init. AREA, HEFF to control vector
o remove SEAICEuseEVP flag o set SEAICE_EVP = 1.E-8 This fixes the 1.E34 gradient, but 1.E-9 would not fix them Question is whether that epsilon is small enough, and why it is so sensitive; Likely culprit is new code in seaice_evp.F: deltaC = SQRT(MAX(deltaC,SEAICE_EPS_SQ)) deltaZ = SQRT(MAX(deltaZ,SEAICE_EPS_SQ))
o in data, remove some default options o in data.seaice, add EVP runtime flags
Initial check-in of lab_sea adjoint (very similar to natl_box_adjoint/code_seaice) * used C-grid and EVP solver * used KPP, GM/Redi * uses inAdExact=.TRUE. * gradient checks fail for both ADM and TLM due to sensitive snow melt code in seaice_grwoth (look for "JZ") f.d. gradients actually seem more problematic than ADM, TLM gradients; issue is related to forward mismatch between g77 an ifort
This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, select a symbolic revision name using the selection box, or choose 'Use Text Field' and enter a numeric revision.
ViewVC Help | |
Powered by ViewVC 1.1.22 |