Changeset 61 for trunk/src/sinobad.d


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

correction du bug du LT avec solsor

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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 
Note: See TracChangeset for help on using the changeset viewer.