code: def rms(array):
      
      #;time should be the last dimension
        
  #    IF sz[ndims-1] EQ 1 THEN RETURN, 0.
      
      nans = np.where(~np.isnan(array))
      cnt=len(nans)
      
      if (cnt == 0) : return np.NaN
      
      nval=np.nansum(~np.isnan(array))
      
      mn = np.nansum(array) / nval
  
      return np.sqrt(np.nansum((array - mn)**2.0)/nval)
docs: ;+
      ; RMS : Returns root-mean-squared deviation of an array
      ; 
      ; @param array {in}{required}{type:NUMERIC} 2-dimensionnal array. Time dimension<br />
      ;   should be the last dimension.
      ;   
      ; @returns RMS array (1 value per time serie)
      ; 
      ; @author Renaud DUSSURGET, LEGOS/CTOH
      ;-