SALE网格设置*ALE_AMBIENT_HYDROSTATIC关键字，计算时空气中出现了水单元？

SALE网格设置*ALE_AMBIENT_HYDROSTATIC关键字，计算时空气中出现了水单元，如图。

最后水单元下落，水面波动然后就停止计算，这会是因为什么？

求助

 ------------ALE网格设置

 *ALE_STRUCTURED_MESH_CONT ROL_POINTS 

$# cpid     unused     unused sfo unused offo

  10001                                   1.0                0.0 

$# n                   x                                     ratio

      1                0.0                                       0.0 

    20              10.0                                       0.0

*ALE_STRUCTURED_MESH_CONT ROL_POINTS 

$# cpid     unused     unused sfo unused offo 

  10002                                   1.0                0.0 

$#      n               x                                    ratio 

          1            0.0                                       0.0 

          2            0.5                                       0.0

        20          19.5                                       0.0    

        21          20.0                                       0.0 

*ALE_STRUCTURED_MESH_CONT ROL_POINTS 

$# cpid     unused     unused sfo unused offo 

  10003                                   1.0                0.0 

$#      n               x                                    ratio

          1           -5.0                                       0.0

        30             5.0                                       0.0 

----------------远场边界，其中set_solid分别是y轴最大最小区域 

----------------mmgbl 1为空气，2为水 

*ALE_AMBIENT_HYDROSTATIC 

$# sid     stype     vecid       grav       pbase     ramptlc 

        2            2           1  -9.8065  101300.0               0 

$# nid     mmgbl 

   1000               1 

   2000               2 

*SET_SOLID_GENERAL 

$ SID 

      2 

$ OPTION   MSHID   XMN   XMX   YMN   YMX   ZMN   ZMX   

  SALECPT            1          1       20         1        2         1       30 

*ALE_AMBIENT_HYDROSTATIC 

$# sid     stype     vecid       grav       pbase     ramptlc 

        3            2           1  -9.8065  101300.0               0

$# nid mmgbl 

   1001          1 

   2001          2 

*SET_SOLID_GENERAL 

$ SID 

      3 

$ OPTION   MSHID   XMN   XMX   YMN   YMX   ZMN   ZMX  

  SALECPT             1         1       20      20        21        1       30 

*NODE 

$# nid      x        y        z    tc    rc 

   1000  0.1      0.0     5.0    0     0 

   1001  0.1    20.0     5.0    0     0 

   2000  0.2    0.05  3.003   0     0

   2001  0.2  19.95  3.003   0     0