Changeset 61


Ignore:
Timestamp:
07/26/11 11:06:10 (13 years ago)
Author:
jbrlod
Message:

correction du bug du LT avec solsor

Location:
trunk/src
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/cost.dat

    r55 r61  
    11 0.0000000000000000e+00 
    2  0.0000000000000000e+00 
    3  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 

  • trunk/src/gcx2.h

    r54 r61  
    1717  } 
    1818  else 
    19 #ifndef K_SOLSORYAO 
    20      YS1=gcx(Yi,Yj); 
     19#ifdef K_SOLSORYAO 
     20   YS1=  x1; 
    2121#else 
    22                YS1=  x1; 
     22               YS1=gcx(Yi,Yj); 
    2323#endif 
    2424  // 
    … …  
    3030backward (YREAL x) 
    3131{ 
     32#ifdef K_SOLSORYAO 
    3233  if(Yt==TU){ 
    3334    YJ1I1=0; 
    … …  
    3839    YJ1I1=1; 
    3940  // 
     41#endif 
    4042} 
    4143 

  • trunk/src/sinobad.d

    r54 r61  
     1//#define FILTER 
    12//#define SOLSORYAO 
    23//#define OPTIMORDER 
    … …  
    56//#Mohamed Berrada 26/02/2009 
    67//#modification (#define) Julien Brajard 06/2011 
     8//#modification (filtre) Mohamed Berrada 19/07/2011 
    79//#============================================================================== 
    810//#DEFVAL ---------------------------------------------DEFVAL 
    … …  
    1517defval TA      50 |__DATE__ | number of time steps 
    1618defval TU       2 | start time step 
     19#ifdef FILTER 
     20defval OFTL      40 
     21defval OFTZ      40 
     22defval OFT      80 
     23#endif 
     24#ifndef FILTER 
     25defval OFT      0 
     26#endif 
    1727//#exec disp_valdef 
     28#ifdef FILTER 
     29defval K_FILTER 
     30#endif 
    1831#ifdef SOLSORYAO 
    1932defval K_SOLSORYAO 
    … …  
    3043option    O_M1QN3 | M2QN1  
    3144//#option O_VARINCR 
    32 //#option    O_GRADTEST 
     45option    O_GRADTEST 
    3346option   O_DBG_NANF 
    3447//#option   O_DBG_TING 
    3548option  O_EXTOBJ "../object/ncutil.o" 
    3649option O_REAL double 
    37 option o_parallel 16 
     50option o_parallel 
    3851 
    3952//#TRAJECTOIRES (DIM) -----------------------------------TRAJ 
    40 traj  Tcst       M   0   1 
    41 traj  Tsbd       M   TU  0   1  TA 
     53#ifdef FILTER 
     54traj  Tcstf      M   0  1 
     55traj  Tlfil      M   1  0    1  OFTL 
     56traj  Tzfil      M   1  OFTL 1  OFTZ 
     57#endif 
     58traj  Tcst       M   0   OFT  1 1 
     59traj  Tsbd       M   TU  OFT  1  TA 
    4260exec disp_traj 
    4361 
     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 
    4467 
    4568space  S0d       M    1           Tcst 
    … …  
    6285modul    sshn_c                 space   S3d inter YA1 YA2       noward            output 1      |target 
    6386 
     87//#modul         ta_c                   space   S3d     input array NPCA                  output 1       
     88#ifndef FILTER 
    6489modul    ta_c                   space   S3d                     noward            output 1      target 
    65  
    66 //#modul         ta_c                   space   S3d     input array NPCA                  output 1       
     90#endif 
    6791 
    6892modul    sa_c                   space   S3d                     noward            output 1      |target 
    6993modul    ua_c                   space   S3d                     noward            output 1      |target 
    7094modul    va_c                   space   S3d                     noward            output 1      |target 
     95 
     96#ifdef FILTER 
     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       
     103#endif 
    71104 
    72105//#------------- 
    … …  
    206239modul   gcx_dynspg_flt          space   S3dt inter YA1 YA2      input 2           output 1      tempo    
    207240#ifndef SOLSORYAO 
    208 modul   solsor_dynspg_flt       space   S3dt inter YA1 YA2      input array 19    output 1      tempo 
     241modul   solsor_dynspg_flt       space   S3dt inter YA1 YA2      input array 19    output 1      tempo hidjac 
     242modul   gcx2                    space   S3dt inter YA1 YA2      input 1           output 1      tempo hidjac 
    209243#endif 
    210244 
    … …  
    214248modul   solsor_flag             space   Ssolsor                 input 4           output 2      tempo 
    215249#define S3dt S3dt_b 
    216 #endif 
    217  
    218  
    219250modul   gcx2                    space   S3dt inter YA1 YA2      input 1           output 1      tempo    
     251#endif 
     252 
     253 
     254 
    220255modul   spguv_dynspg_flt        space   S3dt inter YA1 YA2      input 3           output 2      tempo    
    221256modul   ua                      space   S3dt                    input 3           output 1      tempo 
    … …  
    246281#endif 
    247282 
    248 //exec disp_modul 
     283exec disp_modul 
    249284//############################################################################################# 
    250285 
     286#ifdef FILTER 
     287//#---------->dta_lfil 
     288ctin     dta_lfil       1       from    nu                      1  i    j    k   
     289ctin     dta_lfil       2       from    dta_lfil                1  i    j    k-1 t-1 
     290ctin     dta_lfil       3       from    ztuv_lfil               1  i-1  j    k  t  
     291ctin     dta_lfil       4       from    ztuv_lfil               1  i    j    k  t 
     292ctin     dta_lfil       5       from    ztuv_lfil               2  i    j-1  k  t  
     293ctin     dta_lfil       6       from    ztuv_lfil               2  i    j    k  t  
     294 
     295//#---------->ztuv_lfil 
     296ctin     ztuv_lfil      1       from    dta_lfil                1  i+1  j    k     t-1 
     297ctin     ztuv_lfil      2       from    dta_lfil                1  i    j    k     t-1  
     298ctin     ztuv_lfil      3       from    dta_lfil                1  i    j+1  k     t-1  
     299 
     300//#---------->dta_zfil 
     301ctin     dta_zfil       1       from    dta_lfil                1  i    j    k     t 
     302ctin     dta_zfil       2       from    dta_zfil                1  i    j    k     t-1  
     303ctin     dta_zfil       3       from    zwy_zfil                1  i    j    k+1   t  
     304ctin     dta_zfil       4       from    zwy_zfil                1  i    j    k     t 
     305 
     306//#---------->zwy_zfil 
     307ctin     zwy_zfil       1       from    dta_zfil                1  i    j    k-1   t-1 
     308ctin     zwy_zfil       2       from    dta_zfil                1  i    j    k     t-1  
     309 
     310//#---------->ta_c 
     311ctin     ta_c           1       from    dta_zfil                1  i    j    k-1   t 
     312#endif 
    251313//#ctin  ta_c                   1..NPCA from    pca_ta                  1..NPCA  1 
     314 
     315//exec disp_ct_in 
     316exec disp_modul 
     317 
    252318 
    253319//#---------->sbcgyre 
    … …  
    674740ctin     wa                     1       from    wa                      1  i    j    k+1  t 
    675741ctin     wa                     2       from    hdivn                   1  i    j    k    t 
    676  
    677742//#----------> 
    678743 
    … …  
    811876#endif 
    812877 
     878#ifdef FILTER 
     879// calcul du filtre 
     880 
     881order modinspace S3dtlf 
     882      order YA3   
     883            order YA2 YA1 
     884                  ztuv_lfil 
     885            forder 
     886            order YA2 YA1 
     887                  dta_lfil 
     888            forder 
     889      forder 
     890forder 
     891 
     892order modinspace S3dtzf 
     893    order YA3 YA2 YA1 
     894          zwy_zfil 
     895    forder 
     896    order YA3 YA2 YA1 
     897          dta_zfil 
     898    forder 
     899forder 
     900order modinspace S3d 
     901       order YA3 YA2 YA1 
     902            ta_c 
     903      forder 
     904forder 
     905 
     906#endif    
     907 
    813908#ifndef OPTIMORDER 
    814  
     909   
    815910order modinspace S3dt 
    816911//#---------->sbcgyre 
    … …  
    10071102#endif 
    10081103 
     1104#ifdef FILTER 
     1105 
     1106order spaceintraj Tlfil 
     1107        S3dtlf 
     1108forder 
     1109order spaceintraj Tzfil 
     1110        S3dtzf 
     1111forder 
     1112order spaceintraj Tcst 
     1113        S3d 
     1114forder 
     1115 
     1116#endif 
     1117 
    10091118#ifndef SOLSORYAO 
     1119 
    10101120order spaceintraj Tsbd 
    10111121        S3dt 
    … …  
    10211131insert_fct  arg xchangesavemode 
    10221132 
     1133insert_fct  arg load_eb 
     1134 
    10231135insert_fct arg xwriteout 
    10241136insert_fct arg xwritegrad 

  • trunk/src/sinobad.h

    r55 r61  
    164164void forward_before (int ctrp) 
    165165{       // permet d'intervenir si besoin avant le forward 
     166  kt=Yt-TU; //oceanic time step 
     167  
    166168    if(Yt==TU+1 && neuler==0) 
    167169      r2dt=rdt; 
    … …  
    172174void forward_after (int ctrp) 
    173175{                       // permet d'intervenir si besoin aprÚs le forward 
     176  printf("neuler=%d\n",neuler); 
    174177} 
    175178//_____________________________________________________________________________ 
Note: See TracChangeset for help on using the changeset viewer.