Ignore:
Timestamp:
Feb 12, 2009 2:58:16 PM (15 years ago)
Author:
letzel
Message:
  • NCL scripts in trunk/SCRIPTS/NCL updated
File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/NCL/profiles.ncl

    r219 r236  
    798798   if (prof3d .EQ. 1)then
    799799      t_all = f[:]->time
     800      nt    = dimsizes(t_all)
    800801   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   
    804813   delta_t=t_all(nt-1)/nt
     814   
     815   
     816   print(nt)
     817   print(t_all)
     818   
     819   
     820   ;exit
    805821
    806822   ; ****************************************************       
     
    898914   legend_label=new(nt,double)
    899915   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)         
    905917   end do
    906918   
     
    14161428         if (com) then     
    14171429            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
    14191435               temp_att = f_att->$vNam(varn)$
    14201436               if (log_z .EQ. 1) then
     
    14711487
    14721488         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
    14741494            temp_att = f_att->$vNam(varn)$
    14751495         else
     
    15341554            end do           
    15351555         end if
    1536 
    1537          if (data(varn,0,0) .GT. 1.E35) then
    1538             data(varn,0,:) = data(varn,1,:)
    1539          end if
     1556         
    15401557         if (nof .EQ. 0) then
    15411558            z_(n,:)=z/norm_z
Note: See TracChangeset for help on using the changeset viewer.