Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Sticky Tag: |
Correct an indexing bug in cost_gencost_sshv4.F.
- replace gencost_scalefile, gencost_smooth2Ddiffnbt (now deprecated) with gencost_posproc_c, gencost_posproc_i (gencost_posproc = smooth)
add 2 warnings related to hardcoded (1) num0 and (2) ndaysave to let users know (1) why some fields are set to zero and (2) why entire cost sections are skipped even after passing all tests post ecco_check.
- revised handling of gencost_errfile and gencost_weight.
- remove old prof_etan_mean code
- add includes of PROFILES_SIZE.h
- cost_gencost_bpv4.F, cost_gencost_sshv4.F, ecco_check.F : bring up to date after updated variable naming convention. - ecco_readparms.F : set gencost_is3d automiatically, set mult_gencost default to 1.
- use eccoiter in place of optimcycle. - if autodiff is not compiled then use READ_REC_XY_RL/READ_REC_XYZ_RL instead of active read/write
- cost_gencost_sshv4.F : set gmsl weight to 1. - ecco.h : increase NGENCOST3D to 4. - cost_gencost_bpv4.F : use gencost_outputlevel.
- cost_gencost_boxmean.F, cost_gencost_seaicev4.F : fix gencost for case of undef SST costs (... see devel posting later). - cost_gencost_sshv4.F : remove remaining CPP brackets (ALLOW_GENCOST_SSHV4_OUTPUT, ALLOW_SSH_MEAN_COST_CONTRIBUTION, ALLOW_PSBAR_MEAN) and rely on variables instead (added useEtaMean; used gencost_outputlevel) .
- cost_gencost_boxmean.F, cost_gencost_bpv4.F, ecco_cost_init_fixed.F, ecco_cost_weights.F : treat cases of missing files at run time. - cost_gencost_sshv4.F : rename etabar as etaday to reflect the required averaging frequency, move the test for etaday (now used instead of psbar) definition to ecco_check.F - ecco_check.F : if no gencost term properly defines etaday (e.g. due to missing files) then turn off sshv4 cost all together.
- avoid AD recomputation.
- cost_gencost_customize.F : add automated fill based upon barfile name when it starts with eta, sst, sss, bp, taux, or tauy; rm hard coded gencost_modfld(i,j,bi,bj,1) = empmr(i,j,bi,bj) - cost_gencost_assignperiod.F : revise gencost conditionality - cost_gencost_sshv4.F : rm ALLOW_SSH_COST_CONTRIBUTION etc brackets, add initialization of local arrays, streamline reading of mdt, use gencost_barfile rather than psbarfile, use local variables (e.g. tpob, tpma instead of tpobs, tpmask), only include ecco.h - cost_gencost_sstv4.F : rename igen_amsre as kgen, replace sstbar with gencost_barfile, rm ALLOW_DAILYSST_COST_CONTRIBUTION - cost_gencost_boxmean.F : replace tbarfile with gencost_barfile, rename igen_boxmean as kgen - cost_gencost_sstv4.F, cost_gencost_sshv4.F : rm myiter, mytime args - cost_gencost_bpv4.F (new) : to replace cost_bp using gencost, added to ecco_ad_diff.list
- remove ALLOW_HIGHLAT_ALTIMETRY codes (now in IO routines) - impose minium area for global mean constraint ( slaoffset_sum.GT.0.25 _d 0 * globalArea ) )
- implement slaoffset constraint (gencost_name is 'sshv4-gmsl')
subtract anom_slaobs from pointwise - merge irec loops of PARTS 4 & 5 - compute pointwise at krec = irec + (ndaysave/2) - subtract slaoffset and anom_psMslaobs
subtract time variable global mean in lsc term - remove anom_psMtpobs in PART 4 (is not used in lsc!) - introduce slaoffset then compute it then subtract it from anom_psMslaobs (and anom_slaobs)
revise lsc and pointwise terms wrt time mean removal - rename mean_slaobs2 as mean_slaos_model and mean_slaobs as mean_slaobs_mdt - introduce mean_psMslaobs (and MSK) - use it instead of mean_psMtpobs etc (and MSK) in the computation of lsc and pointwise terms
- fix bug in gencost indexing.
Replace hardcoded mdt year range with runtime param mdt[start,end]date[1,2], rename variable "tpmean" as "mdt" and "topexmeanfile" as "mdtdatfile", add error+stop statement if "topexmeanfile" is used in data.ecco
o move smooth2Ddiffnbt param from pkg/smooth to pkg/ecco o minor bug fix in set gencost default values in ecco_readparms.F
Replace hardcoded value 300 for eccov4 with general runtime param smooth2Ddiffnbt
- refine large scale anomalies cost computation : subsample smoothed field to observed points (at medium day), remove 1/ndaysaveRL factors, remove weekly stepping, output sladiff_subsample
- refine MDT misfit computation : more consistent offset computations, clean up of "doReference or not" logic, mask out mdt points where less than 100 altimeter data points.
- fix bug, reported by O. Wang, in MDT re-referencing, which was always off.
- add compile switches ALLOW_SHALLOW_ALTIMETRY and ALLOW_HIGHLAT_ALTIMETRY.
- introduce gencost_flag(k) that is currently one of 3 values : 0 if cost term k is not used 1 if cost term is fully interactive (i.e. uses gencost_datafile) -1 if ocst term is pre-specified (i.e. uses gencost_name) - introduce using_gencost(k) logical in line with using_topex etc. - ecco_readparms.F : assign using_gencost(k) & gencost_flag(k) based on gencost_datafile(k).NE.' ' (fully interactive cost terms) or on gencost_name(k) (for pre-specified generic cost terms) - throughout : replace "if ( gencost_datafile(num_var) .ne. ' ' )" with "if ( using_gencost(num_var) ) then" - throughout : rm ALLOW_GENCOST_FREEFORM CPP option. The same thing is now achieved dynamically based on gencost_flag. - ecco_check.F : add gencost checks, account for "data"startdate in ecco_check_files, add topexmean check.
- remove #ifdef ALLOW_SMOOTH_CORREL3D brackets. - add more relevant #ifdef ALLOW_SMOOTH ones. - sort out useAtmWind, useSMOOTH, ctrlSmoothCorrel2D.
include ECCO_OPTIONS.h instead of COST_CPPOPTIONS.h
Clean up inclusions of CTRL_SIZE.h
- add run-time switch usePROFILES.
- goals: -- fix startdate.NE.modelstartdate cases and related model crashes. -- alleviate the need for altimeter etc. files (of zeros) for e.g. 1948-1991. - accross the board issues that should be fixed in the updated routines : -- dont try to read a file that does no exist -- take proper account of startdate (e.g. ymod .GE. yday) -- use cal_FullDate to complete date -- if .le. spval (not .lt. spval) then mask out - more specific issues and updates : -- cost_gencost_ssht4.F : hard-code startdate to 19920101 (no run time param avail). -- cost_readscatxfields.F & cost_gencost_sstv4.F : added option for debbuging (if ALLOW_ECCO_DEBUG is defined then to print date info to stdout) -- cost_sla_read_yd.F : added startdate as input parameter -- ecco_summary.F : bug fix; include AD_CONFIG.h was missing.
- dont create ad files unless adjoint run. - bug fix : sla_startdate etc. rather than modelstartdate. - ecco_cost_driver.F : dont do cost_ssh.F when cost_gencost_sshv4.F. - cost_gencost_sshv4.F : only re-reference MDT if altimeter period is covered.
- ifndef ALLOW_PSBAR_MEAN remove sea level mean from psbar - mask out data beyond abs(lat)>66
- cost_gencost_sshv4.F and cost_sshv4.F: compute smoothed 35 day-average fields once a week only -- for speed up. - ecco_cost_driver.F: necessary includes.
Pointwise vs Large-Scale cost function terms for MDT and bin averaged daily SLA (cost_gencost_sshv4.F) and bin averaged daily SST (cost_gencost_sstv4.F) within the gencost framework. note: cost_gencost_sshv4.F v1.1 is the same as cost_sshv4.F v1.5 except that it takes advantage of the gencost framework.
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 |