50 |
C rSurftmp :: free surface r-position that is used to compute hFac_surf |
C rSurftmp :: free surface r-position that is used to compute hFac_surf |
51 |
C adjust_nb_pt :: Nb of grid points where rSurf is adjusted (hFactInf) |
C adjust_nb_pt :: Nb of grid points where rSurf is adjusted (hFactInf) |
52 |
C adjust_volum :: adjustment effect on the volume (domain size) |
C adjust_volum :: adjustment effect on the volume (domain size) |
53 |
INTEGER i,j,k,bi,bj |
INTEGER i,j,bi,bj |
54 |
INTEGER ks |
INTEGER ks |
55 |
_RL hFacInfMOM, Rmin_tmp, hFactmp, adjust_nb_pt, adjust_volum |
_RL hFacInfMOM, Rmin_tmp, hFactmp, adjust_nb_pt, adjust_volum |
56 |
_RL rSurftmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL rSurftmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
73 |
hFac_surfC(i,j,bi,bj) = 0. |
hFac_surfC(i,j,bi,bj) = 0. |
74 |
hFac_surfW(i,j,bi,bj) = 0. |
hFac_surfW(i,j,bi,bj) = 0. |
75 |
hFac_surfS(i,j,bi,bj) = 0. |
hFac_surfS(i,j,bi,bj) = 0. |
76 |
|
PmEpR(i,j,bi,bj) = 0. |
77 |
Rmin_surf(i,j,bi,bj) = Ro_surf(i,j,bi,bj) |
Rmin_surf(i,j,bi,bj) = Ro_surf(i,j,bi,bj) |
78 |
ENDDO |
ENDDO |
79 |
ENDDO |
ENDDO |