Changeset 59


Ignore:
Timestamp:
07/25/11 11:17:12 (13 years ago)
Author:
berrada
Message:

filter and moduls connected at offtime

Location:
branches/branche-mb
Files:
4 added
1 deleted
6 edited

Legend:

Unmodified
Added
Removed

  • branches/branche-mb/scripts/sinobad.i

    r58 r59  
    4848xistate_init 1 ../data_in/file_rest/GYRE_00000400_restart.nc 
    4949load_eb ta_c  // charger l'ebauche ' 
     50setstate nu 0 
     51outoebx nu 1 0 
    5052set_modeltime 0 
    5153print_time OFF 

  • branches/branche-mb/src/cost.dat

    r55 r59  
    22 0.0000000000000000e+00 
    33 0.0000000000000000e+00 
    4  2.9733593674719423e+03 
    5  5.4561410776679122e+02 
    6  2.4475352252509134e+02 
    7  2.6358555318068841e+02 
    8  2.1527711430855038e+02 
    9  2.0204126463990215e+02 
    10  2.1416662011590833e+02 
    11  2.0207208737645496e+02 
    12  2.0211049147089304e+02 
    13  2.0235703823123197e+02 
    14  2.0019969413511828e+02 
    15  2.2286623273509434e+02 
    16  1.9878349444274735e+02 
    17  2.0070292696001030e+02 
    18  2.0014258993979595e+02 
    19  1.9972801663599338e+02 
    20  1.9293537153893590e+02 
    21  2.0814944891652880e+02 
    22  1.8213155756443817e+02 
    23  1.7950672225933866e+02 
    24  1.8080233872884580e+02 
    25  1.7876027971923725e+02 
    26  2.0259497236690703e+02 
    27  1.7881103040469617e+02 
    28  1.7886996041817454e+02 
    29  1.7907554124113310e+02 
    30  1.8045371934564810e+02 
    31  1.7569372841627006e+02 
    32  1.9113648527038990e+02 
    33  1.7516955251891332e+02 
    34  1.9550146370923514e+02 
    35  1.7455342251319595e+02 
    36  2.1077987977826763e+02 
    37  1.7473743190883036e+02 
    38  1.7476659805882767e+02 
    39  1.7491607437730161e+02 
    40  1.7534570209262407e+02 
    41  1.7146968928574532e+02 
    42  1.7918664114219453e+02 
    43  1.7311809756377912e+02 
    44  1.7320818032615955e+02 
    45  1.7322190117805096e+02 
    46  1.7328346667692077e+02 
    47  1.7390681350015083e+02 
    48  1.7091883954920155e+02 
    49  1.7074257389516592e+02 
    50  1.7423942813610816e+02 
    51  1.7069031205380722e+02 
    52  1.7425456057653693e+02 
    53  1.7425631257707835e+02 
    54  1.7425753342598696e+02 

  • branches/branche-mb/src/sinobad.d

    r58 r59  
    1818defval TU       2 | start time step 
    1919#ifdef FILTER 
    20 defval OFT      20 
     20defval OFTL      40 
     21defval OFTZ      40 
     22defval OFT      80 
    2123#endif 
    2224#ifndef FILTER 
    … …  
    5052//#TRAJECTOIRES (DIM) -----------------------------------TRAJ 
    5153#ifdef FILTER 
    52 traj  Tfilh      M   1  0   1  10 
    53 traj  Tfilz      M   1  10   1  OFT 
     54traj  Tcstf      M   0  1 
     55traj  Tlfil      M   1  0    1  OFTL 
     56traj  Tzfil      M   1  OFTL 1  OFTZ 
    5457#endif 
    5558traj  Tcst       M   0   OFT  1 1 
    … …  
    5760exec disp_traj 
    5861 
     62#ifdef FILTER 
     63space  S3df      M    NX NY NZ    Tcstf 
     64space  S3dtlf    M    NX NY NZ    Tlfil 
     65space  S3dtzf    M    NX NY NZ    Tzfil 
     66#endif 
    5967 
    6068space  S0d       M    1           Tcst 
    … …  
    7886 
    7987//#modul         ta_c                   space   S3d     input array NPCA                  output 1       
    80  
    8188#ifndef FILTER 
    8289modul    ta_c                   space   S3d                     noward            output 1      target 
    … …  
    8895 
    8996#ifdef FILTER 
    90 modul    nu                     space   S3d                     noward            output 1      target 
    91 modul    filtre                 space   S0d                     input 0           output 1 
    92 modul    ta_c                   space   S3d                     input 0           output 1 
     97modul    nu                     space   S3df                    noward            output 1      target 
     98modul    dta_lfil               space   S3dtlf                  input 6           output 1      tempo 
     99modul    ztuv_lfil              space   S3dtlf                  input 3           output 2      tempo 
     100modul    dta_zfil               space   S3dtzf                  input 4           output 1      tempo 
     101modul    zwy_zfil               space   S3dtzf                  input 2           output 1      tempo 
     102modul    ta_c                   space   S3d                     input 1           output 1       
    93103#endif 
    94104 
    … …  
    269279#endif 
    270280 
    271 //exec disp_modul 
     281exec disp_modul 
    272282//############################################################################################# 
    273283 
     284#ifdef FILTER 
     285//#---------->dta_lfil 
     286ctin     dta_lfil       1       from    nu                      1  i    j    k   
     287ctin     dta_lfil       2       from    dta_lfil                1  i    j    k-1 t-1 
     288ctin     dta_lfil       3       from    ztuv_lfil               1  i-1  j    k  t  
     289ctin     dta_lfil       4       from    ztuv_lfil               1  i    j    k  t 
     290ctin     dta_lfil       5       from    ztuv_lfil               2  i    j-1  k  t  
     291ctin     dta_lfil       6       from    ztuv_lfil               2  i    j    k  t  
     292 
     293//#---------->ztuv_lfil 
     294ctin     ztuv_lfil      1       from    dta_lfil                1  i+1  j    k     t-1 
     295ctin     ztuv_lfil      2       from    dta_lfil                1  i    j    k     t-1  
     296ctin     ztuv_lfil      3       from    dta_lfil                1  i    j+1  k     t-1  
     297 
     298//#---------->dta_zfil 
     299ctin     dta_zfil       1       from    dta_lfil                1  i    j    k     t 
     300ctin     dta_zfil       2       from    dta_zfil                1  i    j    k     t-1  
     301ctin     dta_zfil       3       from    zwy_zfil                1  i    j    k+1   t  
     302ctin     dta_zfil       4       from    zwy_zfil                1  i    j    k     t 
     303 
     304//#---------->zwy_zfil 
     305ctin     zwy_zfil       1       from    dta_zfil                1  i    j    k-1   t-1 
     306ctin     zwy_zfil       2       from    dta_zfil                1  i    j    k     t-1  
     307 
     308//#---------->ta_c 
     309ctin     ta_c           1       from    dta_zfil                1  i    j    k-1   t 
     310#endif 
    274311//#ctin  ta_c                   1..NPCA from    pca_ta                  1..NPCA  1 
     312 
     313//exec disp_ct_in 
     314exec disp_modul 
     315 
    275316 
    276317//#---------->sbcgyre 
    … …  
    697738ctin     wa                     1       from    wa                      1  i    j    k+1  t 
    698739ctin     wa                     2       from    hdivn                   1  i    j    k    t 
    699  
    700740//#----------> 
    701741 
    … …  
    836876#ifdef FILTER 
    837877// calcul du filtre 
    838 order modinspace S0d 
    839        order YA1 
    840             filtre 
    841       forder 
     878 
     879order modinspace S3dtlf 
     880      order YA3   
     881            order YA2 YA1 
     882                  ztuv_lfil 
     883            forder 
     884            order YA2 YA1 
     885                  dta_lfil 
     886            forder 
     887      forder 
     888forder 
     889 
     890order modinspace S3dtzf 
     891    order YA3 YA2 YA1 
     892          zwy_zfil 
     893    forder 
     894    order YA3 YA2 YA1 
     895          dta_zfil 
     896    forder 
    842897forder 
    843898order modinspace S3d 
    … …  
    846901      forder 
    847902forder 
     903 
    848904#endif    
    849905 
    … …  
    10441100#endif 
    10451101 
     1102#ifdef FILTER 
     1103 
     1104order spaceintraj Tlfil 
     1105        S3dtlf 
     1106forder 
     1107order spaceintraj Tzfil 
     1108        S3dtzf 
     1109forder 
     1110order spaceintraj Tcst 
     1111        S3d 
     1112forder 
     1113 
     1114#endif 
     1115 
    10461116#ifndef SOLSORYAO 
    1047 #ifdef FILTER 
    1048 order spaceintraj Tcst 
    1049         S0d S3d 
    1050 forder 
    1051 #endif 
    10521117 
    10531118order spaceintraj Tsbd 

  • branches/branche-mb/src/sinobad.h

    r58 r59  
    7676  nday_year : current day of the year nyea 
    7777  ndastp    : = nyear*10000 + nmonth*100 + nday*/ 
     78 
     79#ifdef K_FILTER 
     80double pdth_fil=7200.; 
     81double pdtz_fil=7200.; 
     82double fsaht_fil= 1.2e+2;//1000.; 
     83double avt_fil=1.2e-1; 
     84#endif 
    7885 
    7986void cal_ff(); 

  • branches/branche-mb/src/ta_c.h

    r58 r59  
    66//=========================================================================== 
    77//                          methode forward 
    8 forward () 
     8forward (YREAL x1) 
    99{  
    10   YS1 = ta_eb(Yi,Yj,Yk) +dta_c(Yi,Yj,Yk); 
     10  YS1 = ta_eb(Yi,Yj,Yk) +x1; 
    1111 
    1212 
    … …  
    2020//                         methode  backward 
    2121 
    22 backward () 
     22backward (YREAL x1) 
    2323{ 
    24   //gradient se calcule dans le module filtre 
    25    
     24  YJ1I1=1.;  
     25     
    2626  //for(int n=0;n<NPCA;n++)      
    2727  //  Yjac[0][n]    = stdev_ta[n]*shfs_ta(Yi,Yj,Yk,n); 
Note: See TracChangeset for help on using the changeset viewer.