Changeset 228


Ignore:
Timestamp:
01/08/19 11:32:34 (5 years ago)
Author:
dumas
Message:

grounding line flux is not applied on a island point

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/SOURCES/furst_schoof_mod.f90

    r224 r228  
    134134 
    135135 
    136 !afq --                 if (.not.(ilemx(i,j)).and..not.(ilemx(i-1,j))) then 
     136                 if (.not.(ilemx(i,j)).and..not.(ilemx(i-1,j))) then 
    137137                     
    138138                 ! extrapolation pour avoir uxbar(i-1,j) qui sera ensuite prescrit dans call rempli_L2 
    … …  
    159159                 !                 print* 
    160160 
    161 !afq --                 end if 
     161                 end if 
    162162               
    163163              else            !  GL a Est du i staggered,                          o centre, x stag 
    … …  
    167167                 !flux        !     in             in           out  G         out           in 
    168168 
    169 !afq --                 if (.not.(ilemx(i,j)).and..not.(ilemx(i+1,j))) then 
     169                 if (.not.(ilemx(i,j)).and..not.(ilemx(i+1,j))) then 
    170170 
    171171                 ! extrapolation pour avoir uxbar(i,j) qui sera ensuite prescrit dans call rempli_L2 
    … …  
    187187                 gr_line_schoof(i+1,j) = 1 
    188188 
    189 !afq --                 end if 
     189                 end if 
    190190 
    191191              end if 
    … …  
    225225                 !flux        !                    in           out        G  out            in 
    226226 
    227 !afq --                 if (.not.(ilemx(i,j)).and..not.(ilemx(i+1,j))) then 
     227                 if (.not.(ilemx(i,j)).and..not.(ilemx(i+1,j))) then 
    228228                     
    229229                 ! extrapolation pour avoir uxbar(i,j) qui sera ensuite prescrit dans call rempli_L2 
    … …  
    245245                 gr_line_schoof(i+1,j) = 1 
    246246 
    247 !afq --                 end if 
     247                 end if 
    248248               
    249249              else   !  GL a west du i staggered,                          o centre, x stag 
    … …  
    253253                 !flux        !                     in           out  G        out            in 
    254254 
    255 !afq --                 if (.not.(ilemx(i+1,j)).and..not.(ilemx(i+2,j))) then 
     255                 if (.not.(ilemx(i+1,j)).and..not.(ilemx(i+2,j))) then 
    256256                     
    257257                 ! extrapolation pour avoir uxbar(i,j) qui sera ensuite prescrit dans call rempli_L2 
    … …  
    273273                 gr_line_schoof(i+2,j) = 1 
    274274 
    275 !afq --                 end if 
     275                 end if 
    276276               
    277277              end if 
    … …  
    317317                 !flux        !     in            out         G out            in 
    318318 
    319 !afq --                 if (.not.(ilemy(i,j)).and..not.(ilemy(i,j-1))) then 
     319                 if (.not.(ilemy(i,j)).and..not.(ilemy(i,j-1))) then 
    320320 
    321321                 ! extrapolation pour avoir uybar(i,j-1) qui sera ensuite prescrit dans call rempli_L2 
    … …  
    338338 
    339339 
    340 !afq --                 end if 
     340                 end if 
    341341               
    342342              else            !  GL au nord du j staggered,                          o centre, x stag 
    … …  
    346346                 !flux        !     in             in           out  G         out           in 
    347347 
    348 !afq --                 if (.not.(ilemy(i,j)).and..not.(ilemy(i,j+1))) then 
     348                 if (.not.(ilemy(i,j)).and..not.(ilemy(i,j+1))) then 
    349349                     
    350350                 ! extrapolation pour avoir uybar(i,j) qui sera ensuite prescrit dans call rempli_L2 
    … …  
    366366                 gr_line_schoof(i,j+1) = 1 
    367367 
    368 !afq --                 end if 
     368                 end if 
    369369               
    370370              end if 
    … …  
    404404                 !flux        !                    in           out        G  out            in 
    405405 
    406 !afq --                 if (.not.(ilemy(i,j)).and..not.(ilemy(i,j+1))) then 
     406                 if (.not.(ilemy(i,j)).and..not.(ilemy(i,j+1))) then 
    407407 
    408408                 ! extrapolation pour avoir uybar(i,j) qui sera ensuite prescrit dans call rempli_L2 
    … …  
    424424                 gr_line_schoof(i,j+1) = 1 
    425425 
    426 !afq --                 end if 
     426                 end if 
    427427               
    428428              else   !  GL au nord du j staggered,                          o centre, x stag 
    … …  
    432432                 !flux        !                     in           out  G        out            in 
    433433 
    434 !afq --                 if (.not.(ilemy(i,j+1)).and..not.(ilemy(i,j+2))) then 
     434                 if (.not.(ilemy(i,j+1)).and..not.(ilemy(i,j+2))) then 
    435435 
    436436                 ! extrapolation pour avoir uybar(i,j) qui sera ensuite prescrit dans call rempli_L2 
    … …  
    452452                 gr_line_schoof(i,j+2) = 1 
    453453 
    454 !afq --                 end if 
     454                 end if 
    455455 
    456456              end if 
Note: See TracChangeset for help on using the changeset viewer.