Changeset 160 for altifloat/exp_twin_aviso

Timestamp:

07/20/15 11:51:06 (9 years ago)

Author:

jbrlod

Message:

finalize wind effect

Location:

altifloat/exp_twin_aviso

Files:

• assim.i (modified) (1 diff)
• script_prepare_exp_twin.m (modified) (1 diff)

altifloat/exp_twin_aviso/assim.i

@@ -24 +24 @@
 loadstate u_norm 0 ij 0 A 1 ./NU.dat D
 loadstate v_norm 0 ij 0 A 1 ./NV.dat D
+
+#wind coefficient
+#Poulain, JAOT, 2009
+setstate au 0
+setstate av 0
+#setstate bu 0.0061
+#setstate bv -0.0031
+setstate bu 0
+setstate bv 0
+load_wind wind_interp.nc
 
 #Observations :

altifloat/exp_twin_aviso/script_prepare_exp_twin.m

@@ -31 +31 @@
 %% STEP THREE : Prepare filter
 make_filter(meshfile,outmask,20000,rootdir);
+
+%% STEP THREE BIS : Prepare wind
+windfile = [rootdir 'wind_interp.nc'];
+infile = '../../data/wind/wind_20130825-20130930.nc';
+interp_wind(infile,meshfile,windfile);
+
 
 %% STEP FOUR: Prepare init