udated output format

MenakaRevel · Dec 5, 2022 · 289eb19 · 289eb19
1 parent 4914c60
commit 289eb19
Show file tree

Hide file tree

Showing 69 changed files with 93 additions and 57 deletions.
diff --git a/img_code/bin/._Makefile → ESSD_paper/bin/._Makefile b/img_code/bin/._Makefile → ESSD_paper/bin/._Makefile
diff --git a/img_code/bin/Makefile → ESSD_paper/bin/Makefile b/img_code/bin/Makefile → ESSD_paper/bin/Makefile
diff --git a/img_code/bin/print_grid → ESSD_paper/bin/print_grid b/img_code/bin/print_grid → ESSD_paper/bin/print_grid
diff --git a/img_code/bin/print_grid.f90 → ESSD_paper/bin/print_grid.f90 b/img_code/bin/print_grid.f90 → ESSD_paper/bin/print_grid.f90
diff --git a/img_code/bin/print_grid.o → ESSD_paper/bin/print_grid.o b/img_code/bin/print_grid.o → ESSD_paper/bin/print_grid.o
diff --git a/img_code/bin/print_rivvec → ESSD_paper/bin/print_rivvec b/img_code/bin/print_rivvec → ESSD_paper/bin/print_rivvec
diff --git a/img_code/bin/print_rivvec.f90 → ESSD_paper/bin/print_rivvec.f90 b/img_code/bin/print_rivvec.f90 → ESSD_paper/bin/print_rivvec.f90
diff --git a/img_code/bin/print_rivvec.o → ESSD_paper/bin/print_rivvec.o b/img_code/bin/print_rivvec.o → ESSD_paper/bin/print_rivvec.o
diff --git a/img_code/bin/txt_vector → ESSD_paper/bin/txt_vector b/img_code/bin/txt_vector → ESSD_paper/bin/txt_vector
diff --git a/img_code/bin/txt_vector.f90 → ESSD_paper/bin/txt_vector.f90 b/img_code/bin/txt_vector.f90 → ESSD_paper/bin/txt_vector.f90
diff --git a/img_code/bin/txt_vector.o → ESSD_paper/bin/txt_vector.o b/img_code/bin/txt_vector.o → ESSD_paper/bin/txt_vector.o
diff --git a/img_code/f02-along_river_AltiVS.py → ESSD_paper/f02-along_river_AltiVS.py b/img_code/f02-along_river_AltiVS.py → ESSD_paper/f02-along_river_AltiVS.py
diff --git a/img_code/f03-satellite_schematic.py → ESSD_paper/f03-satellite_schematic.py b/img_code/f03-satellite_schematic.py → ESSD_paper/f03-satellite_schematic.py
diff --git a/img_code/f04-flag_map_uparea.py → ESSD_paper/f04-flag_map_uparea.py b/img_code/f04-flag_map_uparea.py → ESSD_paper/f04-flag_map_uparea.py
diff --git a/img_code/f05-unrealistic_map_dist.py → ESSD_paper/f05-unrealistic_map_dist.py b/img_code/f05-unrealistic_map_dist.py → ESSD_paper/f05-unrealistic_map_dist.py
diff --git a/...ode/f06-boxplot_RMSE_disttomouth_slope.py → ...per/f06-boxplot_RMSE_disttomouth_slope.py b/...ode/f06-boxplot_RMSE_disttomouth_slope.py → ...per/f06-boxplot_RMSE_disttomouth_slope.py
diff --git a/img_code/f07-comapre_expert_ordinary.py → ESSD_paper/f07-comapre_expert_ordinary.py b/img_code/f07-comapre_expert_ordinary.py → ESSD_paper/f07-comapre_expert_ordinary.py
@@ -37,7 +37,7 @@
 #
 # from read_patchMS import upstream
 #==============================
-CaMa_dir = "/cluster/data6/menaka/CaMa-Flood_v396a_20200514"
+CaMa_dir = "/cluster/data6/menaka/CaMa-Flood_v4"
 mapname= "glb_06min"
 #==============================
 #=================================

