Changeset 347 for palm/trunk/SCRIPTS/NCL


Ignore:
Timestamp:
Jul 7, 2009 2:21:38 PM (15 years ago)
Author:
heinze
Message:

profiles.ncl now plots initial profiles if available.

File:
1 edited

Legend:

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

    r342 r347  
    408408      vNam_static = new((/dim/),string)
    409409      vNam_temp   = new((/dim/),string)
    410      
     410
    411411      var_char    = stringtocharacter(var)
    412412      no_char     = dimsizes(var_char)
     
    443443
    444444      do j=0,comma-2
    445          vNam_static(j) = vNam_temp(comma - 2 - j)
     445        count_check = 0
     446        do varn=0,dim-1
     447           if(vNam_temp(j) .ne. vNam(varn))
     448             count_check=count_check+1
     449           end if
     450        end do   
     451
     452        if(count_check .eq. dim)
     453            print(" ")
     454            print("The variables 'var="+var+"'" )
     455            print("do not exist on your input file;")
     456            print("be sure to have one comma before and after each variable")
     457            print(" ")
     458            exit
     459        end if
     460
     461        vNam_static(j) = vNam_temp(comma - 2 - j)
    446462      end do
    447463
     
    881897   end if
    882898 
    883    if (prof3d .EQ. 1)then
    884       t_all = f[:]->time
    885       nt    = dimsizes(t_all)
    886    else
    887       t_all = f[:]->time
    888       nt    = dimsizes(t_all)
    889       t_null = False
    890       if (t_all(0) .EQ. 0 .AND. nt .GT. 1)then
    891          delete(t_all)
    892          t_all = f[:]->time(1:)
    893          nt    = dimsizes(t_all)
    894          t_null = True
    895       end if
    896    end if
     899
     900   t_all = f[:]->time
     901   nt    = dimsizes(t_all)
     902
    897903   
    898904   if (nt .EQ. 1)then
     
    11961202   if (log_z .EQ. 1) then
    11971203      data   = new((/dim,(end_time_step-start_time_step)+1,dimz-1/),float)
     1204      data@_FillValue=9.96921e+36
    11981205      data_0 = new((/(end_time_step-start_time_step)+1,dimz-1/),float)
    11991206      data_0 = 0.1
     
    12261233   else
    12271234      data   = new((/dim,(end_time_step-start_time_step)+1,dimz/),float)
     1235      data@_FillValue=9.96921e+36
    12281236      data_0 = new((/(end_time_step-start_time_step)+1,dimz/),float)
    12291237      data_0 = 0.0
     
    14901498
    14911499   res@xyDashPattern = 1*nof
     1500
     1501   
    14921502   
    14931503   do varn = 0,dim-1
     
    15121522         com=isStrSubset(c_var,","+vNam(varn)+"," )     
    15131523         if (com) then     
    1514             if (prof3d .EQ. 0) then
    1515                if (t_null)then
    1516                   temp = f[:]->$vNam(varn)$(1:,:)
    1517                else
    1518                   temp = f[:]->$vNam(varn)$
    1519                end if
     1524            if (prof3d .EQ. 0) then                 
     1525               temp     = f[:]->$vNam(varn)$         
    15201526               temp_att = f_att->$vNam(varn)$
    15211527               if (log_z .EQ. 1) then
     
    15861592         count_var=count_var+1
    15871593
    1588          if (prof3d .EQ. 0) then
    1589             if (t_null)then
    1590                temp = f[:]->$vNam(varn)$(1:,:)
    1591             else
    1592                temp = f[:]->$vNam(varn)$
    1593             end if
     1594         if (prof3d .EQ. 0) then         
     1595            temp = f[:]->$vNam(varn)$
    15941596            temp_att = f_att->$vNam(varn)$
    15951597         else
     
    21962198               end if
    21972199            end if
    2198             if (vNam(varn) .EQ. "wpt") then
     2200            if (vNam(varn) .EQ. "wpt") then       
    21992201               miniwpt=min(data(varn,:,:))-min_value
    22002202               maxiwpt=max(data(varn,:,:))+max_value
Note: See TracChangeset for help on using the changeset viewer.