Re: write matrix.f: Fortran runtime error

Hi

I've solved the issue. The problem was not because of NCL but its because
of my own data.

Thanks!

On Thu, May 31, 2012 at 4:48 PM, Khairunnisa Yahya <kbyahya@ncsu.edu> wrote:

> Hi,
>
> The script I'm using is as follows:
>
> *#!/bin/csh*
> *set timeperiod = 2010050106.2010060105
> set metfile = sim.met.$timeperiod.nc*
> *rm -fv list.usaf list.wban list.row list.col*
> *sed '/USAF.WBAN/d' ../stns.processed/allstn.IJ | awk '{print $1}' | awk
> -F\- '{print $1}' > list.usaf
> sed '/USAF.WBAN/d' ../stns.processed/allstn.IJ | awk '{print $1}' | awk
> -F\- '{print $2}' > list.wban
> sed '/USAF.WBAN/d' ../stns.processed/allstn.IJ | awk '{print $2}' >
> list.row
> sed '/USAF.WBAN/d' ../stns.processed/allstn.IJ | awk '{print $3}' >
> list.col*
> *ncl << EOF
> f = addfile("$metfile","r")
> vars = (/"U10","V10","T2","Q2","PSFC","RAINC","RAINNC"/)
> nvar = dimsizes(vars)
> dat = f->T2
> dims = dimsizes(dat)
> ndims = dimsizes(dims)*
> *lusaf = asciiread("list.usaf", -1, "integer")
> lwban = asciiread("list.wban", -1, "integer")
> lrow = asciiread("list.row", -1, "integer")
> lcol = asciiread("list.col", -1, "integer")
> nsite = dimsizes(lusaf)*
> *do isite = 0, (nsite - 1)*
> *
> print(sprinti("%0.6i",lusaf(isite))+"-"+sprinti("%0.5i",lwban(isite))+"
> "+lrow(isite)+" "+lcol(isite))*
> * xout = new((/dims(0),nvar/), "float")
> do iv = 0, (nvar - 1)
> if (ndims .eq. 3) then
> xout(:,iv) = f->\$vars(iv)\$(:, (lrow(isite) - 1),(lcol(isite) -
> 1))
> end if
> if (ndims .eq. 4) then
> xout(:,iv) = f->\$vars(iv)\$(:,0,(lrow(isite) - 1),(lcol(isite) -
> 1))
> end if
> end do*
> * xopt = True
> **xopt@fout* <xopt@fout>* =
> "sim.met.$timeperiod.extracted/"+sprinti("%0.6i",lusaf(isite))+"-"+sprinti("%0.5i",lwban(isite))
> **xopt@title* <xopt@title>* = "U10 V10 T2 Q2 PSFC RAINC RAINNC"
> xfmt = "f8.2,1x,f8.2,1x,f10.2,1x,e12.4,1x,f10.2,1x,f9.3,1x,f9.3"*
> * write_matrix(xout, xfmt, xopt)*
> *end do*
> *EOF
> *
> Thanks!
>
>
> On Thu, May 31, 2012 at 4:30 PM, Dennis Shea <shea@ucar.edu> wrote:
>
>> You must provide more information than this.
>>
>> Do you have a clean, small script that illustrates the problem?
>>
>>
>> On 5/31/12 12:02 PM, Khairunnisa Yahya wrote:
>>
>>> Hi,
>>> I've recently switched systems from Redhat (3.4.5-2) to Redhat (4.1.2).
>>> Previously I used the NCL version 5.0 and it worked fine. However, when
>>> I switched systems I found out that my script does not work and it gives
>>> this error message:
>>> /At line 100 of file writematrix.f
>>>
>>> Fortran runtime error: No such file or directory
>>> /
>>> I saw the previous posting on the same error and I've tried using NCL
>>> version 6.0.0, 6.0.0-beta and 6.1.0 but they all give the same error
>>> message.
>>> I would appreciate any help, thank you.
>>>
>>>
>>> ______________________________**_________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<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 Fri Jun 1 13:29:41 2012