Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Links to checkpoint65w: | (view) (download) (annotate) |
Sticky Tag: |
rename SEAICE_VECTORIZE_LSR_ZEBRA to SEAICE_LSR_ZEBRA and add more comments
- move ALLOW_COST_ICE, ALLOW_SEAICE_COST_SMR_AREA to SEAICE_OPTIONS.h (note : SEAICE_COST.h is tied to the latter)
since lab_sea is failing anyway, I add two CPP flag SEAICE_VECTORIZE_LSR and SEAICE_LSR_ADJOINT_ITER that avoid many recomputations and change the results also for the forward part. The agreement between adjoint and FD gradient consistenty increases with this change and the grdchk summary reduces from 8.3683902422395E-04 to 8.1651767977887E-04 (big deal (o:)
- remove documentation for retired CPP flags (SEAICE_GROWTH_LEGACY, SEAICE_CAP_HEFF).
- switch to new code (#undef SEAICE_GROWTH_LEGACY).
Replace SEAICE_MULTICATEGORY check with "warning" instead of "error+stop". Update SEAICE_OPTIONS.h for main code and "almost-up-to-date" verification experiments.
comment out SEAICE_MULTICATEGORY, so that Dimitris' new check in seaice_check does not stop the experiment
remove MCPHEE_OCEAN_ICE_HEAT_FLUX MCPHEE_OCEAN_ICE_HEAT_FLUX GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR from SEAICE_OPTIONS.h
- SEAICE_OPTIONS.h : define MCPHEE_OCEAN_ICE_HEAT_FLUX and GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR that only serve the purpose of stopping in seaice_check to aknowledge the change of defaults.
- define SEAICE_DISABLE_HEATCONSFIX, until it confirmed by the group that my "fix" is correct, or otherwise.
- update CPPs and run time params (SEAICE_tempFrz0 = -1.96, SEAICE_dTempFrz_dS = 0., to use constant freezing point). - no result chnages.
- rm retired CPPs.
- remove retired CPPs MCPHEE_OCEAN_ICE_HEAT_FLUX, GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR.
- remove undef SEAICE_AGE.
bring OPTIONS files up-to-date
enable free-drift code
switch off SEAICE_ALLOW_FREEDRIFT (getting too many recomputations after today modifications in seaice_dynsolver.F & seaice_lsr.F)
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 : ----------------------------------------------------------------------
Start a seaice tracer-like code but for now confinded to ice age
define SEAICE_GROWTH_LEGACY and SEAICE_SOLVE4TEMP_LEGACY (except in global_ocean.cs32x15 which does not call seaice_growth)
#define SEAICE_ALLOW_FREEDRIFT to compile and run the free drift code in fwd and adj modes. Note: this has no impact on the result, since SEAICEuseFREEDRIFT has not been set true.
Add EVP solver.
o enable SEAICE_MULTICATEGORY for testing o remove obsolete flag
Exclude the cost_ice_test part of the cost function to produce a "cleaner" adjoint
add missing cvs $Header:$ or $Name:$
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
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 |