#!/bin/csh -f # +-======-+ # Copyright (c) 2003-2007 United States Government as represented by # the Admistrator of the National Aeronautics and Space Administration. # All Rights Reserved. # # THIS OPEN SOURCE AGREEMENT ("AGREEMENT") DEFINES THE RIGHTS OF USE, # REPRODUCTION, DISTRIBUTION, MODIFICATION AND REDISTRIBUTION OF CERTAIN # COMPUTER SOFTWARE ORIGINALLY RELEASED BY THE UNITED STATES GOVERNMENT AS # REPRESENTED BY THE GOVERNMENT AGENCY LISTED BELOW ("GOVERNMENT AGENCY"). # THE UNITED STATES GOVERNMENT, AS REPRESENTED BY GOVERNMENT AGENCY, IS AN # INTENDED THIRD-PARTY BENEFICIARY OF ALL SUBSEQUENT DISTRIBUTIONS OR # REDISTRIBUTIONS OF THE SUBJECT SOFTWARE. ANYONE WHO USES, REPRODUCES, # DISTRIBUTES, MODIFIES OR REDISTRIBUTES THE SUBJECT SOFTWARE, AS DEFINED # HEREIN, OR ANY PART THEREOF, IS, BY THAT ACTION, ACCEPTING IN FULL THE # RESPONSIBILITIES AND OBLIGATIONS CONTAINED IN THIS AGREEMENT. # # Government Agency: National Aeronautics and Space Administration # Government Agency Original Software Designation: GSC-15354-1 # Government Agency Original Software Title: GEOS-5 GCM Modeling Software # User Registration Requested. Please Visit http://opensource.gsfc.nasa.gov # Government Agency Point of Contact for Original Software: # Dale Hithon, SRA Assistant, (301) 286-2691 # # +-======-+ setenv ORIENTATION PORTRAIT @ nmax = $#argv # Parse command line # ------------------ set batch = " " set debug = "no_debug" @ n = 1 while( $n <= $nmax ) if( "$argv[$n]" == "-expid" ) then @ n = $n + 1 set expid = $argv[$n] endif if( "$argv[$n]" == "-source" ) then @ n = $n + 1 set source = $argv[$n] endif if( "$argv[$n]" == "-plot" ) then @ n = $n + 1 set plot = $argv[$n] endif if( "$argv[$n]" == "-output" ) then @ n = $n + 1 set output = $argv[$n] endif if( "$argv[$n]" == "-geosutil" ) then @ n = $n + 1 set geosutil = $argv[$n] endif if( "$argv[$n]" == "-debug" ) then set debug = "debug" endif if( "$argv[$n]" == "-grads" ) then @ n = $n + 1 set grads = $argv[$n] endif if( "$argv[$n]" == "-batch" ) then set batch = "-b" endif if( "$argv[$n]" == "-season" ) then @ n = $n + 1 set seasons = "" set next = $argv[$n] set bit = `echo $next | cut -c1-1` while( "$bit" != "-" ) set seasons = `echo $seasons $next` @ n = $n + 1 if( $n <= $nmax ) then set next = $argv[$n] set bit = `echo $next | cut -c1-1` else set bit = "-" endif end @ n = $n - 1 endif @ n = $n + 1 end if( ${STD_DEV} == 'true' ) setenv DO_STD TRUE set TAYLOR = `echo $TAYLOR | tr "[:lower:]" "[:upper:]"` echo "Portrait Script" echo "---------------" echo " expid: $expid" echo " source: $source" echo " plot: $plot" echo " output: $output" echo "geosutil: $geosutil" echo " debug: $debug" echo " season: $seasons" echo " TAYLOR: $TAYLOR" echo " STD_DEV: $STD_DEV" # Taylor Plots # ------------ if( $TAYLOR == "TRUE" & $plot == "taylor" ) then $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_olr $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_cloud $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_lwp $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_tpw $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_tskin $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_uflux $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_eflux $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_hflux $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog U DYN 850 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog V DYN 850 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog T DYN 850 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog Q MOIST 850 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog U DYN 500 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog V DYN 500 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog T DYN 500 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog Q MOIST 500 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog U DYN 200 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog V DYN 200 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog T DYN 200 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog Q MOIST 200 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog ZLE DYN 300 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog ZLE DYN 500 $expid $output $debug " $grads $batch -p -c "run $GEOSUTIL/plots/taylor/taylor_prog SLP DYN 1000 $expid $output $debug " endif # Loop over Variables for Zonal Mean Comparisons # ---------------------------------------------- set EXPORTS = `echo U V T Q RH2 ZLE OMEGA` set GCS = `echo DYN DYN DYN MOIST MOIST DYN DYN` @ n = 1 foreach EXPORT ($EXPORTS) set GC = $GCS[$n] set this_plot = "progz2_$EXPORT" if( $plot == ${this_plot} ) then $grads $batch -p -c "run $GEOSUTIL/plots/zcmp/progz $EXPORT $GC $expid $output $debug $seasons " endif @ n = $n + 1 end # Compute Horizontal Mean Comparisons for SLP # ------------------------------------------- set EXPORTS = `echo SLP` set GCS = `echo DYN` set LEVS = "1000" @ n = 1 foreach EXPORT ($EXPORTS) set GC = $GCS[$n] set this_plot = "progh_$EXPORT" if( $plot == ${this_plot} ) then foreach LEV ($LEVS) $grads $batch -p -c "run $GEOSUTIL/plots/hcmp/hcmp2 $EXPORT $GC $LEV $expid $output $debug $seasons " end endif @ n = $n + 1 end # Loop over Variables for Horizontal Mean Comparisons # --------------------------------------------------- set EXPORTS = `echo U V T ZLE HE Q RH2 OMEGA` set GCS = `echo DYN DYN DYN DYN DYN MOIST MOIST DYN ` set LEVS = "1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30" @ n = 1 foreach EXPORT ($EXPORTS) set GC = $GCS[$n] set this_plot = "progh_$EXPORT" if( $plot == ${this_plot} ) then foreach LEV ($LEVS) $grads $batch -p -c "run $GEOSUTIL/plots/hcmp/hcmp2 $EXPORT $GC $LEV $expid $output $debug $seasons " end endif @ n = $n + 1 end # Loop over Levels for Quadratics Plots # --------------------------------------------------- if( $plot == hquad ) then set LEVS = "850 500 200 20" foreach LEV ($LEVS) $grads $batch -p -c "run $GEOSUTIL/plots/quad/hquad -EXPID $expid -LEV $LEV -OUTPUT $output -DEBUG $debug -SEASON $seasons " end endif # Surface/2-D # ----------- if( $plot == surf2.misc ) then set EXPORTS = "CN_PRCP:MOIST LS_PRCP:MOIST AN_PRCP:MOIST" set OBS = "TPREC:MOIST" $grads $batch -p -c "run genplot -EXPORT $EXPORTS -OBS $OBS -NAME precip -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 86400 -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT $EXPORTS -OBS $OBS -NAME precip -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 86400 -PREFIX LAND -LAND -SEASON $seasons " $grads $batch -p -c "run $GEOSUTIL/plots/precip/precip $expid $output $debug $seasons " $grads $batch -p -c "run $GEOSUTIL/plots/olr/olr $expid $output $debug $seasons " $grads $batch -p -c "run genplot -EXPORT SH:SURFACE -NAME hflux -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 1 -OCEAN -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT LHFX:SURFACE -NAME eflux -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 1 -OCEAN -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT LWP:MOIST CCWP:MOIST -OBS LWP:MOIST -NAME lwp -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 1000 -OCEAN -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT TQV:AGCM -NAME tpw -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 0.1 -RC SRB -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT TQV:AGCM -NAME tpw -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 0.1 -OCEAN -RC SSMI -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT VENT:SURFACE -NAME speed -EXPID $expid -OUTPUT $output -DEBUG $debug -OCEAN -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT TS:SURFACE -NAME sst -EXPID $expid -OUTPUT $output -DEBUG $debug -OCEAN -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT TS:SURFACE -NAME tskin -EXPID $expid -OUTPUT $output -DEBUG $debug -LAND -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT TAUX:SURFACE -NAME taux -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 100 -OCEAN -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT TAUY:SURFACE -NAME tauy -EXPID $expid -OUTPUT $output -DEBUG $debug -SCALE 100 -OCEAN -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT ALBEDO:SOLAR -NAME albedo -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME NETHEATG -EXPID $expid -EXPORT RADSRF:RADIATION LHFX:SURFACE SH:SURFACE -OBS QNET:RADIATION -OUTPUT $output -DEBUG $debug -SEASON $seasons " endif # Surface/2-D: Clouds # ------------------- if( $plot == surf2.cld ) then $grads $batch -p -c "run genplot -EXPORT CLDLO:SOLAR -NAME cldlo -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT CLDMD:SOLAR -NAME cldmd -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT CLDHI:SOLAR -NAME cldhi -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT CLDTT:SOLAR -NAME cldtt -EXPID $expid -OUTPUT $output -DEBUG $debug -RC SRB -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT CLDTT:SOLAR -NAME cldtt -EXPID $expid -OUTPUT $output -DEBUG $debug -RC ISCCP -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT CLDTT:SOLAR -NAME cldtt -EXPID $expid -OUTPUT $output -DEBUG $debug -OCEAN -RC COADS -SEASON $seasons " endif # Radiation: Shortwave # -------------------- if( $plot == rad.sw ) then $grads $batch -p -c "run genplot -NAME radswt -EXPID $expid -EXPORT SLRTP:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME swgnet -EXPID $expid -EXPORT RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME swgnetc -EXPID $expid -EXPORT RSCS:SOLAR -OBS SLRSFC:SOLAR ALBEDO:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME swgnetcf -EXPID $expid -EXPORT RSCS:SOLAR RSRS:SOLAR -OBS SLRSFC:SOLAR ALBEDO:SOLAR RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME swgdwn -EXPID $expid -EXPORT SLRSF:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME swgdwnc -EXPID $expid -EXPORT SLRSFC:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME swgdwncf -EXPID $expid -EXPORT SLRSFC:SOLAR SLRSF:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME osr -EXPID $expid -EXPORT RSR:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME osrc -EXPID $expid -EXPORT RSC:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME osrcf -EXPID $expid -EXPORT RSC:SOLAR RSR:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME swgup -EXPID $expid -EXPORT SLRSF:SOLAR RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME swgupc -EXPID $expid -EXPORT SLRSFC:SOLAR RSCS:SOLAR -OBS ALBEDO:SOLAR SLRSFC:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -NAME swgupcf -EXPID $expid -EXPORT SLRSFC:SOLAR RSCS:SOLAR SLRSF:SOLAR RSRS:SOLAR -OBS ALBEDO:SOLAR SLRSFC:SOLAR SLRSF:SOLAR RSRS:SOLAR -OUTPUT $output -DEBUG $debug -SEASON $seasons " endif # Radiation: Longwave # ------------------- if( $plot == rad.lw ) then $grads $batch -p -c "run genplot -EXPORT OLC:IRRAD OLR:IRRAD OSRCLR:SOLAR OSR:SOLAR -NAME toacf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT OLR:IRRAD -NAME olr -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT OLC:IRRAD -NAME olrc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT OLC:IRRAD OLR:IRRAD -NAME olrcf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD EMIS:SURFACE -OBS SFCDWN:IRRAD -NAME lwgdwn -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT LCS:IRRAD EMIS:SURFACE -OBS SFCDWNC:IRRAD -NAME lwgdwnc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCS:IRRAD EMIS:SURFACE -OBS SFCDWN:IRRAD SFCDWNC:IRRAD -NAME lwgdwncf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD LWS:IRRAD -OBS SFCUP:IRRAD SFCDWN:IRRAD -NAME lwgnet -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD LCS:IRRAD -OBS SFCUPC:IRRAD SFCDWNC:IRRAD -NAME lwgnetc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCS:IRRAD -OBS SFCUPC:IRRAD SFCDWNC:IRRAD SFCUP:IRRAD SFCDWN:IRRAD -NAME lwgnetcf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD EMIS:SURFACE LWS:IRRAD -OBS SFCUP:IRRAD -NAME lwgup -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT SFCEM:IRRAD EMIS:SURFACE LCS:IRRAD -OBS SFCUPC:IRRAD -NAME lwgupc -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " $grads $batch -p -c "run genplot -EXPORT LWS:IRRAD LCS:IRRAD EMIS:SURFACE -OBS SFCUP:IRRAD SFCUPC:IRRAD -NAME lwgupcf -EXPID $expid -OUTPUT $output -DEBUG $debug -SEASON $seasons " endif exit