Re: issue with dividing the integer by 0

After

   stddev_var_ob = stddev(var_obs)

add

   stddev_var_obs = where(stddev_var_obs.eq.0.0 \
                          ,stddev_var_obs@_FillValue, stddev_var_obs)

On 4/13/12 3:55 PM, Alexander Semenov wrote:
> Hi again,
>
> I'm calculating the normalized root mean square error, and the code
> crashes, because for slp for the particular station equals 0. In my code
> I tried to fix it, but it didn't work. Please have a look and give any
> suggestions on how I can make it work. Thanks.
>
>
> Here's what comes as the program was running:
>
> pacman 1991% ncl rmse-normalized.ncl
>
> Copyright (C) 1995-2011 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 6.0.0
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> (0)------ SELAWIK------
> (0)Calculating for Temperature
> (0)Std dev: 8.30892
> (0)::: 0.825233
> (0)Calculating for Dewpoint
> (0)Std dev: 213.946
> (0)::: 1.01504
> (0)Calculating for WindSpeed
> (0)Std dev: 3.18798
> (0)::: 1.03146
> (0)Calculating for WindDir
> (0)Std dev: 206.94
> (0)::: 0.755161
> (0)Calculating for Sea-Level-Pres
> (0)Std dev: 0
> (0)::: -999
> At line 100 of file writematrix.f
> Fortran runtime error: No such file or directory
>
>
> --
> regards
>
> *******************************************************
> *Alexander Semenov*
> *
> *
> PhD Student - Research Assistant
> International Arctic Research Center
> Department of atmospheric sciences
> University of Alaska Fairbanks
>
> 930 Koyukuk dr, 408c3
> Fairbanks, AK, USA, 99775
> work phone +19074742672
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Apr 16 08:26:03 2012