Changeset 125 for altifloat/src

Timestamp:

06/26/14 09:50:56 (10 years ago)

Author:

jbrlod

Message:

read_bck added

Location:

altifloat/src

Files:

• floater.h (modified) (1 diff)
• floater_delta.d (modified) (1 diff)

altifloat/src/floater.h

@@ -177 +177 @@
 
 }
+
+
+void read_bck (int argv, char *argc[]) {
+/* Read background in a ascii file
+   format
+   Column 1 (%d) : time step (from 1 to jptfl)
+   Column 2 (%d) : lon (Yi) (from 1 to nlon)
+   Column 3 (%d) : lat (Yj) (from 1 to nlat)
+   Column 4 (%f) : u value
+   Column 5 (%f) : v value
+*/
+
+if (argv!=2) {
+  fprintf(stderr,"(read_bck) Error : read_bck with %d (1 nedded)\n",argv-1);
+  fprintf(stderr,"No background read\n");
+  return;
+ }
+ FILE *fid;
+ fid=fopen(argc[1],"r");
+ if (fid==NULL) {
+   fprintf(stderr,"(read_bck) Error : unable to open %s\n",argc[1]);
+   return;
+ }
+ 
+ int it,ilon,ilat,count;
+ count=0;
+ YREAL uval,vval;
+ int status=5;
+ while (status==5) {
+   status=fscanf(fid,"%d %d %d %lf %lf",&it,&ilon,&ilat,&uval,&vval);
+   if (status!=5) break;
+   if (ilon>nlon || ilon<=0 || ilat>nlat || ilat<=0 || it<=0 || it>jptfl) {
+     fprintf(stderr,"(read_bck) Error : bad coordinate values in %s (ilon=%d, ilat=%d, it=%d)\n",argc[1],ilon,ilat,it);
+     fclose(fid);
+     return;
+   }
+   ubck[ilon-1][ilat-1][it-1]=uval;
+   vbck[ilon-1][ilat-1][it-1]=vval;
+   count++;
+ }
+ if (count%(nlon*nlat)!=0) 
+   fprintf(stderr,"(read_bck) Warning : only %d data was loaded(nlon*nlat=%d)\n",count,nlon*nlat);
+ fclose(fid);
+ 
+}
+ 
 void read_obs(int argv, char *argc[]) {
   /* Read observation in a ascii file

altifloat/src/floater_delta.d

@@ -234 +234 @@
 insert_fct arg save_output_uv
 insert_fct arg save_output_rfloat
+insert_fct arg read_bck