mankoff / SICE

Sentinel 3 Snow Products - GEUS
http://snow.geus.dk/
5 stars 5 forks source link

Update pySICEv1.1 #26

Closed BaptisteVandecrux closed 4 years ago

BaptisteVandecrux commented 4 years ago
mankoff commented 4 years ago

When I run sice.py I see divide by 0 errors:

+ parallel --verbose --lb -j 5 'python ./sice.py ./out/2018-08-15/{}' ::: 20180815T110728 20180815T124827 20180815T125127 20180815T125427 20180815T142926 20180815T143226 20180815T161026 20180815T175125
python ./sice.py ./out/2018-08-15/20180815T110728
python ./sice.py ./out/2018-08-15/20180815T124827
python ./sice.py ./out/2018-08-15/20180815T125127
python ./sice.py ./out/2018-08-15/20180815T125427
python ./sice.py ./out/2018-08-15/20180815T142926
End SICE.py ./out/2018-08-15/20180815T110728/ --- 20.892544444 seconds ---
python ./sice.py ./out/2018-08-15/20180815T143226
End SICE.py ./out/2018-08-15/20180815T125427/ --- 188.480614113 seconds ---
/home/kdm/projects/SICE/sice_lib.py:596: RuntimeWarning: divide by zero encountered in log
  p1=np.log(alb_sph[0,:,:])*np.log(alb_sph[0,:,:])
/home/kdm/projects/SICE/sice_lib.py:597: RuntimeWarning: divide by zero encountered in log
  p2=np.log(alb_sph[1,:,:])*np.log(alb_sph[1,:,:])
python ./sice.py ./out/2018-08-15/20180815T161026
End SICE.py ./out/2018-08-15/20180815T143226/ --- 632.308964823 seconds ---
python ./sice.py ./out/2018-08-15/20180815T175125
End SICE.py ./out/2018-08-15/20180815T175125/ --- 113.20127937 seconds ---
End SICE.py ./out/2018-08-15/20180815T161026/ --- 890.873333473 seconds ---
End SICE.py ./out/2018-08-15/20180815T124827/ --- 1362.146538961 seconds ---
/home/kdm/projects/SICE/sice_lib.py:775: RuntimeWarning: divide by zero encountered in true_divide
  rati=r7/r8
/home/kdm/projects/SICE/sice_lib.py:777: RuntimeWarning: divide by zero encountered in true_divide
  an=1./alasta
/home/kdm/projects/SICE/sice_lib.py:778: RuntimeWarning: divide by zero encountered in true_divide
  p   = r7 * np.exp(alam7/alasta)