Changeset 347

Timestamp:

Jul 7, 2009 2:21:38 PM (15 years ago)

Author:

heinze

Message:

profiles.ncl now plots initial profiles if available.

Location:

palm/trunk

Files:

• SCRIPTS/NCL/profiles.ncl (modified) (9 diffs)
• SOURCE/CURRENT_MODIFICATIONS (modified) (1 diff)

palm/trunk/SCRIPTS/NCL/profiles.ncl

@@ -408 +408 @@
       vNam_static = new((/dim/),string)
       vNam_temp   = new((/dim/),string)
-      
+
       var_char    = stringtocharacter(var)
       no_char     = dimsizes(var_char)
@@ -443 +443 @@
 
       do j=0,comma-2
-         vNam_static(j) = vNam_temp(comma - 2 - j)
+        count_check = 0
+        do varn=0,dim-1 
+           if(vNam_temp(j) .ne. vNam(varn)) 
+             count_check=count_check+1
+           end if
+        end do    
+
+        if(count_check .eq. dim)
+            print(" ")
+            print("The variables 'var="+var+"'" )
+            print("do not exist on your input file;")
+            print("be sure to have one comma before and after each variable")
+            print(" ")
+            exit
+        end if
+
+        vNam_static(j) = vNam_temp(comma - 2 - j)
       end do
 
@@ -881 +897 @@
    end if
   
-   if (prof3d .EQ. 1)then
-      t_all = f[:]->time
-      nt    = dimsizes(t_all)
-   else
-      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
+
+   t_all = f[:]->time
+   nt    = dimsizes(t_all)
+
    
    if (nt .EQ. 1)then
@@ -1196 +1202 @@
    if (log_z .EQ. 1) then
       data   = new((/dim,(end_time_step-start_time_step)+1,dimz-1/),float) 
+      data@_FillValue=9.96921e+36
       data_0 = new((/(end_time_step-start_time_step)+1,dimz-1/),float)
       data_0 = 0.1
@@ -1226 +1233 @@
    else
       data   = new((/dim,(end_time_step-start_time_step)+1,dimz/),float) 
+      data@_FillValue=9.96921e+36
       data_0 = new((/(end_time_step-start_time_step)+1,dimz/),float)
       data_0 = 0.0
@@ -1490 +1498 @@
 
    res@xyDashPattern = 1*nof
+
+   
    
    do varn = 0,dim-1 
@@ -1512 +1522 @@
          com=isStrSubset(c_var,","+vNam(varn)+"," )      
          if (com) then      
-            if (prof3d .EQ. 0) then
-               if (t_null)then
-                  temp = f[:]->$vNam(varn)$(1:,:)
-               else
-                  temp = f[:]->$vNam(varn)$
-               end if
+            if (prof3d .EQ. 0) then                  
+               temp     = f[:]->$vNam(varn)$          
                temp_att = f_att->$vNam(varn)$ 
                if (log_z .EQ. 1) then
@@ -1586 +1592 @@
          count_var=count_var+1
 
-         if (prof3d .EQ. 0) then
-            if (t_null)then
-               temp = f[:]->$vNam(varn)$(1:,:) 
-            else
-               temp = f[:]->$vNam(varn)$
-            end if
+         if (prof3d .EQ. 0) then          
+            temp = f[:]->$vNam(varn)$
             temp_att = f_att->$vNam(varn)$ 
          else
@@ -2196 +2198 @@
                end if
             end if
-            if (vNam(varn) .EQ. "wpt") then
+            if (vNam(varn) .EQ. "wpt") then       
                miniwpt=min(data(varn,:,:))-min_value
                maxiwpt=max(data(varn,:,:))+max_value

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -46 +46 @@
 Changed:
 -------
+In case of restart runs without extension, initial profiles are not written 
+to NetCDF-file anymore. (data_output_profiles, modules, read_var_list, write_var_list)
+
 Small change in formatting of the message handling routine concering the output in the
 job protocoll. (message)