/[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.8 by dgoldberg, Wed Oct 10 21:11:46 2012 UTC revision 1.10 by dgoldberg, Thu Mar 7 15:23:19 2013 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4  #include "COST_OPTIONS.h"  #ifdef ALLOW_COST
5    # include "COST_OPTIONS.h"
6    #endif
7  #include "STREAMICE_OPTIONS.h"  #include "STREAMICE_OPTIONS.h"
8    
9    
# Line 27  C     == Global variables === Line 29  C     == Global variables ===
29  # include "STREAMICE.h"  # include "STREAMICE.h"
30  #endif  #endif
31    
32  #include "cost.h"  #ifdef ALLOW_COST
33    # include "cost.h"
34    #endif
35  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
36  # include "tamc.h"  # include "tamc.h"
37  #endif  #endif
# Line 82  C--   Calculate cost function on tile of Line 86  C--   Calculate cost function on tile of
86    
87    
88  !             S_obs(i,j,bi,bj) = 0.0  !             S_obs(i,j,bi,bj) = 0.0
89                 if (streamice_cost_mask(i,j,bi,bj).eq.1.0) THEN
90               cost_func1_streamice (bi,bj) =                cost_func1_streamice (bi,bj) =
91       &       cost_func1_streamice (bi,bj) +       &        cost_func1_streamice (bi,bj) +
92       &       streamice_wgt_surf * (S_obs(i,j,bi,bj)-       &        streamice_wgt_surf * (S_obs(i,j,bi,bj)-
93       &              surf_el_streamice(i,j,bi,bj))**2 / Nx / Ny +       &               surf_el_streamice(i,j,bi,bj))**2 / Nx / Ny +
94       &       streamice_wgt_drift * (H_streamice(i,j,bi,bj)-       &        streamice_wgt_drift * (H_streamice(i,j,bi,bj)-
95       &              H_streamice_prev(i,j,bi,bj))**2 / Nx / Ny       &               H_streamice_prev(i,j,bi,bj))**2 / Nx / Ny
96                 endif
97            
98    
99              end do              end do
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22