Changeset 32
- Timestamp:
- 04/15/11 09:57:32 (13 years ago)
- Location:
- trunk
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/data_out/plot_results.m
r16 r32 1 f=netcdf('exp_T .nc');1 f=netcdf('exp_TS.nc'); 2 2 3 3 addpath('../scripts/matlab_toolbox'); -
trunk/scripts/sinobad.i
r19 r32 64 64 setm_dxmin 1.0e-4 65 65 setm_epsg 1.0e-10 66 setm_ddf1 1 0066 setm_ddf1 1 67 67 runm 68 68 ##SAUVEGARDE -
trunk/src/Yworkdir/Y1sinobad.h
r14 r32 81 81 //~~~~~~~~~~~~~~~~~~~~~~~~~ module ta_c~~~~~~~~~~~~~~~~~~~~~~~~~ 82 82 #define Yid_ta_c 1 83 #define YNBPTARGET_ta_c 084 #define YDEBTARGET_ta_c 085 #define YENDTARGET_ta_c 086 83 /*M3_VALST */ #define Ydelta_ta_c(v, i, j, k, t) Yta_c[i][j][k]->Ydelta[t][v] 87 84 /*M3_VALST */ #define YD_ta_c(v, i, j, k, t) Yta_c[i][j][k]->Ydelta[t][v] … … 101 98 /*M3_VALSGI */ #define YG1_ta_c(i, j, k) Yta_c[i][j][k]->Ygrad[0] 102 99 /*M3_HERESGI*/ #define YG1Y_ta_c Yta_c[Yi][Yj][Yk]->Ygrad[0] 103 /*M3_VALSG */ #define Yepsi_ta_c(v, i, j, k) Yta_c[i][j][k]->Yepsi[v]104 /*M3_VALSG */ #define YE_ta_c(v, i, j, k) Yta_c[i][j][k]->Yepsi[v]105 /*M3_VALSGI */ #define YE1_ta_c(i, j, k) Yta_c[i][j][k]->Yepsi[0]106 100 Yaota_c *Yta_c[YA1_S3d][YA2_S3d][YA3_S3d]; 107 /*M3_VALSG */ #define Ywish_ta_c(v, i, j, k) Yta_c[i][j][k]->Ywish[v]108 /*M3_VALSG */ #define YW_ta_c(v, i, j, k) Yta_c[i][j][k]->Ywish[v]109 /*M3_VALSGI */ #define YW1_ta_c(i, j, k) Yta_c[i][j][k]->Ywish[0]110 /*M3_HERESGI*/ #define YW1Y_ta_c Yta_c[Yi][Yj][Yk]->Ywish[0]111 101 112 102 //~~~~~~~~~~~~~~~~~~~~~~~~~ module sa_c~~~~~~~~~~~~~~~~~~~~~~~~~ 113 103 #define Yid_sa_c 2 104 #define YNBPTARGET_sa_c 0 105 #define YDEBTARGET_sa_c 0 106 #define YENDTARGET_sa_c 0 114 107 /*M3_VALST */ #define Ydelta_sa_c(v, i, j, k, t) Ysa_c[i][j][k]->Ydelta[t][v] 115 108 /*M3_VALST */ #define YD_sa_c(v, i, j, k, t) Ysa_c[i][j][k]->Ydelta[t][v] … … 129 122 /*M3_VALSGI */ #define YG1_sa_c(i, j, k) Ysa_c[i][j][k]->Ygrad[0] 130 123 /*M3_HERESGI*/ #define YG1Y_sa_c Ysa_c[Yi][Yj][Yk]->Ygrad[0] 124 /*M3_VALSG */ #define Yepsi_sa_c(v, i, j, k) Ysa_c[i][j][k]->Yepsi[v] 125 /*M3_VALSG */ #define YE_sa_c(v, i, j, k) Ysa_c[i][j][k]->Yepsi[v] 126 /*M3_VALSGI */ #define YE1_sa_c(i, j, k) Ysa_c[i][j][k]->Yepsi[0] 131 127 Yaosa_c *Ysa_c[YA1_S3d][YA2_S3d][YA3_S3d]; 128 /*M3_VALSG */ #define Ywish_sa_c(v, i, j, k) Ysa_c[i][j][k]->Ywish[v] 129 /*M3_VALSG */ #define YW_sa_c(v, i, j, k) Ysa_c[i][j][k]->Ywish[v] 130 /*M3_VALSGI */ #define YW1_sa_c(i, j, k) Ysa_c[i][j][k]->Ywish[0] 131 /*M3_HERESGI*/ #define YW1Y_sa_c Ysa_c[Yi][Yj][Yk]->Ywish[0] 132 132 133 133 //~~~~~~~~~~~~~~~~~~~~~~~~~ module ua_c~~~~~~~~~~~~~~~~~~~~~~~~~ … … 2403 2403 double Ystate[YNBS_ta_c]; 2404 2404 double Ygrad[YNBS_ta_c]; 2405 double Yepsi[YNBS_ta_c];2406 double Ywish[YNBS_ta_c];2407 2405 double Ydelta[1][YNBS_ta_c]; 2408 2406 … … 2419 2417 double Ystate[YNBS_sa_c]; 2420 2418 double Ygrad[YNBS_sa_c]; 2419 double Yepsi[YNBS_sa_c]; 2420 double Ywish[YNBS_sa_c]; 2421 2421 double Ydelta[1][YNBS_sa_c]; 2422 2422 … … 4051 4051 void Yrazgrad_ta_c(); 4052 4052 void Ysetstate_ta_c(); 4053 void Ysetepsi_ta_c(); 4054 void Yadjust_ta_c(); 4055 void Y3getstate_ta_c(float x[]); 4056 void Y3valstate_ta_c(); 4057 void Y3valgrad_ta_c(float g[]); void Yrazgrad_sa_c(); 4053 void Yrazgrad_sa_c(); 4058 4054 void Ysetstate_sa_c(); 4059 void Yrazgrad_ua_c(); 4055 void Ysetepsi_sa_c(); 4056 void Yadjust_sa_c(); 4057 void Y3getstate_sa_c(float x[]); 4058 void Y3valstate_sa_c(); 4059 void Y3valgrad_sa_c(float g[]); void Yrazgrad_ua_c(); 4060 4060 void Ysetstate_ua_c(); 4061 4061 void Yrazgrad_va_c(); … … 4258 4258 struct Yst_modul YTabMod[87] = { 4259 4259 {"sshn_c", *Ysshn_c, 2, 32, 22, 0, 0, 1, 0, 0, 0, 0, 0, 1.000000e+00, 1.000000e+00, 1.000000e+00, 0, "S3d"}, 4260 {"ta_c", *Yta_c, 3, 32, 22, 31, 0, 1, 0, 0, 1, 0, 0, 1.000000e+00, 1.000000e+00, 1.000000e+00, 0, "S3d"},4261 {"sa_c", *Ysa_c, 3, 32, 22, 31, 0, 1, 0, 0, 0, 0, 0, 1.000000e+00, 1.000000e+00, 1.000000e+00, 0, "S3d"},4260 {"ta_c", *Yta_c, 3, 32, 22, 31, 0, 1, 0, 0, 0, 0, 0, 1.000000e+00, 1.000000e+00, 1.000000e+00, 0, "S3d"}, 4261 {"sa_c", *Ysa_c, 3, 32, 22, 31, 0, 1, 0, 0, 1, 0, 0, 1.000000e+00, 1.000000e+00, 1.000000e+00, 0, "S3d"}, 4262 4262 {"ua_c", *Yua_c, 3, 32, 22, 31, 0, 1, 0, 0, 0, 0, 0, 1.000000e+00, 1.000000e+00, 1.000000e+00, 0, "S3d"}, 4263 4263 {"va_c", *Yva_c, 3, 32, 22, 31, 0, 1, 0, 0, 0, 0, 0, 1.000000e+00, 1.000000e+00, 1.000000e+00, 0, "S3d"}, -
trunk/src/Yworkdir/Y2sinobad.h
r14 r32 111 111 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){ 112 112 if ((Ysshn_c(Yw1, Yw2)= new Yaosshn_c ())==NULL) Ycreate_errnew();}} 113 }/*M3_ADJUS*/114 void Yadjust_ta_c()115 {116 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1)117 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2)118 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3)119 for (int Yws=0; Yws<YNBS_ta_c; ++Yws)120 { Ystate_ta_c(Yws, Yw1, Yw2, Yw3) -= Yepsi_ta_c(Yws, Yw1, Yw2, Yw3) * Ygrad_ta_c(Yws, Yw1, Yw2, Yw3);121 }122 }/*M3_GCTOTB */123 void YgradCTOtab_ta_c(YREAL tab[])124 {125 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){126 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){127 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){128 for (int Yws=0; Yws<YNBS_ta_c; ++Yws){129 tab[Y3windice] += YG_ta_c(Yws, Yw1, Yw2, Yw3); ++Y3windice;}}}}130 }/*M3_VSTA*/131 void Y3valstate_ta_c()132 {133 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){134 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){135 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){136 for (int Yws=0; Yws<YNBS_ta_c; ++Yws){137 Y3x[Y3windice] = Yta_c(Yw1, Yw2, Yw3)->Ystate[Yws]; ++Y3windice;}}}}138 }/*M3_GSTA*/139 void Y3getstate_ta_c(float x[])140 {141 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){142 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){143 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){144 for (int Yws=0; Yws<YNBS_ta_c; ++Yws){145 Yta_c(Yw1, Yw2, Yw3)->Ystate[Yws] = x[Y3windice]; ++Y3windice;}}}}146 }/*M3_VGRAD*/147 void Y3valgrad_ta_c(float g[])148 {149 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){150 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){151 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){152 for (int Yws=0; Yws<YNBS_ta_c; ++Yws){153 g[Y3windice] = Yta_c(Yw1, Yw2, Yw3)->Ygrad[Yws]; ++Y3windice;}}}}154 113 }/*M3_TBTOG */ 155 114 void YtabTOgrad_ta_c(YREAL tab[]) … … 265 224 void Yio_ta_c(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val) 266 225 { Yiorwf ("ta_c", Yws, Yw1, Yw2, Yw3, Ywt, &YS_ta_c(Yws, Yw1, Yw2, Yw3), val); 267 }/*M0_OUTOB*/268 int Youtoobs_ta_c(YioKind yiokind, int Yws, int Ywt, int arbpdt)269 {270 Yaobs.time = arbpdt;271 Yaobs.imod = Yimod("ta_c");272 Yaobs.smod = Yws;/*M3_OUTOB*/273 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1)274 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2)275 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3)276 { if(select_io(yiokind, "ta_c", Yws, Yw1, Yw2, Yw3, Ywt, &Ystate_ta_c(Yws, Yw1, Yw2, Yw3)))277 { Yaobs.iaxe = Yw1; Yaobs.jaxe = Yw2; Yaobs.kaxe = Yw3;278 Yaobs.vobs = Ystate_ta_c(Yws, Yw1, Yw2, Yw3);279 if(!Yobs_insert (&Yaobs)) return(0);280 }281 }/*M9_OUTOB*/282 return(1);283 }/*M3_SEPS*/284 void Ysetepsi_ta_c(double val)285 {286 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){287 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){288 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){289 for (int Yws=0; Yws<YNBS_ta_c; ++Yws){290 Yta_c(Yw1, Yw2, Yw3)->Yepsi[Yws] = val;}}}}291 226 }/*M3_CREY*/ 292 227 void Ycreate_ta_c() … … 296 231 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){ 297 232 if ((Yta_c(Yw1, Yw2, Yw3) = new Yaota_c ())==NULL) Ycreate_errnew();}}} 298 }/*M3_SWISH*/ 299 void Ysetwish_ta_c(double val) 233 }/*M3_ADJUS*/ 234 void Yadjust_sa_c() 235 { 236 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1) 237 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2) 238 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3) 239 for (int Yws=0; Yws<YNBS_sa_c; ++Yws) 240 { Ystate_sa_c(Yws, Yw1, Yw2, Yw3) -= Yepsi_sa_c(Yws, Yw1, Yw2, Yw3) * Ygrad_sa_c(Yws, Yw1, Yw2, Yw3); 241 } 242 }/*M3_GCTOTB */ 243 void YgradCTOtab_sa_c(YREAL tab[]) 300 244 { 301 245 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){ 302 246 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){ 303 247 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){ 304 for (int Yws=0; Yws<YNBS_ta_c; ++Yws){ 305 Ywish_ta_c(Yws, Yw1, Yw2, Yw3) = val;}}}} 248 for (int Yws=0; Yws<YNBS_sa_c; ++Yws){ 249 tab[Y3windice] += YG_sa_c(Yws, Yw1, Yw2, Yw3); ++Y3windice;}}}} 250 }/*M3_VSTA*/ 251 void Y3valstate_sa_c() 252 { 253 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){ 254 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){ 255 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){ 256 for (int Yws=0; Yws<YNBS_sa_c; ++Yws){ 257 Y3x[Y3windice] = Ysa_c(Yw1, Yw2, Yw3)->Ystate[Yws]; ++Y3windice;}}}} 258 }/*M3_GSTA*/ 259 void Y3getstate_sa_c(float x[]) 260 { 261 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){ 262 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){ 263 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){ 264 for (int Yws=0; Yws<YNBS_sa_c; ++Yws){ 265 Ysa_c(Yw1, Yw2, Yw3)->Ystate[Yws] = x[Y3windice]; ++Y3windice;}}}} 266 }/*M3_VGRAD*/ 267 void Y3valgrad_sa_c(float g[]) 268 { 269 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){ 270 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){ 271 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){ 272 for (int Yws=0; Yws<YNBS_sa_c; ++Yws){ 273 g[Y3windice] = Ysa_c(Yw1, Yw2, Yw3)->Ygrad[Yws]; ++Y3windice;}}}} 306 274 }/*M3_TBTOG */ 307 275 void YtabTOgrad_sa_c(YREAL tab[]) … … 417 385 void Yio_sa_c(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val) 418 386 { Yiorwf ("sa_c", Yws, Yw1, Yw2, Yw3, Ywt, &YS_sa_c(Yws, Yw1, Yw2, Yw3), val); 387 }/*M0_OUTOB*/ 388 int Youtoobs_sa_c(YioKind yiokind, int Yws, int Ywt, int arbpdt) 389 { 390 Yaobs.time = arbpdt; 391 Yaobs.imod = Yimod("sa_c"); 392 Yaobs.smod = Yws;/*M3_OUTOB*/ 393 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1) 394 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2) 395 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3) 396 { if(select_io(yiokind, "sa_c", Yws, Yw1, Yw2, Yw3, Ywt, &Ystate_sa_c(Yws, Yw1, Yw2, Yw3))) 397 { Yaobs.iaxe = Yw1; Yaobs.jaxe = Yw2; Yaobs.kaxe = Yw3; 398 Yaobs.vobs = Ystate_sa_c(Yws, Yw1, Yw2, Yw3); 399 if(!Yobs_insert (&Yaobs)) return(0); 400 } 401 }/*M9_OUTOB*/ 402 return(1); 403 }/*M3_SEPS*/ 404 void Ysetepsi_sa_c(double val) 405 { 406 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){ 407 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){ 408 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){ 409 for (int Yws=0; Yws<YNBS_sa_c; ++Yws){ 410 Ysa_c(Yw1, Yw2, Yw3)->Yepsi[Yws] = val;}}}} 419 411 }/*M3_CREY*/ 420 412 void Ycreate_sa_c() … … 424 416 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){ 425 417 if ((Ysa_c(Yw1, Yw2, Yw3) = new Yaosa_c ())==NULL) Ycreate_errnew();}}} 418 }/*M3_SWISH*/ 419 void Ysetwish_sa_c(double val) 420 { 421 for (int Yw1=0; Yw1<YA1_S3d; ++Yw1){ 422 for (int Yw2=0; Yw2<YA2_S3d; ++Yw2){ 423 for (int Yw3=0; Yw3<YA3_S3d; ++Yw3){ 424 for (int Yws=0; Yws<YNBS_sa_c; ++Yws){ 425 Ywish_sa_c(Yws, Yw1, Yw2, Yw3) = val;}}}} 426 426 }/*M3_TBTOG */ 427 427 void YtabTOgrad_ua_c(YREAL tab[]) … … 13367 13367 { 13368 13368 if(1==0); 13369 else if (strcmp(nmmod, " ta_c") == 0)13370 YW_ ta_c(Yws, Yw1, Yw2, Yw3) = YS_ta_c(Yws, Yw1, Yw2, Yw3) - vobs;13369 else if (strcmp(nmmod, "sa_c") == 0) 13370 YW_sa_c(Yws, Yw1, Yw2, Yw3) = YS_sa_c(Yws, Yw1, Yw2, Yw3) - vobs; 13371 13371 else if (strcmp(nmmod, "tb") == 0) 13372 13372 YW_tb(Yws, Yw1, Yw2, Yw3) = YS_tb(Yws, Yw1, Yw2, Yw3, Ywt) - vobs; … … 13383 13383 { 13384 13384 if(1==0); 13385 else if (strcmp(nmmod, " ta_c") == 0)13386 YG_ ta_c(Yws, Yw1, Yw2, Yw3) += Ycostdiff(YW_ta_c(Yws, Yw1, Yw2, Yw3), YW_ta_c(Yws, Yw1, Yw2, Yw3));13385 else if (strcmp(nmmod, "sa_c") == 0) 13386 YG_sa_c(Yws, Yw1, Yw2, Yw3) += Ycostdiff(YW_sa_c(Yws, Yw1, Yw2, Yw3), YW_sa_c(Yws, Yw1, Yw2, Yw3)); 13387 13387 else if (strcmp(nmmod, "tb") == 0) 13388 13388 YG_tb(Yws, Yw1, Yw2, Yw3, Ywt) += Ycostdiff(YW_tb(Yws, Yw1, Yw2, Yw3), YW_tb(Yws, Yw1, Yw2, Yw3)); … … 13399 13399 { 13400 13400 if(1==0); 13401 else if (strcmp(nmmod, " ta_c") == 0) Youtoobs_ta_c(yiokind, numout, pdt, arbpdt);13401 else if (strcmp(nmmod, "sa_c") == 0) Youtoobs_sa_c(yiokind, numout, pdt, arbpdt); 13402 13402 else if (strcmp(nmmod, "tb") == 0) Youtoobs_tb(yiokind, numout, pdt, arbpdt); 13403 13403 else if (strcmp(nmmod, "sb") == 0) Youtoobs_sb(yiokind, numout, pdt, arbpdt); … … 13408 13408 void Yadjust_all () 13409 13409 { 13410 Yadjust_ ta_c ();13410 Yadjust_sa_c (); 13411 13411 } 13412 13412 void Yrazgrad_all () … … 13689 13689 void Ysetwish_mod (int imod, YREAL val) 13690 13690 { if(1==0); 13691 else if (imod== 1) Ysetwish_ta_c(val);13691 else if (imod==2) Ysetwish_sa_c(val); 13692 13692 else if (imod==38) Ysetwish_tb(val); 13693 13693 else if (imod==39) Ysetwish_sb(val); … … 13789 13789 void YgradCTOtab_target (YREAL tab[]) 13790 13790 { //Y3windice = 0; 13791 YgradCTOtab_ ta_c (tab);13791 YgradCTOtab_sa_c (tab); 13792 13792 } 13793 13793 void YtabTOgrad_target (YREAL tab[]) 13794 13794 { //Y3windice = 0; 13795 YtabTOgrad_ ta_c (tab);13795 YtabTOgrad_sa_c (tab); 13796 13796 } 13797 13797 void Y3getstate_all (float x[]) 13798 13798 { Y3windice = 0; 13799 Y3getstate_ ta_c (x);13799 Y3getstate_sa_c (x); 13800 13800 } 13801 13801 void Y3valstate_all () 13802 13802 { Y3windice = 0; 13803 Y3valstate_ ta_c ();13803 Y3valstate_sa_c (); 13804 13804 } 13805 13805 void Y3valgrad_all (float g[]) 13806 13806 { Y3windice = 0; 13807 Y3valgrad_ ta_c (g);13807 Y3valgrad_sa_c (g); 13808 13808 } 13809 13809 void YdeltaEQPstate_traj(int itraj, char *norkmod, int frompdt, int topdt, double pfact) … … 13814 13814 if ( !strcmp(norkmod, "sshn_c") || !strcmp(norkmod, "Y#A") ) 13815 13815 YdeltaEQPstate_sshn_c (0, 1, pfact); 13816 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T"))13816 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") ) 13817 13817 YdeltaEQPstate_ta_c (0, 1, pfact); 13818 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") )13818 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") ) 13819 13819 YdeltaEQPstate_sa_c (0, 1, pfact); 13820 13820 if ( !strcmp(norkmod, "ua_c") || !strcmp(norkmod, "Y#A") ) … … 13993 13993 void YdeltaEQPstate_target(double pfact) 13994 13994 { 13995 YdeltaEQPstate_ ta_c (0, 1, pfact);13995 YdeltaEQPstate_sa_c (0, 1, pfact); 13996 13996 } 13997 13997 void YdeltaEQPgrad_target(double pfact) 13998 13998 { 13999 YdeltaEQPgrad_ ta_c (0, 1, pfact);13999 YdeltaEQPgrad_sa_c (0, 1, pfact); 14000 14000 } 14001 14001 void YstateEQPdelta_traj(int itraj, char *norkmod, int frompdt, int topdt, double pfact) … … 14006 14006 if ( !strcmp(norkmod, "sshn_c") || !strcmp(norkmod, "Y#A") ) 14007 14007 YstateEQPdelta_sshn_c (0, 1, pfact); 14008 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T"))14008 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") ) 14009 14009 YstateEQPdelta_ta_c (0, 1, pfact); 14010 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") )14010 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") ) 14011 14011 YstateEQPdelta_sa_c (0, 1, pfact); 14012 14012 if ( !strcmp(norkmod, "ua_c") || !strcmp(norkmod, "Y#A") ) … … 14185 14185 void YstateEQPdelta_target(double pfact) 14186 14186 { 14187 YstateEQPdelta_ ta_c (0, 1, pfact);14187 YstateEQPdelta_sa_c (0, 1, pfact); 14188 14188 } 14189 14189 void YstateEQPOdelta_target(double pfact, char *codop) 14190 14190 { 14191 YstateEQPOdelta_ ta_c (0, 1, pfact, codop);14191 YstateEQPOdelta_sa_c (0, 1, pfact, codop); 14192 14192 } 14193 14193 void YstateEQAPTdelta_target(double pfact, YREAL tab[]) 14194 14194 { 14195 YstateEQAPTdelta_ ta_c (0, 1, pfact, tab);14195 YstateEQAPTdelta_sa_c (0, 1, pfact, tab); 14196 14196 } 14197 14197 void YgradEQPstate_traj(int itraj, char *norkmod, int frompdt, int topdt, double pfact) … … 14202 14202 if ( !strcmp(norkmod, "sshn_c") || !strcmp(norkmod, "Y#A") ) 14203 14203 YgradEQPstate_sshn_c(pfact); 14204 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T"))14204 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") ) 14205 14205 YgradEQPstate_ta_c(pfact); 14206 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") )14206 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") ) 14207 14207 YgradEQPstate_sa_c(pfact); 14208 14208 if ( !strcmp(norkmod, "ua_c") || !strcmp(norkmod, "Y#A") ) … … 14381 14381 void YgradEQPstate_target(double pfact) 14382 14382 { 14383 YgradEQPstate_ ta_c (pfact);14383 YgradEQPstate_sa_c (pfact); 14384 14384 } 14385 14385 void YgradEQPOstate_target(double pfact, char *codop) 14386 14386 { 14387 YgradEQPOstate_ ta_c (pfact, codop);14387 YgradEQPOstate_sa_c (pfact, codop); 14388 14388 } 14389 14389 void YstateTOtab_traj (int itraj, char *norkmod, int frompdt, int topdt, YREAL tab[]) … … 14394 14394 if ( !strcmp(norkmod, "sshn_c") || !strcmp(norkmod, "Y#A") ) 14395 14395 YstateTOtab_sshn_c (tab); 14396 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T"))14396 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") ) 14397 14397 YstateTOtab_ta_c (tab); 14398 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") )14398 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") ) 14399 14399 YstateTOtab_sa_c (tab); 14400 14400 if ( !strcmp(norkmod, "ua_c") || !strcmp(norkmod, "Y#A") ) … … 14573 14573 void YstateTOtab_target (YREAL tab[]) 14574 14574 { //Y3windice = 0; 14575 YstateTOtab_ ta_c (tab);14575 YstateTOtab_sa_c (tab); 14576 14576 } 14577 14577 void YgradTOtab_traj (int itraj, char *norkmod, int frompdt, int topdt, YREAL tab[]) … … 14582 14582 if ( !strcmp(norkmod, "sshn_c") || !strcmp(norkmod, "Y#A") ) 14583 14583 YgradTOtab_sshn_c (tab); 14584 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T"))14584 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") ) 14585 14585 YgradTOtab_ta_c (tab); 14586 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") )14586 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") ) 14587 14587 YgradTOtab_sa_c (tab); 14588 14588 if ( !strcmp(norkmod, "ua_c") || !strcmp(norkmod, "Y#A") ) … … 14761 14761 void YgradTOtab_target (YREAL tab[]) 14762 14762 { //Y3windice = 0; 14763 YgradTOtab_ ta_c (tab);14763 YgradTOtab_sa_c (tab); 14764 14764 } 14765 14765 void YtabTOgrad_traj (int itraj, char *norkmod, int frompdt, int topdt, YREAL tab[]) … … 14770 14770 if ( !strcmp(norkmod, "sshn_c") || !strcmp(norkmod, "Y#A") ) 14771 14771 YtabTOgrad_sshn_c (tab); 14772 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T"))14772 if ( !strcmp(norkmod, "ta_c") || !strcmp(norkmod, "Y#A") ) 14773 14773 YtabTOgrad_ta_c (tab); 14774 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") )14774 if ( !strcmp(norkmod, "sa_c") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") ) 14775 14775 YtabTOgrad_sa_c (tab); 14776 14776 if ( !strcmp(norkmod, "ua_c") || !strcmp(norkmod, "Y#A") ) … … 15127 15127 { 15128 15128 if (!strcmp(argv[0], "setepsi_all") || !strcmp(argv[0], "SETEPSI_ALL")){all=1; val=atof(argv[1]);} else val=atof(argv[2]); 15129 if (!strcmp(argv[1], " ta_c") || all)15130 {Ysetepsi_ ta_c(val);codret=1;}15129 if (!strcmp(argv[1], "sa_c") || all) 15130 {Ysetepsi_sa_c(val);codret=1;} 15131 15131 } 15132 15132 else codret=0; -
trunk/src/cost.dat
r14 r32 2 2 0.0000000000000000e+00 3 3 0.0000000000000000e+00 4 2.9733593674719423e+03 5 2.9046564674634410e+03 6 2.8844078750994972e+03 7 2.1533662170138909e+03 8 3.7645353988317231e+02 9 2.1743622020262393e+02 10 1.3325376402313574e+02 11 1.0089472892137428e+02 12 9.3476389696007189e+01 13 8.8719928304913211e+01 14 9.0834543682672248e+01 15 8.8130698029267435e+01 16 9.0170152793152795e+01 17 8.8265477839190027e+01 18 8.8281377590830218e+01 19 8.8710953227195915e+01 20 8.9177210220475857e+01 21 9.0343293430849968e+01 22 9.0156125898119669e+01 23 8.7443014766499118e+01 24 8.7378898959154185e+01 25 9.0197838294519727e+01 26 9.0231123456027703e+01 27 9.0251747653212306e+01 28 9.0266310888836969e+01 29 8.7354177212549217e+01 30 8.7352312614036819e+01 31 9.0268596517281566e+01 32 9.0270225161585003e+01 4 7.9500011475282708e+00 5 6.9584567722863468e+01 6 6.3544557715163670e+00 7 4.5195261866713539e+00 8 4.1368024754435542e+00 9 4.1173925937788400e+00 10 4.4705004404786965e+00 11 4.4772946613443398e+00 12 4.3418385150596306e+00 13 4.3353102149522806e+00 14 4.2495486911837279e+00 15 4.2535273008791705e+00 16 4.2345336365761872e+00 17 4.2197520449371364e+00 18 4.2092179852297607e+00 19 4.1123363902328558e+00 20 4.1113762968592953e+00 21 4.1107720985439284e+00 22 4.1103175097543936e+00 23 4.2097782534273280e+00 -
trunk/src/sinobad.d
r14 r32 47 47 modul sshn_c space S3d inter YA1 YA2 noward output 1 |target 48 48 49 modul ta_c space S3d noward output 1 target49 modul ta_c space S3d noward output 1 |target 50 50 51 51 #modul ta_c space S3d input array NPCA output 1 52 52 53 modul sa_c space S3d noward output 1 |target53 modul sa_c space S3d noward output 1 target 54 54 modul ua_c space S3d noward output 1 |target 55 55 modul va_c space S3d noward output 1 |target
Note: See TracChangeset
for help on using the changeset viewer.