Changeset 236 for palm/trunk/SCRIPTS/NCL/profiles.ncl
- Timestamp:
- Feb 12, 2009 2:58:16 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SCRIPTS/NCL/profiles.ncl
r219 r236 798 798 if (prof3d .EQ. 1)then 799 799 t_all = f[:]->time 800 nt = dimsizes(t_all) 800 801 else 801 t_all = f[:]->time(1:) 802 end if 803 nt = dimsizes(t_all) 802 t_all = f[:]->time 803 nt = dimsizes(t_all) 804 t_null = False 805 if (t_all(0) .EQ. 0 .AND. nt .GT. 1)then 806 delete(t_all) 807 t_all = f[:]->time(1:) 808 nt = dimsizes(t_all) 809 t_null = True 810 end if 811 end if 812 804 813 delta_t=t_all(nt-1)/nt 814 815 816 print(nt) 817 print(t_all) 818 819 820 ;exit 805 821 806 822 ; **************************************************** … … 898 914 legend_label=new(nt,double) 899 915 do p=start_time_step,end_time_step 900 if (t_all(p)/3600 .LT. 1) then 901 legend_label(p-start_time_step)=decimalPlaces(t_all(p)/3600,3,True) 902 else 903 legend_label(p-start_time_step)=decimalPlaces(t_all(p)/3600,0,True) 904 end if 916 legend_label(p-start_time_step)=decimalPlaces(t_all(p)/3600,3,True) 905 917 end do 906 918 … … 1416 1428 if (com) then 1417 1429 if (prof3d .EQ. 0) then 1418 temp = f[:]->$vNam(varn)$(1:,:) 1430 if (t_null)then 1431 temp = f[:]->$vNam(varn)$(1:,:) 1432 else 1433 temp = f[:]->$vNam(varn)$ 1434 end if 1419 1435 temp_att = f_att->$vNam(varn)$ 1420 1436 if (log_z .EQ. 1) then … … 1471 1487 1472 1488 if (prof3d .EQ. 0) then 1473 temp = f[:]->$vNam(varn)$(1:,:) 1489 if (t_null)then 1490 temp = f[:]->$vNam(varn)$(1:,:) 1491 else 1492 temp = f[:]->$vNam(varn)$ 1493 end if 1474 1494 temp_att = f_att->$vNam(varn)$ 1475 1495 else … … 1534 1554 end do 1535 1555 end if 1536 1537 if (data(varn,0,0) .GT. 1.E35) then 1538 data(varn,0,:) = data(varn,1,:) 1539 end if 1556 1540 1557 if (nof .EQ. 0) then 1541 1558 z_(n,:)=z/norm_z
Note: See TracChangeset
for help on using the changeset viewer.