/[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.6 by dgoldberg, Thu Oct 4 15:40:16 2012 UTC revision 1.8 by dgoldberg, Wed Oct 10 21:11:46 2012 UTC
# Line 57  C     == Local variables Line 57  C     == Local variables
57        ithi = mybxhi(mythid)        ithi = mybxhi(mythid)
58    
59        WRITE(suff,'(I10.10)') myIter        WRITE(suff,'(I10.10)') myIter
60        suff = suff(1:ILNBLNK(suff))//'.data'        suff = suff(1:ILNBLNK(suff))//'.bin'
61    
62        CALL READ_FLD_XY_RS( 'land_ice_surf.'//suff, ' ',        CALL READ_FLD_XY_RS( 'land_ice_surf.'//suff, ' ',
63       &      S_obs, 0, myThid )       &      S_obs, 0, myThid )
# Line 70  CADJ STORE surf_el_streamice  = comlev1, Line 70  CADJ STORE surf_el_streamice  = comlev1,
70  CADJ &     kind = isbyte  CADJ &     kind = isbyte
71  CADJ STORE H_streamice_prev  = comlev1, key = ikey_dynamics,  CADJ STORE H_streamice_prev  = comlev1, key = ikey_dynamics,
72  CADJ &     kind = isbyte  CADJ &     kind = isbyte
73    CADJ STORE H_streamice  = comlev1, key = ikey_dynamics,
74    CADJ &     kind = isbyte
75  #endif  #endif
76    
77  C--   Calculate cost function on tile of this instance  C--   Calculate cost function on tile of this instance
# Line 83  C--   Calculate cost function on tile of Line 85  C--   Calculate cost function on tile of
85    
86               cost_func1_streamice (bi,bj) =               cost_func1_streamice (bi,bj) =
87       &       cost_func1_streamice (bi,bj) +       &       cost_func1_streamice (bi,bj) +
88       &       0.5 * (S_obs(i,j,bi,bj)-       &       streamice_wgt_surf * (S_obs(i,j,bi,bj)-
89       &              surf_el_streamice(i,j,bi,bj))**2 +       &              surf_el_streamice(i,j,bi,bj))**2 / Nx / Ny +
90       &       0.5 * (H_streamice(i,j,bi,bj)-       &       streamice_wgt_drift * (H_streamice(i,j,bi,bj)-
91       &              H_streamice_prev(i,j,bi,bj))**2       &              H_streamice_prev(i,j,bi,bj))**2 / Nx / Ny
92        
93    
94              end do              end do
95            end do            end do
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22