Changeset 178 for trunk/SOURCES/dragging_param_beta_sedim_mod.f90

Timestamp:

01/30/18 11:33:39 (6 years ago)

Author:

aquiquet

Message:

Parametrised beta is the same when using sediments or not (for no sediment condition)

File:

• trunk/SOURCES/dragging_param_beta_sedim_mod.f90 (modified) (1 diff)

trunk/SOURCES/dragging_param_beta_sedim_mod.f90

@@ -209 +209 @@
 
 ! new parametrisation of beta on Neff:
-betamx(:,:)= beta_slope*(neffmx(:,:)**beta_expo)/neffmin
-betamy(:,:)= beta_slope*(neffmy(:,:)**beta_expo)/neffmin
+betamx(:,:)= beta_slope*(neffmx(:,:)**beta_expo) !/neffmin
+betamy(:,:)= beta_slope*(neffmy(:,:)**beta_expo) !/neffmin
 
 where (ilemx(:,:)) betamx(:,:) = betamx(:,:) * coef_ile
 where (ilemy(:,:)) betamy(:,:) = betamy(:,:) * coef_ile
 
-!where (h_sedimmx(:,:).gt.seuil_sedim) betamx(:,:) = betamx(:,:) * coef_sedim
-!where (h_sedimmy(:,:).gt.seuil_sedim) betamy(:,:) = betamy(:,:) * coef_sedim
-where (h_sedimmx(:,:).gt.seuil_sedim) betamx(:,:) = beta_slope*(neffmin**beta_expo)/neffmin
-where (h_sedimmy(:,:).gt.seuil_sedim) betamy(:,:) = beta_slope*(neffmin**beta_expo)/neffmin
+where (h_sedimmx(:,:).gt.seuil_sedim) betamx(:,:) = betamx(:,:) * coef_sedim
+where (h_sedimmy(:,:).gt.seuil_sedim) betamy(:,:) = betamy(:,:) * coef_sedim
+!where (h_sedimmx(:,:).gt.seuil_sedim) betamx(:,:) = beta_slope*(neffmin**beta_expo)/neffmin
+!where (h_sedimmy(:,:).gt.seuil_sedim) betamy(:,:) = beta_slope*(neffmin**beta_expo)/neffmin
 
 betamx(:,:)=max(betamx(:,:),betamin)