diff --git a/img_code/fe05-uparea_flag.py → ESSD_paper/fe05-uparea_flag.py b/img_code/fe05-uparea_flag.py → ESSD_paper/fe05-uparea_flag.py
diff --git a/img_code/fe09-boxplot_expert_ordinary.py → ESSD_paper/fe09-boxplot_expert_ordinary.py b/img_code/fe09-boxplot_expert_ordinary.py → ESSD_paper/fe09-boxplot_expert_ordinary.py
diff --git a/img_code/mk_fig_mool.sh → ESSD_paper/mk_fig_mool.sh b/img_code/mk_fig_mool.sh → ESSD_paper/mk_fig_mool.sh
@@ -24,7 +24,7 @@ export OMP_NUM_THREADS=$NCPUS
 
 # got to working dirctory
 # cd $PBS_O_WORKDIR
-cd "/cluster/data6/menaka/AltiMaP/img_code"
+cd "/cluster/data6/menaka/AltiMaP/ESSD_paper"
 
 
 mkdir -p fig

diff --git a/img_code/my_colorbar.py → ESSD_paper/my_colorbar.py b/img_code/my_colorbar.py → ESSD_paper/my_colorbar.py
diff --git a/README.md b/README.md
@@ -1,3 +1,29 @@
+### Dataset Varaibles
+# VS metadata
+ID: Identification number of VS
+station         :VS name
+dataname        :dataset name
+lon 	        :longitude [degrees east]
+lat	            :latitude [degrees north]
+satellite	    :name of the satellite
+# MERIT Hydro-related	
+flag	        :allocation flag
+elevation	    :elevation at VS location on MERIT Hydro [m]
+dist_to_mouth	:distance to the unit-catchment mouth [km]
+kx1	            :best x-coordinate with respect to the 10°×10° higher resolution tile
+ky1	            :best y-coordinate with respect to the 10°×10° higher resolution tile
+kx2	            :second-best option of x-coordinate with respect to the 10°×10° high-resolution tile
+ky2	            :second-best option of y-coordinate with respect to the 10°×10° high-resolution tile
+dist1	        :distance from the second-best location to the VS [km]
+dist2	        :distance from the second-best location to the VS [km]
+rivwth	        :River width of the allocated location [m]
+# Coarse-resolution river network-related	
+ix	            :x-coordinate with respect to coarse resolution
+iy	            :y-coordinate with respect to coarse resolution
+EGM08	        :EGM 2008 datum elevation[m]
+EGM96	        :EGM 1996 datum elevation[m]
+
+
 ### Altimetry Virtual Stations on CaMa-Flood River Network
 
 1. Allocate virtual stations to the CaMa-Flood map

