Changeset 236 for palm/trunk/SCRIPTS/NCL/profiles.ncl

Timestamp:

Feb 12, 2009 2:58:16 PM (14 years ago)

Author:

letzel

Message:

• NCL scripts in trunk/SCRIPTS/NCL updated

File:

• palm/trunk/SCRIPTS/NCL/profiles.ncl (modified) (5 diffs)

palm/trunk/SCRIPTS/NCL/profiles.ncl

@@ -798 +798 @@
    if (prof3d .EQ. 1)then
       t_all = f[:]->time
+      nt    = dimsizes(t_all)
    else
-      t_all = f[:]->time(1:) 
-   end if
-   nt    = dimsizes(t_all)
+      t_all = f[:]->time
+      nt    = dimsizes(t_all)
+      t_null = False
+      if (t_all(0) .EQ. 0 .AND. nt .GT. 1)then
+         delete(t_all)
+         t_all = f[:]->time(1:)
+         nt    = dimsizes(t_all)
+         t_null = True
+      end if
+   end if
+   
    delta_t=t_all(nt-1)/nt
+   
+   
+   print(nt)
+   print(t_all)
+   
+   
+   ;exit
 
    ; ****************************************************       
@@ -898 +914 @@
    legend_label=new(nt,double)
    do p=start_time_step,end_time_step
-      if (t_all(p)/3600 .LT. 1) then
-         legend_label(p-start_time_step)=decimalPlaces(t_all(p)/3600,3,True)
-      else
-         legend_label(p-start_time_step)=decimalPlaces(t_all(p)/3600,0,True)
-      end if
+      legend_label(p-start_time_step)=decimalPlaces(t_all(p)/3600,3,True)         
    end do
    
@@ -1416 +1428 @@
          if (com) then      
             if (prof3d .EQ. 0) then
-               temp = f[:]->$vNam(varn)$(1:,:)
+               if (t_null)then
+                  temp = f[:]->$vNam(varn)$(1:,:)
+               else
+                  temp = f[:]->$vNam(varn)$
+               end if
                temp_att = f_att->$vNam(varn)$ 
                if (log_z .EQ. 1) then
@@ -1471 +1487 @@
 
          if (prof3d .EQ. 0) then
-            temp = f[:]->$vNam(varn)$(1:,:)
+            if (t_null)then
+               temp = f[:]->$vNam(varn)$(1:,:)
+            else
+               temp = f[:]->$vNam(varn)$
+            end if
             temp_att = f_att->$vNam(varn)$ 
          else
@@ -1534 +1554 @@
             end do            
          end if
-
-         if (data(varn,0,0) .GT. 1.E35) then
-            data(varn,0,:) = data(varn,1,:)
-         end if
+         
          if (nof .EQ. 0) then
             z_(n,:)=z/norm_z