Changeset 329 for palm/trunk/SCRIPTS/NCL

Timestamp:

Jun 2, 2009 2:16:11 PM (15 years ago)

Author:

heinze

Message:

Small change in formatting in message.f90. Bugfix in cross_sections.ncl in case of normalizing.

File:

• palm/trunk/SCRIPTS/NCL/cross_sections.ncl (modified) (12 diffs)

palm/trunk/SCRIPTS/NCL/cross_sections.ncl

@@ -1019 +1019 @@
 
    if (xyc .EQ. 1) then
-      d=(ye-ys+1)/(major_ticks_y-1)
+      d= (ye-ys+1)/(major_ticks_y-1)
       e=(xe-xs+1)/(major_ticks_x-1)
       array_yl       =new(major_ticks_y,integer)
@@ -1036 +1036 @@
          if (ar .GT. 0)
             do min_ar=0,3
-               array_minor_yl(4*(ar-1)+min_ar)= int2dble(array_yl(ar-1))+int2dble(array_yl(ar)-array_yl(ar-1))/5*(min_ar+1)
+                array_minor_yl(4*(ar-1)+min_ar)= y_d(array_yl(ar-1))+(y_d(array_yl(ar))-y_d(array_yl(ar-1)))/5*(min_ar+1)
             end do 
          end if  
@@ -1045 +1045 @@
          if (br .GT. 0)
             do min_br=0,3
-               array_minor_xb(4*(br-1)+min_br)= int2dble(array_xb(br-1))+int2dble(array_xb(br)-array_xb(br-1))/5*(min_br+1)
+               array_minor_xb(4*(br-1)+min_br)= x_d(array_xb(br-1))+(x_d(array_xb(br))-x_d(array_xb(br-1)))/5*(min_br+1)
             end do 
          end if 
@@ -1069 +1069 @@
          if (ar .GT. 0)
             do min_ar=0,3
-               array_minor_yl(4*(ar-1)+min_ar)= int2dble(array_yl(ar-1))+int2dble(array_yl(ar)-array_yl(ar-1))/5*(min_ar+1)
+               array_minor_yl(4*(ar-1)+min_ar)= z_d(array_yl(ar-1))+(z_d(array_yl(ar))-z_d(array_yl(ar-1)))/5*(min_ar+1)
             end do 
          end if  
@@ -1078 +1078 @@
          if (br .GT. 0)
             do min_br=0,3
-               array_minor_xb(4*(br-1)+min_br)= int2dble(array_xb(br-1))+int2dble(array_xb(br)-array_xb(br-1))/5*(min_br+1)
+               array_minor_xb(4*(br-1)+min_br)= x_d(array_xb(br-1))+(x_d(array_xb(br))-x_d(array_xb(br-1)))/5*(min_br+1)
             end do 
          end if 
@@ -1099 +1099 @@
       do ar=1,major_ticks_y-1
          array_yl(ar)=d*(ar-1)+d-1
-         array_yl_labels(ar) = y_d(array_yl(ar))
+         array_yl_labels(ar) = z_d(array_yl(ar))
          if (ar .GT. 0)
             do min_ar=0,3
-               array_minor_yl(4*(ar-1)+min_ar)= int2dble(array_yl(ar-1))+int2dble(array_yl(ar)-array_yl(ar-1))/5*(min_ar+1)
+               array_minor_yl(4*(ar-1)+min_ar)= z_d(array_yl(ar-1))+(z_d(array_yl(ar))-z_d(array_yl(ar-1)))/5*(min_ar+1)
             end do 
          end if  
@@ -1108 +1108 @@
       do br=1,major_ticks_x-1
          array_xb(br)=e*(br-1)+e-1 
-         array_xb_labels(br) = x_d(array_xb(br)) 
+         array_xb_labels(br) = y_d(array_xb(br)) 
          if (br .GT. 0)
             do min_br=0,3
-               array_minor_xb(4*(br-1)+min_br)= int2dble(array_xb(br-1))+int2dble(array_xb(br)-array_xb(br-1))/5*(min_br+1)
+               array_minor_xb(4*(br-1)+min_br)= y_d(array_xb(br-1))+(y_d(array_xb(br))-y_d(array_xb(br-1)))/5*(min_br+1)
             end do 
          end if 
       end do
    end if
-
+ 
    if (axes_explicit .EQ. 1)then
       cs_res@tmYLMode = "Explicit"
       cs_res@tmXBMode = "Explicit"
-      cs_res@tmYLValues = array_yl
-      if (xyc .EQ. 1)then
+      if (xyc .EQ. 1)then 
+         cs_res@tmYLValues = y_d(array_yl) 
+         cs_res@tmXBValues = x_d(array_xb)      
          cs_res@tmYLLabels = array_yl_labels/norm_y
          cs_res@tmXBLabels = array_xb_labels/norm_x
@@ -1136 +1137 @@
       end if
       if (xzc .EQ. 1)then
+         cs_res@tmYLValues = z_d(array_yl) 
+         cs_res@tmXBValues = x_d(array_xb)  
          cs_res@tmYLLabels = array_yl_labels/norm_z
          cs_res@tmXBLabels = array_xb_labels/norm_x
@@ -1150 +1153 @@
       end if
       if (yzc .EQ. 1)then
+         cs_res@tmYLValues = z_d(array_yl) 
+         cs_res@tmXBValues = y_d(array_xb)  
          cs_res@tmYLLabels = array_yl_labels/norm_z
          cs_res@tmXBLabels = array_xb_labels/norm_y
@@ -1163 +1168 @@
          end if
       end if
-      cs_res@tmXBValues = array_xb      
       cs_res@tmYLMinorValues = array_minor_yl
       cs_res@tmXBMinorValues = array_minor_xb
@@ -1301 +1305 @@
             temp = f[:]->$vNam(varn)$
             data_att = f_att->$vNam(varn)$
-            data(varn,:,:,:,:)=temp(:,zs:ze,ys:ye,xs:xe)                
+            if(vNam(varn) .eq. "ts_xy" .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy")
+              ;these variables depend von zu1_xy and that's why they have only one z-layer
+              data(varn,:,0,:,:)=temp(:,0,ys:ye,xs:xe) 
+            else
+              data(varn,:,:,:,:)=temp(:,zs:ze,ys:ye,xs:xe) 
+            end if
+            delete(temp)                
          end if
          if ( xzc .eq. 1 ) then
@@ -1330 +1340 @@
          
          unit(varn) = data_att@units
+         delete(data_att)
 
       end if