diff --git a/f01-uparea_criteria.py → img_code/f01-uparea_criteria.py b/f01-uparea_criteria.py → img_code/f01-uparea_criteria.py
diff --git a/f02-slope_criteria.py → img_code/f02-slope_criteria.py b/f02-slope_criteria.py → img_code/f02-slope_criteria.py
diff --git a/f03-flag_criteria.py → img_code/f03-flag_criteria.py b/f03-flag_criteria.py → img_code/f03-flag_criteria.py
diff --git a/f04-ele_diff_hist.py → img_code/f04-ele_diff_hist.py b/f04-ele_diff_hist.py → img_code/f04-ele_diff_hist.py
diff --git a/f05-ele_diff_braided.py → img_code/f05-ele_diff_braided.py b/f05-ele_diff_braided.py → img_code/f05-ele_diff_braided.py
diff --git a/f06-ele_diff_map.py → img_code/f06-ele_diff_map.py b/f06-ele_diff_map.py → img_code/f06-ele_diff_map.py
diff --git a/f07-ele_diff_elevation_uparea.py → img_code/f07-ele_diff_elevation_uparea.py b/f07-ele_diff_elevation_uparea.py → img_code/f07-ele_diff_elevation_uparea.py
diff --git a/f09-flag_map.py → img_code/f09-flag_map.py b/f09-flag_map.py → img_code/f09-flag_map.py
diff --git a/f10-high_resolution_map.py → img_code/f10-high_resolution_map.py b/f10-high_resolution_map.py → img_code/f10-high_resolution_map.py
diff --git a/f11-down_dist_map.py → img_code/f11-down_dist_map.py b/f11-down_dist_map.py → img_code/f11-down_dist_map.py
diff --git a/f12-compare_cmf_data.py → img_code/f12-compare_cmf_data.py b/f12-compare_cmf_data.py → img_code/f12-compare_cmf_data.py
diff --git a/f13-boxplot_bias_cmf_obs.py → img_code/f13-boxplot_bias_cmf_obs.py b/f13-boxplot_bias_cmf_obs.py → img_code/f13-boxplot_bias_cmf_obs.py
diff --git a/f14-boxplot_amplitude_diff_cmf_obs.py → ...ode/f14-boxplot_amplitude_diff_cmf_obs.py b/f14-boxplot_amplitude_diff_cmf_obs.py → ...ode/f14-boxplot_amplitude_diff_cmf_obs.py
diff --git a/f15-scatter_mean_bias_vs_ele_diff.py → ...code/f15-scatter_mean_bias_vs_ele_diff.py b/f15-scatter_mean_bias_vs_ele_diff.py → ...code/f15-scatter_mean_bias_vs_ele_diff.py
diff --git a/f16-down_dist_hist.py → img_code/f16-down_dist_hist.py b/f16-down_dist_hist.py → img_code/f16-down_dist_hist.py
diff --git a/f17-dist_to_mouth_vs_mean_bias.py → img_code/f17-dist_to_mouth_vs_mean_bias.py b/f17-dist_to_mouth_vs_mean_bias.py → img_code/f17-dist_to_mouth_vs_mean_bias.py
diff --git a/f18-boxplots_bias_height_distom.py → img_code/f18-boxplots_bias_height_distom.py b/f18-boxplots_bias_height_distom.py → img_code/f18-boxplots_bias_height_distom.py
diff --git a/f19-boxplot_RMSE_elediff_distom.py → img_code/f19-boxplot_RMSE_elediff_distom.py b/f19-boxplot_RMSE_elediff_distom.py → img_code/f19-boxplot_RMSE_elediff_distom.py
diff --git a/f20-boxplot_RMSE_cmf_obs.py → img_code/f20-boxplot_RMSE_cmf_obs.py b/f20-boxplot_RMSE_cmf_obs.py → img_code/f20-boxplot_RMSE_cmf_obs.py
diff --git a/f21-scatter_elediff_disttom_RMSE.py → img_code/f21-scatter_elediff_disttom_RMSE.py b/f21-scatter_elediff_disttom_RMSE.py → img_code/f21-scatter_elediff_disttom_RMSE.py
diff --git a/f22-compare_cmf_interplotaion_data.py → ...ode/f22-compare_cmf_interplotaion_data.py b/f22-compare_cmf_interplotaion_data.py → ...ode/f22-compare_cmf_interplotaion_data.py
diff --git a/f23-compare_cmf_elediff.py → img_code/f23-compare_cmf_elediff.py b/f23-compare_cmf_elediff.py → img_code/f23-compare_cmf_elediff.py
diff --git a/f24-compare_cmf_all.py → img_code/f24-compare_cmf_all.py b/f24-compare_cmf_all.py → img_code/f24-compare_cmf_all.py
diff --git a/f25-compare_all_map.py → img_code/f25-compare_all_map.py b/f25-compare_all_map.py → img_code/f25-compare_all_map.py
diff --git a/f26-higher_resolution_AltiVS.py → img_code/f26-higher_resolution_AltiVS.py b/f26-higher_resolution_AltiVS.py → img_code/f26-higher_resolution_AltiVS.py
diff --git a/f27-along_river_AltiVS.py → img_code/f27-along_river_AltiVS.py b/f27-along_river_AltiVS.py → img_code/f27-along_river_AltiVS.py
diff --git a/f28-high-res_river_profile.py → img_code/f28-high-res_river_profile.py b/f28-high-res_river_profile.py → img_code/f28-high-res_river_profile.py
diff --git a/f29-unrealstic_map.py → img_code/f29-unrealstic_map.py b/f29-unrealstic_map.py → img_code/f29-unrealstic_map.py
diff --git a/f30-allocation_flag_map.py → img_code/f30-allocation_flag_map.py b/f30-allocation_flag_map.py → img_code/f30-allocation_flag_map.py
diff --git a/f31-high-res_map.py → img_code/f31-high-res_map.py b/f31-high-res_map.py → img_code/f31-high-res_map.py
diff --git a/f32-alloc_high-res.py → img_code/f32-alloc_high-res.py b/f32-alloc_high-res.py → img_code/f32-alloc_high-res.py
diff --git a/f33-flag_map_example.py → img_code/f33-flag_map_example.py b/f33-flag_map_example.py → img_code/f33-flag_map_example.py
diff --git a/f34-RMSE_flags.py → img_code/f34-RMSE_flags.py b/f34-RMSE_flags.py → img_code/f34-RMSE_flags.py
diff --git a/f35-unreal_map_flag.py → img_code/f35-unreal_map_flag.py b/f35-unreal_map_flag.py → img_code/f35-unreal_map_flag.py
diff --git a/f36-best_alloc_vs_near_alloc.py → img_code/f36-best_alloc_vs_near_alloc.py b/f36-best_alloc_vs_near_alloc.py → img_code/f36-best_alloc_vs_near_alloc.py
diff --git a/f37-pdf_disttomouth.py → img_code/f37-pdf_disttomouth.py b/f37-pdf_disttomouth.py → img_code/f37-pdf_disttomouth.py
diff --git a/f38-pdf_elevation_diff.py → img_code/f38-pdf_elevation_diff.py b/f38-pdf_elevation_diff.py → img_code/f38-pdf_elevation_diff.py
diff --git a/f39-boxplot_RMSE_slope.py → img_code/f39-boxplot_RMSE_slope.py b/f39-boxplot_RMSE_slope.py → img_code/f39-boxplot_RMSE_slope.py
diff --git a/f40-uparea_flag.py → img_code/f40-uparea_flag.py b/f40-uparea_flag.py → img_code/f40-uparea_flag.py
diff --git a/f41-sensitivity_unreal.py → img_code/f41-sensitivity_unreal.py b/f41-sensitivity_unreal.py → img_code/f41-sensitivity_unreal.py
diff --git a/f42-unrealVS_rivwth.py → img_code/f42-unrealVS_rivwth.py b/f42-unrealVS_rivwth.py → img_code/f42-unrealVS_rivwth.py
diff --git a/f43-AltiMap_Ordinary.py → img_code/f43-AltiMap_Ordinary.py b/f43-AltiMap_Ordinary.py → img_code/f43-AltiMap_Ordinary.py
diff --git a/f44-unrealVS_bias_geoid.py → img_code/f44-unrealVS_bias_geoid.py b/f44-unrealVS_bias_geoid.py → img_code/f44-unrealVS_bias_geoid.py
diff --git a/s01-allocate_VS.sh b/s01-allocate_VS.sh
@@ -1,7 +1,7 @@
 #! /bin/bash
 
 ### SET "mool PBS" @ IIS U-Tokyo
