66 |
rv_old_si(i,j,bi,bj) = 0. _d 0 |
rv_old_si(i,j,bi,bj) = 0. _d 0 |
67 |
zu_old_si(i,j,bi,bj) = 0. _d 0 |
zu_old_si(i,j,bi,bj) = 0. _d 0 |
68 |
zv_old_si(i,j,bi,bj) = 0. _d 0 |
zv_old_si(i,j,bi,bj) = 0. _d 0 |
69 |
h_after_uflux_si(i,j,bi,bj) = 0. _d 0 |
! h_after_uflux_si(i,j,bi,bj) = 0. _d 0 |
70 |
#ifdef STREAMICE_HYBRID_STRESS |
#ifdef STREAMICE_HYBRID_STRESS |
71 |
streamice_taubx (i,j,bi,bj) = 0. _d 0 |
streamice_taubx (i,j,bi,bj) = 0. _d 0 |
72 |
streamice_tauby (i,j,bi,bj) = 0. _d 0 |
streamice_tauby (i,j,bi,bj) = 0. _d 0 |
99 |
!#endif |
!#endif |
100 |
! NEW DIRECTIVES - DNG |
! NEW DIRECTIVES - DNG |
101 |
|
|
102 |
|
#ifdef ALLOW_GENTIM2D_CONTROL |
103 |
CALL CTRL_MAP_GENTIM2D (myTime, myIter, myThid) |
CALL CTRL_MAP_GENTIM2D (myTime, myIter, myThid) |
104 |
|
#endif |
105 |
|
|
106 |
if (myIter.eq.0) then |
if (myIter.eq.0) then |
107 |
CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid ) |
CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid ) |
108 |
|
CALL WRITE_FLD_XY_RL |
109 |
|
& ("surf_el_init","",surf_el_streamice,0,myThid) |
110 |
endif |
endif |
111 |
CALL STREAMICE_VELMASK_UPD (myThid) |
CALL STREAMICE_VELMASK_UPD (myThid) |
112 |
|
|
113 |
|
#ifdef ALLOW_STREAMICE_FLUX_CONTROL |
114 |
|
CALL STREAMICE_APPLY_FLUX_CTRL ( myTime, myIter, myThid ) |
115 |
|
#endif |
116 |
|
|
117 |
CALL STREAMICE_VEL_SOLVE( myThid ) |
CALL STREAMICE_VEL_SOLVE( myThid ) |
118 |
|
|
119 |
if(.not.STREAMICE_diagnostic_only) THEN |
if(.not.STREAMICE_diagnostic_only) THEN |
120 |
|
|
121 |
CALL STREAMICE_ADVECT_THICKNESS ( myThid, deltaT ) |
CALL STREAMICE_ADVECT_THICKNESS ( myThid, myIter, deltaT ) |
122 |
|
|
123 |
endif |
endif |
124 |
|
|