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

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

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

revision 1.14 by dgoldberg, Sat Jun 8 22:15:34 2013 UTC revision 1.16 by dgoldberg, Wed Jun 12 20:48:08 2013 UTC
# Line 61  C     ZERO OUT FLOATING POINT ARRAYS Line 61  C     ZERO OUT FLOATING POINT ARRAYS
61            area_shelf_streamice(i,j,bi,bj) = 0. _d 0            area_shelf_streamice(i,j,bi,bj) = 0. _d 0
62            mass_ice_streamice(i,j,bi,bj) = 0. _d 0            mass_ice_streamice(i,j,bi,bj) = 0. _d 0
63            BDOT_streamice(i,j,bi,bj) = 0. _d 0            BDOT_streamice(i,j,bi,bj) = 0. _d 0
64            ADOT_streamice(i,j,bi,bj) = 0. _d 0            ADOT_streamice(i,j,bi,bj) = streamice_adot_uniform
65            C_basal_friction(i,j,bi,bj) = C_basal_fric_const            C_basal_friction(i,j,bi,bj) = C_basal_fric_const
66            B_glen(i,j,bi,bj) = B_glen_isothermal            B_glen(i,j,bi,bj) = B_glen_isothermal
67            H_streamice_prev(i,j,bi,bj) = 0. _d 0            H_streamice_prev(i,j,bi,bj) = 0. _d 0
# Line 357  C The 0 is the "iteration" argument. The Line 357  C The 0 is the "iteration" argument. The
357                            
358         lenx = sNx*nSx*nPx*delX(1)         lenx = sNx*nSx*nPx*delX(1)
359         leny = sNy*nSy*nPy*delY(1)         leny = sNy*nSy*nPy*delY(1)
        print *, 'lenx', lenx  
        print *, 'leny', leny  
360         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
361          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
362           DO j=1,sNy           DO j=1,sNy
# Line 377  C The 0 is the "iteration" argument. The Line 375  C The 0 is the "iteration" argument. The
375        ELSE IF (STREAMICEbasalTracConfig.EQ.'1DPERIODIC' ) THEN        ELSE IF (STREAMICEbasalTracConfig.EQ.'1DPERIODIC' ) THEN
376                            
377         lenx = sNx*nSx*nPx*delX(1)         lenx = sNx*nSx*nPx*delX(1)
        print *, 'lenx', lenx  
378         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
379          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
380           DO j=1,sNy           DO j=1,sNy
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22