-#PBS -q F40
+#PBS -q E40
 #PBS -l select=1:ncpus=40:mem=100gb
 #PBS -l place=scatter
 #PBS -j oe
@@ -17,7 +17,7 @@ export OMP_NUM_THREADS=$NCPUS
 
 # got to working dirctory
 # cd $PBS_O_WORKDIR
-cd "/cluster/data6/menaka/Altimetry"
+cd "/cluster/data6/menaka/AltiMaP"
 
 #CaMA-Flood directory
 # CaMa_dir="/cluster/data6/menaka/CaMa-Flood_v396a_20200514"
@@ -50,11 +50,12 @@ fi
 
 echo "starting calculations........"
 # echo "            ID                                      station            dataname         lon       lat       ix      iy     ele_diff     EGM08     EGM96        satellite" > tmp.txt
-printf '%30s%67s%12s%10s%10s%8s%12s%10s%17s%6s%15s%12s%10s%8s%8s%8s%12s%12s%10s%8s%12s%10s\n' ID station dataname lon lat ix iy elevation EGM08 EGM96 flag satellite dist_to_mouth kx1 ky1 kx2 ky2 dist1 dist2 rivwth > tmp.txt
+# printf '%30s%67s%12s%10s%10s%8s%12s%10s%17s%6s%15s%12s%10s%8s%8s%8s%12s%12s%10s%8s%12s%10s\n' ID station dataname lon lat ix iy elevation EGM08 EGM96 flag satellite dist_to_mouth kx1 ky1 kx2 ky2 dist1 dist2 rivwth > tmp.txt
+printf '%13s%64s%12s%12s%10s%17s%6s%12s%15s%10s%8s%8s%8s%14s%12s%12s%10s%8s%12s%10s\n' ID station dataname lon lat satellite flag elevation dist_to_mouth kx1 ky1 kx2 ky2 dist1 dist2 rivwth ix iy EGM08 EGM96 > tmp.txt
 SOUTH=-60
 while [ $SOUTH -lt 90 ];
 do
-  WEST=-180
+  WEST=-80 #-180
   while [ $WEST -lt 180 ];
   do
     CNAME=`./src/set_name $WEST $SOUTH`

