Changeset 10


Ignore:
Timestamp:
03/16/11 14:02:55 (13 years ago)
Author:
jbrlod
Message:

design a reference version

Location:
trunk
Files:
159 added
1 deleted
14 edited

Legend:

Unmodified
Added
Removed

  • trunk/scripts/sinobad.i

    r9 r10  
    88#paramÚres de contrÃŽle : sshn_c, ta_c, sa_c, ua_c, va_c 
    99#obs (3D a t=48 Yt=50) : sshb, tb, sb,ub,vb 
     10#ATTENTION : la minimisation ne se passe pas bien encore 
    1011 
    1112##INITIALISATION 
    … …  
    5253savestate va_c  1   ijk   0    A       1       ../data_out/va_c_init.dat 
    5354 
     55##TEST DE LA FONCTION OBJECTIVE 
     56goto FINTEST 
     57testof 1.e+00 0.5e-01 2. 10 1 
     58goto FINRUN 
     59FINTEST 
    5460##DEBUT DE L'ASSIMILATION' 
    5561 
    … …  
    6369setm_dxmin   1.0e-4 
    6470setm_epsg    1.0e-10 
    65 setm_ddf1    1 
     71setm_ddf1    100 
    6672runm 
    6773##SAUVEGARDE 

  • trunk/src/cost.dat

    r9 r10  
    33 0.0000000000000000e+00 
    44 2.9733593674719423e+03 
     5 2.9095139883673278e+03 
     6 2.4330312362758295e+03 
     7 2.0642529387596464e+03 
     8 1.9211236976402699e+03 
     9 1.7958968029352288e+03 
     10 1.6298045122109095e+03 
     11 1.5076021103427668e+03 
     12 1.5000076375290523e+03 
     13 1.4276397352270269e+03 
     14 1.3110968128873515e+03 
     15 1.2801012445744680e+03 
     16 1.3104813902608603e+03 
     17 1.2829552646049826e+03 
     18 1.2833299063134032e+03 
     19 1.2847663831547602e+03 
     20 1.2891132528206504e+03 
     21 1.2818428187933409e+03 
     22 1.2764982097924344e+03 
     23 1.2809727315336945e+03 
     24 1.2417113477527846e+03 
     25 1.6709113736713609e+03 
     26 1.2583333872621713e+03 
     27 1.2582520427303459e+03 
     28 1.2550104229729704e+03 
     29 1.2499876871051058e+03 
     30 1.2656188241658340e+03 
     31 1.2643832627133945e+03 
     32 1.2554231947828837e+03 
     33 1.2476636563166230e+03 
     34 1.2489695361098022e+03 
     35 1.2530459381768840e+03 
     36 1.2484393921698424e+03 
     37 1.2470596346544980e+03 
     38 1.2477431325355756e+03 
     39 1.2511259156307392e+03 
     40 1.2407190501617613e+03 
     41 1.2406815136035698e+03 

  • trunk/src/sa.h

    r1 r10  
    4141    YJ1I5=1.; 
    4242#endif 
     43if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1) YJ1I5=0; 
    4344  } 
    4445  else{ 

  • trunk/src/ta.h

    r1 r10  
    4242    //if(Yi==28 && Yj==12)  
    4343    YJ1I5=1.; 
     44    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1) YJ1I5=0; 
    4445    //#endif 
    4546  } 
    46   else{ 
     47   else{ 
    4748    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1){ 
    4849      YJ1I1=0.;     YJ1I2=0.;   YJ1I3=0.;     YJ1I4=0.;   YJ1I5=0.; 

  • trunk/src/ta_de_trazdf_imp.h

    r1 r10  
    1515  /*     4.384959433157585e-05 
    1616         2.134224909062430e-08*/ 
    17   if(Yt==TU) 
     17 if (Yi==30 && Yj==1 && Yt==16) 
     18  { 
     19      printf("Yi=%i,Yj=%i,Yk=%i,Yt=%i,,x1=%f,x2=%f,x3=%f,x4=%f\n",Yk,x1,x2,x3,x4); 
     20      }   
     21 
     22if(Yt==TU) 
    1823    YS1=0.; 
    1924  else{ 

  • trunk/src/ua.h

    r1 r10  
    5050#ifdef YE_ua_c 
    5151    YJ1I3=1.; 
     52    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1) YJ1I3=0; 
    5253#endif 
    5354    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1  ){ 

  • trunk/src/ua_dyn_hpg_zco.h

    r1 r10  
    1616  /*     A 1.579814851804421e-04 
    1717         R 1.270370352774595e-03*/ 
     18 
     19  /* if (Yi==30 && Yj==1 && Yk==0 && Yt==16) 
     20  { 
     21      printf("ua_dyn_hpg_zco:Yk=%i,x1=%f,x2=%f\n",Yk,x1,x2); 
     22      }   
     23  */ 
    1824  if(Yt==TU) 
    1925    YS1=0.; 

  • trunk/src/ua_dyn_ldf_lap.h

    r1 r10  
    1313        x4      from    hdivn                   1  i    j    k    t-1 
    1414        x5      from    hdivn                   1  i+1  j    k    t-1*/ 
     15 
     16 /*if (Yi==30 && Yj==1 && Yk==0 && Yt<13) 
     17  { 
     18    printf("ua_dyn_ldf_lap:Yt=%i,x1=%f,x2=%f,x3=%f, x4=%f, x5=%f\n\n",Yt,x1,x2,x3,x4,x5); 
     19      }   
     20 */ 
    1521 
    1622  if(Yt==TU) 

  • trunk/src/ua_vor_ene.h

    r1 r10  
    1818  /*     1.258248811116506e-12 
    1919         1.604676681128015e-08*/ 
    20   if(Yt==TU) 
     20  /*  if (Yi==30 && Yj==1 && Yk==0 && Yt<13) 
     21  { 
     22    printf("ua_vor_ene:Yt=%i,x1=%f,x2=%f,x3=%f, x4=%f, x5=%f ,x6=%f, x7=%f\n",Yt,x1,x2,x3,x4,x5,x6,x7); 
     23    }    */ 
     24 
     25if(Yt==TU) 
    2126    YS1=0.; 
    2227  else{ 

  • trunk/src/va.h

    r1 r10  
    1616  /*     A 1.966672497070418e-07 
    1717         R 1.869538462034785e-07*/ 
    18  
     18  /*if (Yi>25 && Yj<3 && Yt<13 && Yk==0) 
     19    { 
     20      printf("va:Yi=%i, Yj=%i, Yt=%i,x1=%f,x2=%f,x3=%f\n",Yi,Yj,Yt,x1,x2,x3); 
     21      }*/ 
    1922 
    2023  if(Yt==TU){ 
    … …  
    4144#ifdef YE_va_c 
    4245    YJ1I3=1.; 
     46    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1) YJ1I3=0; 
    4347#endif 
    4448  } 

  • trunk/src/zrhs_ua_dynzdf_imp.h

    r1 r10  
    1414  /*     2.060938356583834e-07 
    1515         2.819052929015263e-10*/ 
     16  /* if (Yi==30 && Yj==1 && Yk==0 && Yt==16) 
     17  { 
     18      printf("zrhs_ua_dynzdf_imp:Yk=%i,x1=%f,x2=%f,x3=%f\n",Yk,x1,x2,x3); 
     19      }  */ 
    1620 
    1721  if(Yt==TU) 

  • trunk/src/zwt_trazdf_imp.h

    r1 r10  
    1515  /*      0 
    1616          0*/ 
     17 
     18  /* if (Yi==30 && Yj==1 && Yt==16) 
     19  { 
     20      printf("Yk=%i,x1=%f,x2=%f,x3=%f,x4=%f\n",Yk,x1,x2,x3,x4); 
     21      }  */ 
    1722  if(Yt==TU) 
    1823    YS1=0.; 

  • trunk/src/zwxyz_vor_ene.h

    r1 r10  
    1313        x4      from    va                      1  i+1  j    k    t-1*/ 
    1414 
     15  /* if (Yi>25 && Yj<3 && Yt<13 && Yk==0) 
     16    { 
     17      printf("zwxyz_vor_ene:Yi=%i, Yj=%i, Yt=%i,x1=%f,x2=%f,x3=%f, x4=%f, YS2=%f\n",Yi,Yj,Yt,x1,x2,x3,x4,e1v(Yi,Yj)*x3); 
     18      }*/ 
     19 
    1520  if(Yt==TU) 
    1621    YS1=0.; 
    … …  
    2126      YS1=e2u(Yi,Yj)*x1; 
    2227      YS2=e1v(Yi,Yj)*x3; 
     28 
    2329      //pour ntot==5---> 
    2430      double y31=(x4+x3)*(e2v(Yi+1,Yj)-e2v(Yi,Yj)); 
Note: See TracChangeset for help on using the changeset viewer.