| 1 | !!!! user defined modules used in cfdXZ code !!!! |
|---|
| 2 | |
|---|
| 3 | !***************************************************************** |
|---|
| 4 | module stress_parameters |
|---|
| 5 | ! container for parameters of stress boundary |
|---|
| 6 | real :: Patch_Radius !! [m] |
|---|
| 7 | real :: Patch_Decay_Scale !! [m] |
|---|
| 8 | real :: delta |
|---|
| 9 | real :: tau_0 !! [kg/(ms2)] |
|---|
| 10 | real :: Shear !! [1/s] |
|---|
| 11 | end module stress_parameters |
|---|
| 12 | !***************************************************************** |
|---|
| 13 | |
|---|
| 14 | |
|---|
| 15 | !***************************************************************** |
|---|
| 16 | module user_parameters |
|---|
| 17 | !******* ice stress problem |
|---|
| 18 | ! x_0,x_1,x_2 positions of ice/water interface [m] |
|---|
| 19 | ! surface_V_of_x surface BC for v, applied in [0,x_1] and |
|---|
| 20 | ! [x_2,Lx] |
|---|
| 21 | ! z_N center position for stratified layer [m] |
|---|
| 22 | ! z_N_scale vertical scale of stratified layer [m] |
|---|
| 23 | ! HV_RATIO nu_horiz/nu where nu is taken as nu_vertical |
|---|
| 24 | ! [1] |
|---|
| 25 | ! (invoked only when HV_METHOD=1) |
|---|
| 26 | ! Smax,Smin max/min Salinity values for ambient profile |
|---|
| 27 | ! [ppt] |
|---|
| 28 | ! Tmax,Tmin max/min Temp values used for ambient profile |
|---|
| 29 | ! [deg C] |
|---|
| 30 | ! Km_of_z eddy viscosity for momentum [m2/s] |
|---|
| 31 | ! Kt_of_z eddy diffusivity for temperature [m2/s] |
|---|
| 32 | ! Kh horizontal eddy diffusivity for both momentum |
|---|
| 33 | ! and temperature when HV_METHOD=2 [m2/s] |
|---|
| 34 | !********************************* |
|---|
| 35 | real :: x_0,x_1,x_2,V_0 |
|---|
| 36 | real, allocatable :: surface_V_of_x(:) |
|---|
| 37 | real, allocatable :: Km_of_z(:),Kt_of_z(:) |
|---|
| 38 | real :: z_N, z_N_scale |
|---|
| 39 | real :: Smax,Smin,Tmax,Tmin,dTdz |
|---|
| 40 | real :: HV_RATIO |
|---|
| 41 | integer :: HV_METHOD |
|---|
| 42 | character(80) :: LAPACK_FLAG |
|---|
| 43 | real :: alpha_rho |
|---|
| 44 | character(len=80) :: pomme_filename,pomme_data_flag,ts_profil |
|---|
| 45 | end module user_parameters |
|---|
| 46 | !***************************************************************** |
|---|
| 47 | |
|---|
| 48 | |
|---|
| 49 | !***************************************************************** |
|---|
| 50 | module immersed_boundary |
|---|
| 51 | integer :: n_j,n_left,n_right |
|---|
| 52 | integer,allocatable :: ibdry_owner(:),local_region(:,:) |
|---|
| 53 | real,allocatable :: ibdry_locations(:,:),ibdry_normals(:,:) |
|---|
| 54 | real,allocatable :: ibdry_data(:,:,:),ibdry_integrals(:,:,:) |
|---|
| 55 | real,allocatable :: ibdry_rhs(:,:), z_limits(:,:) |
|---|
| 56 | real,allocatable :: mask1(:,:), mask2(:,:),mask3(:,:) |
|---|
| 57 | real :: delta_1, delta_2, delta_3 |
|---|
| 58 | real :: xalpha,xbeta,x_left,x_right,z_left,z_right |
|---|
| 59 | end module immersed_boundary |
|---|
| 60 | !***************************************************************** |
|---|
| 61 | |
|---|
| 62 | |
|---|