diff --git a/s03-unreal_obs.sh b/s03-unreal_obs.sh
@@ -24,7 +24,7 @@ export OMP_NUM_THREADS=$NCPUS
 
 # got to working dirctory
 # cd $PBS_O_WORKDIR
-cd "/cluster/data6/menaka/Altimetry"
+cd "/cluster/data6/menaka/AltiMaP"
 
 #CaMA-Flood directory
 # CaMa_dir="/cluster/data6/menaka/CaMa-Flood_v396a_20200514"
@@ -42,12 +42,14 @@ map="glb_06min"
 dataname="HydroWeb"
 
 # observation list
-obstxt="/cluster/data6/menaka/Altimetry/out/altimetry_"$map"_20220725.txt"
+obstxt="/cluster/data6/menaka/AltiMaP/out/altimetry_"$map"_20221205.txt"
 
 # out dir
 outdir="./out"
 mkdir -p $outdir
 
+outname="biased_removed_altimetry_"$map"_"${day}".txt"
+
 # threshold for finding outliers
 threshold=10.0 #m
 
@@ -57,7 +59,10 @@ threshold=10.0 #m
 # mkdir -p $figdir
 day=$(date +"%Y%m%d")
 
-python src/unreal_obs.py $dataname $map $CaMa_dir $TAG $obstxt $threshold > "$outdir/unreal_obs_${day}.txt" #& #> /dev/null 2>&1 & 
+printf '%13s%64s%12s%12s%10s%17s%6s%12s%15s%10s%8s%8s%8s%14s%12s%12s%10s%8s%12s%10s\n' ID station dataname lon lat satellite flag elevation dist_to_mouth kx1 ky1 kx2 ky2 dist1 dist2 rivwth ix iy EGM08 EGM96 > tmp.txt
+python src/unreal_obs.py $dataname $map $CaMa_dir $TAG $obstxt $threshold >> tmp.txt  #& #> /dev/null 2>&1 & 
+
+mv tmp.txt "$outdir/$outname"
 
 wait
 

diff --git a/src/allocate_VS.f90 b/src/allocate_VS.f90
@@ -526,11 +526,14 @@ program SET_MAP
     ! print*, trim(station), lon0, lat0, ix, iy, flag, dist1, dist2
     ! print*, "================================================"
     if (iXX > 0 .or. iYY > 0) then
-        print '(a30,2x,a65,2x,a10,2x,2f10.2,2x,2i8.0,2x,3f10.2,2x,a15,2x,f13.2,2x,i4.0,2x,4i8.0,2f12.2,2x,f13.2)',& 
-        &trim(adjustl(id)), trim(station), trim(dataname), lon0, lat0, iXX, iYY, ele1m(kx,ky),&
-        &egm08, egm96, trim(sat), diffdist*1e-3, flag, kx1, ky1, kx2, ky2, dist1, dist2, riv1m(kx,ky)       ! elevtn(iXX,iYY)-
-    ! else
-    !     print*, "no data"
+        ! print '(a30,2x,a65,2x,a10,2x,2f10.2,2x,2i8.0,2x,3f10.2,2x,a15,2x,f13.2,2x,i4.0,2x,4i8.0,2f12.2,2x,f13.2)',& 
+        ! &trim(adjustl(id)), trim(station), trim(dataname), lon0, lat0, iXX, iYY, ele1m(kx,ky),&
+        ! &egm08, egm96, trim(sat), diffdist*1e-3, flag, kx1, ky1, kx2, ky2, dist1, dist2, riv1m(kx,ky)       ! elevtn(iXX,iYY)-
+
+        print '(a13,4x,a60,2x,a10,2x,2f10.2,2x,a15,2x,i4.0,2x,f10.2,2x,f13.2,2x,4i8.0,2x,3f12.2,2x,2i8.0,2x,2f10.2)',& 
+        &trim(adjustl(id)), trim(station), trim(dataname), lon0, lat0, trim(sat),& 
+        &flag, ele1m(kx,ky),diffdist*1e-3, kx1, ky1, kx2, ky2, dist1, dist2, riv1m(kx,ky),&
+        &iXX, iYY, egm08, egm96
     end if
         ! write(27,'(a14,2x,a40,2x,a10,2x,2f10.2,2x,2i8.0,2x,3f10.2,2x,a15)') trim(adjustl(id)),& 
         ! &trim(station),trim(dataname), lon0, lat0,iXX, iYY,elevtn(iXX,iYY)-ele1m(kx,ky),&

