Changeset 101


Ignore:
Timestamp:
06/06/13 17:51:29 (11 years ago)
Author:
jbrlod
Message:

add filter (with a bug)

Location:
altifloat
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • altifloat/obs_float/uoptim1.dat

    r99 r101  
    5375371 10 15   0.000000000000000e+00 
    5385381 10 16   0.000000000000000e+00 
    539 1 10 17   3.395363688468933e-01 
    540 1 10 18   8.488409221172333e-02 
     5391 10 17   0.000000000000000e+00 
     5401 10 18   0.000000000000000e+00 
    5415411 10 19   0.000000000000000e+00 
    5425421 10 20   0.000000000000000e+00 
    … …  
    5955951 11 15   0.000000000000000e+00 
    5965961 11 16   0.000000000000000e+00 
    597 1 11 17   1.455155909061432e-01 
    598 1 11 18   3.637889772653580e-02 
     5971 11 17   0.000000000000000e+00 
     5981 11 18   0.000000000000000e+00 
    5995991 11 19   0.000000000000000e+00 
    6006001 11 20   0.000000000000000e+00 

  • altifloat/obs_float/uoptim2.dat

    r100 r101  
    5375371 10 15   0.000000000000000e+00 
    5385381 10 16   0.000000000000000e+00 
    539 1 10 17   3.881162367761135e-01 
    540 1 10 18   1.005681008100510e-01 
     5391 10 17   0.000000000000000e+00 
     5401 10 18   0.000000000000000e+00 
    5415411 10 19   0.000000000000000e+00 
    5425421 10 20   0.000000000000000e+00 
    … …  
    5955951 11 15   0.000000000000000e+00 
    5965961 11 16   0.000000000000000e+00 
    597 1 11 17   4.893954098224640e-01 
    598 1 11 18   2.525196895003319e-01 
     5971 11 17   0.000000000000000e+00 
     5981 11 18   0.000000000000000e+00 
    5995991 11 19   0.000000000000000e+00 
    6006001 11 20   0.000000000000000e+00 
    … …  
    6536531 12 15   0.000000000000000e+00 
    6546541 12 16   0.000000000000000e+00 
    655 1 12 17   1.327397972345352e-01 
    656 1 12 18   1.024922579526901e-01 
     6551 12 17   0.000000000000000e+00 
     6561 12 18   0.000000000000000e+00 
    6576571 12 19   0.000000000000000e+00 
    6586581 12 20   0.000000000000000e+00 

  • altifloat/obs_float/voptim1.dat

    r99 r101  
    5375371 10 15   0.000000000000000e+00 
    5385381 10 16   0.000000000000000e+00 
    539 1 10 17   7.136353105306625e-02 
    540 1 10 18   1.784088276326656e-02 
     5391 10 17   0.000000000000000e+00 
     5401 10 18   0.000000000000000e+00 
    5415411 10 19   0.000000000000000e+00 
    5425421 10 20   0.000000000000000e+00 
    … …  
    5955951 11 15   0.000000000000000e+00 
    5965961 11 16   0.000000000000000e+00 
    597 1 11 17   3.058437071740627e-02 
    598 1 11 18   7.646092679351568e-03 
     5971 11 17   0.000000000000000e+00 
     5981 11 18   0.000000000000000e+00 
    5995991 11 19   0.000000000000000e+00 
    6006001 11 20   0.000000000000000e+00 

  • altifloat/obs_float/voptim2.dat

    r100 r101  
    5375371 10 15   0.000000000000000e+00 
    5385381 10 16   0.000000000000000e+00 
    539 1 10 17   7.240399427246302e-02 
    540 1 10 18   1.842492126161233e-02 
     5391 10 17   0.000000000000000e+00 
     5401 10 18   0.000000000000000e+00 
    5415411 10 19   0.000000000000000e+00 
    5425421 10 20   0.000000000000000e+00 
    … …  
    5955951 11 15   0.000000000000000e+00 
    5965961 11 16   0.000000000000000e+00 
    597 1 11 17   9.550390206277370e-02 
    598 1 11 18   5.933662271127105e-02 
     5971 11 17   0.000000000000000e+00 
     5981 11 18   0.000000000000000e+00 
    5995991 11 19   0.000000000000000e+00 
    6006001 11 20   0.000000000000000e+00 
    … …  
    6536531 12 15   0.000000000000000e+00 
    6546541 12 16   0.000000000000000e+00 
    655 1 12 17   4.170881211757660e-02 
    656 1 12 18   3.267061710357666e-02 
     6551 12 17   0.000000000000000e+00 
     6561 12 18   0.000000000000000e+00 
    6576571 12 19   0.000000000000000e+00 
    6586581 12 20   0.000000000000000e+00 

  • altifloat/src/floater.h

    r100 r101  
    77YREAL vr_dy=0; 
    88 
    9  
     9//FILTER COEFFICIENTS 
     10#ifdef K_FILTER  
     11YREAL c1=0.5; 
     12YREAL c2,c3,c4,c5; 
     13#endif 
    1014 
    1115//Options du run incrŽmental 
    … …  
    1721//Declaration  
    1822 
    19 void appli_start(int argc, char *argv[]){} 
     23void appli_start(int argc, char *argv[]){ 
     24#ifdef K_FILTER 
     25c2=c3=c4=c5=0.125; 
     26#endif 
     27} 
    2028void before_it(int nit){} 
    2129void cost_function(int pdt){} 
    … …  
    6068  for (j=0;j<nlon;j++) 
    6169        for (k=0;k<nlat;k++) { 
    62                 YS1_u_d(j,k)=0; 
    63                 YS1_v_d(j,k)=0; 
     70#ifdef K_FILTER 
     71          YS1_uc_d(j,k)=0; 
     72          YS1_vc_d(j,k)=0; 
     73#else      
     74          YS1_u_d(j,k)=0; 
     75          YS1_v_d(j,k)=0; 
     76#endif 
    6477} 
    6578} 

  • altifloat/src/floater_delta.d

    r100 r101  
    1 //#define FILTER 
     1#define FILTER 
    22 
    33defval jtlag 12//nombre de pas de temps Delta t/ delta t 
    … …  
    7373modul lapfilv space  S_euler_df_t clonol lapfilu 
    7474modul u_d      space S_euler_d  input 1 output 1 
    75 modul v_d      space S_euler_d  input 1 output 1 
     75modul v_d      space S_euler_d  clonol u_d 
    7676 
    7777#else 
    … …  
    126126ctin lapfilv    6      from  lapfilv  1  i   j+1 t-1 
    127127 
    128 ctin uc_d       1      from  lapfilu   1  i   j   t 
    129 ctin vc_d       1      from  lapfilv   1  i   j   t 
     128ctin u_d       1      from  lapfilu   1  i   j   t 
     129ctin v_d       1      from  lapfilv   1  i   j   t 
    130130#endif 
    131131 
    … …  
    172172#ifdef FILTER 
    173173order spaceintraj Tfil 
    174         S_euler_df 
     174        S_euler_df_t 
    175175forder 
    176176#endif 
Note: See TracChangeset for help on using the changeset viewer.