Re: Problem using conditional statements (if then else) in NCL

Hi,

Thanks for the help!

I've tryed to do as you both told me but it still didn't work. the "any"
function caused NCL to make all values equal (=100) and spliting the “if “
statement in two cause an "end if error"!
Maybe it can be a problem in my NCL version (5.1.1)?

Rosa

2009/7/24 Dennis Shea <shea_at_ucar.edu>

> Maybe, the "any" function
> http://www.ncl.ucar.edu/Document/Functions/Built-in/any.shtml
>
> if (.not.all(ismissing(umrl)) .and. any(umrl .gt. 100)) then
>
> Also, rather than
>
> > ;******************* criar as novas variaveis TS e ur ****************
> > TS = new(dimsizes(temp),float)
> > TS!0 = "time"
> > TS!1 = "lev"
> > TS!2 = "lat"
> > TS!3 = "lon"
> > TS&time = time
> > TS&lev = lev
> > TS&lat = lat
> > TS&lon = lon
> > TS_at_long_name = "temperatura do ar a superficie"
> > TS_at_units = "oC"
> > TS_at_grid_type = "gaussian"
> > TS = temp - 273.15
>
> I would suggest, either of the following for 'cleaner' code
>
> TS = temp ; create a copy including all meta data
> TS = temp - 273.15
> TS_at_long_name = "temperatura do ar a superficie"
> TS_at_units = "oC"
> TS_at_grid_type = "gaussian"
>
> or
> TS = temp - 273.15
> copy_VarCoords(temp, TS) ; copy coordinate information
> TS_at_long_name = "temperatura do ar a superficie"
> TS_at_units = "oC"
> TS_at_grid_type = "gaussian"
>
> The same for 'ur'
>
> Rosa Maria Nascimento dos Santos wrote:
>
>> Hi,
>>
>> I have a problem when trying to use "if then else" statments i my script
>> it gets the following error mensage: "fatal:Conditional statements (if and
>> do while) require SCALAR logical values, see all and any functions". I've
>> tryed to use the functions "all" and "ismissing" to solve it but it didn't
>> work. What can i do to solve this?! Could anyone help me?!
>>
>> Thanks!
>>
>> Rosa
>> P.S.: I'm sending a part of the script bellow!
>> ====================================================
>> begin
>>
>> file_in = addfile("/home/rosa/CCSM3/dados_paulo/
>> BRAMS.1997_mes06_dia18_hora18.nc <http://BRAMS.1997_mes06_dia18_hora18.nc
>> >","r")
>>
>>
>> att_names = getvaratts(file_in) ; get cdf_file file's global
>> attributes
>> psnm = file_in->psnm
>> uvel = file_in->uvel
>> vvel = file_in->vvel
>> zgeo = file_in->zgeo
>> temp = file_in->temp
>> umrl = file_in->umrl
>> time = file_in->time
>> lev = file_in->lev
>> lat = file_in->lat
>> lon = file_in->lon
>>
>> ;******************* criar as novas variaveis TS e ur ****************
>> TS = new(dimsizes(temp),float)
>> TS!0 = "time"
>> TS!1 = "lev"
>> TS!2 = "lat"
>> TS!3 = "lon"
>> TS&time = time
>> TS&lev = lev
>> TS&lat = lat
>> TS&lon = lon
>> TS_at_long_name = "temperatura do ar a superficie"
>> TS_at_units = "oC"
>> TS_at_grid_type = "gaussian"
>> TS = temp - 273.15
>> ;============================================================
>> ur = new(dimsizes(umrl),float)
>> ur!0 = "time"
>> ur!1 = "lev"
>> ur!2 = "lat"
>> ur!3 = "lon"
>> ur&time = time
>> ur&lev = lev
>> ur&lat = lat
>> ur&lon = lon
>> ur_at_long_name = "temperatura do ar a superficie"
>> ur_at_units = "oC"
>> ur_at_grid_type = "gaussian"
>>
>> if (.not.all(ismissing(umrl)) .and. umrl .gt. 100) then
>> ur = 100
>> else
>> ur = umrl
>> end if
>>
>> print(ur)
>> ==================================================
>>
>> --
>> Rosa Maria Nascimento Santos
>> Profra. Dra. do Curso de Meteorologia
>> EST / UEA
>> Av. Darcy Vargas, 1200 - Parque 10
>> CEP 690050-020, Manaus – AM – BR
>> Fone: +55 92 32365573 Ramal 22
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>

-- 
Rosa Maria Nascimento Santos
Profra. Dra. do Curso de Meteorologia
EST / UEA
Av. Darcy Vargas, 1200  -  Parque 10
CEP 690050-020, Manaus – AM – BR
Fone: +55 92 32365573 Ramal 22

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jul 24 2009 - 08:33:56 MDT