Changeset 122 for altifloat/src/floater.h

Timestamp:

06/23/14 17:55:10 (10 years ago)

Author:

jbrlod

Message:

correct a bug in spaceorder / add save_output_rfloat

File:

• altifloat/src/floater.h (modified) (3 diffs)

altifloat/src/floater.h

@@ -69 +69 @@
       pjobs[i][j]=0;
       piret[i][j]=0;
-      pjret[i][j]=0;
+      pjret[i][j]=0; 
       pmask[i][j]=0;
       }
@@ -248 +248 @@
 
 void save_uv(int iti, int itn) {
-  /* save u,v fields in umod,vmod */
+  /* save u,v fields in umod,vmod 
+   and save piret,pjret*/
   int j,k,it;
-  for (it=iti;it<itn;it++)
+  int Yifloat=0;
+  for (it=iti;it<itn;it++) {
+    Yifloat=0;
     for (j=0;j<nlon;j++)
       for (k=0;k<nlat;k++) {
@@ -260 +263 @@
 #endif
       }
-}
+    for (j=0;j<jpnfl;j++) 
+      if (pmask[iti][j]==1 && is_activ(j,iti)>0 && it>iti) {
+        pjret[it][j]=YS1_r_float(Yifloat,it-iti);
+        piret[it][j]=YS2_r_float(Yifloat,it-iti);
+        Yifloat++;
+      }
+  }
+}
+
+void save_output_rfloat (int argc, char *argv[]) {
+  FILE *fid;
+  fid=fopen(argv[1],"w");
+  if (fid==NULL) {
+    printf("\nfailed to open %s",argv[1]);
+    exit(3);
+  }
+  int j,it;
+  for (j=0;j<jpnfl;j++)
+    for (it=0;it<jptfl;it++)
+      
+      fprintf(fid,"%d %d %f %f\n",it,j,piret[it][j],pjret[it][j]);
+        
+                                             
+  fclose(fid);
+
+}
+
 void save_output_uv (int argc, char *argv[]) {
   FILE *fid;