Changeset 347 for palm/trunk/SCRIPTS
 Timestamp:
 Jul 7, 2009 2:21:38 PM (14 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SCRIPTS/NCL/profiles.ncl
r342 r347 408 408 vNam_static = new((/dim/),string) 409 409 vNam_temp = new((/dim/),string) 410 410 411 411 var_char = stringtocharacter(var) 412 412 no_char = dimsizes(var_char) … … 443 443 444 444 do j=0,comma2 445 vNam_static(j) = vNam_temp(comma  2  j) 445 count_check = 0 446 do varn=0,dim1 447 if(vNam_temp(j) .ne. vNam(varn)) 448 count_check=count_check+1 449 end if 450 end do 451 452 if(count_check .eq. dim) 453 print(" ") 454 print("The variables 'var="+var+"'" ) 455 print("do not exist on your input file;") 456 print("be sure to have one comma before and after each variable") 457 print(" ") 458 exit 459 end if 460 461 vNam_static(j) = vNam_temp(comma  2  j) 446 462 end do 447 463 … … 881 897 end if 882 898 883 if (prof3d .EQ. 1)then 884 t_all = f[:]>time 885 nt = dimsizes(t_all) 886 else 887 t_all = f[:]>time 888 nt = dimsizes(t_all) 889 t_null = False 890 if (t_all(0) .EQ. 0 .AND. nt .GT. 1)then 891 delete(t_all) 892 t_all = f[:]>time(1:) 893 nt = dimsizes(t_all) 894 t_null = True 895 end if 896 end if 899 900 t_all = f[:]>time 901 nt = dimsizes(t_all) 902 897 903 898 904 if (nt .EQ. 1)then … … 1196 1202 if (log_z .EQ. 1) then 1197 1203 data = new((/dim,(end_time_stepstart_time_step)+1,dimz1/),float) 1204 data@_FillValue=9.96921e+36 1198 1205 data_0 = new((/(end_time_stepstart_time_step)+1,dimz1/),float) 1199 1206 data_0 = 0.1 … … 1226 1233 else 1227 1234 data = new((/dim,(end_time_stepstart_time_step)+1,dimz/),float) 1235 data@_FillValue=9.96921e+36 1228 1236 data_0 = new((/(end_time_stepstart_time_step)+1,dimz/),float) 1229 1237 data_0 = 0.0 … … 1490 1498 1491 1499 res@xyDashPattern = 1*nof 1500 1501 1492 1502 1493 1503 do varn = 0,dim1 … … 1512 1522 com=isStrSubset(c_var,","+vNam(varn)+"," ) 1513 1523 if (com) then 1514 if (prof3d .EQ. 0) then 1515 if (t_null)then 1516 temp = f[:]>$vNam(varn)$(1:,:) 1517 else 1518 temp = f[:]>$vNam(varn)$ 1519 end if 1524 if (prof3d .EQ. 0) then 1525 temp = f[:]>$vNam(varn)$ 1520 1526 temp_att = f_att>$vNam(varn)$ 1521 1527 if (log_z .EQ. 1) then … … 1586 1592 count_var=count_var+1 1587 1593 1588 if (prof3d .EQ. 0) then 1589 if (t_null)then 1590 temp = f[:]>$vNam(varn)$(1:,:) 1591 else 1592 temp = f[:]>$vNam(varn)$ 1593 end if 1594 if (prof3d .EQ. 0) then 1595 temp = f[:]>$vNam(varn)$ 1594 1596 temp_att = f_att>$vNam(varn)$ 1595 1597 else … … 2196 2198 end if 2197 2199 end if 2198 if (vNam(varn) .EQ. "wpt") then 2200 if (vNam(varn) .EQ. "wpt") then 2199 2201 miniwpt=min(data(varn,:,:))min_value 2200 2202 maxiwpt=max(data(varn,:,:))+max_value
Note: See TracChangeset
for help on using the changeset viewer.