source: altifloat/src/floater.h @ 108

//constantes
YREAL rdtflo=2*3600;
//YREAL rdtflo=1;
YREAL ur_dx=0;
YREAL ur_dy=0;
YREAL vr_dx=0;
YREAL vr_dy=0;

//FILTER COEFFICIENTS
#ifdef K_FILTER 
//YREAL c1;
//YREAL c2,c3,c4,c5;
//YREAL dtfil=26042000;
YREAL dtfil=28125000;

#endif

//Options du run incrŽmental
int nb_extiter=16;
short inc_save=1;
char dirsave[50]="../obs_float/";
char savename[50]="optim25float-filtre15-o";
extern float Y3ddf1;
//Declaration 

void appli_start(int argc, char *argv[]){
#ifdef K_FILTER
  //c1=0.5;
  //c2=c3=c4=c5=0.125;
#endif
}
void before_it(int nit){}
void cost_function(int pdt){}
void adjust_target(){}
void after_it(int nit){
        
        //printf("u_d(9,16)=%f %f %f %f\n",YS1_u_d(9,16),YS1_u_d(10,16),YS1_u_d(10,17),YS1_u_d(9,17));
        }
void forward_before(int ctrp){
}
void forward_after(int ctrp){
}

void backward_before(int ctrp){
}

void backward_after(int ctrp){
        }

short select_io(int indic, char *nmod, int sortie, int iaxe,
                int jaxe, int kaxe, int pdt, YREAL *val){
        //if(indic==YIO_LOADSTATE && !strcmp(nmod,"u")) {
        //printf("s=%d t=%d Yi=%d Yj=%d val=%f\n",sortie,pdt,iaxe,jaxe,*val);   
        //}
        
        return(1);
}

//in case we want to see how it s loading the data we do this above

void update_uv(){
        int j,k;
        for (j=0;j<nlon;j++)
        for (k=0;k<nlat;k++) {
                YS1_u(j,k)+=YS1_u_d(j,k);
                YS1_v(j,k)+=YS1_v_d(j,k);
                //printf("u=%f delu=%f \n",YS1_u(j,k),YS1_u_d(j,k));    

}
}

void erase_udvd() {
  int j,k;
  for (j=0;j<nlon;j++)
        for (k=0;k<nlat;k++) {
#ifdef K_FILTER
          YS1_uc_d(j,k)=0;
          YS1_vc_d(j,k)=0;
#else     
          YS1_u_d(j,k)=0;
          YS1_v_d(j,k)=0;
#endif
}
}
//void load_aviso(int argc, char *argv[]){
        
//      }

void run_inc(int argv, char *argc[]) {
//Run the incremental optimisation
//1 forward of the complete model
//2 initialize du et dv to zero
//3 minimize the incremental cost function
//4 update u and v

/*
Options
nb_extiter : number of extern loop (default=6)
inc_save : save option
        0  : no save
        1 (default): save after each extern loop r_float,u and v
dirsave : directory in case inc_save>0 (default="../obs_float/")
savename : basename of save fils (default "optim")
*/
char filesave[100];
int i;

char sactiv[20]="activ";
char sMD[20]="M";
char sonly[20]="only";
char *liste[3];
liste[0]=sactiv;
liste[1]=sMD;
liste[2]=sonly;
/*Check M1QN3 config*/
 if (YNbItRun<=0)
                {        printf("runm(2): number of run iteration not seted; use set_nbiter please\n");
                         //return(0);
                }
        if (Y3ddf1<=0)
                {       printf("runm(2): expected positive fcost decrease missed; use setm_ddf1\n");
                         //return(0);
                }
        if (YioInsertObsCtr<0)
                   printf("runm(2): warning : oh oh, run with no obs !!! \n");
                
        YTypeAdjust = ADJUST_M1QN3; //d'office avec M1QN3
        
//Savestate config
 YioModulot = OFF; //On sauvegarde tous les pas de temps
    YioWrite = ON;
    YioRead= OFF;
   // Yio_savestate(cdes[1], cdes[3], 0, cdes[7]);
    YioState=0 ;//Save all states
    YioBin=OFF; //ascii output
    YioAscii=ON; //ascii output
    YioAxes=ON; //save axe numbers

for (i=0;i<nb_extiter;i++) {
     strcpy(liste[1],"M");
         Yactraj(3, liste);
         Yset_modeltime(0);
         before_it(1);
         //printf("---forward(i=%d)---\n",i);
         Yforward(-1, 0);
         
         if (inc_save==1 && i>0) {
          //save rfloat
          YioTime=ON;
          sprintf(filesave,"%sr_float_%s_%d.dat",dirsave,savename,i);
          Yio_savestate("r_float","i",0,filesave);
         }
             strcpy(liste[1],"D");
         Yactraj(3,liste);
         Yset_modeltime(0);
         
        
        Y3run ('0');
#ifdef K_FILTER
        //compute values for u_d and v_d
        strcpy(liste[1],"Tfil");
        Yactraj(3,liste);
        Yset_modeltime(0);
        before_it(1);
        Yforward(-1, 0);
        strcpy(liste[1],"T_euler_d");
        Yactraj(3,liste);
        Yset_modeltime(0);
        before_it(1);
        Yforward(-1, 0);
#endif


    update_uv();
    erase_udvd();

    if (inc_save==1) {
   
    
    //save u
    YioTime=OFF;
    sprintf(filesave,"%su%s%d.dat",dirsave,savename,i+1);
    Yio_savestate("u","ij",0,filesave);
    
    //save v
    sprintf(filesave,"%sv%s%d.dat",dirsave,savename,i+1);
    Yio_savestate("v","ij",0,filesave);
    
    }
    }
    /*End of the optimization, we run another time the model to compute the final value of r_float
    */
     strcpy(liste[1],"M");
    Yactraj(3, liste);
         Yset_modeltime(0);
         before_it(1);
         Yforward(-1, 0);
 if (inc_save==1) {
          //save rfloat
          YioTime=ON;
          sprintf(filesave,"%sr_float_%s_%d.dat",dirsave,savename,nb_extiter);
          Yio_savestate("r_float","i",0,filesave);
         }
    
}