Changeset 232

Timestamp:

08/26/08 09:04:08

Author:

dan

Message:

--

Files:

• scout/trunk/cdl_fixed_point.pl (modified) (22 diffs)

scout/trunk/cdl_fixed_point.pl

@@ -1 +1 @@
-sub fixed_point() {
+
-
-  # top timestamp
@@ -842 +843 @@
-      for my $j (0..$#this_wind_gust_data) {
-        push @{$wind_gust[$i]{'data'}}, $this_wind_gust_data[$j];
+        print( "$this_wind_gust_data $j: $this_wind_gust_data[$j] wind_gust $i: $wind_gust[$i]{'data'} \n" );
-      }
-    }
@@ -943 +945 @@
-      
-      #DWR 4/5/2008
-
+
+
+
-      
-      # we need to loop through the attributes, so find out how many there are
@@ -972 +976 @@
-        #also canceling conversion if air_pressure > 900
-        #print "$this_air_pressure_data[$j]\n";
-
+
+
+
+
+
+
+
-print STATION_ID_SQLFILE "-- this_air_pressure_data[j]=".$this_air_pressure_data[$j]." missing_value_value=".$missing_value_value." Fill_value_value=".$Fill_value_value."\n";
-        push @{$air_pressure[$i]{'data'}},
@@ -1051 +1061 @@
-      $varget = NetCDF::varget($ncid, $sea_surface_eastward_current[$i]{'var_id'},
-        (0), $time_dim{'dim_size'}, \@this_sea_surface_eastward_current_data);
+
-      if ($varget < 0) {die "ABORT! Cannot get $sea_surface_eastward_current[$i]{'var_name'} data.\n";}
-      # get all the attributes for this variable
@@ -1072 +1083 @@
-
-      #DWR 4/5/2008
-$units_value
+'m_s-1'
-
-      # we need to loop through the attributes, so find out how many there are
@@ -1098 +1109 @@
-      # add a NULL where missing value
-      for my $j (0..$#this_sea_surface_eastward_current_data) {
-
-
+    
+    
-        push @{$sea_surface_eastward_current[$i]{'data'}},
-          ($this_sea_surface_eastward_current_data[$j] * $this_slope + $this_intercept);
@@ -1491 +1502 @@
-                $this_val = sprintf("%.2f",$water_level[$i]{'data'}[$j]);
-                
-water_level[$i]{'reference'}
+this_val
-                
-                if (length $water_level[$i]{'reference_to_mllw'} > 0) {
@@ -1499 +1510 @@
-                  $this_val = 'NULL';
-                }
-
-
+DataVal + $
+
-                if (length $water_level[$i]{'reference_to_msl'} > 0) {
-                  $this_val = sprintf("%.2f",$water_level[$i]{'reference_to_msl'});
@@ -1517 +1528 @@
-                #DWR 4/5/2008
-                $NavD88DataVal = $this_val;
+=cut                
-            }
-            obsKMLSubRoutines::KMLAddObsToHash( 'water_level', 
@@ -1526 +1538 @@
-                                              'm(MLLW)',
-                                              $rObsHash );
+=comment
-            obsKMLSubRoutines::KMLAddObsToHash( 'water_level', 
-                                              $KMLTimeStamp[$j],
@@ -1542 +1555 @@
-                                              'm(NAVD88)',
-                                              $rObsHash );
+=cut                                              
-          }                                                      
-        }
@@ -2008 +2022 @@
-
-              if ($wind_gust[$j]{'data'}[$k] == $missing_value_value
-
-
+
+
+
-                print WIND_SQLFILE ',NULL';
-              }
@@ -2089 +2104 @@
-
-                  if ($wind_gust[$j]{'data'}[$k] != $missing_value_value
-
-&& $wind_gust[$j]{'data'}[$k] != '') 
+ )
+#&& $wind_gust[$j]{'data'}[$k] != '') DWR v1.1.2.0
-                  {
-                    $WindGstDataVal = sprintf("%.2f",$wind_gust[$j]{'data'}[$k]);
@@ -2372 +2387 @@
-            }           
-            my $strUnits;
+            #DWR v1.1.1.0
+            #Added UnitsStringConversionfucntion.
+            $strUnits = obsKMLSubRoutines::UnitsStringConversion( $salinity[$i]{'units'}, $XMLControlFile ); 
+            
-            if( lc( $salinity[$i]{'units'} ) eq 'ppt' )
-            {
@@ -2510 +2529 @@
-              my $ECurrentDataVal = 'NULL';
-              my $NCurrentDataVal = 'NULL';
-
+
+
+
+
+
+
+
+
+
+
+
-              if( $this_time_stamp_sec > $oldest_ok_timestamp )
-              {             
@@ -2524 +2553 @@
-                  {
-                    $Height = sprintf("%.2f",$sea_surface_eastward_current[$i]{'height'});                
-
+ 
-                  $ECurrentDataVal = sprintf("%.2f",$sea_surface_eastward_current[$i]{'data'}[$k]);
-                  
-                  $NCurrentDataVal = sprintf("%.2f",$sea_surface_northward_current[$j]{'data'}[$k]);
-
+
-                }
-
-
-
-
-
-
-
-
-
-
+
-                                                    $KMLTimeStamp[$k],
-$ECurrentDataVal
+@VectorDataVal[1]
-                                                    1,
-                                                    $strPlatformID,
-                                                    $Height,
-
+
+
+
+
+
+
+
+
+
+
+
-                                                    $rObsHash );
-                $strUnits = '';
@@ -2631 +2661 @@
-            print WAVE_SQLFILE ','.'\''.$height_var{'positive'}.'\'';
-            if ($significant_wave_height[$i]{'data'}[$j] == $missing_value_value
-
-
+
+
+
-              print WAVE_SQLFILE ',NULL';
-            }
@@ -2642 +2673 @@
-            }
-            if ($dominant_wave_period[$i]{'data'}[$j] == $missing_value_value
-
-
+
+
+
-              print WAVE_SQLFILE ',NULL';
-            }
@@ -2689 +2721 @@
-              }
-              if ($significant_wave_height[$i]{'data'}[$j] != $missing_value_value
-
-&& $significant_wave_height[$i]{'data'}[$j] != '') 
+ )
+#&& $significant_wave_height[$i]{'data'}[$j] != '')  DWR v1.1.2.0
-              {
-                $SigWaveHgtDataVal = sprintf("%.2f",$significant_wave_height[$i]{'data'}[$j]);
-              }
-              if ($dominant_wave_period[$i]{'data'}[$j] != $missing_value_value
-
-&& $dominant_wave_period[$i]{'data'}[$j] != '') 
+ )
+#&& $dominant_wave_period[$i]{'data'}[$j] != '') DWR v1.1.2.0
-              {
-                $DomWaveDataVal = sprintf("%.2f",$dominant_wave_period[$i]{'data'}[$j]);
@@ -2749 +2781 @@
-      
-      my $strXMLPath; 
-Y-%m-%dT%H-%M-
+H%M
-      #chomp( $strDate );      
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-      print( "XMLFilePath: $strXMLPath\n" );
-      obsKMLSubRoutines::BuildKMLFile( \%ObsHash, $strXMLPath );
@@ -2767 +2816 @@
-  }
-}
+
+sub get_mag_and_dir {
+
+my ($x, $y, $scale) = @_;
+my ($mag,$angle);
+
+$mag = sprintf("%d",$scale*sqrt($x*$x+$y*$y));
+#print "$mag\n";
+
+$angle = atan2($y,$x);
+$angle = sprintf("%d",180/3.1416*$angle);
+$angle = 90 - $angle;
+#only return positive degrees
+if ($angle < 0) 
+{ 
+  $angle = 360 + $angle; 
+}
+
+my @result = ($mag, $angle);
+
+return (@result);
+
+}