/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_cost_surf_accum.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_cost_surf_accum.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by heimbach, Fri Sep 28 02:55:40 2012 UTC revision 1.4 by heimbach, Tue Oct 2 00:43:21 2012 UTC
# Line 28  C     == Global variables === Line 28  C     == Global variables ===
28  #endif  #endif
29    
30  #include "cost.h"  #include "cost.h"
31    #ifdef ALLOW_AUTODIFF_TAMC
32    # include "tamc.h"
33    #endif
34    
35  C     == Routine arguments ==  C     == Routine arguments ==
36  C     myThid - Thread number for this instance of the routine.  C     myThid - Thread number for this instance of the routine.
# Line 61  C     == Local variables Line 64  C     == Local variables
64  C--   Calculate mask for tracer cells  (0 => land, 1 => water)  C--   Calculate mask for tracer cells  (0 => land, 1 => water)
65  !       k=1  !       k=1
66    
67    #ifdef ALLOW_AUTODIFF_TAMC
68    CADJ STORE surf_el_streamice  = comlev1, key = ikey_dynamics,
69    CADJ &     kind = isbyte
70    #endif
71    
72  C--   Calculate cost function on tile of this instance  C--   Calculate cost function on tile of this instance
73        do bj = jtlo,jthi        do bj = jtlo,jthi
74          do bi = itlo,ithi          do bi = itlo,ithi
# Line 70  C--   Calculate cost function on tile of Line 78  C--   Calculate cost function on tile of
78               cost_func1_streamice (bi,bj) =               cost_func1_streamice (bi,bj) =
79       &       cost_func1_streamice (bi,bj) +       &       cost_func1_streamice (bi,bj) +
80       &       0.5 * (S_obs(i,j,bi,bj)-       &       0.5 * (S_obs(i,j,bi,bj)-
81                       surf_el_streamice(i,j,bi,bj))**2       &              surf_el_streamice(i,j,bi,bj))**2
82    
83              end do              end do
84            end do            end do
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22