Changeset 183 for palm/trunk/SCRIPTS

Timestamp:

Aug 4, 2008 3:39:12 PM (16 years ago)

Author:

letzel

Message:

• mrun: (optionally) link local input
• diffusion_s: bugfix for calculation of horizontal fluxes at vertical walls
• NCL scripts in trunk/SCRIPTS/NCL updated

Location:

palm/trunk/SCRIPTS

Files:

• NCL/spectra.ncl (modified) (5 diffs)
• mrun (modified) (3 diffs)

palm/trunk/SCRIPTS/NCL/spectra.ncl

@@ -13 +13 @@
       delete(parameter@_FillValue)
    else
-      if (isfilepresent("~/palm/current_version/trunk/SRIPTS/NCL/.ncl_preferences")) then
-         parameter = asciiread("~/palm/current_version/trunk/SRIPTS/NCL/.ncl_preferences",129,"string")
+      if (isfilepresent("~/palm/current_version/trunk/SCRIPTS/NCL/.ncl_preferences")) then
+         parameter = asciiread("~/palm/current_version/trunk/SCRIPTS/NCL/.ncl_preferences",129,"string")
          delete(parameter@_FillValue)
       else
          print(" ")
-         print("'.ncl_preferences' does not exist in '~/palm/current_version/trunk/SRIPTS/NCL/'")
+         print("'.ncl_preferences' does not exist in '~/palm/current_version/trunk/SCRIPTS/NCL/'")
          print(" ")
          exit
@@ -372 +372 @@
    if (sort .EQ. "time")
       plot  = new(dim*dimz,graphic)
-      plot_ = new(dim*dimz,graphic)
       np=dimt
       res@lgTitleString = "Time [h]"
    else
       plot  = new(dim*dimt,graphic)
-      plot_ = new(dim*dimt,graphic)
       np=dimz
       res@lgTitleString = "Height [m]"  
@@ -455 +453 @@
                res@trXMaxF = max(x_axis)
                res@gsnLeftString      = vNam(varn)
-               res@gsnRightString     = "Height: "+z(p)+"m"
+               res@gsnRightString     = "Height = "+z(p)+"m"
                if (norm .NE. 1)then
                   res@tiYAxisString      = vNam(varn)+" / "+norm
@@ -482 +480 @@
                   res@trXMaxF = max(x_axis)
                   res@gsnLeftString      = vNam(varn)
-                  res@gsnRightString     = "Time: "+legend_label(p)+"h"
+                  res@gsnRightString     = "Time = "+legend_label(p)+"h"
                   if (norm .NE. 1)then
                      res@tiYAxisString      = vNam(varn)+" / "+norm
@@ -521 +519 @@
    resP@txFontHeightF          = 0.014
 
-   do m=0,n-1
-      plot_(m)=plot(n-1-m)
-   end do
-
    if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
-      gsn_panel(wks,plot_,(/n,1/),resP)
+      gsn_panel(wks,plot,(/n,1/),resP)
    else    
       do i = 0,n-1, no_lines*no_columns
          if( (i+no_lines*no_columns) .gt. (n-1)) then
-            gsn_panel(wks,plot_(i:n-1),(/no_lines,no_columns/),resP)
+            gsn_panel(wks,plot(i:n-1),(/no_lines,no_columns/),resP)
          else
-           gsn_panel(wks,plot_(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP)
+           gsn_panel(wks,plot(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP)
          end if
       end do

palm/trunk/SCRIPTS/mrun

@@ -147 +147 @@
      #                     adjustments for lcxt4 (Bergen Center for Computational
      #                     Science)
+     # 22/05/08 - Marcus - If environment variable link_local_input is set to 
+     #                     true, mrun tries "ln -f" on local input and resorts 
+     #                     to "cp" or "cp -r" on error
      # 27/05/08 - Siggi  - PATH is set to PALM_BIN everywhere (missing so far)
      # 14/07/08 - Siggi  - adjustments for lcsgih
-
  
     # VARIABLENVEREINBARUNGEN + DEFAULTWERTE
@@ -196 +198 @@
  job_on_file=""
  keep_data_from_previous_run=false
+ link_local_input=false
  link_local_output=false
  localhost_realname=$(hostname)
@@ -2539 +2542 @@
                 printf "\n      providing $numprocs files for the respective processors"
                 mkdir  ${localin[$i]}
-                cp -r  ${absnamein[$i]}/*  ${localin[$i]}
+                if [[ $link_local_input = true ]]
+                    then
+                    printf "      using ln -f\n"
+                    cd ${absnamein[$i]}
+                    for file in $(ls *)
+                      do
+                      ln -f $file  ${localin[$i]}
+                    done
+                    cd $TEMPDIR
+                fi
+                # If "ln -f" fails of if "$link_local_input = false" do a normal "cp -r"
+                if [[ ! -f "${localin[$i]}/_0000" ]]
+                    then
+                    if [[ $link_local_input = true ]]
+                        then
+                        printf "      ln failed for .../_0000, using cp...\n"
+                    fi
+                    cp -r  ${absnamein[$i]}/*  ${localin[$i]}
+                fi
 
              else
                    # BEREITSTELLUNG AUF EINPROZESSORRECHNERN
-                cp  ${absnamein[$i]}  ${localin[$i]}
+                if [[ $link_local_input = true ]]
+                then
+                    printf "      using ln -f\n"
+                    ln -f  ${absnamein[$i]}  ${localin[$i]}
+                fi
+                # If "ln -f" fails of if "$link_local_input = false" do a normal "cp"
+                if [[ ! -f "${localin[$i]}" ]]
+                then
+                    if [[ $link_local_input = true ]]
+                    then
+                        printf "      ln failed, using cp...\n"
+                    fi
+                    cp  ${absnamein[$i]}  ${localin[$i]}
+                fi
              fi
           fi