ATN 15.Sep.2019

ASTE release 1:

code:
MITgcm: checkpoint c65q
code_ASTE_R1
  - include logdifkr code
  - include sshv4 that allows for monthly inputs
  - include sea ice thickness and velocity costs (bug fix included)
  - update bpv4 with run-time input for smoothing length-scale
  - update to latest read/write subroutines

  -- NOT yet working for cost_gencost_boxmean.
     (At the time I got the cost_boxmean to work, the code
      changed drastically with many fields being renamed/moved...
      In my attempt to "update" the code to bring it into
      consistency with what Gael had changed, there's still
      a mistake somewhere, so this cost is not working.)

  -- NOT yet working for any costs related to transports
     across defined transects, for the exact same reason above.

  If there's desire to have these two costs working, I have "older"
  version that works (results for my Nguyen et al. 2017 and up-coming
  Bering Strait paper).

Lastly, I've made several attempts to update the code to more
present versions, however, not only does the adjoint run significantly
slower but after some checkpoint the adjoint is no longer stable and
ASTE would consistently crash.  As a result, I've stayed with ckpt65q
and only implement updates where I needed.

=======================================
In adjoint mode: 

code_ASTE_R1/SIZE.h
OL[x,y]=4

build_ASTE_R1:
-mods ../code_ASTE_R1

NAMELISTS_ASTE_R1:
cp data_lfs data
 tempAdvScheme=30,
 saltAdvScheme=30,

 nonlinFreeSurf=0,
 useRealFreshWaterFlux=.FALSE.,
 linFSconserveTr=.FALSE.,

# 2002.01.01 to 2017.12.31 + 1 day: 366*4+365*12=5844days*24*3600/1200 + 1ts
 nTimeSteps=420769,
 deltaT = 1200.,

=======================================

In forward mode, using non linear free surface, real freshwater flux, advection scheme 7:

code_ASTE_R1/SIZE.h
OL[x,y]=8

build_ASTE_R1:
'-mods ../code_ASTE_R1 ../code_adv7_bypass_tamc'

NAMELISTS_ASTE_R1:
cp data_nlfs_adv7 data
 tempAdvScheme=7,
 saltAdvScheme=7,

 select_rStar=2,
 nonlinFreeSurf=4,
 useRealFreshWaterFlux=.TRUE.,
 linFSconserveTr=.FALSE.,

# 2002.01.01 to 2017.12.31 + 1 day: 366*4+365*12=5844days*24*3600/600 + 1ts
 nTimeSteps=841537,
 deltaT = 600.,

========================================

Model binary inputs can be found in:
input_binaries
input_weight_ASTER1.tgz

Raw JRA55 atmospheric forcing:
https://ecco.jpl.nasa.gov/drive/files/ECCO2/atmos/jra55


========================================
To enable momentum diagnostic output for closed budget analyses, use the follow additional code and recompile: 

'-mods ../../code_ASTE_R1 ../code_adv7_bypass_tamc ../code_ASTE_R1_mombudg'

Then copy NAMELISTS_MOMBUDG/data.diagnostics to your namelists directory for only momentum budget diagnostics, or merge that data.diagnostics with NAMELISTS_ASTE_R1/data.diagnostics to also output state/trsp variables in addition to momentum budget terms.