Changeset 160 for altifloat/src

Timestamp:

07/20/15 11:51:06 (9 years ago)

Author:

jbrlod

Message:

finalize wind effect

Location:

altifloat/src

Files:

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

altifloat/src/floater.h

@@ -25 +25 @@
 YREAL uwind[nlon][nlat][jptfl];
 YREAL vwind[nlon][nlat][jptfl];
+YREAL utotal[nlon][nlat][jptfl];
+YREAL vtotal[nlon][nlat][jptfl];
 
 //Background
@@ -88 +90 @@
         ubck[i][j][k]=0;
         vbck[i][j][k]=0;
+        utotal[i][j][k]=0;
+        vtotal[i][j][k]=0;
       }
 }
@@ -277 +281 @@
   
   int time[Ntime];
-  float u10[Ntime][nlat][nlon];
-  float v10[Ntime][nlat][nlon];
-  
+  float (*u10)[nlat][nlon] = new float[Ntime][nlat][nlon];
+  float (*v10)[nlat][nlon] = new float[Ntime][nlat][nlon];
+
+  //  float u10[Ntime][Nlat][Nlon];
+  //float v10[Ntime][Nlat][Nlon];
+ 
+  if ((retval = nc_get_var_float(ncid,u10id,&u10[0][0][0])))
+    ERR(retval);
+  
+  if ((retval = nc_get_var_float(ncid,v10id,&v10[0][0][0])))
+    ERR(retval);
+
+ 
   if ((retval = nc_get_var_int(ncid,timeid,&time[0])))
     ERR(retval);
   
-  if ((retval = nc_get_var_float(ncid,u10id,&u10[0][0][0])))
-    ERR(retval);
-  
-  if ((retval = nc_get_var_float(ncid,v10id,&v10[0][0][0])))
-    ERR(retval);
 
   //Load time dimensions
@@ -500 +509 @@
 void update_uv(){
   int j,k,t;
-  for (t=0;t<jtlag+1;t++)
+  for (t=0;t<jtlag;t++)
     for (j=0;j<nlon;j++)
       for (k=0;k<nlat;k++) {
-        YS1_u(j,k,t)+=YS1_u_d(j,k);
-        YS1_v(j,k,t)+=YS1_v_d(j,k);
+        YS1_u(j,k,t+1)+=YS1_u_d(j,k);
+        YS1_v(j,k,t+1)+=YS1_v_d(j,k);
         //printf("u=%f delu=%f \n",YS1_u(j,k),YS1_u_d(j,k));    
 
@@ -564 +573 @@
     for (j=0;j<nlon;j++)
       for (k=0;k<nlat;k++) {
-        umod[j][k][it]=YS1_u(j,k,it-iti);
-        vmod[j][k][it]=YS1_v(j,k,it-iti);
+        umod[j][k][it]=YS1_u(j,k,it-iti+1);
+        vmod[j][k][it]=YS1_v(j,k,it-iti+1);
+        utotal[j][k][it]=YS1_utot(j,k,it-iti+1);
+        vtotal[j][k][it]=YS1_vtot(j,k,it-iti+1);
 #ifdef UPDATE_BCK
-        ubck[j][k][it]=YS1_u(j,k,it-iti);
-        vbck[j][k][it]=YS1_v(j,k,it-iti);
+        ubck[j][k][it]=YS1_u(j,k,it-iti+1);
+        vbck[j][k][it]=YS1_v(j,k,it-iti+1);
 #endif
       }
@@ -630 +641 @@
         for (k=0;k<nlat;k++) {
           fprintf(fid,"%d %d %d %f %f\n",it,j,k,umod[j][k][it],vmod[j][k][it]);
+        }
+                                             
+  fclose(fid);
+
+}
+
+void save_output_uvtot (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,k,it;
+    for (it=0;it<jptfl;it++)
+      for (j=0;j<nlon;j++)
+        for (k=0;k<nlat;k++) {
+          fprintf(fid,"%d %d %d %f %f\n",it,j,k,utotal[j][k][it],vtotal[j][k][it]);
         }
                                              

altifloat/src/floater_delta.d

@@ -317 +317 @@
 insert_fct     myforward
 insert_fct arg load_wind
+insert_fct arg save_output_uvtot