- Timestamp:
- 05/19/17 17:44:41 (7 years ago)
- Location:
- trunk/SOURCES
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SOURCES/3D-physique-gen_mod.f90
r102 r109 229 229 integer,dimension(nx,ny) :: FRONTFACEX!< type de front sur les faces x 230 230 integer,dimension(nx,ny) :: FRONTFACEY!< type de front sur les faces y 231 integer,dimension(nx,ny) :: gr_line_schoof ! points ou on impose le flux de schoof (pour sorties) 231 232 integer,dimension(nx,ny) :: ICE !< presence de glace au point considere, 232 233 !< seuil 0 si pose, 1 si flottant -
trunk/SOURCES/furst_schoof_mod.f90
r104 r109 83 83 countx(:,:)=0. 84 84 county(:,:)=0. 85 85 gr_line_schoof(:,:)=0 86 86 87 archimtab(:,:) = Bsoc(:,:)+H(:,:)*alpha_flot -sealevel 87 88 … … 127 128 countx (i-1,j) = countx(i-1,j)+1 128 129 imx_diag (i-1,j) = 1 129 130 gr_line_schoof(i-1,j) = 1 131 130 132 ! extrapolation pour avoir uxbar(i,j) 131 133 … … 136 138 countx (i,j) = countx(i,j)+1 137 139 imx_diag (i,j) = 1 140 gr_line_schoof(i,j) = 1 138 141 ! print*,'schoof uxbar(i,j)=',uxbar(i,j) 139 142 ! print*,'Hgl',Hgl … … 156 159 countx (i,j) = countx(i,j)+1 157 160 imx_diag (i,j) = 1 158 161 gr_line_schoof(i,j) = 1 162 159 163 ! extrapolation pour avoir uxbar(i+1,j) 160 164 … … 164 168 countx (i+1,j) = countx(i+1,j)+1 165 169 imx_diag (i+1,j) = 1 170 gr_line_schoof(i+1,j) = 1 166 171 end if 167 172 … … 193 198 countx (i,j) = countx(i,j)+1 194 199 imx_diag (i,j) = 1 195 200 gr_line_schoof(i,j) = 1 201 196 202 ! extrapolation pour avoir uxbar(i+1,j) 197 203 … … 201 207 countx (i+1,j) = countx(i+1,j)+1 202 208 imx_diag (i+1,j) = 1 203 209 gr_line_schoof(i+1,j) = 1 210 204 211 else ! GL a west du i staggered, o centre, x stag 205 212 … … 216 223 countx (i+1,j) = countx(i+1,j)+1 217 224 imx_diag (i+1,j) = 1 225 gr_line_schoof(i+1,j) = 1 218 226 219 227 ! extrapolation pour avoir uxbar(i+1,j) … … 224 232 countx (i+2,j) = countx(i+2,j)+1 225 233 imx_diag (i+2,j) = 1 234 gr_line_schoof(i+2,j) = 1 226 235 end if 227 236 … … 259 268 county (i,j-1) = county(i,j-1)+1 260 269 imy_diag (i,j-1) = 1 270 gr_line_schoof(i,j-1) = 1 261 271 262 272 ! extrapolation pour avoir uybar(i,j) … … 267 277 county (i,j) = county(i,j)+1 268 278 imy_diag (i,j) = 1 279 gr_line_schoof(i,j) = 1 269 280 270 281 else ! GL au nord du j staggered, o centre, x stag … … 282 293 county (i,j) = county(i,j)+1 283 294 imy_diag(i,j) = 1 295 gr_line_schoof(i,j) = 1 284 296 285 297 ! extrapolation pour avoir uybar(i,j+1) … … 290 302 county (i,j+1) = county(i,j+1)+1 291 303 imy_diag (i,j+1) = 1 304 gr_line_schoof(i,j+1) = 1 292 305 end if 293 306 … … 319 332 county (i,j) = county(i,j)+1 320 333 imy_diag (i,j) = 1 334 gr_line_schoof(i,j) = 1 321 335 322 336 ! extrapolation pour avoir uybar(i,j+1) … … 327 341 county (i,j+1) = county(i,j+1)+1 328 342 imy_diag (i,j+1) = 1 343 gr_line_schoof(i,j+1) = 1 329 344 330 345 else ! GL au nord du j staggered, o centre, x stag … … 342 357 county (i,j+1) = county(i,j+1)+1 343 358 imy_diag (i,j+1) = 1 344 359 gr_line_schoof(i,j+1) = 1 360 345 361 ! extrapolation pour avoir uybar(i,j+1) 346 362 … … 350 366 county (i,j+2) = county(i,j+2)+1 351 367 imy_diag (i,j+2) = 1 368 gr_line_schoof(i,j+2) = 1 369 352 370 end if 353 371
Note: See TracChangeset
for help on using the changeset viewer.