diff --git a/src/unreal_obs.py b/src/unreal_obs.py
@@ -1,6 +1,8 @@
 #!/opt/local/bin/python
 # -*- coding: utf-8 -*-
-
+'''
+Located biased VS
+'''
 import datetime
 import numpy as np
 from numpy import ma
@@ -251,50 +253,49 @@ def meanHydroWeb(station,egm08=0.0,egm96=0.0):
     line2   = re.split("_",station)
     riv     = line2[1]
     stream  = line2[2]
+    dataname= line[2]
     lon     = float(line[3])
     lat     = float(line[4])
-    ix      = int(line[5])-1
-    iy      = int(line[6])-1
+    sat     = line[5]
+    #--
+    flag    = int(line[6])
     elev    = float(line[7])
-    EGM08   = float(line[8])
-    EGM96   = float(line[9])
-    sat     = line[10].strip()
-    dist    = float(line[11])
-    flag    = int(line[12])
-    kx      = int(line[13])
-    ky      = int(line[14])
-    # rivwth  = float(line[19])
-    #-----------------------
-    # print (riv,station,kx,ky)
-    # if riv != rivername0:
-    #     continue
-    # if stream != stream0:
-    #     continue
-    nums.append(num)
-    river.append(riv)
-    pname.append(station)
-    lons.append(lon)
-    lats.append(lat)
-    xlist.append(ix)
-    ylist.append(iy)
-    lelev.append(elev)
-    egm08.append(EGM08)
-    egm96.append(EGM96)
-    llsat.append(sat)
-    ldtom.append(dist)
-    lflag.append(flag)
-    kxlst.append(kx)
-    kylst.append(ky)
-    # print (riv,station)
-#=============================
-pnum=len(pname)
-#=====================================
-for point in np.arange(pnum):
-    meanW, stdW = meanHydroWeb(pname[point],egm96=egm96[point],egm08=egm08[point])
-    iXX = xlist[point]
-    iYY = ylist[point]
-    elev= lelev[point]
-    sat = llsat[point]
+    dist    = float(line[8])
+    kx1     = int(line[9])
+    ky1     = int(line[10])
+    kx2     = int(line[11])
+    ky2     = int(line[12])
+    dist1   = float(line[13])
+    dist2   = float(line[14])
+    rivwth  = float(line[15])
+    #--
+    ix      = int(line[16])
+    iy      = int(line[17])
+    EGM08   = float(line[18])
+    EGM96   = float(line[19])
+
+    # nums.append(num)
+    # river.append(riv)
+    # pname.append(station)
+    # lons.append(lon)
+    # lats.append(lat)
+    # xlist.append(ix)
+    # ylist.append(iy)
+    # lelev.append(elev)
+    # egm08.append(EGM08)
+    # egm96.append(EGM96)
+    # llsat.append(sat)
+    # ldtom.append(dist)
+    # lflag.append(flag)
+    # kxlst.append(kx)
+    # kylst.append(ky)
+    meanW, stdW = meanHydroWeb(station,egm96=EGM96,egm08=EGM08)
+    # iXX = xlist[point]
+    # iYY = ylist[point]
+    # elev= lelev[point]
+    # sat = llsat[point]
     # if meanW > elevtn[iYY,iXX] + upthr or meanW < elevtn[iYY,iXX] - dwthr:
     if meanW > elev + upthr or meanW < elev - dwthr:
-        print ("%69s%4d%12.3f%12.3f%12.3f%15s%12.3f")%(pname[point],lflag[point],elev,meanW,meanWSE_VICBC[iYY+dYY,iXX+dXX],sat,meanW - elev)
+        flag=flag+100
+        # print ("%69s%4d%12.3f%12.3f%12.3f%15s%12.3f")%(pname[point],lflag[point],elev,meanW,meanWSE_VICBC[iYY+dYY,iXX+dXX],sat,meanW - elev)
+    print ("%13s%64s%12s%12.2f%12.2f%17s%6d%12.2f%15.2f%10d%8d%8d%8d%14.2f%12.2f%12.2f%10d%8d%12.2f%10.2f")%(num,station,dataname,lon,lat,sat,flag,elev,dist,kx1,ky1,kx2,ky2,dist1,dist2,rivwth,ix,iy,EGM08,EGM96)